Noninvasive 4D pressure difference mapping derived from 4D flow MRI in patients with repaired aortic coarctation: comparison with young healthy volunteers

To assess spatial and temporal pressure characteristics in patients with repaired aortic coarctation compared to young healthy volunteers using time-resolved velocity-encoded three-dimensional phase-contrast magnetic resonance imaging (4D flow MRI) and derived 4D pressure difference maps. After in vitro validation against invasive catheterization as gold standard, 4D flow MRI of the thoracic aorta was performed at 1.5T in 13 consecutive patients after aortic coarctation repair without recoarctation and 13 healthy volunteers. Using in-house developed processing software, 4D pressure difference maps were computed based on the Navier-Stokes equation. Pressure difference amplitudes, maximum slope of pressure amplitudes and spatial pressure range at mid systole were retrospectively measured by three readers, and twice by one reader to assess inter- and intraobserver agreement. In vitro, pressure differences derived from 4D flow MRI showed excellent agreement to invasive catheter measurements. In vivo, pressure difference amplitudes, maximum slope of pressure difference amplitudes and spatial pressure range at mid systole were significantly increased in patients compared to volunteers in the aortic arch, the proximal descending and the distal descending thoracic aorta (p < 0.05). Greatest differences occurred in the proximal descending aorta with values of the three parameters for patients versus volunteers being 19.7 ± 7.5 versus 10.0 ± 2.0 (p < 0.001), 10.9 ± 10.4 versus 1.9 ± 0.4 (p = 0.002), and 8.7 ± 6.3 versus 1.6 ± 0.9 (p < 0.001). Inter- and intraobserver agreements were excellent (p < 0.001). Noninvasive 4D pressure difference mapping derived from 4D flow MRI enables detection of altered intraluminal aortic pressures and showed significant spatial and temporal changes in patients with repaired aortic coarctation.

Aerodynamic Optimization of High-Speed Trains Nose using a Genetic Algorithm and Artificial Neural Network

An aerodynamic optimization of the train aerodynamic characteristics in term of front wind action sensitivity is carried out in this paper. In particular, a genetic algorithm (GA) is used to perform a shape optimization study of a high-speed train nose. The nose is parametrically defined via Bézier Curves, including a wider range of geometries in the design space as possible optimal solutions. Using a GA, the main disadvantage to deal with is the large number of evaluations need before finding such optimal. Here it is proposed the use of metamodels to replace Navier-Stokes solver. Among all the posibilities, Rsponse Surface Models and Artificial Neural Networks (ANN) are considered. Best results of prediction and generalization are obtained with ANN and those are applied in GA code. The paper shows the feasibility of using GA in combination with ANN for this problem, and solutions achieved are included.

A non-perturbative renormalization group study of the stochastic NavierStokes equation

We study the renormalization group flow of the average action of the stochastic Navier-Stokes equation with power-law forcing. Using Galilean invariance, we introduce a nonperturbative approximation adapted to the zero-frequency sector of the theory in the parametric range of the Hölder exponent 4−2 ɛ of the forcing where real-space local interactions are relevant. In any spatial dimension d, we observe the convergence of the resulting renormalization group flow to a unique fixed point which yields a kinetic energy spectrum scaling in agreement with canonical dimension analysis. Kolmogorov's −5/3 law is, thus, recovered for ɛ=2 as also predicted by perturbative renormalization. At variance with the perturbative prediction, the −5/3 law emerges in the presence of a saturation in the ɛ dependence of the scaling dimension of the eddy diffusivity at ɛ=3/2 when, according to perturbative renormalization, the velocity field becomes infrared relevant.

Estudio teórico experimental de la estabilidad lateral en vehículos cisterna. Metodología para la determinación del umbral de vuelco.

El oleaje de la masa de líquido contenida en el interior del tanque de los camiones cisterna puede afectar significativamente a su dinámica de marcha y su estabilidad. Plateando como objetivo, el aumento de la seguridad, cualquier mejora en el diseño del vehículo o en la geometría de la virola, redundará en un aumento de la misma. En la actualidad, el marco normativo exige a los vehículos un umbral de vuelco mínimo, verificable mediante ensayo o cálculo en condiciones cuasi-estáticas. La experiencia en su aplicación ha demostrado que, éste último resulta extremadamente conservador, limitando las posibilidades de diseño del fabricante; además, el grado de equivalencia con las condiciones dinámicas a las que se puede ver sometido el vehículo durante su vida útil, no es claro. En esta tesis, se proponen una serie de mejoras en el método de cálculo simplificado, que aumentan la precisión en la estimación del umbral de vuelco, junto con el empleo de un coeficiente de ponderación del ángulo requerido, para tener en cuenta el efecto dinámico del oleaje. Además, se ha desarrollado un método de cálculo alternativo, que disminuye el error de la predicción, aunque requiere conocer un mayor número de parámetros del vehículo. Para cuantificar el efecto de la carga líquida, se ha realizado un análisis experimental, con tanques a escala, de las geometrías más frecuentes, junto con la propuesta de una metodología de simulación de la interacción fluido-vehículo, basada en el acoplamiento entre un modelo obtenido mediante una analogía mecánica que reproduce el efecto del oleaje, con un modelo del vehículo, empleando la técnica de multi-cuerpo, que ha sido validada experimentalmente. Para la obtención de los parámetros que la definen, se propone dos alternativas, la experimentación y estimación parámetros o el empleo la dinámica computacional de fluidos, basada en las ecuaciones de Navier-Stokes, incorporando la técnica de volumen de fluido para capturar la superficie libre. Se ofrecen resultados de ambas y su comparación. Finalmente, se han analizado las diferentes tipologías de maniobras existentes para la evaluación de la estabilidad lateral de los vehículos con ensayos dinámicos, calificando y cuantificando su equivalencia.

Aerodynamic optimization of the ICE2 high-speed train nose using a genetic algorithm and metamodels

An aerodynamic optimization of the ICE 2 high-speed train nose in term of front wind action sensitivity is carried out in this paper. The nose is parametrically defined by Be?zier Curves, and a three-dimensional representation of the nose is obtained using thirty one design variables. This implies a more complete parametrization, allowing the representation of a real model. In order to perform this study a genetic algorithm (GA) is used. Using a GA involves a large number of evaluations before finding such optimal. Hence it is proposed the use of metamodels or surrogate models to replace Navier-Stokes solver and speed up the optimization process. Adaptive sampling is considered to optimize surrogate model fitting and minimize computational cost when dealing with a very large number of design parameters. The paper introduces the feasi- bility of using GA in combination with metamodels for real high-speed train geometry optimization.

Comparing linear and non linear wind flow models

Assessing wind conditions on complex terrain has become a hard task as terrain complexity increases. That is why there is a need to extrapolate in a reliable manner some wind parameters that determine wind farms viability such as annual average wind speed at all hub heights as well as turbulence intensities. The development of these tasks began in the early 90´s with the widely used linear model WAsP and WAsP Engineering especially designed for simple terrain with remarkable results on them but not so good on complex orographies. Simultaneously non-linearized Navier Stokes solvers have been rapidly developed in the last decade through CFD (Computational Fluid Dynamics) codes allowing simulating atmospheric boundary layer flows over steep complex terrain more accurately reducing uncertainties. This paper describes the features of these models by validating them through meteorological masts installed in a highly complex terrain. The study compares the results of the mentioned models in terms of wind speed and turbulence intensity.

RANS simulations of wind flow at the Bolund experiment

As part of their development, the predictions of numerical wind flow models must be compared with measurements in order to estimate the uncertainty related to their use. Of course, the most rigorous such comparison is under blind conditions. The following paper includes a detailed description of three different wind flow models, all based on a Reynolds-averaged Navier-Stokes approach and two-equation k-ε closure, that were tested as part of the Bolund blind comparison (itself based on the Bolund experiment which measured the wind around a small coastal island). The models are evaluated in terms of predicted normalized wind speed and turbulent kinetic energy at 2 m and 5 m above ground level for a westerly wind direction. Results show that all models predict the mean velocity reasonably well; however accurate prediction of the turbulent kinetic energy remains achallenge.

Simulación de componentes no estacionarias de viento en orografías complejas

Realistic operation of helicopter flight simulators in complex topographies (such as urban environments) requires appropriate prediction of the incoming wind, and this prediction should be made in real time. Unfortunately, the wind topology around complex topographies shows time-dependent, fully nonlinear, turbulent patterns (i.e., wakes) whose simulation cannot be made using computationally inexpensive tools based on corrected potential approximations. Instead, the full Navier-Stokes plus some kind of turbulent modeling is necessary, which is quite computationally expensive. The complete unsteady flow depends on two parameters, namely the velocity and orientation of the free stream flow. The aim of this MSc thesis is to develop a methodology for the real time simulation of these complex flows. For simplicity, the flow around a single building (20 mx20 m cross section and 100 m height) is considered, with free stream velocity in the range 5-25 m/s. Because of the square cross section, the problem shows two reflection symmetries, which allows for restricting the orientations to the range 0° < a. < 45°. The methodology includes an offline preprocess and the online operation. The preprocess consists in three steps: An appropriate, unstructured mesh is selected in which the flow is sim¬ulated using OpenFOAM, and this is done for 33 combinations of 3 free stream intensities and 11 orientations. For each of these, the simulation proceeds for a sufficiently large time as to eliminate transients. This step is quite computationally expensive. Each flow field is post-processed using a combination of proper orthogonal decomposition, fast Fourier transform, and a convenient optimization tool, which identifies the relevant frequencies (namely, both the basic frequencies and their harmonics) and modes in the computational mesh. This combination includes several new ingredients to filter errors out and identify the relevant spatio-temporal patterns. Note that, in principle, the basic frequencies depend on both the intensity and the orientation of the free stream flow. The outcome of this step is a set of modes (vectors containing the three velocity components at all mesh points) for the various Fourier components, intensities, and orientations, which can be organized as a third order tensor. This step is fairly computationally inexpensive. The above mentioned tensor is treated using a combination of truncated high order singular value, decomposition and appropriate one-dimensional interpolation (as in Lorente, Velazquez, Vega, J. Aircraft, 45 (2008) 1779-1788). The outcome is a tensor representation of both the relevant fre¬quencies and the associated Fourier modes for a given pair of values of the free stream flow intensity and orientation. This step is fairly compu¬tationally inexpensive. The online, operation requires just reconstructing the time-dependent flow field from its Fourier representation, which is extremely computationally inex¬pensive. The whole method is quite robust.

Computation of nonlinear streaky structures in boundary layer flow

In this work, the Reduced Navier Stokes (RNS) are numerically integrated, and used to calculate nonlinear finite amplitude streaks. These structures are interesting since they can have a stabilizing effect and delay the transition to the turbulent regime. RNS formulation is also used to compute the family of nonlinear intrinsic streaks that emerge from the leading edge in absence of any external perturbation. Finally, this formulation is generalized to include the possibility of having a curved bottom wall

Physics of Vibrating Airfoils at Low Reduced Frequency

The unsteady aerodynamics of low pressure turbine vibrating airfoils in flap mode is studied in detail using a frequency domain linearized Navier-Stokes solver. Both the travelling-wave and influence coefficient formulations of the problem are used to highlight key aspects of the physics and understand different trends such as the effect of reduced frequency and Mach number. The study is focused in the low-reduced frequency regime which is of paramount relevance for the design of aeronautical low-pressure turbines and compressors. It is concluded that the effect of the Mach number on the unsteady pressure phase can be neglected in first approximation and that the unsteadiness of the vibrating and adjacent airfoils is driven by vortex shedding mechanisms. Finally a simple model to estimate the work-per-cycle as a function of the reduced frequency and Mach Number is provided. The edge-wise and torsion modes are presented in less detail but it is shown that acoustic waves are essential to explain its behaviour. The non-dimensional work-per-cycle of the edge-wise mode shows a large dependence with the Mach number while in the torsion mode a large number of airfoils is needed to reconstruct the work-per-cycle departing from the influence coefficients.

Estudio matemático del comportamiento de un fluido incompresible en flujo laminar

El objetivo de este trabajo es realizar un estudio del comportamiento de un fluido en flujo laminar. En primer lugar se realizará una breve introducción Histórica, con el fin de situar al lector en la época donde suceden los acontecimientos más importantes relativos a la mecánica de fluidos. Se presentarán a las figuras más importantes de la mecánica de fluidos moderna, como Sir Gabriel Stokes o Osborne Reynolds, así como los acontecimientos más importantes. A continuación se expondran las ecuaciones que rigen el comportamiento de un fluido, las ecuaciones de Navier-Stokes y la ecuación de Reynolds, con el objetivo de ayudar al lector a entender los análisis posteriores que se realizarán sobre dichas ecuaciones. En tercer lugar se analizará la ecuación de Reynolds, la existencia y unicidad de soluciones, para seguidamente, realizar una simulación del problema. Dicha simulación se ha realizado en un script, bajo la herramienta Matlab, se explicará como se ha realizado la simulación y se expondrán ejemplos de los casos de un fluido compresible y uno incompresible para diversas ecuaciones. Seguidamente se describe el fenómeno de la capa límite y se explica con un ejemplo en un fluido compresible. Para finalizar se exponen las conclusiones obtenidas en la realización de este trabajo.

The efficiency of propulsion by a rotating flagellum

[At very low Reynolds number, the regime in which fluid dynamics is governed by Stokes equations, a helix that translates along its axis under an external force but without an external torque will necessarily rotate. By the linearity of the Stokes equations, the same helix that is caused to rotate due to an external torque will necessarily translate. This is the physics that underlies the mechanism of flagellar propulsion employed by many microorganisms. Here, I examine the linear relationships between forces and torques and translational and angular velocities of helical objects to understand the nature of flagellar propulsion.]

Active control of waves in a cochlear model with subpartitions.

Multiscale asymptotic methods developed previously to study macromechanical wave propagation in cochlear models are generalized here to include active control of a cochlear partition having three subpartitions, the basilar membrane, the reticular lamina, and the tectorial membrane. Activation of outer hair cells by stereocilia displacement and/or by lateral wall stretching result in a frequency-dependent force acting between the reticular lamina and basilar membrane. Wavelength-dependent fluid loads are estimated by using the unsteady Stokes' equations, except in the narrow gap between the tectorial membrane and reticular lamina, where lubrication theory is appropriate. The local wavenumber and subpartition amplitude ratios are determined from the zeroth order equations of motion. A solvability relation for the first order equations of motion determines the subpartition amplitudes. The main findings are as follows: The reticular lamina and tectorial membrane move in unison with essentially no squeezing of the gap; an active force level consistent with measurements on isolated outer hair cells can provide a 35-dB amplification and sharpening of subpartition waveforms by delaying dissipation and allowing a greater structural resonance to occur before the wave is cut off; however, previously postulated activity mechanisms for single partition models cannot achieve sharp enough tuning in subpartitioned models.

Numerical simulation of magnetite segregation in a dense medium cyclone

Numerical simulations of turbulent driven flow in a dense medium cyclone with magnetite medium have been conducted using Fluent. The predicted air core shape and diameter were found to be close to the experimental results measured by gamma ray tomography. It is possible that the Large eddy simulation (LES) turbulence model with Mixture multi-phase model can be used to predict the air/slurry interface accurately although the LES may need a finer grid. Multi-phase simulations (air/water/medium) are showing appropriate medium segregation effects but are over-predicting the level of segregation compared to that measured by gamma-ray tomography in particular with over prediction of medium concentrations near the wall. Further, investigated the accurate prediction of axial segregation of magnetite using the LES turbulence model together with the multi-phase mixture model and viscosity corrections according to the feed particle loading factor. Addition of lift forces and viscosity correction improved the predictions especially near the wall. Predicted density profiles are very close to gamma ray tomography data showing a clear density drop near the wall. The effect of size distribution of the magnetite has been fully studied. It is interesting to note that the ultra-fine magnetite sizes (i.e. 2 and 7 mu m) are distributed uniformly throughout the cyclone. As the size of magnetite increases, more segregation of magnetite occurs close to the wall. The cut-density (d(50)) of the magnetite segregation is 32 gm, which is expected with superfine magnetite feed size distribution. At higher feed densities the agreement between the [Dungilson, 1999; Wood, J.C., 1990. A performance model for coal-washing dense medium cyclones, Ph.D. Thesis, JKMRC, University of Queensland] correlations and the CFD are reasonably good, but the overflow density is lower than the model predictions. It is believed that the excessive underflow volumetric flow rates are responsible for under prediction of the overflow density. (c) 2006 Elsevier Ltd. All rights reserved.

Study of flow behaviour in a three-product cyclone using computational fluid dynamics

Simplicity in design and minimal floor space requirements render the hydrocyclone the preferred classifier in mineral processing plants. Empirical models have been developed for design and process optimisation but due to the complexity of the flow behaviour in the hydrocyclone these do not provide information on the internal separation mechanisms. To study the interaction of design variables, the flow behaviour needs to be considered, especially when modelling the new three-product cyclone. Computational fluid dynamics (CFD) was used to model the three-product cyclone, in particular the influence of the dual vortex finder arrangement on flow behaviour. From experimental work performed on the UG2 platinum ore, significant differences in the classification performance of the three-product cyclone were noticed with variations in the inner vortex finder length. Because of this simulations were performed for a range of inner vortex finder lengths. Simulations were also conducted on a conventional hydrocyclone of the same size to enable a direct comparison of the flow behaviour between the two cyclone designs. Significantly, high velocities were observed for the three-product cyclone with an inner vortex finder extended deep into the conical section of the cyclone. CFD studies revealed that in the three-product cyclone, a cylindrical shaped air-core is observed similar to conventional hydrocyclones. A constant diameter air-core was observed throughout the inner vortex finder length, while no air-core was present in the annulus. (c) 2006 Elsevier Ltd. All rights reserved.