High-order incompressible large-eddy simulation of fully inhomogeneous turbulent flows

The subgrid-scale (SGS) eddy-viscosity model developed by Vreman [Phys. Fluids 16 (2004) 3670] and its dynamic version [Phys. Fluids 19 (2007) 065110] are tested in large-eddy simulations (LES) of the turbulent flow in an Re = 12,000 lid-driven cubical cavity by comparison to the direct numerical simulation (DNS) data of Leriche and Gavrilakis [Phys. Fluids 12 (2000) 1363]. This appears to be the first test of this class of model to flows without any homogeneous flow directions, which is typical of flows in complex geometries. Additional LES predictions at Re = 18,000 and Re = 22,000 are compared to the DNS data of Leriche [J. Sci. Comp. 27 (2006)]. The new LES framework yielded excellent agreement for both the mean velocity and Reynolds stress profiles and matches DNS data better than the more traditional Smagorinsky-based SGS models.     

Scalar Statistics along Inertial Particle Trajectory in Isotropic Turbulence

The statistics of a passive scalar along inertial particle trajectory in homogeneous isotropic turbulence with a mean scalar gradient is investigated by using direct numerical simulation. We are interested in the influence of particle inertia on such statistics, which is crucial for further understanding and development of models in non-isothermal gas-particle flows. The results show that the scalar variance along particle trajectory decreases with the increasing particle inertia firstly; when the particle＇s Stokes number St is less than 1.0, it reaches the minimal value when St is around 1.0, then it increases if St increases further. However, the scalar dissipation rate along the particle trajectory shows completely contrasting behavior in comparison with the scalar variance. The mechanical-to-thermal time scale ratios averaged along particle, （r）p, are approximately two times smaller than that computed in the Eulerian frame r, and stay at nearly 1.77 with a weak dependence on particle inertia. In addition, the correlations between scalar dissipation and flow structure characteristics along particle trajectories, such as strain and vorticity, are also computed, and they reach their maximum and minimum, 0.31 and 0.25, respectively, when St is around 1.0.     

Vorticity,backscatter and counter-gradient transport predictions using two-level simulation of turbulent flows

The two-level simulation (TLS) method evolves both the large-and the small-scale fields in a two-scale approach and has shown good predictive capabilities in both isotropic and wall-bounded high Reynolds number (Re) turbulent flows in the past. Sensitivity and ability of this modelling approach to predict fundamental features (such as backscatter, counter-gradient turbulent transport, small-scale vorticity, etc.) seen in high Re turbulent flows is assessed here by using two direct numerical simulation (DNS) datasets corresponding to a forced isotropic turbulence at Taylor’s microscale-based Reynolds number Re_λ ≈ 433 and a fully developed turbulent flow in a periodic channel at friction Reynolds number Re_τ ≈ 1000. It is shown that TLS captures the dynamics of local co-/counter-gradient transport and backscatter at the requisite scales of interest. These observations are further confirmed through a posteriori investigation of the flow in a periodic channel at Re_τ = 2000. The results reveal that the TLS method can capture both the large- and the small-scale flow physics in a consistent manner, and at a reduced overall cost when compared to the estimated DNS or wall-resolved LES cost.     

On mean flow universality of turbulent wall flows. I. High Reynolds number flow analysis

ABSTRACT

The universality and mathematical physical structure of wall-bounded turbulent flows is a topic of discussions over many decades. There is no agreement about questions like what is the physical mean flow structure, how universal is it, and how universal are theoretical concepts for local and global flow variations. These questions are addressed by using latest direct numerical simulation (DNS) data at moderate Reynolds numbers Re and experimental data up to extreme Re. The mean flow structure is explained by analytical models for three canonical wall-bounded turbulent flows (channel flow, pipe flow, and the zero-pressure gradient turbulent boundary layer). Thorough comparisons with DNS and experimental data provide support for the validity of models. Criteria for veritable physics derived from observations are suggested. It is shown that the models presented satisfy these criteria. A probabilistic interpretation of the mean flow structure shows that the physical constraints of equal entropies and equally likely mean velocity values in a region unaffected by boundary effects impose a universal log-law structure. The structure of wall-bounded turbulent flows is much more universal than previously expected. There is no discrepancy between local logarithmic velocity variations and global friction law and bulk velocity variations. Flow effects are limited to the minimum: the difference of having a bounded or unbounded domain, and the variation range of mean velocity values allowed by the geometry.     

A one-equation turbulence model for recirculating flows

A one-equation turbulence model which relies on the turbulent kinetic energy transport equation has been developed to predict the flow properties of the recirculating flows. The turbulent eddy-viscosity coefficient is computed from a recalibrated Bradshaw’s assumption that the constant a₁ = 0.31 is recalibrated to a function based on a set of direct numerical simulation (DNS) data. The values of dissipation of turbulent kinetic energy consist of the near-wall part and isotropic part, and the isotropic part involves the von Karman length scale as the turbulent length scale. The performance of the new model is evaluated by the results from DNS for fully developed turbulence channel flow with a wide range of Reynolds numbers. However, the computed result of the recirculating flow at the separated bubble of NACA4412 demonstrates that an increase is needed on the turbulent dissipation, and this leads to an advanced tuning on the self-adjusted function. The improved model predicts better results in both the non-equilibrium and equilibrium flows, e.g. channel flows, backward-facing step flow and hump in a channel.     

Anisotropic Reynolds stress tensor representation in shear flows using DNS and experimental data

Tensorial decompositions and projections are used to study the performance of algebraic non-linear models and predict the anisotropy of the Reynolds stresses. Direct numerical simulation (DNS) data for plane channel flows at friction Reynolds number (Re_τ = 180, 395, 590, 1000), and for the boundary layer using both DNS (Re_τ = 359, 830, 1271) and experimental data (Re_τ = 2680, 3891, 4941, 7164) are used to build and evaluate the models. These data are projected into tensorial basis formed from the symmetric part of mean velocity gradient and non-persistence-of-straining tensor. Six models are proposed and their performances are investigated. The scalar coefficients for these six different levels of approximations of the Reynolds stress tensor are derived, and made dimensionless using the classical turbulent scales, the kinetic turbulent energy (κ) and its dissipation rate (ε). The dimensionless coefficients are then coupled with classical wall functions. One model is selected by comparing the predicted Reynolds stress components with experimental and DNS data, presenting a good prediction for the shear and normal Reynolds stresses.     

Conditional velocity statistics in the double scalar mixing layer – A mapping closure approach

In this work we use 3D direct numerical simulations (DNS) to investigate the average velocity conditioned on a conserved scalar in a double scalar mixing layer (DSML). The DSML is a canonical multistream flow designed as a model problem for the extensively studied piloted diffusion flames. The conditional mean velocity appears as an unclosed term in advanced Eulerian models of turbulent non-premixed combustion, like the conditional moment closure and transported probability density function (PDF) methods. Here it accounts for inhomogeneous effects that have been found significant in flames with relatively low Damköhler numbers. Today there are only a few simple models available for the conditional mean velocity and these are discussed with reference to the DNS results. We find that both the linear model of Kutznetzov and the Li and Bilger model are unsuitable for multi stream flows, whereas the gradient diffusion model of Pope shows very close agreement with DNS over the whole range of the DSML. The gradient diffusion model relies on a model for the conserved scalar PDF and here we have used a presumed mapping function PDF, that is known to give an excellent representation of the DNS. A new model for the conditional mean velocity is suggested by arguing that the Gaussian reference field represents the velocity field, a statement that is evidenced by a near perfect agreement with DNS. The model still suffers from an inconsistency with the unconditional flux of conserved scalar variance, though, and a strategy for developing fully consistent models is suggested.     

The effect of subgrid-scale models on the entrainment of a passive scalar in a turbulent planar jet

Classical large-eddy simulation (LES) modelling assumes that the passive subgrid-scale (SGS) models do not influence large-scale quantities, even though there is now ample evidence of this in many flows. In this work, direct numerical simulation (DNS) and large-eddy simulations of turbulent planar jets at Reynolds number Re_H = 6000 including a passive scalar with Schmidt number Sc = 0.7 are used to study the effect of several SGS models on the flow integral quantities e.g. velocity and scalar jet spreading rates. The models analysed are theSmagorinsky, dynamic Smagorinsky, shear-improved Smagorinsky and the Vreman. Detailed analysis of the thin layer bounding the turbulent and non-turbulent regions – the so-called turbulent/non-turbulent interface (TNTI) – shows that this region raises new challenges for classical SGS models. The small scales are far from equilibrium and contain a high fraction of the total kinetic energy and scalar variance, but the situation is worse for the scalar than for the velocity field. Both a-priori and a-posteriori (LES) tests show that the dynamic Smagorinsky and shear-improved models give the best results because they are able to accurately capture the correct statistics of the velocity and passive scalar fluctuations near the TNTI. The results also suggest the existence of a critical resolution Δ_x, of the order of the Taylor scale λ, which is needed for the scalar field. Coarser passive scalar LES i.e. Δ_x ≥ λ results in dramatic changes in the integral quantities. This fact is explained by the dynamics of the small scales near the jet interface.     

Scaling in large Prandtl number turbulent thermal convection

We study the scaling properties of heat transfer Nu in turbulent thermal convection at large Prandtl number Pr using a quasi-linear theory. We show that two regimes arise, depending on the Reynolds number Re. At low Reynolds number, NuPr ^-1/2 and Re are a function of RaPr ^-3/2. At large Reynolds number NuPr ^1/3 and RePr are function only of RaPr ^2/3 (within logarithmic corrections). In practice, since Nu is always close to Ra ^1/3, this corresponds to a much weaker dependence of the heat transfer in the Prandtl number at low Reynolds number than at large Reynolds number. This difference may solve an existing controversy between measurements in SF₆ (large Re) and in alcohol/water (lower Re). We link these regimes with a possible global bifurcation in the turbulent mean flow. We further show how a scaling theory could be used to describe these two regimes through a single universal function. This function presents a bimodal character for intermediate range of Reynolds number. We explain this bimodality in term of two dissipation regimes, one in which fluctuation dominate, and one in which mean flow dominates. Altogether, our results provide a six parameters fit of the curve Nu(Ra, Pr) which may be used to describe all measurements at Pr≥0.7. Received 27 February 2002 / Received in final form 29 May 2002 Published online 31 July 2002     

Local multifractal thermodynamics of 3D turbulence

Using results of a direct numerical simulation (DNS) of 3D turbulence we show that the observed generalized scaling (i.e. scaling moments versus moments of different orders) is consistent with a lognormal-like distribution of turbulent energy dissipation fluctuations with moderate amplitudes for all space scales available in this DNS (beginning from the molecular viscosity scale up to largest ones). Local multifractal thermodynamics has been developed to interpret the data obtained using the generalized scaling, and a new interval of space scales with inverse cascade of generalized energy has been found between dissipative and inertial intervals of scales for sufficiently large values of the Reynolds number. Received 21 July 2000     

Large scale correlations for energy injection mechanisms in swirling turbulent flows

We report an experimental study of large scale correlations in the power injected in turbulent swirling flows generated in the gap between two coaxial rotating disks. We measure the pressure fluctuations on the blades of one disk, as well as the pressure drop between the leading and the trailing edges of the rotating blades, i.e. the local drag force. Measurements at different positions on one blade and on two successive blades display a correlation length much larger than the ones usually expected in turbulent flows. The time lag for which the correlation between two points is maximum, strongly depends on the global flow configuration. These results help us to understand the statistical properties of the injected power fluctuations in turbulent swirling flows. Received 2 September 1999     

Supersonic Turbulent Boundary Layer： DNS and RANS

We assess the performance of a few turbulence models for Reynolds averaged Navier-Stokes （RANS） simulation of supersonic boundary layers, compared to the direct numerical simulations （DNS） of supersonic flat-plate turbulent boundary layers, carried out by Gao et al. [Chin. Phys. Lett. 22（2005）1709] and Huang et al. [Sci. Chin. 48 （2005） 614], as well as some available experimental data. The assessment is made for two test cases, with incoming Mach numbers and Reynolds numbers M = 2.25, Re = 365, 000//in, and M = 4.5, Re = 1.7 × 10^7/m, respectively. It is found that in the first case the prediction of RANS models agrees well with the DNS and the experimental data, while for the second case the agreement of the DNS models with experiment is less satisfactory. The compressibility effect on the RANS models is discussed.     

Density-PDFs and Lagrangian statistics of highly compressible turbulence

In isothermal, highly compressible turbulent flows, density fluctuations follow a log-normal distribution. We establish a connection between these density fluctuations and the probability-density-functions (PDF) of Lagrangian tracer particles advected with the flow. Our predicted particle statistics is tested against large scale numerical simulations, which were performed with 512³ collocation points and 2 million tracer particles integrated over several dynamical times.     

On mean flow universality of turbulent wall flows. II. Asymptotic flow analysis

Understanding of the structure of turbulent flows at extreme Reynolds numbers (Re) is relevant because of several reasons: almost all turbulence theories are only valid in the high Re limit, and most turbulent flows of practical relevance are characterized by very high Re. Specific questions about wall-bounded turbulent flows at extreme Re concern the asymptotic laws of the mean velocity and turbulence statistics, their universality, the convergence of statistics towards their asymptotic profiles, and the overall physical flow organization. In extension of recent studies focusing on the mean flow at moderate and relatively high Re, the latter questions are addressed with respect to three canonical wall-bounded flows (channel flow, pipe flow, and the zero-pressure gradient turbulent boundary layer). Main results reported here are the asymptotic logarithmic law for the mean velocity and corresponding scale-separation laws for bulk flow properties, the Reynolds shear stress, the turbulence production and turbulent viscosity. A scaling analysis indicates that the establishment of a self-similar turbulence state is the condition for the development of a strict logarithmic velocity profile. The resulting overall physical flow structure at extreme Re is discussed.     

Collisionless magnetohydrodynamic turbulence in two dimensions

In this paper we give a formulation of two-dimensional (2D) collisionless magnetohydrodynamic (MHD) turbulence that includes the effects of both electron inertia and electron pressure (or parallel electron compressibility) and is applicable to strongly magnetized collisionless plasmas. We place particular emphasis on the departures from the 2D classical MHD turbulence results produced by the collisionless MHD effects. We investigate the fractal/multi-fractal aspects of spatial intermittency. The fractal model for intermittent collisionless MHD turbulence appears to be able to describe the observed k⁻¹ spectrum in the solar wind. Multi-fractal scaling behaviors in the inertial range are first deduced, and are then extrapolated down to the dissipative microscales. We then consider a parabolic-profile model for the singularity spectrum f (α), as an explicit example of a multi-fractal scenario. These considerations provide considerable insights into the basic mechanisms underlying spatial intermittency in 2D fully developed collisionless MHD turbulence.     

On minimum aspect ratio for duct flow facilities and the role of side walls in generating secondary flows

To the surprise of some of our colleagues, we recently recommended aspect ratios of at least 24 (instead of accepted values over last few decades ranging from 5 to 12) to minimise effects of side walls in turbulent duct flow experiments, in order to approximate the two-dimensional channel flow. Here we compile available results from hydraulics and civil engineering literature, where this was already documented in the 1980s. This is of great importance due to the large amount of computational studies (mainly direct numerical simulations, DNSs) for spanwise-periodic turbulent channel flows, and the extreme complexity of constructing a fully developed duct flow facility with aspect ratio of 24 for high Reynolds numbers with adequate probe resolution. Results from this non-traditional literature for the turbulence community are compared to our recent database of DNS of turbulent duct flows with aspect ratios ranging from 1 to 18 at Re_{τ, c} values of 180 and 330, leading to very good agreement between their experimental and our computational results at these low Reynolds numbers. The DNS results also reveal the complexity of a multitude of streamwise vortical structures in addition to the secondary corner flows (which extend up to z ? 5h). These time-dependent and meandering streamwise structures are located at the core of the duct and scale with its half-height. Comparisons of these structures with the vortical motions found in spanwise-periodic channels reveal similitudes in their time-averages and the same rate of decay of their mean kinetic energy ～ T^{? 1}_A, with T_A being the averaging time. However, differences between the two flows are identified and ideas for their future analysis are proposed.     

Investigation of the ignition processes of a multi-injection flame in a Diesel engine environment using the flamelet model

In this work, the first flamelet analysis is conducted of a highly resolved DNS of a multi-injection flame with both auto-ignition and ignition induced by flame-flame interaction. A novel method is proposed to identify the different combustion modes of ignition processes using generalized flamelet equations. The state-of-the-art DNS database generated by Rieth et al. (US National Combustion Meeting, 2019) for a multi-injection flame in a Diesel engine environment is investigated. Three-dimensional flamelets are extracted from the DNS at different time instants with a focus on auto-ignition and interaction-ignition processes. The influences of mixture field interactions and the scalar dissipation rate on the ignition process are investigated by varying the species composition boundary conditions of the transient flamelet equations. Budget analyses of the generalized flamelet equations show that the transport along the mixture fraction iso-surface is insignificant during the auto-ignition process, but becomes important when interaction-ignition occurs, which is further confirmed through a flamelet regime classification method.     

Turbulence/flame/wall interactions in non-premixed inclined slot-jet flames impinging at a wall using direct numerical simulation

In the present work, three-dimensional turbulent non-premixed oblique slot-jet flames impinging at a wall were investigated using direct numerical simulation (DNS). Two cases are considered with the Damköhler number (Da) of case A being twice that of case B. A 17 species and 73-step mechanism for methane combustion was employed in the simulations. It was found that flame extinction in case B is more prominent compared to case A. Reignition in the lower branch of combustion for case A occurs when the scalar dissipation rate relaxes, while no reignition occurs in the lower branch for case B due to excessive scalar dissipation rate. A method was proposed to identify the flame quenching edges of turbulent non-premixed flames in wall-bounded flows based on the intersections of mixture fraction and OH mass fraction iso-surfaces. The flame/wall interactions were examined in terms of the quenching distance and the wall heat flux along the quenching edges. There is essentially no flame/wall interaction in case B due to the extinction caused by excessive turbulent mixing. In contrast, significant interactions between flames and the wall are observed in case A. The quenching distance is found to be negatively correlated with wall heat flux as previously reported in turbulent premixed flames. The influence of chemical reactions and wall on flow topologies was identified. The FS/U and FC/U topologies are found near flame edges, and the NNN/U topology appears when reignition occurs. The vortex-dominant topologies, FC/U and FS/S, play an increasingly important role as the jet turbulence develops.     

An a priori assessment of the Partially Stirred Reactor (PaSR) model for MILD combustion

Moderate or Intense Low-oxygen Dilution (MILD) combustion has drawn increasing attention as it allows to avoid the thermo-chemical conditions prone to the formation of pollutant species while ensuring high energy efficiency and fuel flexibility. MILD combustion is characterized by a strong competition between turbulent mixing and chemical kinetics so that turbulence-chemistry interactions are naturally strengthened and unsteady phenomena such as local extinction and re-ignition may occur. The underlying physical mechanisms are not fully understood yet and the validation of combustion models featuring enhanced predictive capabilities is required. Within this context, high-fidelity data from Direct Numerical Simulation (DNS) represent a great opportunity for the assessment and the validation of combustion closure formulations. In this study, the performance of the Partially Stirred Reactor (PaSR) combustion model in MILD conditions is a priori assessed on Direct Numerical Simulations (DNS) of turbulent combustion of MILD mixtures in a cubical domain. Modeled quantities of interest, such as heat release rate and reaction rates of major and minor species, are compared to the corresponding filtered quantities extracted from the DNS. Different submodels for the key model parameters, i.e., the chemical time scale τ_c and the mixing time scale τ_mix, are considered and their influence on the results is evaluated. The results show that the mixing time scale is the leading scale in the investigated cases. The best agreement with the DNS data regarding the prediction of heat release rate and chemical source terms is achieved by the PaSR model that employs a local dynamic approach for the estimation of the mixing time scale. An overestimation of the OH species source terms occurs in limited zones of the computational domain, characterized by low heat release rates.