Direct numerical simulation of passive scalar in decaying compressible turbulence

1IntroductionDirectnumericalsimulation(DNS)becomesanimportanttoolinrecentresearchofturbulence[1].DNSofcompressibleturbulenceismoredifficultthanthatoftheincompressibleturbulence.WhentheturbulentMachnumberisgreaterthan0.3theshockletsmayappearinthecompressibleturbulentflowfields.Thereasonandmechanismofshockletsexistencearenotclearyet.TheturbulentMachnumberinDNScannotbeveryhighwiththepresentexistingnumericalmethodsandcomputerresource.Fortheproblemofcompressibleisotropicturbulencewiththeinitia…     

Numerical demonstration of fluctuation dynamo at low magnetic Prandtl numbers

Direct numerical simulations of incompressible nonhelical randomly forced MHD turbulence are used to demonstrate for the first time that the fluctuation dynamo exists in the limit of large magnetic Reynolds number Rm>1 and small magnetic Prandtl number Pm<1. The dependence of the critical Rmc for dynamo on the hydrodynamic Reynolds number Re is obtained for 1 less than or similar Re less than or similar 6700. In the limit Pm<1, Rmc is about 3 times larger than for the previously well-established dynamo at large and moderate Prandtl numbers: Rmc less than or similar 200 for Re greater than or similar 6000 compared to Rmc approximately 60 for Pm>or=1. It is not yet possible to determine numerically whether the growth rate of the magnetic energy is proportional, Rm1/2 in the limit Rm-->infinity, as it should be if the dynamo is driven by the inertial-range motions at the resistive scale.     

Magnetic-field generation in Kolmogorov turbulence

We analyze the initial, kinematic stage of magnetic field evolution in an isotropic and homogeneous turbulent conducting fluid with a rough velocity field, v(l) approximately l(alpha), alpha<1. This regime is relevant to the problem of magnetic field generation in fluids with small magnetic Prandtl number, i.e., with Ohmic resistivity much larger than viscosity. We propose that the smaller the roughness exponent alpha, the larger the magnetic Reynolds number that is needed to excite magnetic fluctuations. This implies that numerical or experimental investigations of magnetohydrodynamic turbulence with small Prandtl numbers need to achieve extremely high resolution in order to describe magnetic phenomena adequately.     

Particle subgrid scale modelling in large-eddy simulations of particle-laden turbulence

This study is concerned with particle subgrid scale (SGS) modelling in large-eddy simulations (LESs) of particle-laden turbulence. Although many particle-laden LES studies have neglected the effect of the SGS on the particles, several particle SGS models have been proposed in the literature. In this research, the approximate deconvolution method (ADM) and the stochastic models of Fukagata et al. (Dynamics of Brownian particles in a turbulent channel flow, Heat Mass Transf. 40 (2004), 715–726) Shotorban and Mashayek (A stochastic model for particle motion in large-eddy simulation, J. Turbul. 7 (2006), 1–13) and Berrouk et al. (Stochastic modelling of inertial particle dispersion by subgrid motion for LES of high Reynolds number pipe flow, J. Turbul. 8 (2007), pp. 1–20) are analysed. The particle SGS models are assessed using both a priori and a posteriori simulations of inertial particles in a periodic box of decaying, homogeneous and isotropic turbulence with an initial Reynolds number of Re_λ = 74. The model results are compared with particle statistics from a direct numerical simulation (DNS). Particles with a large range of Stokes numbers are tested using various filter sizes and stochastic model constant values. Simulations with and without gravity are performed to evaluate the ability of the models to account for the crossing trajectory and continuity effects. The results show that ADM improves results but is only capable of recovering a portion of the SGS turbulent kinetic energy. Conversely, the stochastic models are able to recover sufficient SGS energy, but show a large range of results dependent on the Stokes number and filter size. The stochastic models generally perform best at small Stokes numbers, but are unable to predict preferential concentration.     

A LES-CMC formulation for premixed flames including differential diffusion

A finite volume large eddy simulation–conditional moment closure (LES-CMC) numerical framework for premixed combustion developed in a previous studyhas been extended to account for differential diffusion. The non-unity Lewis number CMC transport equation has an additional convective term in sample space proportional to the conditional diffusion of the progress variable, that in turn accounts for diffusion normal to the flame front and curvature-induced effects. Planar laminar simulations are first performed using a spatially homogeneous non-unity Lewis number CMC formulation and validated against physical-space fully resolved reference solutions. The same CMC formulation is subsequently used to numerically investigate the effects of curvature for laminar flames having different effective Lewis numbers: a lean methane–air flame with Le_eff = 0.99 and a lean hydrogen–air flame with Le_eff = 0.33. Results suggest that curvature does not affect the conditional heat release if the effective Lewis number tends to unity, so that curvature-induced transport may be neglected. Finally, the effect of turbulence on the flame structure is qualitatively analysed using LES-CMC simulations with and without differential diffusion for a turbulent premixed bluff body methane–air flame exhibiting local extinction behaviour. Overall, both the unity and the non-unity computations predict the characteristic M-shaped flame observed experimentally, although some minor differences are identified. The findings suggest that for the high Karlovitz number (from 1 to 10) flame considered, turbulent mixing within the flame weakens the differential transport contribution by reducing the conditional scalar dissipation rate and accordingly the conditional diffusion of the progress variable.     

A multi agent model for the limit order book dynamics

By means of the Howard-Busse method of the optimum theory of turbulence we investigate numerically the effect of strong rotation on the upper bound on the convective heat transport in a horizontal fluid layer of infinite Prandtl number Pr. We discuss the case of fields with one wave number for regions of Rayleigh and Taylor numbers R and Ta where no analytical asymptotic bounds on the Nusselt number Nu can be derived by the Howard-Busse method. Nevertheless we observe that when R > 10⁸ and Ta is large enough the wave number of the optimum fields comes close to the analytical asymptotic result α₁ = (R/5)^1/4. We detect formation of a nonlinear structure similar to the nonlinear vortex discussed by Bassom and Chang [Geophys. Astrophys. Fluid Dyn. 76, 223 (1994)]. In addition we obtain evidence for a reshaping of the horizontal structure of the optimum fields for large values of Rayleigh and Taylor numbers.     

Flame front analysis of high-pressure turbulent lean premixed methane–air flames

An experimental study on lean turbulent premixed methane–air flames at high pressure is conducted by using a turbulent Bunsen flame configuration. A single equivalence ratio flame at Φ = 0.6 is explored for pressures ranging from atmospheric pressure to 0.9 MPa. LDA measurements of the cold flow indicate that turbulence intensities and the integral length scale are not sensitive to pressure. Due to the decreased kinematic viscosity with increasing pressure, the turbulent Reynolds numbers increase, and isotropic turbulence scaling relations indicate a large decrease of the smallest turbulence scales. Available experimental results and PREMIX code computations indicate a decrease in laminar flame propagation velocities with increasing pressure, essentially between the atmospheric pressure and 0.5 MPa. The u′/S_L ratio increases therefore accordingly. Instantaneous flame images are obtained by Mie scattering tomography. The images and their analysis show that pressure increase generates small scale flame structures. In an attempt to generalize these results, the variance of the flamelet curvatures, the standard deviation of the flamelet orientation angle, and the flamelet crossing lengths have been plotted against which is proportional to the ratio between the integral and Taylor length scales, and which increases with pressure. These three parameters vary linearly with the ratio between large and small turbulence scales and clearly indicate the strong effect of this parameter on premixed turbulent flame dynamics and structure. An obvious consequence is the increase in flame surface density and hence burning rate with pressure, as confirmed by its direct determination from 2D tomographic images.     

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     

Large-eddy simulations of temporally accelerating turbulent channel flow

The unsteady turbulent channel flow subject to the temporal acceleration is considered in this study. Large-eddy simulations were performed to study the response of the turbulent flow to the temporal acceleration. The simulations were started with the fully developed turbulent channel flow at an initial Reynolds number of Re₀ = 3500 (based on the channel half-height and the bulk-mean velocity), and then a constant temporal acceleration was applied. During the acceleration, the Reynolds number of the channel flow increased linearly from the initial Reynolds number to the final Reynolds number of Re₁ = 22,600. The effect of grid resolution, domain size, time step size on the simulation results was assessed in a preliminary study using simulations of the accelerating turbulent flow as well as simulations of the steady turbulent channel flow at various Reynolds numbers. Simulation parameters were carefully chosen from the preliminary study to ascertain the accuracy of the simulation. From the accelerating turbulent flow simulations, the delays in the response of various flow properties to the temporal acceleration were measured. The distinctive features of the delays responsible for turbulence production, energy redistribution, and radial propagation were identified. Detailed turbulence statistics including the wall shear stress response during the acceleration were examined. The results reveal the changes in the near-wall structures during the acceleration. A self-sustaining mechanism of turbulence is proposed to explain the response of the turbulent flow to the temporal acceleration. Although the overall flow characteristics are similar between the channel and pipe flows, some differences were observed between the two flows.     

Vortex dynamics of a trapezoidal bluff body placed inside a circular pipe

The influence of Reynolds number and blockage ratio on the vortex dynamics of a trapezoidal bluff body placed inside a circular pipe is studied experimentally and numerically. Low aspect ratio, high blockage ratio, curved end conditions (junction of pipe and bluff body), axisymmetric upstream flow with shear and turbulence are some of the intrinsic features of this class of bluff body flows which have been scarcely addressed in the literature. A large range (200:200,000) of Reynolds number (Re_D) is covered in this study, encompassing all the three pipe flow regimes (laminar, transition and turbulent). Four different flow regimes are defined based on the distinct features of Strouhal number (St)–Re_D relation: steady, laminar irregular, transition and turbulent. The wake in the steady regime is stationary with no oscillations in the shear layer. The laminar regime is termed as irregular owing to irregular vortex shedding. The vortex shedding in this regime is observed to be symmetric. The emergence of separation bubble downstream of the bluff body on either side is another interesting feature of this regime, which is further observed to be symmetric. Two pairs of mean streamwise vortices are noticed in the near-wake regime, which are termed as reverse dipole-type wake topology. Beyond the irregular laminar regime, the Strouhal number falls gradually and vortex shedding becomes more periodic. This regime is named transition and occurs close to the Reynolds number at which transition to turbulence takes place in a fully developed pipe. The turbulent regime is characterised by a nearly constant Strouhal number. Typical Karman-type vortex shedding is noticed in this regime. The convection velocity, wake width formation length and irrecoverable pressure loss are quantified to highlight the influence of blockage ratio. These results will be useful to develop basic understanding of vortex dynamics of confined bluff body flow for several practical applications.     

Three-Dimensional Stability of Burgers Vortices

Burgers vortices are explicit stationary solutions of the Navier-Stokes equations which are often used to describe the vortex tubes observed in numerical simulations of three-dimensional turbulence. In this model, the velocity field is a two-dimensional perturbation of a linear straining flow with axial symmetry. The only free parameter is the Reynolds number Re = Γ/ν, where Γ is the total circulation of the vortex and ν is the kinematic viscosity. The purpose of this paper is to show that Burgers vortices are asymptotically stable with respect to small three-dimensional perturbations, for all values of the Reynolds number. This general result subsumes earlier studies by various authors, which were either restricted to small Reynolds numbers or to two-dimensional perturbations. Our proof relies on the fact that the linearized operator at Burgers vortex has a simple and very specific dependence upon the axial variable. This allows to reduce the full linearized equations to a vectorial two-dimensional problem, which can be treated using an extension of the techniques developed in earlier works. Although Burgers vortices are found to be stable for all Reynolds numbers, the proof indicates that perturbations may undergo an important transient amplification if Re is large, a phenomenon that was indeed observed in numerical simulations.     

Fractal modelling of turbulent mixing

The aim of this work is to propose a new model for turbulent flows, called the fractal model (FM), applicable both in a Reynolds averaged Navier–Stokes (RANS) and a large-eddy simulation (LES) formulation, with the ultimate goal of applying it to simulate turbulent combustion irrelevant of its mode (premixed or non-premixed). The model is able to turn itself off in the laminar zones of the flow, and in particular near walls. It is based on the fractal theory. It describes the physics of the smaller spatial scales and therefore represents a small-scales model.

FM describes the physics of the small scales of turbulence based on the phenomenological concept of vortex cascade and on the self-similar behaviour of turbulence in the inertial range. Such a model is used in each cell of a numerical calculation. A characteristic length Δ is associated to each cell, and the local energy u ³ _Δ/Δ is distributed over a certain number of eddies, which depends on the local Reynolds number Re _Δ. Each vortex of the cascade generates N _c vortices; the recursive process of vortex generation terminates at the dissipative scale level, i.e. when the eddy Reynolds number is equal to one. FM is also able to estimate the volume fraction occupied by the dissipative fine structures of turbulence; this quantity is critical in reactive turbulent flows.

The physics of small scales is summarized by a turbulent ‘viscosity’ μ_t, to be added to the molecular one. μ_t is zero where the flow is laminar and, in particular, goes to zero at solid walls. Assuming μ_t to be isotropic, FM is applicable in a RANS formulation (IFM, isotropic fractal model). The model can be extended to the anisotropic case (AFM, anisotropic fractal model) and therefore used to close the transport equations in an LES approach. In the present paper, the model (IFM) is used in a RANS approach and is validated through a test case studied experimentally by Johnson and Bennett, and numerically (with LES) by Akselvoll and Moin. The results obtained are in good agreement both with the experimental and the numerical ones. Other tests are being performed.     

Preferential concentration of heavy particles in turbulence

Particle-laden flows are of relevant interest in many industrial and natural systems. When the carrier flow is turbulent, a striking feature is the phenomenon called preferential concentration: particles denser than the fluid have the tendency to inhomogeneously distribute in space, forming clusters and depleted regions. We present an investigation of clustering of small water droplets in homogeneous and isotropic active-grid-generated turbulence. We investigate the effect of Reynolds number (R_λ) and Stokes number (St) on particles clustering in the range R_λ ～ 200?400 and St ～ 2?10. Using Voronoï diagrams, we characterise clustering level and cluster properties (geometry, typical dimension and fractality). The exact same Voronoï analysis is then applied to investigate clustering properties of specific topological points of the velocity field of homogeneous isotropic turbulence obtained from direct numerical simulations at R_λ ～ 220 and 300. The goal is to compare clustering properties of actual particles with those of such points in order to explore the relevance of possible clustering mechanisms, including centrifugal effects (heavy particles sampling preferentially low-vorticity regions) and sweep-stick mechanisms (heavy particles preferentially sticking to low-acceleration points). Our study points towards a leading role of zero-acceleration points and sweep-stick effects, at least for the experimental conditions considered in this study.     

Simulation of homogeneous turbulent shear flows at higher Reynolds numbers: numerical challenges and a remedy

In a recent study, Isaza and Collins [J. Fluid Mech., 637 (2009), pp. 213–239] found the asymptotic state of homogeneous turbulent shear flows (HTSFs) to be sensitively dependent on the initial shear parameter (), and yet be almost independent of the initial Reynolds number (R _λ≡qλ/ν). The stringent resolution criteria they employed, however, restricted their studies to relatively low Reynolds numbers. In this paper, we present higher resolution direct numerical simulations of HTSFs over a wider range of Reynolds numbers, aided in part by an improved parallelisation scheme that utilises two-dimensional domain decomposition. We maximise the time-window for our simulations by determining appropriate settings for the initial energy spectrum, viscosity and domain configuration, thereby ensuring that we attain the highest possible asymptotic Reynolds number at the chosen grid resolution. In the course of our study, we find that the pseudo-spectral method suffers from Gibbs oscillations while resolving the thin vortical structures that tend to form in HTSFs. The nonlinear growth of these oscillations leads to spurious energy buildup in the high-wavenumber region of the spectrum, and contaminates the flow field. Consequently, the growth of the integral length scale is found to be numerically stunted, well before the intended final Reynolds number is attained. The issue is rectified by the application of exponential-type spectral filters, which stabilise the simulations and extend the runtime window, permitting attainment of larger asymptotic Reynolds numbers. Various large-scale flow statistics are then studied, and their dependence on the initial value of the shear parameter and Reynolds number corroborates the findings of Isaza and Collins.     

Numerical study of influence of molecular diffusion in the Mild combustion regime

In this paper, the importance of molecular diffusion versus turbulent transport in the moderate or intense low-oxygen dilution (Mild) combustion mode has been numerically studied. The experimental conditions of Dally et al. [Proc. Combust. Inst. 29 (2002) 1147–1154] were used for modelling. The EDC model was used to describe the turbulence–chemistry interaction. The DRM-22 reduced mechanism and the GRI 2.11 full mechanism were used to represent the chemical reactions of an H₂/methane jet flame. The importance of molecular diffusion for various O₂ levels, jet Reynolds numbers and H₂ fuel contents was investigated. Results show that the molecular diffusion in Mild combustion cannot be ignored in comparison with the turbulent transport. Also, the method of inclusion of molecular diffusion in combustion modelling has a considerable effect on the accuracy of numerical modelling of Mild combustion. By decreasing the jet Reynolds number, decreasing the oxygen concentration in the airflow or increasing H₂ in the fuel mixture, the influence of molecular diffusion on Mild combustion increases.     

Ginzburg–Landau description of laminar-turbulent oblique band formation in transitional plane Couette flow

Plane Couette flow, the flow between two parallel planes moving in opposite directions, is an example of wall-bounded flow experiencing a transition to turbulence with an ordered coexistence of turbulent and laminar domains in some range of Reynolds numbers [R _g, R _t] . When the aspect-ratio is sufficiently large, this coexistence occurs in the form of alternately turbulent and laminar oblique bands. As R goes up trough the upper threshold R _t, the bands disappear progressively to leave room to a uniform regime of featureless turbulence. This continuous transition is studied here by means of under-resolved numerical simulations understood as a modelling approach adapted to the long time, large aspect-ratio limit. The state of the system is quantitatively characterised using standard observables (turbulent fraction and turbulence intensity inside the bands). A pair of complex order parameters is defined for the pattern which is further analysed within a standard Ginzburg–Landau formalism. Coefficients of the model turn out to be comparable to those experimentally determined for cylindrical Couette flow.     

Generalised scale-by-scale energy-budget equations and large-eddy simulations of anisotropic scalar turbulence at various Schmidt numbers

A necessary condition for the accurate prediction of turbulent flows using large-eddy simulation (LES) is the correct representation of energy transfer between the different scales of turbulence in the LES. For scalar turbulence, transfer of energy between turbulent length scales is described by a transport equation for the second moment of the scalar increment. For homogeneous isotropic turbulence, the underlying equation is the well-known Yaglom equation. In the present work, we study the turbulent mixing of a passive scalar with an imposed mean gradient by homogeneous isotropic turbulence. Both direct numerical simulations (DNS) and LES are performed for this configuration at various Schmidt numbers, ranging from 0.11 to 5.56. As the assumptions made in the derivation of the Yaglom equation are violated for the case considered here, a generalised Yaglom equation accounting for anisotropic effects, induced by the mean gradient, is derived in this work. This equation can be interpreted as a scale-by-scale energy-budget equation, as it relates at a certain scale r terms representing the production, turbulent transport, diffusive transport and dissipation of scalar energy. The equation is evaluated for the conducted DNS, followed by a discussion of physical effects present at different scales for various Schmidt numbers. For an analysis of the energy transfer in LES, a generalised Yaglom equation for the second moment of the filtered scalar increment is derived. In this equation, new terms appear due to the interaction between resolved and unresolved scales. In an a-priori test, this filtered energy-budget equation is evaluated by means of explicitly filtered DNS data. In addition, LES calculations of the same configuration are performed, and the energy budget as well as the different terms are thereby analysed in an a-posteriori test. It is shown that LES using an eddy viscosity model is able to fulfil the generalised filtered Yaglom equation for the present configuration. Further, the dependence of the terms appearing in the filtered energy-budget equation on varying Schmidt numbers is discussed.     

Dynamos at extreme magnetic Prandtl numbers: insights from shell models

We present a magnetohydrodynamic (MHD) shell model suitable for computation of various energy fluxes of MHD turbulence for very small and very large magnetic Prandtl numbers Pm; such computations are inaccessible to direct numerical simulations. For small Pm, we observe that both kinetic and magnetic energy spectra scale as k^?5/3 in the inertial range, but the dissipative magnetic energy scales as k^?11/3exp?(? k/k_η). Here the kinetic energy at large length scale feeds the large-scale magnetic field that cascades to small-scale magnetic field, which gets dissipated by Joule heating. The large-Pm dynamo has a similar behaviour except that the dissipative kinetic energy scales as k^?13/3. For this case, the large-scale velocity field transfers energy to the large-scale magnetic field, which gets transferred to small-scale velocity and magnetic fields; the energy of the small-scale magnetic field also gets transferred to the small-scale velocity field, and the energy thus accumulated is dissipated by the viscous force.