A Kármán–Howarth–Monin equation for variable-density turbulence

2018 ◽  
Vol 843 ◽  
pp. 382-418 ◽  
Author(s):  
Chris C. K. Lai ◽  
John J. Charonko ◽  
Katherine Prestridge
Keyword(s):  
Energy Transfer ◽  
Initial Density ◽  
Density Ratio ◽  
Variable Density ◽  
Navier Stokes ◽  
Mean Flow ◽  
Development Region ◽  
The Mean ◽  
Variable Density Turbulence ◽  
Transfer Term

We present a generalisation of the Kármán–Howarth–Monin (K–H–M) equation to include variable-density (VD) effects. The derived equation (i) reduces to the original K–H–M equation when density is a constant and (ii) leads to a VD analogue of the $4/5$-law with the same value of constant ($=4/5$) appearing as the prefactor of the dissipation rate. The equation is employed to understand negative turbulent kinetic energy production in a $\text{SF}_{6}$ turbulent round jet with an initial density ratio of 4.2. From a Reynolds-averaged Navier–Stokes (RANS) perspective, negative production means that the mean flow is strengthened at the expense of the energy of turbulent fluctuations. We show that the associated energy transfer is accomplished by the deformation of smaller turbulent eddies into large ones in the development region of the jet and is captured by the linear scale-by-scale energy transfer term in the VD K–H–M equation. The nonlinear transfer term of the VD K–H–M equation depicts a conventional forward cascade for all eddies having a size less than the Eulerian integral length scale, regardless of their orientation. The net effect is a retarded energy cascade in the non-Boussinesq jet that has not been accounted for by existing turbulence theories. Implications of this observation for turbulence modelling are discussed.

Measurement of the Mean Flow Field in a Smooth Rotating Channel With Coriolis and Buoyancy Effects

2017 ◽  
Author(s):  
Ruquan You ◽  
Haiwang Li ◽  
Zhi Tao ◽  
Kuan Wei
Keyword(s):  
Flow Field ◽  
Coriolis Force ◽  
Indium Tin Oxide ◽  
Buoyancy Force ◽  
Flow Direction ◽  
Density Ratio ◽  
Mean Flow ◽  
The Mean ◽  
Static Conditions ◽  
Rotating Channel

The mean flow field in a smooth rotating channel was measured by particle image velocimetry under the effect of buoyancy force. In the experiments, the Reynolds number, based on the channel hydraulic diameter (D) and the bulk mean velocity (Um), is 10000, and the rotation numbers are 0, 0.13, 0.26, 0.39, 0.52, respectively. The four channel walls are heated with Indium Tin Oxide (ITO) heater glass, making the density ratio (d.r.) about 0.1 and the maximum value of buoyancy number up to 0.27. The mean flow field was simulated on a 3D reconstruction at the position of 3.5<X/D<6.5, where X is along the mean flow direction. The effect of Coriolis force and buoyancy force on the mean flow was taken into consideration in the current work. The results show that the Coriolis force pushes the mean flow to the trailing side, making the asymmetry of the mean flow with that in the static conditions. On the leading surface, due to the effect of buoyancy force, the mean flow field changes considerably. Comparing with the case without buoyancy force, separated flow was captured by PIV on the leading side in the case with buoyancy force. More details of the flow field will be presented in this work.

scholarly journals Self-consistent triple decomposition of the turbulent flow over a backward-facing step under finite amplitude harmonic forcing

2019 ◽  
Vol 475 (2225) ◽  
pp. 20190018
Author(s):  
E. Yim ◽  
P. Meliga ◽  
F. Gallaire
Keyword(s):  
Turbulent Flow ◽  
Finite Amplitude ◽  
Navier Stokes ◽  
Reynolds Stresses ◽  
Mean Flow ◽  
Backward Facing Step ◽  
Eddy Viscosity Model ◽  
Coherent Interaction ◽  
Self Consistent ◽  
The Mean

We investigate the saturation of harmonically forced disturbances in the turbulent flow over a backward-facing step subjected to a finite amplitude forcing. The analysis relies on a triple decomposition of the unsteady flow into mean, coherent and incoherent components. The coherent–incoherent interaction is lumped into a Reynolds averaged Navier–Stokes (RANS) eddy viscosity model, and the mean–coherent interaction is analysed via a semi-linear resolvent analysis building on the laminar approach by Mantič-Lugo & Gallaire (2016 J. Fluid Mech. 793 , 777–797. ( doi:10.1017/jfm.2016.109 )). This provides a self-consistent modelling of the interaction between all three components, in the sense that the coherent perturbation structures selected by the resolvent analysis are those whose Reynolds stresses force the mean flow in such a way that the mean flow generates exactly the aforementioned perturbations, while also accounting for the effect of the incoherent scale. The model does not require any input from numerical or experimental data, and accurately predicts the saturation of the forced coherent disturbances, as established from comparison to time-averages of unsteady RANS simulation data.

Improved κ-ɛ Modelling of Impeller Flow Performance of a Mixed-Flow Pump under off-Design Operating States

1993 ◽  
Vol 207 (4) ◽  
pp. 219-229 ◽  
Author(s):  
S M Fraser ◽  
Y Zhang
Keyword(s):  
Stokes Equations ◽  
Three Dimensional ◽  
Careful Analysis ◽  
Navier Stokes ◽  
Mean Flow ◽  
Mixed Flow ◽  
Navier Stokes Equations ◽  
Flow Through ◽  
The Mean ◽  
Flow Pump

Three-dimensional turbulent flow through the impeller passage of a model mixed-flow pump has been simulated by solving the Navier-Stokes equations with an improved κ-ɛ model. The standard κ-ɛ model was found to be unsatisfactory for solving the off-design impeller flow and a converged solution could not be obtained at 49 per cent design flowrate. After careful analysis, it was decided to modify the standard κ-ɛ model by including the extra rates of strain due to the acceleration of impeller rotation and geometrical curvature and removing the mathematical ill-posedness between the mean flow turbulence modelling and the logarithmic wall function.

On the interactions between large-scale structure and fine-grained turbulence in a free shear flow I. The development of temporal interactions in the mean

1976 ◽  
Vol 352 (1669) ◽  
pp. 213-247 ◽  
Keyword(s):  
Energy Transfer ◽  
Kinetic Energy ◽  
Shear Layer ◽  
Turbulent Kinetic Energy ◽  
Viscous Dissipation ◽  
Mean Flow ◽  
Fine Grained ◽  
Large Eddy ◽  
The Mean ◽  
Three Components

In this problem a mean turbulent shear layer originally exists, homogeneous in the streamwise direction, formed perhaps by previous instabilities, but in equilibrium with the fine-grained turbulence. At a given time, a large eddy of a fixed horizontal wavenumber is initiated. We study the subsequent time development of the non-equilibrium interactions between the three components of flow as they adjust towards ultimate simultaneous equilibrium, using the integrated energy-balance conservation equations to derive the amplitude equations. This necessarily involves the usual averaging procedure and a conditional or phase-averaging procedure by which the large structure motion is educed from the total fluctuations. In general, the mean flow growth is due to the energy transfer to both fluctuating components, the large eddy gains energy from the mean motion and exchanges energy with the fine-grained turbulence, while the fine-grained turbulence gains energy from the mean flow and exchanges with the large eddy and converts its energy to heat through viscous dissipation of the smallest scales. The closure problem is obtained via the shape assumptions which enter into the interaction integrals. The situation in which the fine-grained turbulent kinetic energy production and viscous dissipation are in local balance is considered, the displacement from equilibrium being due only to the energy transfer from the large eddy. The large eddy shape is taken to be two-dimensional, instability-wavelike, with its vorticity axis perpendicular to the direction of the mean outer stream. Prior to averaging, detailed but approximate calculations of the wave-induced turbulent Reynolds stresses are obtained; the product of these stresses with the appropriate large-eddy rates of strain give the energy transfer mechanism between the two disparate scales of fluctuations. Coupled, nonlinear amplitude or energy density equations for the three components of motion are obtained, the coefficients of which are the interaction integrals guided by the shape assumptions. It is found that for the special case of parallel flow, the energy of the large eddy first undergoes a hydrodynamic-instability type of amplification but eventually decays due to the energy transfer to the fine-grained turbulence, while the turbulent kinetic energy is displaced from an original level of equilibrium to a new one because of the ability of the large eddy to negotiate an indirect energy transfer from the mean flow. For the growing shear layer, approximate considerations show that if the mechanism of energy transfer from the large to the small scale is eventually weakened by the shear layer growth compared to the large-eddy production mechanism so that the amplification and decay process repeats, ‘bursts’ of the remnant of the same large eddy will occur repeatedly until an ultimate equilibrium is reached among the three interacting components of motion. However, for the large eddy whose wavenumber corresponds to that of the initially most amplified case, the ‘bursting’ phenomenon is much less pronounced and equilibrium is very nearly reached at the end of the very first ‘burst’.

Oscillatory forcing of flow through porous media. Part 1. Steady flow

2002 ◽  
Vol 465 ◽  
pp. 213-235 ◽  
Author(s):  
D. R. GRAHAM ◽  
J. J. L. HIGDON
Keyword(s):  
Porous Media ◽  
Flow Rate ◽  
Stokes Equations ◽  
Nonlinear Flow ◽  
Navier Stokes ◽  
Full Time ◽  
Inertial Effects ◽  
Mean Flow ◽  
Navier Stokes Equations ◽  
The Mean

Oscillatory forcing of a porous medium may have a dramatic effect on the mean flow rate produced by a steady applied pressure gradient. The oscillatory forcing may excite nonlinear inertial effects leading to either enhancement or retardation of the mean flow. Here, in Part 1, we consider the effects of non-zero inertial forces on steady flows in porous media, and investigate the changes in the flow character arising from changes in both the strength of the inertial terms and the geometry of the medium. The steady-state Navier–Stokes equations are solved via a Galerkin finite element method to determine the velocity fields for simple two-dimensional models of porous media. Two geometric models are considered based on constricted channels and periodic arrays of circular cylinders. For both geometries, we observe solution multiplicity yielding both symmetric and asymmetric flow patterns. For the cylinder arrays, we demonstrate that inertial effects lead to anisotropy in the effective permeability, with the direction of minimum resistance dependent on the solid volume fraction. We identify nonlinear flow phenomena which might be exploited by oscillatory forcing to yield a net increase in the mean flow rate. In Part 2, we take up the subject of unsteady flows governed by the full time-dependent Navier–Stokes equations.

scholarly journals Dynamics of a contracting fluid compound filament with a variable density ratio

2021 ◽  
Vol 24 (2) ◽  
pp. first
Author(s):  
Truong V. Vu ◽  
Vinh T. Nguyen ◽  
Phan H. Nguyen ◽  
Nang X. Ho ◽  
Binh D. Pham ◽  
...  
Keyword(s):  
Interfacial Tension ◽  
Stokes Equations ◽  
Weighting Function ◽  
Density Ratio ◽  
Variable Density ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Front Tracking Method ◽  
Fixed Grid ◽  
Eulerian Grid

Introduction: Compound fluid filaments appear in many applications, e.g., drug delivery and processing or microfluidic systems. This paper focuses on the numerical simulation of an incompressible, immiscible, and Newtonian fluid for the contraction process of a fluid compound filament by solving the Navier-Stokes equations. The front-tracking method is used to solve this problem, which uses connected segments (Lagrangian grid) that move on a fixed grid (Eulerian grid) to represent the interface between the liquids. Methods: The interface points are advected by the velocity interpolated from those of the fixed grid using the area weighting function. The coordinates of the interface points are used to construct the indicators specifying the different fluids and compute the interfacial tension force. Results: The simulation results show that under the effects of the interfacial tension, the capsuleshaped filament can transform into a spherical compound droplet (i.e., non-breakup) or can break up into smaller spherical compound and simple droplets (i.e., breakup). When the density ratio of the outer to middle fluids increases, the filament changes from non-breakup to breakup upon contraction. Conclusion: Increasing the density ratio enhances the breakup of the compound filament during contraction. The breakup is also promoted by increasing the initial length of the filament.

A Numerical Study of Vortex Breakdown in Turbulent Swirling Flows

1999 ◽  
Vol 122 (1) ◽  
pp. 179-183 ◽  
Author(s):  
Robert E. Spall ◽  
Blake M. Ashby
Keyword(s):  
Reynolds Stress ◽  
Stokes Equations ◽  
Vortex Breakdown ◽  
Numerical Study ◽  
Reynolds Stress Model ◽  
Navier Stokes ◽  
Stress Model ◽  
Mean Flow ◽  
Navier Stokes Equations ◽  
The Mean

Solutions to the incompressible Reynolds-averaged Navier–Stokes equations have been obtained for turbulent vortex breakdown within a slightly diverging tube. Inlet boundary conditions were derived from available experimental data for the mean flow and turbulence kinetic energy. The performance of both two-equation and full differential Reynolds stress models was evaluated. Axisymmetric results revealed that the initiation of vortex breakdown was reasonably well predicted by the differential Reynolds stress model. However, the standard K-ε model failed to predict the occurrence of breakdown. The differential Reynolds stress model also predicted satisfactorily the mean azimuthal and axial velocity profiles downstream of the breakdown, whereas results using the K-ε model were unsatisfactory. [S0098-2202(00)01601-1]

The hybrid Euler–Lagrange procedure using an extension of Moffatt's method

2010 ◽  
Vol 661 ◽  
pp. 45-72 ◽  
Author(s):  
A. M. SOWARD ◽  
P. H. ROBERTS
Keyword(s):  
Fluid Dynamic ◽  
Fluid Motion ◽  
Large Eddy Simulations ◽  
Navier Stokes ◽  
Constant Density ◽  
Mean Flow ◽  
Tensor Calculus ◽  
Inviscid Incompressible Fluid ◽  
Large Eddy ◽  
The Mean

The hybrid Euler–Lagrange (HEL) description of fluid mechanics, pioneered largely by Andrews & McIntyre (J. Fluid Mech., vol. 89, 1978, pp. 609–646), has had to face the fact, in common with all Lagrangian descriptions of fluid motion, that the variables used do not describe conditions at the coordinate x, upon which they depend, but conditions elsewhere at some displaced position xL(x, t) = x + ξ(x, t), generally dependent on time t. To address this issue, we employ ‘Lie dragging’ techniques of general tensor calculus to extend a method introduced by Moffatt (J. Fluid Mech., vol. 166, 1986, pp. 359–378) in the fluid dynamic context, whereby the point x is dragged to xL(x, t) by a ‘fictitious steady flow’ η*(x, t) in a unit of ‘fictitious time’. Whereas ξ(x, t) is a Lagrangian concept intimately linked to the location xL(x, t), the ‘dragging velocity’ η*(x, t) has an essentially Eulerian character, because it describes the fictitious velocity at x itself. For the case of constant-density fluids, we show, using solenoidal η*(x, t) instead of solenoidal ξ(x, t), how the HEL theory can be cast into Eulerian form. A useful aspect of this Eulerian development is that the mean flow itself remains solenoidal, a feature that traditional HEL theories lack. Our method realizes the objective sought by Holm (Physica D, vol. 170, 2002, pp. 253–286) in his derivation of the Navier–Stokes–α equation, which is the basis of one of the methods currently employed to represent the sub-grid scales in large-eddy simulations. His derivation, based on expansion to second order in ξ, contained an error which, when corrected, implied a violation of Kelvin's theorem on the constancy of circulation in inviscid incompressible fluid. We show that this is rectified when the expansion is in η* rather than ξ, Kelvin's theorem then being satisfied to all orders for which the expansion converges. We discuss the implications of our approach using η* for the Navier–Stokes–α theory.

scholarly journals Vortex Dynamics in the Wake of Three Generic Types of Freestream Turbines

2017 ◽  
Vol 140 (2) ◽  
Author(s):  
Matthieu Boudreau ◽  
Guy Dumas
Keyword(s):  
Vortex Dynamics ◽  
Energy Recovery ◽  
Three Dimensional ◽  
Axial Flow ◽  
Navier Stokes ◽  
Mean Flow ◽  
Tip Vortices ◽  
The Mean ◽  
Optimal Efficiency ◽  
Oscillating Foil

An analysis of the vortex dynamics in the wake of three different freestream turbine concepts is conducted to gain a better understanding of the main processes affecting the energy recovery in their wakes. The turbine technologies considered are the axial-flow turbine (AFT), the crossflow turbine (CFT), also known as the H-Darrieus turbine, and the oscillating-foil turbine (OFT). The analysis is performed on single turbines facing a uniform oncoming flow and operating near their optimal efficiency conditions at a Reynolds number of 107. Three-dimensional (3D) delayed detached-eddy simulations (DDES) are carried out using a commercial finite volume Navier–Stokes solver. It is found that the wake dynamics of the AFT is significantly affected by the triggering of an instability, while that of the CFT and the OFT are mainly governed by the mean flow field stemming from the tip vortices' induction.

Validation of Particle-Laden Turbulent Flow Simulations Including Turbulence Modulation

2015 ◽  
Vol 137 (7) ◽  
Author(s):  
Yvonne Reinhardt ◽  
Leonhard Kleiser
Keyword(s):  
Numerical Study ◽  
Density Ratio ◽  
Particle Diameter ◽  
Stokes Number ◽  
Navier Stokes ◽  
Mean Flow ◽  
Turbulent Length Scale ◽  
Turbulence Modulation ◽  
Flow Simulations ◽  
Wide Range

The objective of the present numerical study is the validation of wall-bounded, turbulent particle-laden air flow simulations for a wide range of flow and particle parameters (i.e., flow and particle Reynolds numbers, Stokes number, particle-to-fluid density ratio, ratio of particle diameter to turbulent length scale) covering the one-, two- and four-way coupling regimes. The applied computational fluid dynamics (CFD) model follows the Eulerian two-fluid approach in a Reynolds-Averaged Navier–Stokes (RANS) context and is based on the kinetic theory of granular flow (KTGF) for closures concerning the particulate phase. The fluid turbulence is modeled applying a low-Reynolds-number k–epsilon turbulence model. The main focus is put on the modeling of turbulence coupling between the fluid and the particle phase. Different from common practice, the choice of a model accounting for turbulence modulation is made dependent on the prevailing coupling regime. For the case of four-way coupling, a new modulation model is suggested that well predicts turbulence augmentation and attenuation. The predictive capabilities of the present approach are evaluated by comparing simulation results to experimental benchmark data of various pipe and channel flows. Very good agreement with reference data is obtained for the mean flow and turbulence profiles of both phases.

Sign in / Sign up

Export Citation Format

Share Document