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.     

一种改进后的亚格子特征尺度表达式

针对大涡模拟涡粘性亚格子模式中使用网格尺度为特征尺度存在的不足,提出一种改进的亚格子特征尺度表达式,利用结合Sagaut混合尺度模式获得的改进模式与Smagorinsky模式、Germano模式和Sagaut模式等进行对比研究.时间发展混合层和中性大气边界层的大涡模拟结果表明,改进的有效亚格子特征尺度能合理反映亚格子脉动的时空分布特性,同时扩展了适用范围.在对耗散程度、流动演化过程、可解湍流强度和可解雷诺应力等方面的模拟,改进模式优于基准模式,表明亚格子脉动有效特征尺度在大涡模拟中具有一定的理论基础和广泛的应用价值.     

A new mixed subgrid-scale model for large eddy simulation of turbulent drag-reducing flows of viscoelastic fluids

A mixed subgrid-scale(SGS) model based on coherent structures and temporal approximate deconvolution(MCT) is proposed for turbulent drag-reducing flows of viscoelastic fluids. The main idea of the MCT SGS model is to perform spatial filtering for the momentum equation and temporal filtering for the conformation tensor transport equation of turbulent flow of viscoelastic fluid, respectively. The MCT model is suitable for large eddy simulation(LES) of turbulent dragreducing flows of viscoelastic fluids in engineering applications since the model parameters can be easily obtained. The LES of forced homogeneous isotropic turbulence(FHIT) with polymer additives and turbulent channel flow with surfactant additives based on MCT SGS model shows excellent agreements with direct numerical simulation(DNS) results. Compared with the LES results using the temporal approximate deconvolution model(TADM) for FHIT with polymer additives, this mixed SGS model MCT behaves better, regarding the enhancement of calculating parameters such as the Reynolds number.For scientific and engineering research, turbulent flows at high Reynolds numbers are expected, so the MCT model can be a more suitable model for the LES of turbulent drag-reducing flows of viscoelastic fluid with polymer or surfactant additives.     

An improved dynamic subgrid-scale model and its application to large eddy simulation of rotating channel flows

Rotating turbulence occurs extensively in nature and engineering circumstances. Meanwhile, understanding physical mechanisms of the rotating turbulence is important to the fundamental research of turbulence. The turbulent flow in rotating frames undergoes two kinds of Coriolis force effects. First, a secondary flow is induced in the case that there is a mean vorticity component perpendicular to the rotating axis. Second, there are augmenting or suppressing effects on the turbulence if there i…     

Predictions of Conjugate Heat Transfer in Turbulent Channel Flow Using Advanced Wall-Modeled Large Eddy Simulation Techniques

In this paper, advanced wall-modeled large eddy simulation (LES) techniques are used to predict conjugate heat transfer processes in turbulent channel flow. Thereby, the thermal energy transfer process involves an interaction of conduction within a solid body and convection from the solid surface by fluid motion. The approaches comprise a two-layer RANS–LES approach (zonal LES), a hybrid RANS–LES representative, the so-called improved delayed detached eddy simulation method (IDDES) and a non-equilibrium wall function model (WFLES), respectively. The results obtained are evaluated in comparison with direct numerical simulation (DNS) data and wall-resolved LES including thermal cases of large Reynolds numbers where DNS data are not available in the literature. It turns out that zonal LES, IDDES and WFLES are able to predict heat and fluid flow statistics along with wall shear stresses and Nusselt numbers accurately and that are physically consistent. Furthermore, it is found that IDDES, WFLES and zonal LES exhibit significantly lower computational costs than wall-resolved LES. Since IDDES and especially zonal LES require considerable extra work to generate numerical grids, this study indicates in particular that WFLES offers a promising near-wall modeling strategy for LES of conjugated heat transfer problems. Finally, an entropy generation analysis using the various models showed that the viscous entropy production is zero inside the solid region, peaks at the solid–fluid interface and decreases rapidly with increasing wall distance within the fluid region. Except inside the solid region, where steep temperature gradients lead to high (thermal) entropy generation rates, a similar behavior is monitored for the entropy generation by heat transfer process.     

不可压反应流动直接模拟和亚网格模型的检验

本文用谱方法对三维槽道不可压湍流反应流动进行了直接模拟,用直接模拟数据对大涡模拟亚网格质量流和燃烧模型进行了检验,结果发现,引入壁面阻尼修正的模型与精确值的符合比较好.     

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.     

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.     

Relevance of the subgrid-scales for large eddy simulations of turbulence-radiation interactions in a turbulent plane jet

An a priori study based on direct numerical simulation (DNS) of a non-isothermal turbulent plane jet has been carried out in order to analyse the role of the small-scales of turbulence on thermal radiation. Filtered DNS and large eddy simulation (LES) without subgrid-scale (SGS) model have been estimated for the radiative heat transfer. The comparison of the results highlights the subgrid-scale influence over the filtered radiation quantities, such as the radiative intensity and the radiative emission. The influence of the optical thickness is also studied. It is shown that the subgrid-scales are not significant near the centerline of the jet, where the radiative heat transfer is more important, and therefore that the SGS can be neglected in this configuration. However, when the optical thickness increases, the SGS become relevant and SGS modeling may be needed.     

Multi-scale properties of large eddy simulations: correlations between resolved-scale velocity-field increments and subgrid-scale quantities

We provide analytical and numerical results concerning multi-scale correlations between the resolved velocity field and the subgrid-scale (SGS) stress-tensor in large eddy simulations (LES). Following previous studies for Navier–Stokes equations, we derive the exact hierarchy of LES equations governing the spatio-temporal evolution of velocity structure functions of any order. The aim is to assess the influence of the subgrid model on the inertial range intermittency. We provide a series of predictions, within the multifractal theory, for the scaling of correlation involving the SGS stress and we compare them against numerical results from high-resolution Smagorinsky LES and from a-priori filtered data generated from direct numerical simulations (DNS). We find that LES data generally agree very well with filtered DNS results and with the multifractal prediction for all leading terms in the balance equations. Discrepancies are measured for some of the sub-leading terms involving cross-correlation between resolved velocity increments and the SGS tensor or the SGS energy transfer, suggesting that there must be room to improve the SGS modelisation to further extend the inertial range properties for any fixed LES resolution.     

A priori analysis of sub-grid variance of a reactive scalar using DNS data of high Ka flames

Direct numerical simulations (DNS) of low and high Karlovitz number (Ka) flames are analysed to investigate the behaviour of the reactive scalar sub-grid scale (SGS) variance in premixed combustion under a wide range of combustion conditions (regimes). An order of magnitude analysis is performed to assess the importance of various terms in the variance evolution equation and the analysis is validated using the DNS results. This analysis sheds light on the relative behaviour among turbulent transport and production, scalar dissipation and chemical processes involved in the evolution of the SGS variance at different Ka. The common expectation is that the variance equation shifts from a reaction-dissipation balance at low Ka to a production–dissipation balance at high Ka with diminishing reaction contribution. However, in large eddy simulation (LES), a high Ka alone does not make the reaction term negligible, as the relative importance of the reaction term has a concurrent increase with filter size. The filter size can be relatively large compared with the Kolmogorov length scale in practical LES of high Ka flames, and as a consequence a reaction–production–dissipation balance may prevail in the variance equation even in a high Ka configuration, and this possibility is quantified using the DNS analysis in this work. This has implications from modelling perspectives, and therefore two commonly used closures in LES for the SGS scalar dissipation rate are investigated a priori to estimate the importance of the above balance in LES modelling. The results are explained to highlight the interplay among turbulence, chemistry and dissipation processes as a function of Ka.     

A priori evaluation of large eddy simulation subgrid-scale scalar flux models in isotropic passive-scalar and anisotropic buoyancy-driven homogeneous turbulence

Direct numerical simulation (DNS) of passive (non-buoyant) and active (buoyant) scalar homogeneous turbulence is carried out using a standard pseudo-spectral numerical method. The flow settings simulated include stationary forced and decaying passive-scalar turbulence, as well as decaying anisotropic active-scalar turbulence. The Schmidt number is unity in all cases. The results are compared with, and are found to be in very good agreement with, previous similar DNS studies. The well-validated DNS data are divided into 19 sets, and are employed to study different large eddy simulation (LES) subgrid-scale (SGS) models for the SGS scalar flux. The models examined include three eddy-viscosity-type models (Smagorinsky, Vreman and Sigma with a constant SGS Schmidt number), a Dynamic Structure model and two versions of the Gradient (Gradient and Modulated Gradient) model. The models are investigated with respect to their ability to predict the orientation, and the magnitude, of the SGS scalar flux. Eddy-viscosity models are found to predict the magnitude of the SGS scalar flux accurately, but are poor at predicting the orientation of the SGS scalar flux. The Dynamic Structure and Gradient models are better than eddy-viscosity models at predicting both the magnitude and direction. However, neither of them can be realised in an actual LES, without carrying additional transport equations. Based on these observations, four new models are proposed – combining directions from Dynamic Structure and Gradient models, and magnitudes from Smagorinsky and Vreman eddy-viscosity models. These models are expected to be better than eddy-viscosity and Modulated Gradient models, and this is confirmed by preliminary a posteriori tests.     

Large-eddy simulations of the Richtmyer-Meshkov instability of rectangular interfaces accelerated by shock waves

A parallel algorithm and code MVFT(multi-viscous-fluid and turbulence) of large-eddy simulation(LES) is developed from our MVPPM(multi-viscous-fluid piecewise parabolic method),and performed to solve the multi compressible Navier-Stokes(N-S) equations.The effect of the unresolved subgrid-scale(SGS) motions on the large scales is represented by different SGS stress models in LES.A Richtmyer-Meshkov instability experiment of the evolution of a rectangular block of SF6,which occupies half of the height of the ...     

Influences of subgrid scale dynamics on resolvable scale statistics in large-eddy simulations

Recently, the -expansion and recursive renormalization group (RNG) theories as well as approximate inertial manifolds (AIM) have been exploited as means of systematically modeling subgrid scales in large-eddy simulations (LES). Although these theoretical approaches are rather complicated mathematically, their key approximations can be investigated using direct numerical simulations (DNS). In fact, the differences among these theories can be traced to whether they retain or neglect interactions between the subgrid-subgrid and subgrid-resolvable scales. In this paper, we focus on the influence of these two interactions on the evolution of the resolvable scales in LES: the effect^A which keeps only the interactions between the small and large scales; and, the effect^B which, on the other hand, keeps only the interactions among the subgrid-subgrid scales. The performance of these models is analyzed using the velocity fields of the direct numerical simulations. Specifically, our comparison is based on the analysis of the energy and enstrophy spectra, as well as higher-order statistics of the velocity and velocity derivatives. We found that the energy spectrum and higher-order statistics for the simulations with the effect^A (referred to, hereafter, as model^A) are in very good agreement with the filtered DNS. The comparison between the computations with effect^B (referred to, hereafter, as model^B) and the filtered DNS, however, is not satisfactory. Moreover, the decorrelation between the filtered DNS and model^A is much slower than that of the filtered DNS and model^B. Therefore, we conclude that the model^A, taking into acciunt the interactions between the subgrid and resolvable scales, is a faithful subgrid model for LES for the range of Reynolds numbers considered.     

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.     

Lagrangian view of time irreversibility of fluid turbulence

A turbulent flow is maintained by an external supply of kinetic energy, which is eventually dissipated into heat at steep velocity gradients. The scale at which energy is supplied greatly differs from the scale at which energy is dissipated, the more so as the turbulent intensity(the Reynolds number) is larger. The resulting energy flux over the range of scales, intermediate between energy injection and dissipation, acts as a source of time irreversibility. As it is now possible to follow accurately fluid particles in a turbulent flow field, both from laboratory experiments and from numerical simulations, a natural question arises: how do we detect time irreversibility from these Lagrangian data? Here we discuss recent results concerning this problem. For Lagrangian statistics involving more than one fluid particle, the distance between fluid particles introduces an intrinsic length scale into the problem. The evolution of quantities dependent on the relative motion between these fluid particles, including the kinetic energy in the relative motion, or the configuration of an initially isotropic structure can be related to the equal-time correlation functions of the velocity field, and is therefore sensitive to the energy flux through scales, hence to the irreversibility of the flow. In contrast, for singleparticle Lagrangian statistics, the most often studied velocity structure functions cannot distinguish the "arrow of time". Recent observations from experimental and numerical simulation data, however, show that the change of kinetic energy following the particle motion, is sensitive to time-reversal. We end the survey with a brief discussion of the implication of this line of work.     

Particle transport behavior in air channel flow with multi-group Lagrangian tracking

The particle motions of dispersion and transport in air channel flow are investigated using a large eddy simulation(LES) and Lagrangian trajectory method. The mean and fluctuating velocities of the fluids and particles are obtained,and the results are in good agreement with the data in the literature. Particle clustering is observed in the near-wall and low-speed regions. To reveal the evolution process and mechanism of particle dispersion and transport in the turbulent boundary layer, a multi-group Lagrangian tracking is applied when the two-phase flow has become fully developed: the fluid fields are classified into four sub-regions based on the flow characteristics, and particles in the turbulent region are divided accordingly into four groups when the gas–particle flow is fully developed. The spatiotemporal transport of the four groups of particles is then tracked and analyzed. The detailed relationship between particle dispersion and turbulent motion is investigated and discussed.     

Secondary vortices over surfaces with spanwise varying drag

A spanwise heterogeneity of streamwise drag is known to lead to the formation of large secondary motions of Prandtl's second kind. Based on the data sets extracted from direct numerical simulations (DNS) of fully developed turbulent channel flow where streamwise stripes of free-slip surface with varying spanwise extension are introduced, we investigate the topological structure of the secondary motions. We find a complex restructuring of the secondary motion with increasing extent of free-slip/no-slip region where the width of the free-slip region in viscous units appears to be one important governing parameter for the vortex formation. The most striking feature of this restructuring is a change in the rotational direction of the major vortex pair such that the related high- and low-momentum pathways are found at different locations. The present results reveal that the spanwise inhomogeneity of the Reynolds stress distribution is strongly related to the observed change of rotational direction. In addition, it is shown that the vorticity source remains largely unchanged and mainly restricted to a rather small region close to the discontinuity in the boundary condition, despite the fact that the topology of secondary motions substantially changes with variation of the spanwise length scale. This suggests a complex interplay between the vortices that are generated at the surface discontinuities and the surrounding flow.     

Proposed stochastic parameterisations of subgrid turbulence in large eddy simulations of turbulent channel flow

Stochastic and deterministic subgrid parameterisations are developed for the large eddy simulation (LES) of a turbulent channel flow with friction-velocity-based Reynolds number of Re_τ = 950 and centreline-based Reynolds number of Re₀ = 20,580. The subgrid model coefficients (eddy viscosities) are determined from the statistics of truncated reference direct numerical simulations (DNSs). The stochastic subgrid model consists of a mean-field shift, a drain eddy viscosity acting on the resolved field and a stochastic backscatter force of variance proportional to the backscatter eddy viscosity. The deterministic variant consists of a net eddy viscosity acting on the resolved field, which represents the net effect of the drain and backscatter. LES adopting the stochastic and deterministic models is shown to reproduce the time-averaged kinetic energy spectra of the DNS within the resolved scales.     

Rotation effect on near-wall turbulence statistics and flow structures

Turbulent flow in an axially rotating pipe, involving complicated physical mecha- nism of turbulence, is a typical problem for the study of rotating turbulent flow. The pipe rotation induces two effects on the flow. One is the stabilizing effect due to the centrifu- gal and Coriolis forces, which accounts for the relaminarization of the turbulence[1—3] and the reduction of the friction coefficient at the pipe wall. The behavior is also related to the wall streaks inclining to the azimuthal di…