On the role of global flow instability analysis in closed loop flow control

Control of linear flow instabilities has been demonstrated to be an effective theoretical flow control methodology, capable of modifying transitional flow on canonical geometries such as the plane channel and the flat-plate boundary layer.

Global Stability Analysis of a Compressible Turbulent Flow around a High-Lift Configuration

The purpose of this work is to analyze a complex high lift configuration for which significant regions of separated flow are present. Current state of the art methods have some diffculty to predict the origin and the progression of this separated flow when increasing the angle of attack. The mechanisms responsible for the maximum lift limit on multi-element wing con?gurations are not clear; this stability analysis could help to understand the physics behind the phenomenon and to find a relation between the flow separation and the instability onset. The methodology presented herein consists in the computation of a steady base flow solution based on a finite volume discretization and a proposal of the solution for a generalized eigenvalue problem corresponding to the perturbed and linearized problem. The eigenvalue problem has been solved with the Arnoldi iterative method, one of the Krylov subspace projection methods. The described methodology was applied to the NACA0012 test case in subsonic and in transonic conditions and, finally, for the first time to the authors knowledge, on an industrial multi-component geometry, such as the A310 airfoil, in order to identify low frequency instabilities related to the separation. One important conclusion is that for all the analyzed geometries, one unstable mode related to flow separation appears for an angle of attack greater than the one correspondent to the maximum lift coe?cient condition. Finally, an adjoint study was carried out in order to evaluate the receptivity and the structural sensitivity of the geometries, giving an indication of the domain region that could be modified resulting in the biggest change of the flowfield.

Analysis of air flow past and through the 2415-3S airfoil for an unmanned aerial vehicle with internal propulsion system

This paper deals with the prediction of velocity fields on the 2415-3S airfoil which will be used for an unmanned aerial vehicle with internal propulsion system and in this way analyze the air flow through an internal duct of the airfoil using computational fluid dynamics. The main objective is to evaluate the effect of the internal air flow past the airfoil and how this affects the aerodynamic performance by means of lift and drag forces. For this purpose, three different designs of the internal duct were studied; starting from the base 2415-3S airfoil developed in previous investigation, basing on the hypothesis of decreasing the flow separation produced when the propulsive airflow merges the external flow, and in this way obtaining the best configuration. For that purpose, an exhaustive study of the mesh sensitivity was performed. It was used a non-structured mesh since the computational domain is three-dimensional and complex. The selected mesh contains approximately 12.5 million elements. Both the computational domain and the numerical solution were made with commercial CAD and CFD software, respectively. Air, incompressible and steady was analyzed. The boundary conditions are in concordance with experimental setup in the AF 6109 wind tunnel. The k-e model is utilized to describe the turbulent flow process as followed in references. Results allowed obtaining velocity contours as well as lift and drag coefficients and also the location of separation and reattachment regions in some cases for zero degrees of angle of attack on the internal and external surfaces of the airfoil. Finally, the selection of the configuration with the best aerodynamic performance was made, selecting the option without curved baffles.

Numerical analysis of air flow past the 2415-3S airfoil for an unmanned aerial vehicle with internal propulsion system

This paper deals with the prediction of pressure and velocity fields on the 2415-3S airfoil which will be used for and unmanned aerial vehicle with internal propulsion system and in this way analyze the air flow through an internal duct of the airfoil using computational fluid dynamics. The main objective is to evaluate the effect of the internal air flow past the airfoil and how this affects the aerodynamic performance by means of lift and drag forces. For this purpose, three different designs of the internal duct were studied; starting from the base 2415-3S airfoil developed in previous investigation, basing on the hypothesis of decreasing the flow separation produced when the propulsive airflow merges the external flow, and in this way obtaining the best configuration. For that purpose, an exhaustive study of the mesh sensitivity was performed. It was used a non-structured mesh since the computational domain is tridimensional and complex. The selected mesh contains approximately 12.5 million elements. Both the computational domain and the numerical solution were made with commercial CAD and CFD software respectively. Air, incompressible and steady was analyzed. The boundary conditions are in concordance with experimental setup in the AF 6109 wind tunnel. The k-ε model is utilized to describe the turbulent flow process as followed in references. Results allowed obtaining pressure and velocity contours as well as lift and drag coefficients and also the location of separation and reattachment regions in some cases for zero degrees of angle of attack on the internal and external surfaces of the airfoil. Finally, the selection of the configuration with the best aerodynamic performance was made, selecting the option without curved baffles.

Pulsatile blood flow interaction with arterial walls of aorta : autoregulation and impedance pressure boundary condition and its biomedical applications

Para las decisiones urgentes sobre intervenciones quirúrgicas en el sistema cardiovascular se necesitan simulaciones computacionales con resultados fiables y que consuman un tiempo de cálculo razonable. Durante años los investigadores han trabajado en diversos métodos numéricos de cálculo que resulten atractivos para los cirujanos. Estos métodos, precisos pero costosos desde el punto de vista del coste computacional, crean un desajuste entre la oferta de los ingenieros que realizan las simulaciones y los médicos que operan en el quirófano. Por otra parte, los métodos de cálculo más simplificados reducen el tiempo de cálculo pero pueden proporcionar resultados no realistas. El objetivo de esta tesis es combinar los conceptos de autorregulación e impedancia del sistema circulatorio, la interacción flujo sanguíneo-pared arterial y modelos geométricos idealizados tridimensionales de las arterias pero sin pérdida de realismo, con objeto de proponer una metodología de simulación que proporcione resultados correctos y completos, con tiempos de cálculo moderados. En las simulaciones numéricas, las condiciones de contorno basadas en historias de presión presentan inconvenientes por ser difícil conocerlas con detalle, y porque los resultados son muy sensibles ante pequeñas variaciones de dichas historias. La metodología propuesta se basa en los conceptos de autorregulación, para imponer la demanda de flujo aguas abajo del modelo en el ciclo cardiaco, y la impedancia, para representar el efecto que ejerce el flujo en el resto del sistema circulatorio sobre las arterias modeladas. De este modo las historias de presión en el contorno son resultados del cálculo, que se obtienen de manera iterativa. El método propuesto se aplica en una geometría idealizada del arco aórtico sin patologías y en otra geometría correspondiente a una disección Stanford de tipo A, considerando la interacción del flujo pulsátil con las paredes arteriales. El efecto de los tejidos circundantes también se incorpora en los modelos. También se hacen aplicaciones considerando la interacción en una geometría especifica de un paciente anciano que proviene de una tomografía computarizada. Finalmente se analiza una disección Stanford tipo B con tres modelos que incluyen la fenestración del saco. Clinicians demand fast and reliable numerical results of cardiovascular biomechanic simulations for their urgent pre-surgery decissions. Researchers during many years have work on different numerical methods in order to attract the clinicians' confidence to their colorful contours. Though precise but expensive and time-consuming methodologies create a gap between numerical biomechanics and hospital personnel. On the other hand, simulation simplifications with the aim of reduction in computational time may cause in production of unrealistic outcomes. The main objective of the current investigation is to combine ideas such as autoregulation, impedance, fluid-solid interaction and idealized geometries in order to propose a computationally cheap methodology without excessive or unrealistic simplifications. The pressure boundary conditions are critical and polemic in numerical simulations of cardiovascular system, in which a specific arterial site is of interest and the rest of the netwrok is neglected but represented by a boundary condition. The proposed methodology is a pressure boundary condition which takes advantage of numerical simplicity of application of an imposed pressure boundary condition on outlets, while it includes more sophisticated concepts such as autoregulation and impedance to gain more realistic results. Incorporation of autoregulation and impedance converts the pressure boundary conditions to an active and dynamic boundary conditions, receiving feedback from the results during the numerical calculations and comparing them with the physiological requirements. On the other hand, the impedance boundary condition defines the shapes of the pressure history curves applied at outlets. The applications of the proposed method are seen on idealized geometry of the healthy arotic arch as well as idealized Stanford type A dissection, considering the interaction of the arterial walls with the pulsatile blood flow. The effect of surrounding tissues is incorporated and studied in the models. The simulations continue with FSI analysis of a patient-specific CT scanned geometry of an old individual. Finally, inspiring of the statistic results of mortality rates in Stanford type B dissection, three models of fenestrated dissection sac is studied and discussed. Applying the developed boundary condition, an alternative hypothesis is proposed by the author with respect to the decrease in mortality rates in patients with fenestrations.

Effect of the pulse trajectory on ultrasonic fluid velocity measurement

In a general situation a non-uniform velocity field gives rise to a shift of the otherwise straight acoustic pulse trajectory between the transmitter and receiver transducers of a sonic anemometer. The aim of this paper is to determine the effects of trajectory shifts on the velocity as measured by the sonic anemometer. This determination has been accomplished by developing a mathematical model of the measuring process carried out by sonic anemometers; a model which includes the non-straight trajectory effect. The problem is solved by small perturbation techniques, based on the relevant small parameter of the problem, the Mach number of the reference flow, M. As part of the solution, a general analytical expression for the deviations of the computed measured speed from the nominal speed has been obtained. The correction terms of both the transit time and of the measured speed are of M 2 order in rotational velocity field. The method has been applied to three simple, paradigmatic flows: one-directional horizontal and vertical shear flows, and mixed with a uniform horizontal flow.

Discrete linear stability and sensitivity analysis of fluid flows

Este trabajo presenta un método discreto para el cálculo de estabilidad hidrodinámica y análisis de sensibilidad a perturbaciones externas para ecuaciones diferenciales y en particular para las ecuaciones de Navier-Stokes compressible. Se utiliza una aproximación con variable compleja para obtener una precisión analítica en la evaluación de la matriz Jacobiana. Además, mapas de sensibilidad para la sensibilidad a las modificaciones del flujo de base y a una fuerza constante permiten identificar las regiones del campo fluido donde una modificacin (ej. fuerza puntual) tiene un efecto estabilizador del flujo. Se presentan cuatro casos de prueba: (1) un caso analítico para comprobar la derivación discreta, (2) una cavidad cerrada a bajo Reynolds para mostrar la mayor precisión en el cálculo de los valores propios con la aproximación de paso complejo, (3) flujo 2D en un cilindro circular para validar la metodología, y (4) flujo en un cavidad abierta, presentado para validar el método en casos de inestabilidades convectivamente inestables. Los tres últimos casos mencionados (2-4) se resolvieron con las ecuaciones de Navier-Stokes compresibles, utilizando un método Discontinuous Galerkin Spectral Element Method. Se obtuvo una buena concordancia para el caso de validación (3), cuando se comparó el nuevo método con resultados de la literatura. Además, este trabajo muestra que para el cálculo de los modos propios directos y adjuntos, así como para los mapas de sensibilidad, el uso de variables complejas es de suprema importancia para obtener una predicción precisa. El método descrito es aplicado al análisis para la estabilización de la estela generada por un disco actuador, que representa un modelo sencillo para hélices, rotores de helicópteros o turbinas eólicas. Se explora la primera bifurcación del flujo para un disco actuador, y se sugiere que está asociada a una inestabilidad de tipo Kelvin-Helmholtz, cuya estabilidad se controla con en el número de Reynolds y en la resistencia del disco actuador (o fuerza resistente). En primer lugar, se verifica que la disminución de la resistencia del disco tiene un efecto estabilizador parecido a una disminución del Reynolds. En segundo lugar, el análisis hidrodinmico discreto identifica dos regiones para la colocación de una fuerza puntual que controle las inestabilidades, una cerca del disco y otra en una zona aguas abajo. En tercer lugar, se muestra que la inclusión de un forzamiento localizado cerca del actuador produce una estabilización más eficiente que al forzar aguas abajo. El análisis de los campos de flujo controlados confirma que modificando el gradiente de velocidad cerca del actuador es más eficiente para estabilizar la estela. Estos resultados podrían proporcionar nuevas directrices para la estabilización de la estela de turbinas de viento o de marea cuando estén instaladas en un parque eólico y minimizar las interacciones no estacionarias entre turbinas. ABSTRACT A discrete framework for computing the global stability and sensitivity analysis to external perturbations for any set of partial differential equations is presented. In particular, a complex-step approximation is used to achieve near analytical accuracy for the evaluation of the Jacobian matrix. Sensitivity maps for the sensitivity to base flow modifications and to a steady force are computed to identify regions of the flow field where an input could have a stabilising effect. Four test cases are presented: (1) an analytical test case to prove the theory of the discrete framework, (2) a lid-driven cavity at low Reynolds case to show the improved accuracy in the calculation of the eigenvalues when using the complex-step approximation, (3) the 2D flow past a circular cylinder at just below the critical Reynolds number is used to validate the methodology, and finally, (4) the flow past an open cavity is presented to give an example of the discrete method applied to a convectively unstable case. The latter three (2–4) of the aforementioned cases were solved with the 2D compressible Navier–Stokes equations using a Discontinuous Galerkin Spectral Element Method. Good agreement was obtained for the validation test case, (3), with appropriate results in the literature. Furthermore, it is shown that for the calculation of the direct and adjoint eigenmodes and their sensitivity maps to external perturbations, the use of complex variables is paramount for obtaining an accurate prediction. An analysis for stabilising the wake past an actuator disc, which represents a simple model for propellers, helicopter rotors or wind turbines is also presented. We explore the first flow bifurcation for an actuator disc and it suggests that it is associated to a Kelvin- Helmholtz type instability whose stability relies on the Reynolds number and the flow resistance applied through the disc (or actuator forcing). First, we report that decreasing the disc resistance has a similar stabilising effect to an decrease in the Reynolds number. Second, a discrete sensitivity analysis identifies two regions for suitable placement of flow control forcing, one close to the disc and one far downstream where the instability originates. Third, we show that adding a localised forcing close to the actuator provides more stabilisation that forcing far downstream. The analysis of the controlled flow fields, confirms that modifying the velocity gradient close to the actuator is more efficient to stabilise the wake than controlling the sheared flow far downstream. An interesting application of these results is to provide guidelines for stabilising the wake of wind or tidal turbines when placed in an energy farm to minimise unsteady interactions.

Fluid-solid interaction in arteries incorporating the autoregulation concept in boundary conditions

In pre-surgery decisions in hospital emergency cases, fast and reliable results of the solid and fluid mechanics problems are of great interest to clinicians. In the current investigation, an iterative process based on a pressure-type boundary condition is proposed in order to reduce the computational costs of blood flow simulations in arteries, without losing control of the important clinical parameters. The incorporation of cardiovascular autoregulation, together with the well-known impedance boundary condition, forms the basis of the proposed methodology. With autoregulation, the instabilities associated with conventional pressure-type or impedance boundary conditions are avoided without an excessive increase in computational costs. The general behaviour of pulsatile blood flow in arteries, which is important from the clinical point of view, is well reproduced through this new methodology. In addition, the interaction between the blood and the arterial walls occurs via a modified weak coupling, which makes the simulation more stable and computationally efficient. Based on in vitro experiments, the hyperelastic behaviour of the wall is characterised and modelled. The applications and benefits of the proposed pressure-type boundary condition are shown in a model of an idealised aortic arch with and without an ascending aorta dissection, which is a common cardiovascular disorder.

Methodology for instantaneous average exhaust gas mass flow rate measurement

This paper presents a new methodology for measurement of the instantaneous average exhaust mass flow rate in reciprocating internal combustion engines to be used to determinate real driving emissions on light duty vehicles, as part of a Portable Emission Measurement System (PEMS). Firstly a flow meter, named MIVECO flow meter, was designed based on a Pitot tube adapted to exhaust gases which are characterized by moisture and particle content, rapid changes in flow rate and chemical composition, pulsating and reverse flow at very low engine speed. Then, an off-line methodology was developed to calculate the instantaneous average flow, considering the ?square root error? phenomenon. The paper includes the theoretical fundamentals, the developed flow meter specifications, the calibration tests, the description of the proposed off-line methodology and the results of the validation test carried out in a chassis dynamometer, where the validity of the mass flow meter and the methodology developed are demonstrated.

Vascular MADs: Two novel MAD-related genes selectively inducible by flow in human vascular endothelium

Vascular endothelium is an important transducer and integrator of both humoral and biomechanical stimuli within the cardiovascular system. Utilizing a differential display approach, we have identified two genes, Smad6 and Smad7, encoding members of the MAD-related family of molecules, selectively induced in cultured human vascular endothelial cells by steady laminar shear stress, a physiologic fluid mechanical stimulus. MAD-related proteins are a recently identified family of intracellular proteins that are thought to be essential components in the signaling pathways of the serine/threonine kinase receptors of the transforming growth factor β superfamily. Smad6 and Smad7 possess unique structural features (compared with previously described MADs), and they can physically interact with each other, and, in the case of Smad6, with other known human MAD species, in endothelial cells. Transient expression of Smad6 or Smad7 in vascular endothelial cells inhibits the activation of a transfected reporter gene in response to both TGF-β and fluid mechanical stimulation. Both Smad6 and Smad7 exhibit a selective pattern of expression in human vascular endothelium in vivo as detected by immunohistochemistry and in situ hybridization. Thus, Smad6 and Smad7 constitute a novel class of MAD-related proteins, termed vascular MADs, that are induced by fluid mechanical forces and can modulate gene expression in response to both humoral and biomechanical stimulation in vascular endothelium.

The leukocyte response to fluid stress

Leukocyte migration from a hemopoietic pool across marrow endothelium requires active pseudopod formation and adhesion. Leukocytes rarely show pseudopod formation while in circulation. At question then is the mechanism that serves to minimize leukocyte pseudopod formation in the circulation. We tested the hypothesis that fluid shear stress acts to prevent pseudopod formation. When individual human leukocytes (neutrophils, monocytes) spreading on glass surfaces in vitro were subjected to fluid shear stress (≈1 dyn/cm2), an instantaneous retraction of pseudopods was observed. Removal of the fluid shear stress in turn led to the return of pseudopod projection and cell spreading. When steady shear stress was prolonged over several minutes, leukocyte swelling occurs together with an enhanced random motion of cytoplasmic granules and a reduction of cytoplasmic stiffness. The response to shear stress could be suppressed by K+ channel blockers and chelation of external Ca2+. In rat mesentery microvessels after occlusion, circulating leukocytes project pseudopods in free suspension or when attached to the endothelium, even though immediately after occlusion only few pseudopods were present. When flow was restored, pseudopods on adhering leukocytes were retracted and then the cells began to roll and detach from the endothelium. In conclusion, plasma shear stress in the circulation serves to reduce pseudopod projection and adhesion of circulating leukocytes and vice versa reduction of shear stress leads to pseudopod projection and spreading of leukocytes on the endothelium.

Behavior of several two-dimensional fluid equations in singular scenarios

We give conditions that rule out formation of sharp fronts for certain two-dimensional incompressible flows. We show that a necessary condition of having a sharp front is that the flow has to have uncontrolled velocity growth. In the case of the quasi-geostrophic equation and two-dimensional Euler equation, we obtain estimates on the formation of semi-uniform fronts.

Fluid shear stress inhibits TNF-α activation of JNK but not ERK1/2 or p38 in human umbilical vein endothelial cells: Inhibitory crosstalk among MAPK family members

Atherosclerosis preferentially occurs in areas of turbulent flow and low fluid shear stress, whereas laminar flow and high shear stress are atheroprotective. Inflammatory cytokines, such as tumor necrosis factor-α (TNF-α) and IL-1 stimulate expression of endothelial cell (EC) genes that may promote atherosclerosis. TNF-α and IL-1 regulate gene expression in ECs, in part, by stimulating mitogen-activated protein kinases (MAPK), which phosphorylate transcription factors. We hypothesized that steady laminar flow inhibits cytokine-mediated activation of MAPK in EC. To test this hypothesis, we determined the effects of flow (shear stress = 12 dynes/cm2) on TNF-α and IL-1-stimulated activity of three MAPK in human umbilical vein ECs (HUVEC): extracellular signal-regulated kinase (ERK1/2), p38, and c-Jun N-terminal kinase (JNK). Flow alone stimulated ERK1/2 and p38 activity but decreased JNK activity compared with static controls. TNF-α or IL-1 alone activated ERK1/2, p38, and JNK maximally at 15 min in HUVEC. Preexposing HUVEC for 10 min to flow inhibited TNF-α and IL-1 activation of JNK by 46% and 49%, respectively, but had no significant effect on ERK1/2 or p38 activation. Incubation of HUVEC with PD98059, which inhibits flow-mediated ERK1/2 activation, prevented flow from inhibiting cytokine activation of JNK. Phorbol 12-myristate 13-acetate, which strongly activates ERK1/2, also inhibited TNF-α activation of JNK. These findings indicate that fluid shear stress inhibits TNF-α-mediated signaling events in HUVEC via the activation of the ERK1/2 signaling pathway. Inhibition of TNF-α signal transduction represents a mechanism by which steady laminar flow may exert atheroprotective effects on the endothelium.

Nitric oxide mediates the increase in local cerebral blood flow during focal seizures.

The role of nitric oxide (NO) in the increase in local cerebral blood flow (LCBF) elicited by focal cortical epileptic seizures was investigated in anesthetized adult rats. Seizures were induced by topical bicuculline methiodide applied through two cranial windows drilled over homotopic sites of the frontal cortex, and LCBF was measured by quantitative autoradiography by using 4-iodo[N-methyl-14C]antipyrine. Superfusion of an inhibitor of NO synthase, N omega-nitro-L-arginine (NA; 1 mM), for 45 min abolished the increase of LCBF induced by topical bicuculline methiodide (10 mM) [164 +/- 18 ml/100 g per min in the artificial cerebrospinal fluid (aCSF)-superfused side and 104 +/- 12 ml/100 g per ml in the NA-superfused side; P < 0.005]. This effect was reversed by coapplication of an excess of L-arginine substrate (10 mM) (218 +/- 22 ml/100 g per min in the aCSF-superfused side and 183 +/- 31 ml/100 g per min in the NA + L-Arg-superfused side) but not by 10 mM D-arginine, a stereoisomer with poor affinity for NO synthase (193 +/- 17 ml/100 g per min in the aCSF-superfused side and 139 +/- 21 ml/100 g per min in the NA + D-Arg-superfused side; P < 0.005). Superfusion of the guanylyl cyclase inhibitor methylene blue attenuated the LCBF increase elicited by topical bicuculline methiodide by 25% +/- 16% (P < 0.05). The present findings suggest that NO is the mediator of the vasodilation in response to focal epileptic seizures.

Investigating the Role of cAMP-Signaling in the Regulation of Vascular Endothelial Cell Adaptive Responses to Fluid Shear Stress

The circulating blood exerts a force on the vascular endothelium, termed fluid shear stress (FSS), which directly impacts numerous vascular endothelial cell (VEC) functions. For example, high rates of linear and undisturbed (i.e. laminar) blood flow maintains a protective and quiescent VEC phenotype. Meanwhile, deviations in blood flow, which can occur at vascular branchpoints and large curvatures, create areas of low, and/or oscillatory FSS, and promote a pro-inflammatory, pro-thrombotic and hyperpermeable phenotype. Indeed, it is known that these areas are prone to the development of atherosclerotic lesions. Herein, we show that cyclic nucleotide phosphodiesterase (PDE) 4D (PDE4D) activity is increased by FSS in human arterial endothelial cells (HAECs) and that this activation regulates the activity of cAMP-effector protein, Exchange Protein-activated by cAMP-1 (EPAC1), in these cells. Importantly, we also show that these events directly and critically impact HAEC responses to FSS, especially when FSS levels are low. Both morphological events induced by FSS, as measured by changes in cell alignment and elongation in the direction of FSS, and the expression of critical FSS-regulated genes, including Krüppel-like factor 2 (KLF2), endothelial nitric oxide synthase (eNOS) and thrombomodlin (TM), are mediated by EPAC1/PDE4D signaling. At a mechanistic level, we show that EPAC1/PDE4D acts through the vascular endothelial-cadherin (VECAD)/ platelet-cell adhesion molecule-1 (PECAM1)/vascular endothelial growth factor receptor 2 (VEGFR2) mechanosensor to activate downstream signaling though Akt. Given the critical role of PDE4D in mediating these effects, we also investigated the impact of various patterns of FSS on the expression of individual PDE genes in HAECs. Notably, PDE2A was significantly upregulated in response to high, laminar FSS, while PDE3A was upregulated under low, oscillatory FSS conditions only. These data may provide novel therapeutic targets to limit FSS-dependent endothelial cell dysfunction (ECD) and atherosclerotic development.