Validation of large eddy simulation method for study of flatness and skewness of decaying compressible magnetohydrodynamic turbulence

In the present work we study potential applicability of large eddy simulation (LES) method for prediction of flatness and skewness of compressible magnetohydrodynamic (MHD) turbulence. The knowledge of these quantities characterizes non-Gaussian properties of turbulence and can be used for verification of hypothesis on Gaussianity for the turbulent flow under consideration. Prediction accuracy of these quantities by means of LES method directly determines efficiency of reconstruction of probability density function (PDF) that depends on used subgrid-scale (SGS) parameterizations. Applicability of LES approach for studying of PDF properties of turbulent compressible magnetic fluid flow is investigated and potential feasibilities of five SGS parameterizations by means of comparison with direct numerical simulation results are explored. The skewness and the flatness of the velocity and the magnetic field components under various hydrodynamic Reynolds numbers, sonic Mach numbers, and magnetic Reynolds numbers are studied. It is shown that various SGS closures demonstrate the best results depending on change of similarity numbers of turbulent MHD flow. The case without any subgrid modeling yields sufficiently good results as well. This indicates that the energy pile-up at the small scales that is characteristic for the model without any subgrid closure, does not significantly influence on determination of PDF. It is shown that, among the subgrid models, the best results for studying of the flatness and the skewness of velocity and magnetic field components are demonstrated by the Smagorinsky model for MHD turbulence and the model based on cross-helicity for MHD case. It is visible from the numerical results that the influence of a choice subgrid parametrization for the flatness and the skewness of velocity is more essential than for the same characteristics of magnetic field.     

Assessment of Subgrid-Scale Models for Decaying Compressible MHD Turbulence

Large eddy simulation method is formulated for study of compressible magnetohydrodynamic turbulence and assessment of different subgrid-scale models as applied for decaying case is performed. The filtered equations of compressible magnetopause using the mass-weighted filtering operation are obtained. Mass-weighted filtered equations for large-scale turbulent component comprise subgrid-scale terms and five models for closure of the subgrid-scale terms are suggested. In present paper the obtained results of numerical computations for large eddy simulation are compared with the results of direct numerical simulation of three-dimensional compressible magnetohydrodynamic turbulence. Assessment of five subgrid-scale models of large eddy simulation for MHD flows is fulfilled. The comparisons between large eddy simulation and direct numerical simulation are carried out regarding the temporal evolution of the global quantities kinetic and magnetic energy, cross helicity and the spectra of kinetic and magnetic energy.     

Simulating turbulence in complex geometries

We first review the state-of-the art in direct numerical simulation and present a new class of spectral methods on unstructured grids for handling complex-geometry domains. Subsequently, we concentrate on the classical problem of the turbulent wake behind a circular cylinder and compare the accuracy of spectral DNS versus other LES results available in the literature. We find that DNS provides consistent agreement with the experimental results, but that LES predictions are inconsistent and depend strongly on the interaction between numerical discretization and the subgrid model. We also demonstrate via a simple vorticity-based analysis of the turbulent near-wake that eddy-viscosity models are inappropriate for sudgrid modeling. In contrast, preliminary a priori tests suggest that scale-similarity models may be a good candidate. We close the paper by forecasting the use of dynamic DNS and comment on its role in simulating turbulence in complex geometries.     

Filtered particle tracking in isotropic turbulence and stochastic modeling of subgrid-scale dispersion

A numerical study based on the Eulerian–Lagrangian formulation is performed for dispersed phase motion in a turbulent flow. The effect of spatial filtering, commonly employed in large-eddy simulations, and the role of the subgrid scale turbulence on the statistics of heavy particles, including preferential concentration, are studied through a priori analysis of DNS of particle-laden forced isotropic turbulence. In simulations where the subgrid scale kinetic energy attains 30–35% of the total we observe the impact of residual fluid motions on particles of a smaller inertia. It is shown that neglecting the influence of subgrid scale fluctuations has a significant effect on the preferential concentration of those particles. A stochastic Langevin model is proposed to reconstruct the residual (or subgrid scale) fluid velocity along particle trajectories. The computation results for a selection of particle inertia parameters are performed to appraise the model through comparisons of particle turbulent kinetic energy and the statistics of preferential concentrations.     

Large Eddy Simulation of a Propagating Turbulent Premixed Flame

A large eddy simulation of a turbulent premixed flame propagatingthrough a chamber containing a square obstruction is presented anddiscussed. The governing equations for compressible, reacting flowsare Favre filtered and turbulence closure is achieved using thedynamic Smagorinsky subgrid model. A simple flame surface densitymodel based on the flamelet concept is employed for the subgrid scalereaction rate. The simulation gives very good agreement with experimentalresults for the speed and the shape of the flame as it propagates throughthe chamber. The peak pressures, however, are underpredicted by20–30%. Reasons for this are discussed and it is concluded that amore sophisticated combustion model is required for complex flowssuch as this one, if a more accurate prediction of the pressureis to be achieved.     

LES of spatially developing 3D compressible mixing layer

The spatial development of a turbulent compressible mixing layer is investigated by means of large eddy simulation (LES). The subgrid viscosity is represented by the so-called mixed-scale model, adapted to compressible flows. Two different shock capturing schemes and three sets of inlet white-noise perturbations are investigated. The comparison between numerical and experimental results gives an overall good agreement.     

Optimized group velocity control scheme and DNS of decaying compressible turbulence of relative high turbulent Mach number

To overcome the difficulty in the DNS of compressible turbulence at high turbulent Mach number, a new difference scheme called GVC8 is developed. We have succeeded in the direct numerical simulation of decaying compressible turbulence up to turbulent Mach number 0.95. The statistical quantities thus obtained at lower turbulent Mach number agree well with those from previous authors with the same initial conditions, but they are limited to simulate at lower turbulent Mach numbers due to the so‐called start‐up problem. The energy spectrum and coherent structure of compressible turbulent flow are analysed. The scaling law of compressible turbulence is studied. The computed results indicate that the extended self‐similarity holds in decaying compressible turbulence despite the occurrence of shocklets, and compressibility has little effects on relative scaling exponents when turbulent Mach number is not very high. Copyright © 2005 John Wiley & Sons, Ltd.     

Pr数对湍流被动标量输运的影响

采用直接数值模拟的方法,研究分子Pг数对湍流被动标量输运的影响,并提供充分的证据证明,湍流Pг数明显依赖于分子Pг数.在算例中,湍流雷诺平均PгT数与分子Pг数的倒数呈线性关系;湍流亚格子Pгt数与分子Pг数的关系较为复杂,在分子Pг数为1附近时,湍流亚格子Pгt数出现极小值.     

A thermodynamic and dynamic subgrid closure model for two‐material cells

A general and robust subgrid closure model for two‐material cells is proposed. The conservative quantities of the entire cell are apportioned between two materials, and then, pressure and velocity are fully or partially equilibrated by modeling subgrid wave interactions. An unconditionally stable and entropy‐satisfying solution of the processes has been successfully found. The solution is valid for arbitrary level of relaxation. The model is numerically designed with care for general materials and is computationally efficient without recourse to subgrid iterations or subcycling in time. The model is implemented and tested in the Lagrange‐remap framework. Two interesting results are observed in 1D tests. First, on the basis of the closure model without any pressure and velocity relaxation, a material interface can be resolved without creating numerical oscillations and/or large nonphysical jumps in the problem of the modified Sod shock tube. Second, the overheating problem seen near the wall surface can be solved by the present entropy‐satisfying closure model. The generality, robustness, and efficiency of the model make it useful in principle in algorithms, such as ALE methods, volume of fluid methods, and even some mixture models, for compressible two‐phase flow computations. Copyright © 2013 John Wiley & Sons, Ltd.     

Corrélations de vitesse lagrangienne et échelle intégrale temporelle en simulation des grandes échellesLagrangian velocity correlations and timescales by a large-eddy simulation

Velocity correlations and Lagrangian timescales are studied numerically by means of a direct numerical simulation (DNS) and a large-eddy simulation (LES) coupled with a subgrid Lagrangian stochastic model, in the case of a homogeneous and isotropic turbulence. A Langevin model is used to determine the subgrid component of the velocity of fluid particles. Numerical results of Lagrangian velocity correlations and timescales are presented. These quantities play an important role in turbulent mixing and scalar dispersion. To cite this article: G. Wei et al., C. R. Mecanique 334 (2006).     

Numerical simulation of vortical and coherent structures in compressible jet flows

The compressible flows of plane free jets and jets of the intake-stroke of a rectangular piston-engine model are investigated by numerical simulations. The observed vortical structures appear to be the well-known coherent structures of turbulent shear layers. The simulated structures are compared to experimental data by means of density fields and turbulent statistics taken from different authors. The computed flow depends on physical as well as on numerical parameters. The good agreement with the experimental data is obtained by direct simulation without any turbulence model.     

可压缩自由剪切流混合转捩大涡模拟

针对湍流气动光学效应与冲压发动机气体混合机理问题,开展了可压缩混合层流动空间模式大涡模拟和时间模式直接数值模拟研究.通过对流场(包含亚/亚混合、超/亚混合两种情况)失稳、转捩直至完全湍流的空间发展过程的研究表明,对流Mach数0.4状态下流场失稳以二维最不稳定扰动为主;非线性发展中,基频涡对并及展向涡撕裂主控流动转捩,流场发生混合转捩;转捩后脉动流场基本达到局部各向同性,此时,湍流Mach数低于0.3,流动压缩性可近似忽略.     

Lagrangian evaluations of viscous models for velocity gradient dynamics in non-stationary turbulence

Simple autonomous dynamical models of velocity gradients are found to be useful in understanding the essential physics of non-linear turbulent processes. Such models can also be employed as closure models for the Lagrangian PDF methods of turbulence computations. The pressure Hessian and the viscous processes incumbent in the exact velocity gradient evolution equation are non-local in nature. Several models have been proposed for these processes. In this work, we focus specifically on two models meant for the incumbent viscous process: the linear Lagrangian diffusion model (LLDM) and the recent fluid deformation closure model (RFDM). Performance of both the models have indeed been examined earlier, but most evaluations have been restricted to statistical stationary flow fields. In this work, we subject these models to further scrutiny. Our evaluation procedure (i) uses direct numerical simulation data of decaying isotropic (non-stationary) turbulence, (ii) follows identified fluid particles (the so-called Lagrangian evolution), (iii) uses both compressible and nearly incompressible flow fields. In nearly incompressible regime, the RFD model is found to be satisfactory, while the LLDM model overestimates viscous effects at late times. In the compressible regime, both the models show inadequacies. For compressible flows, we propose an alternative modelling strategy which shows improvement over both LLD and RFD models.     

Interaction of isotropic turbulence with a normal shock wave

The evolution of incompressible and compressible isotropic 2-d turbulent fields interacting with a normal shock wave up to Mach numbers of 2.4 was investigated by means of direct numerical simulation using an ENO scheme. A comparison of statistics with linear analysis results is presented. Vorticity amplification in the DNS agrees well with the linear theory. Energy spectra are enhanced more in the small scales than in the large scales for incoming incompressible turbulence. The amplification rate for initially compressible turbulence is comparatively small.     

Validation of Lagrangian Two-Way Coupled Point-Particle Models in Large-Eddy Simulations

The accuracy of point-particle models with two-way coupling for particles of Kolmogorov-length-scale size is assessed. Turbulent kinetic energy budgets are analyzed in physical and in spectral space. It is shown that the force projection of the two-way coupling consistently models the direct transfer of kinetic energy on the particle surfaces and the enhanced viscous dissipation in the vicinity of the particles. Direct and large-eddy simulations of particle-laden flows in isotropic decaying turbulence are conducted and compared with direct-particle fluid simulations, where the particle-fluid interaction is fully resolved. An analysis in spectral space shows that turbulence modulation by particles mainly occurs at larger scales, although the momentum transfer takes place at the smallest scales. Therefore, the turbulent kinetic energy cascade of the single phase dominates in particle-laden flows. It is shown that point-particle models do not interfere with subgrid scale models, which usually act on the smallest scale. Consequently, point-particle models predict sufficiently accurate the turbulence modulation in direct numerical simulations and even when a subgrid scale model is used. The resolution of the LES does not affect the accuracy of the point-particle model, when the subgrid kinetic energy is negligible.     

An unsteady incompressible Navier–Stokes solver for large eddy simulation of turbulent flows

An unsteady incompressible Navier–Stokes solver that uses a dual time stepping method combined with spatially high‐order‐accurate finite differences, is developed for large eddy simulation (LES) of turbulent flows. The present solver uses a primitive variable formulation that is based on the artificial compressibility method and various convergence–acceleration techniques are incorporated to efficiently simulate unsteady flows. A localized dynamic subgrid model, which is formulated using the subgrid kinetic energy, is employed for subgrid turbulence modeling. To evaluate the accuracy and the efficiency of the new solver, a posteriori tests for various turbulent flows are carried out and the resulting turbulence statistics are compared with existing experimental and direct numerical simulation (DNS) data. Copyright © 1999 John Wiley & Sons, Ltd.     

Assessment of LES Subgrid-scale Models and Investigation of Hydrodynamic Behaviour for an Axisymmetrical Bluff Body Flow

This work is concerned with the investigation of fluid-mechanical behaviour and the performance of different subgrid-scale models for LES in the numerical prediction of a confined axisymmetrical bluff-body flow. Four subgrid-scale turbulence models comprising the Smagorinsky model, Dynamic Smagorinsky model, WALE model and subgrid turbulent kinetic energy model, are validated and compared directly against the experimental data. Two different mesh counts are used for the LES studies, one with a higher mesh resolution in the shear layer than the other. It is found that increasing the mesh resolution improves the time-averaged fluctuating velocity profiles, but has less effect on the time-averaged filtered velocity profiles. A comparison against experiment shows that the recirculation zone length is well predicted using LES. The accuracy of the four different subgrid scale models is then assessed by comparing the LES results using the dense mesh with the experiment. Comparisons with the time-averaged axial and radial velocity profiles demonstrate that LES displays good agreement with the experimental data, with the essential flow features captured both qualitative and quantitatively. The subgrid velocity also matches well with the experimental results, but a slight underprediction of the inner shear layer is observed for all subgrid models. In general, it is found that the Smagorinsky and WALE models are more dissipative than the Dynamic Smagorinsky model and subgrid TKE model. Comparison of the spectra against the experiment shows that LES can capture dominant features of the turbulent flow with reasonable accuracy, and weak spectral peaks related to the Kevin-Helmholtz instability and helical vortex shedding are present.     

Implicit LES of free and wall‐bounded turbulent flows based on the discontinuous Galerkin/symmetric interior penalty method

This paper presents the second validation step of a compressible discontinuous Galerkin solver with symmetric interior penalty (DGM/SIP) for the direct numerical simulation (DNS) and the large eddy simulation (LES) of complex flows. The method has already been successfully validated for DNS of an academic flow and has been applied to flows around complex geometries (e.g. airfoils and turbomachinery blades). During these studies, the advantages of the dissipation properties of the method have been highlighted, showing a natural tendency to dissipate only the under‐resolved scales (i.e the smallest scales present on the mesh), leaving the larger scales unaffected. This phenomenon is further enhanced as the polynomial order is increased. Indeed, the order increases the dissipation at the largest wave numbers, while its range of impact is reduced. These properties are spectrally compatible with a subgrid‐scale model, and hence DGM may be well suited to be used for an implicit LES (ILES) approach. A validation of this DGM/ILES approach is here investigated on canonical flows, allowing to study the impact of the discretisation on the turbulence for under‐resolved computations. The first test case is the LES of decaying homogeneous isotropic turbulence (HIT) at very high Reynolds number. This benchmark allows to assess the spectral behaviour of the method for implicit LES. The results are in agreement with theory and are even slightly more accurate than other numerical results from literature, obtained using a pseudo‐spectral (PS) method with a state‐of‐the‐art subgrid‐scale model. The second benchmark is the LES of the channel flow. Three Reynolds numbers are considered: Re_τ=395, 590 and 950. The results are compared with DNS of Moser et al. and Hoyas et al., also using PS methods. Both averaged velocity and fluctuations are globally in good agreement with the reference, showing the ability of the method to predict equilibrium wall‐bounded flow turbulence. To show that the method is able to perform accurate DNS, a DNS of HIT at Re_λ=64 and a DNS of the channel flow at Re_τ=180 are also performed. The effects of the grid refinement are investigated on the channel flow at Re_τ=395, highlighting the improvement of the results when refining the mesh in the spanwise direction. Finally, the modification of the ILES parameters, that is the Riemann solver and of the SIP coefficient, is studied on both cases, showing a significant influence on the choice of the Riemann solver. Copyright © 2015 John Wiley & Sons, Ltd.     

磁流体湍流及数值模拟研究综述

自从20世纪90年代初,AJAX计划的提出掀起了磁流体技术的研究热潮.本文总结了国内外关于磁流体湍流研究的成果,包括两方面:磁流体湍流的特点和数值模型.相比于普通湍流,磁流体湍流表现出电磁诱导的各向异性、线性焦耳耗散、湍流重层流化等特性,本文同时对磁流体湍流的结构和脉动抑制现象进行了描述:数值模型研究从雷诺平均方程和大涡模拟两方面作了介绍.当前研究集中在不可压缩磁流体湍流,距离能够应用于航空航天问题还有很长的路要走.     

Modeling air entrainment and transport in a hydraulic jump using two-fluid RANS and DES turbulence models

Both RaNS (Reynolds-averaged Navier–Stokes) and DES (Detached Eddy Simulation) type turbulence models were used in conjunction with a two-fluid model of bubbly flow and a new subgrid air entrainment model to predict air entrainment and transport in a hydraulic jump. It was found that the void fraction profiles predicted by both methods are in agreement with the experimental data in the lower shear layer region, which contains the air bubbles entrained at the so-called toe of the hydraulic jump. In contrast, in the upper roller region behind the toe, the averaged results of the DES turbulence model gives accurate predictions while a RaNS turbulence model does not. This is because the DES turbulence model successfully captures the strong fluctuations on the free surface which allows it to entrain air near the top of the roller region. In contrast, RaNS type turbulence model results in a steady, smooth interface which fails to capture the wave-induced bubble sources in that region. To our knowledge, this study is the first successful quantitative numerical simulation of the overall void fraction profiles in a hydraulic jump.