A Reynolds stress model for turbulent flows of viscoelastic fluids

A second-order closure is developed for predicting turbulent flows of viscoelastic fluids described by a modified generalised Newtonian fluid model incorporating a nonlinear viscosity that depends on a strain-hardening Trouton ratio as a means to handle some of the effects of viscoelasticity upon turbulent flows. Its performance is assessed by comparing its predictions for fully developed turbulent pipe flow with experimental data for four different dilute polymeric solutions and also with two sets of direct numerical simulation data for fluids theoretically described by the finitely extensible nonlinear elastic – Peterlin model. The model is based on a Newtonian Reynolds stress closure to predict Newtonian fluid flows, which incorporates low Reynolds number damping functions to properly deal with wall effects and to provide the capability to handle fluid viscoelasticity more effectively. This new turbulence model was able to capture well the drag reduction of various viscoelastic fluids over a wide range of Reynolds numbers and performed better than previously developed models for the same type of constitutive equation, even if the streamwise and wall-normal turbulence intensities were underpredicted.     

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…     

A fractal flame-wrinkling large eddy simulation model for premixed turbulent combustion

Combustion plays an important role in a wide variety of industrial applications, such as gas-turbines, furnaces, spark-ignition engines, and various air-breathing engines. The ability to predict and understand the behavior of reacting flows in practical devices is fundamental to improved combustors with higher efficiency and reduced levels of emissions. At present, large eddy simulation is considered the most promising approach for premixed combustion modeling since the large-scale energy containing flow structures are resolved on the grid. However, the typically thin reaction zone cannot be resolved. To overcome this difficulty flamelet models, in which the reaction is assumed to take place in thin layers, wrinkled by the turbulence can sometimes be used. In these models, the turbulent flame speed can be represented as the product of the laminar flame speed, S_u, corrected for the effects of stretch (strain and curvature) and the flame-wrinkling, Ξ. In this study, we propose to model Ξ using fractal theory. This model requires sub-models for the fractal dimension, and the inner and outer cut-offs—the latter being set by the grid. A model is proposed for the inner cut-off, whereas an empirical parameterization is used to provide the fractal dimension. The proposed model is applied to flame kernel growth in homogeneous isotropic turbulence in a fan-stirred bomb and to a lean premixed flame in a plane symmetric dump combustor. Good qualitative and quantitative agreement with experimental data were obtained for the proposed model in both cases. Comparison with other well-known turbulent flame speed closure models shows that the proposed model behaves at least as good, or even better, than the reference models.     

A restricted nonlinear large eddy simulation model for high Reynolds number flows

Recent numerical studies of the restricted nonlinear (RNL) model have demonstrated its ability to reproduce important features of wall turbulence despite a severe reduction in the number of degrees of freedom. In these prior studies, the RNL model included full resolution of the viscous term. In this work, we extend the RNL model to arbitrarily high Reynolds numbers by developing a RNL large eddy simulation (LES) framework along with a method to systematically identify an appropriate streamwise wavenumber support based on spectral properties of wall turbulence. This method leads to a band-limited RNL–LES system which is successful in reproducing some of the most important statistical features captured in previous low to moderate Reynolds number simulations, e.g. the mean velocity and second-order moment profiles. The RNL–LES framework offers a new approach to understanding the connection between coherent structures and the momentum transfer mechanisms of wall turbulence at arbitrarily high Reynolds numbers, where resolution of the viscous terms can become computationally expensive even with the relatively low computational complexity afforded through the dynamical restriction of the RNL model.     

A-priori evaluations of subgrid-scale terms for large-eddy simulation of compressible turbulent flows

The subgrid-scale terms for different formulations of the energy equation are evaluated from a-priori tests using the direct numerical simulation (DNS) data of a compressible mixing layer at a moderate Mach number of M = 0.65. To extend the generality of the results, the simulations were performed with three different initial conditions for the velocity fields. To examine the impact of strong temperature variations on the subgrid scales, a non-isothermal mixing layer with lower to upper free-stream temperature ratio of 3 is also considered. For cold simulations, with equal free-stream temperatures, the total energy equation is shown to be the best choice in view of the accuracy and the subgrid-scale modelling requirements. For hot simulations, with the free-stream temperature ratio equal to 3, the total enthalpy equation is found to be the best formulation for the energy equation. Furthermore, it is shown that the subgrid-scale pressure dilatation term, which has been largely neglected so far, is of the same order of the subgrid-scale heat flux. Based on the present results, the contribution of the subgrid-scale pressure dilatation can be up to 46% of the total sugbrid-scale activity. Moreover, the time evolutions of the volume-average mean kinetic energy, turbulent kinetic energy, production, dissipation, and pressure dilatation terms are considered. Unlike the subgrid-scale pressure dilatation term, the volume-average pressure dilatation terms are negligible, and compressibility does not affect the large-scale evolutions of the mean and turbulent kinetic energies.     

Burst detection in turbulent channel flows based on large eddy simulation databases

As one of the important coherent structures in the near-wall region, turbulent burst is responsible for the production and transport of major turbulent kinetic energy and Rey- nolds stress[1]. Nearly half of turbulent kinetic energy or Reynolds stress is produced in the near-wall region, and 80% flows in outer region only contribute 20% of them. Both ejection and sweeping events contribute 60―70% of the turbulent shear stress respec- tively[2]. Recently, turbulent burst process has been foun…     

On the effect of an anisotropy-resolving subgrid-scale model on large eddy simulation predictions of turbulent open channel flow with wall roughness

A large eddy simulation (LES) was conducted of turbulent flow in a channel with a rough wall on one side and a free surface on the other by adopting an anisotropy-resolving subgrid-scale (SGS) model. A shear Reynolds number of Re_τ = 395 was used based on the mean friction velocity and channel height. To investigate the grid dependency of the LES results caused by the SGS model, three grid resolutions were tested under the same definition of a roughness shape by using the immersed boundary method. The results obtained were compared with direct numerical simulation data with and without the wall roughness and those without the extra anisotropic term. The primary focus was on how the present anisotropic SGS model with coarser grid resolutions can properly provide the effects of roughness on the mean velocity and turbulent stresses, leading to a considerable reduction of the computational cost of LES.     

Calibration of a new very large eddy simulation (VLES) methodology for turbulent flow simulation

Science China Physics, Mechanics & Astronomy - Following the idea of Speziale’s Very Large Eddy Simulation (VLES) method, a new unified hybrid simulation approach was proposed which can...     

Grid and subgrid-scale interactions in viscoelastic turbulent flow and implications for modelling

Using direct numerical simulations of turbulent plane channel flow of homogeneous polymer solutions, described by the Finitely Extensible Nonlinear Elastic-Peterlin (FENE-P) rheological constitutive model, a-priori analyses of the filtered momentum and FENE-P constitutive equations are performed. The influence of the polymer additives on the subgrid-scale (SGS) energy is evaluated by comparing the Newtonian and the viscoelastic flows, and a severe suppression of SGS stresses and energy is observed in the viscoelastic flow. All the terms of the transport equation of the SGS kinetic energy for FENE-P fluids are analysed, and an approximated version of this equation for use in future large eddy simulation closures is suggested. The terms responsible for kinetic energy transfer between grid-scale (GS) and SGS energy (split into forward/backward energy transfer) are evaluated in the presence of polymers. It is observed that the probability and intensity of forward scatter events tend to decrease in the presence of polymers.     

Evaluation of numerical strategies for large eddy simulation of particulate two-phase recirculating flows

Predicting particle dispersion in recirculating two-phase flows is a key issue for reacting flows and a potential application of large eddy simulation (LES) methods. In this study, Euler/Euler and Euler/Lagrange LES approaches are compared in the bluff body configuration from Borée et al. [J. Borée, T. Ishima, I. Flour, The effect of mass loading and inter-particle collisions on the development of the polydispersed two-phase flow downstream of a confined bluff body, J. Fluid Mech. 443 (2001) 129–165] where glass beads are injected into a complex recirculating flow. These tests are performed for non-reacting, non-evaporating sprays but are mandatory validations before computing realistic combustion chambers. Two different codes (one explicit and compressible and the other implicit and incompressible) are also tested on the same configuration. Results show that the gas flow is well predicted by both codes. The dispersed phase is also well predicted by both codes but the Lagrangian approach predicts root-mean-square values more accurately than the Eulerian approach. The effects of mesh, solvers and numerical schemes are discussed for each method.     

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 simulation of orientation and rotation of ellipsoidal particles in isotropic turbulent flows

The rotational motion and orientational distribution of ellipsoidal particles in turbulent flows are of significance in environmental and engineering applications. Whereas the translational motion of an ellipsoidal particle is controlled by the turbulent motions at large scales, its rotational motion is determined by the fluid velocity gradient tensor at small scales, which raises a challenge when predicting the rotational dispersion of ellipsoidal particles using large eddy simulation (LES) method due to the lack of subgrid scale (SGS) fluid motions. We report the effects of the SGS fluid motions on the orientational and rotational statistics, such as the alignment between the long axis of ellipsoidal particles and the vorticity, the mean rotational energy at various aspect ratios against those obtained with direct numerical simulation (DNS) and filtered DNS. The performances of a stochastic differential equation (SDE) model for the SGS velocity gradient seen by the particles and the approximate deconvolution method (ADM) for LES are investigated. It is found that the missing SGS fluid motions in LES flow fields have significant effects on the rotational statistics of ellipsoidal particles. Alignment between the particles and the vorticity is weakened; and the rotational energy of the particles is reduced in LES. The SGS-SDE model leads to a large error in predicting the alignment between the particles and the vorticity and over-predicts the rotational energy of rod-like particles. The ADM significantly improves the rotational energy prediction of particles in LES.     

A new phase-screen method for electromagnetic wave propagation in turbulent flows using large-eddy simulation

Use of large-eddy simulation (LES) data in electromagnetic wave propagation modeling is not very common because of the high computational cost involved. A new phase-screen method is proposed to model radio wave propagation, in the atmospheric turbulence, using the resolved scale refractivity obtained from LES. The proposed method offers the same level of accuracy, as the one already existing in the literature, at much cheaper cost.     

Modelling of the subgrid scale wrinkling factor for large eddy simulation of turbulent premixed combustion

We propose a model for assessing the unresolved wrinkling factor in the large eddy simulation of turbulent premixed combustion. It relies essentially on a power-law dependence of the wrinkling factor on the filter size and an original expression for the ‘active’ corrugating strain rate. The latter is written as the turbulent strain multiplied by an efficiency function that accounts for viscous effects and the kinematic constraint of Peters. This yields functional expressions for the fractal dimension and the inner cut-off length scale, the latter being (i) filter-size independent and (ii) consistent with the Damköhler asymptotic behaviours at both large and small Karlovitz numbers. A new expression for the wrinkling factor that incorporates finite Reynolds number effects is further proposed. Finally, the model is successfully assessed on an experimental filtered database.     

Modeling of filtered heat release for large eddy simulation of compressible infinitely fast reacting flows

A priori and a posteriori studies for large eddy simulation of the compressible turbulent infinitely fast reacting shear layer are presented. The filtered heat release appearing in the energy equation is unclosed and the accuracy of different models for the filtered scalar dissipation rate and the conditional filtered scalar dissipation rate of the mixture fraction in closing this term is analyzed. The effect of different closures of the subgrid transport of momentum, energy and scalars on the modeling of the filtered heat release via the resolved fields is also considered. Three explicit models of these subgrid fluxes are explored, each with an increasing level of reconstruction and all of them regularized by a Smagorinsky-type term. It is observed that a major part of the error in the prediction of the conditional filtered scalar dissipation comes from the unsatisfactory modeling of the filtered dissipation itself. The error can be substantial in the turbulent fluctuation (rms) of the dissipation fields. It is encouraging that all models give good predictions of the mean and rms density in a posteriori LES of this flow with realistic heat release corresponding to large density change. Although a posteriori results show a small sensitivity to subgrid modeling errors in the current problem, extinction–reignition phenomena involving finite-rate chemistry would demand more accurate modeling of the dissipation rates. A posteriori results also show that the resolved fields obtained with the approximate reconstruction using moments (ARM) agree better with the filtered direct numerical simulation since the level of reconstruction in the modeled subfilter fluxes is increased.     

Large eddy simulation of turbulent combustion systems

This paper reviews recent and ongoing work on numerical models for turbulent combustion systems based on a classical LES approach. The work is confined to single-phase reacting flows. First, important physico-chemical features of combustion-LES are discussed along with several aspects of overall LES models. Subsequently, some numerical issues, in particular questions associated with the reliability of LES results, are outlined. The details of chemistry, its reduction, and tabulation are not addressed here. Second, two illustrative applications dealing with non-premixed and premixed flame configurations are presented. The results show that combustion-LES is able to provide predictions very close to measured data for configurations where the flow is governed by large turbulent structures. To meet the future demands, new key challenges in specific modelling areas are suggested, and opportunities for advancements in combustion-LES techniques are highlighted. From a predictive point of view, the main target must be to provide a reliable method to aid combustion safety studies and the design of combustion systems of practical importance.