Vortex rings interacting with a wall

Accurate numerical simulations of vortex rings impinging on flat boundaries have revealed the same features observed in the experiment of Walkeret al. (1987). They observed atRe _v >1250 the formation of azimuthal instabilities in the secondary ring during its compression within the primary ring. In the present numerical simulation the number of waves agrees very well with those observed in the experiment. The distributions of the vorticities together with the distribution of each term in the vorticity equations give insights on the formation of azimuthal instabilities.     

Three-dimensional vortex visualization in stratified spin-up

Abstract  

We present the results of three-dimensional time-dependent numerical simulations of incremental spin-up of a thermally stratified fluid. The fluid inside a vertical cylindrical container of radius R and height 2H is water characterized by the kinematic viscosity ν and thermal diffusivity κ. Initially, its density (temperature) varies linearly with height and is characterized by a constant buoyancy frequency N, which is proportional to the density gradient. The system undergoes an abrupt change in the rotation rate from its initial value Ω_i, when the fluid is in a solid-body rotation state, to the final value Ω_f. The aim of this contribution is to show the formation of columnar vortices in a high Rossby number spin-up flow.     

Numerical simulation of the non-Newtonian blood flow through a mechanical aortic valve

This work focuses on the comparison between Newtonian and non-Newtonian blood flows through a bileaflet mechanical heart valve in the aortic root. The blood, in fact, is a concentrated suspension of cells, mainly red blood cells, in a Newtonian matrix, the plasma, and consequently its overall behavior is that of a non-Newtonian fluid owing to the action of the cells’ membrane on the fluid part. The common practice, however, assumes the blood in large vessels as a Newtonian fluid since the shear rate is generally high and the effective viscosity becomes independent of the former. In this paper, we show that this is not always the case even in the aorta, the largest artery of the systemic circulation, owing to the pulsatile and transitional nature of the flow. Unexpectedly, for most of the pulsating cycle and in a large part of the fluid volume, the shear rate is smaller than the threshold level for the blood to display a constant effective viscosity and its shear thinning character might affect the system dynamics. A direct inspection of the various flow features has shown that the valve dynamics, the transvalvular pressure drop and the large-scale features of the flow are very similar for the Newtonian and non-Newtonian fluid models. On the other hand, the mechanical damage of the red blood cells (hemolysis), induced by the altered stress values in the flow, is larger for the non-Newtonian fluid model than for the Newtonian one.     

Turbulent Flow and Dispersion of Inertial Particles in a Confined Jet Issued by a Long Cylindrical Pipe

In this work we examine first the flow field of a confined jet produced by a turbulent flow in a long cylindrical pipe issuing in an abrupt angle diffuser. Second, we examine the dispersion of inertial micro-particles entrained by the turbulent flow. Specifically, we examine how the particle dispersion field evolves in the multiscale flow generated by the interactions between the large-scale structures, which are geometry dependent, with the smaller turbulent scales issued by the pipe which are advected downstream. We use Large-Eddy-Simulation (LES) for the flow field and Lagrangian tracking for particle dispersion. The complex shape of the domain is modelled using the immersed-boundaries method. Fully developed turbulence inlet conditions are derived from an independent LES of a spatially periodic cylindrical pipe flow. The flow field is analyzed in terms of local velocity signals to determine spatial coherence and decay rate of the coherent K–H vortices and to make quantitative comparisons with experimental data on free jets. Particle dispersion is analyzed in terms of statistical quantities and also with reference to the dynamics of the coherent structures. Results show that the particle dynamics is initially dominated by the Kelvin–Helmholtz (K–H) rolls which form at the expansion and only eventually by the advected smaller turbulence scales.     

Normal and oblique collisions of a vortex ring with a wall

The oblique collision of a vortex ring with a solid wall, atRe=/=1389, has been analysed by the direct simulation of the Navier-Stokes equations in Cartesian coordinates. In accordance with a previous experimental study [1], the secondary vorticity produced at the wall is organized into a loop-like vortex in the region of the ring furthest away from the wall. As the ring approaches the wall, the region closest is subjected to a high rate of stretching which increases the vorticity in the core. The vorticity gradients along the core generate bi-helical vortex lines continually displaced towards the region of the ring furthest away from the wall. The analysis of the vorticity and straining fields revealed that the pressure gradient along the core is responsible for the convective motion that displaces these vortex lines and accumulates secondary vorticity in the region far from the wall. This vorticity rolls up and forms a secondary structure which by self-induction moves away from the wall.The fundamental role of the differential stretching has been demonstrated by comparing the case of oblique collision with that of normal collision and with the collision of a two-dimensional vortex pair with an oblique wall.

---

Sommario L'interazione di un vortice ad anello con una parete obliqua, aRe=1389, è stata analizzata mediante la simulazione diretta delle equazioni di Navier-Stokes in coordinate cartesiane. In accordo con un precedente esperimento [1] è stato evidenziato che la vorticità secondaria, prodotta alla parete, si organizza in una strutura vorticosa a loop nella regione dell'anello più lontana dalla parete. Quando il vortice si avvicina alla parete, la parte più vicina è soggetta ad un'elevata deformazione che aumenta il valore della vorticità nel core. La distribuzione non uniforme di vorticità lungo il core del vortice genera delle linee di vorticità elicoidali che vengono transportate verso la regione dell'anello più lontana dalla parete. L'analisi dei campi di vorticità e di deformazione ha rivelato che il gradiente di pressione, dovuto al campo di deformazione non uniforme lungo il core del vortice, è responsabile di un moto convettivo che trasporta le linee di vorticità ed accumula la vorticità secondaria nella regione del vortice più lontana dalla parete, dove la struttura secondaria viene generata.Il ruolo fondamentale della deformazione non uniforme è stato evidenziato mediante il confronto della collisione obliqua coni casi di collisione normale e di collisione di una coppia di vortici bidimensionali con una parete obliqua.

     

A Non-Adiabatic Flamelet Progress�CVariable Approach for LES of Turbulent Premixed Flames

A progress variable/flame surface density/probability density function method has been employed for a Large Eddy Simulation of a CH₄/Air turbulent premixed bluff body flame. In particular, both mean and variance of the progress variable are transported and subgrid spatially filtered gradient contributes to model the flame surface density (that introduces the effect of the subgrid flame reaction zone) and to presume a probability density function (that introduces the effect of subgrid fluctuations on chemistry). Chemistry is preliminarly tabulated in terms of laminar premixed flames and enthalpy is included as a new coordinate in their tabulation to take into account heat losses in the flowfield. Then, the PDF is used to build a turbulent flamelet library. The filtered mass, momentum, enthalpy and scalar equations mentioned above are integrated by an explicit scheme using finite differences, 2nd?Corder accurate in space and third order in time, over a cylindrical non-uniform grid using a staggered mesh. The bluff-body geometry is modelled by using the Immersed Boundary Method. The numerical predictions are compared with the available experimental data.     

Identification of zones in a free evolving vortex ring

A direct numerical simulation of a three-dimensional free evolving vortex ring has been performed to understand some features of the dynamics of vortex rings. In previous experiments, a persistent wake has been observed and the direct simulation has clarified the causes of this wake formation. A kinematic alignment analysis of the vorticity vector with the eigenvectors of the strain-rate tensor, together with vorticity visualisations, have been employed to explain the growth of the Widnall instabilities. The analysis has shown the presence of two different zones, the first characterised by intense vorticity and low strain rate and the second with opposite features. In the former, there is a ring core which maintains its coherence even in the long term, in the latter there are elongated structures contributing to the deformation of the ring.