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.     

A velocity-estimation subgrid model constrained by subgrid scale dissipation

Purely dissipative eddy-viscosity subgrid models have proven very successful in large-eddy simulations (LES) at moderate resolution. Simulations at coarse resolutions where the underlying assumption of small-scale universality is not valid, warrant more advanced models. However, non-eddy viscosity models are often unstable due to the lack of sufficient dissipation. This paper proposes a simple modeling approach which incorporates the dissipative nature of existing eddy viscosity models into more physically appealing non-eddy viscosity SGS models. The key idea is to impose the SGS dissipation of the eddy viscosity model as a constraint on the non-eddy viscosity model when determining the coefficients in the non-eddy viscosity model. We propose a new subgrid scale model (RSEM), which is based on estimation of the unresolved velocity field. RSEM is developed in physical space and does not require the use of finer grids to estimate the subgrid velocity field. The model coefficient is determined such that total SGS dissipation matches that from a target SGS model in the mean or least-squares sense. The dynamic Smagorinsky model is used to provide the target dissipation. Results are shown for LES of decaying isotropic turbulence and turbulent channel flow. For isotropic turbulence, RSEM displays some level of backward dissipation, while yielding as good results as the dynamic Smagorinsky model. For channel flow, the results from RSEM are better than those from the dynamic Smagorinsky model for both statistics and instantaneous flow structures.     

格子Boltzmann亚格子模型的研究

为了将格子Boltzmann法应用于大雷诺数流动的模拟,本文将Smagorinsky亚格子模型和LBGK模型相结合,并对该亚格子LBM模型进行了研究。利用该亚格子LBM模型,对二维顶盖驱动流进行了模拟,得到了若干大雷诺数下流线图和方腔中心线上无量纲速度分布。计算结果与基准解进行比较,两者相互吻合。     

Lie-symmetry group and modeling in non-isothermal fluid mechanics

The symmetry group of the non-isothermal Navier–Stokes equations is used to develop physics-preserving turbulence models for the subgrid stress tensor and the subgrid heat flux. The Reynolds analogy is not used. The theoretical properties of the models are investigated. In particular, their compatibility with the scaling laws of the flow is proven. A numerical test, in the configuration of an air flow in a ventilated and differentially heated room is presented.     

Large eddy simulation of city micro-atmospheric environment

Air quality is one of the important conditions for a better residence life in the populated urban area and it is closed related to the micro-atmospheric environment. Atmospheric environment is controlled by air motion with multi-scales in the city, while air flows in the residence area are of micro-scale atmospheric motion. This paper introduces a modern numerical simulation method, i.e. large eddy simulation (LES), for studying micro-atmospheric flows in the city residence area. For the complex flow features in the residence area, the proper application of LES is studied and various numerical methods are compared in order to investigate their effects on the prediction accuracy of micro-atmospheric flows, for instance, roughness elements and immersed boundary method for complex terrain, different subgrid models and so on. The wind field (including turbulence properties) and contaminant dispersion are computed by the proposed method for a model and a realistic residence area, and the numerical results are in good agreement with the experimental measurements. Supported by the National Natural Science Foundation of China (Grant No. 10572073) and Foundation for Development of Science and Technology in Macao (Grant No. 022/2006/A)     

双大涡模拟方法模拟提升管内气固两相流动

基于Smagorinsky涡黏模型以及颗粒动理学理论,建立了气固两相流双大涡模拟模型。考虑大涡模拟中过滤尺度的影响,给出颗粒相亚格子压力和热传导系数计算模型。考虑颗粒聚团对两相作用的影响,给出了考虑颗粒聚团作用的气固两相多尺度曳力系数模型。数值模拟了提升管内气固两相流动特性,合理地预测出了提升管内气固两相环-核流动结构。模拟结果与Knowlton等实测结果相吻合。     

Lattice Boltzmann method for three-dimensional moving particles in a Newtonian fluid

A lattice Boltzmann method is developed to simulate three-dimensional solid particle motions in fluids. In the present model, a uniform grid is used and the exact spatial location of the physical boundary of the suspended particles is determined using an interpolation scheme. The numerical accuracy and efficiency of the proposed lattice Boltzmann method is demonstrated by simulating the sedimentation of a single sphere in a square cylinder. Highly accurate simulation results can be achieved with few meshes, compared with the previous lattice Boltzmann methods. The present method is expected to find applications on the flow systems with moving boundaries, such as the blood flow in distensible vessels, the particle－flow interaction and the solidification of alloys.     

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.     

多孔介质中流体流动及扩散的耦合格子Boltzmann模型

多孔介质中高Péclet数和大黏性比下混溶流体的流动和扩散广泛存在于二氧化碳驱油、化工生产等工业过程中.用数值方法对该问题进行研究时,关键在于如何正确描述高Péclet数和大黏性比下多孔介质内流体的行为.为此,提出了一种基于多松弛模型和格子动理模型的耦合格子Boltzmann模型.通过Chapman-Enskog分析,证明该模型能有效求解不可压Navier-Stokes方程和对流扩散方程.数值结果表明,该模型不仅具有二阶精度和良好的稳健性,而且对于高Péclet数和大黏性比的问题具有良好的数值稳定性,为模拟此类问题提供了有效工具.     

二维卡门涡街的格子Boltzmann仿真

格子气自动机和格子Boltzmann方法的迅速发展提供了一类求解流体力学问题的新的方法。本文中，我们介绍了Boltzmann方法，解决了格子气方法中的缺点，通过选择适当平衡分布及其参数，导出了Navier－Stokers方程，并得到了声速和粘性系数。最后在微机上模拟了在无限长平板流动问题及绕单分离板的流动问题，得到了卡门涡街。结果表该模型有格子气方法及其它的数值方法所没有的优点，计算更精确、更直观、更有效。     

Lattice Boltzmann Simulation of Shear-Thinning Fluids

It is shown how shear-thinning flow can be simulated without the need for numerical differentiation by following a lattice Boltzmann approach. The basic idea of is to combine the Cross model of viscosity with a 3D multiple relaxation time lattice Boltzmann method and to extract the required velocity derivatives from intrinsic quantities of the lattice Boltzmann scheme. Computational results are presented for a simple benchmark and for the simulation of liquid composite moulding.     

Filter size definition in anisotropic subgrid models for large eddy simulation on irregular grids

The definition of the characteristic filter size to be used for subgrid scales models in large eddy simulation using irregular grids is still an unclosed problem. We investigate some different approaches to the definition of the filter length for anisotropic subgrid scale models and we propose a tensorial formulation based on the inertial ellipsoid of the grid element. The results demonstrate an improvement in the prediction of several key features of the flow when the anisotropicity of the grid is explicitly taken into account with the tensorial filter size.     

基于总能形式的耦合的双分布函数热晶格玻尔兹曼数值方法

双分布函数热晶格玻尔兹曼数值方法在微尺度热流动系统中得到广泛的应用. 本文基于晶格玻尔兹曼平衡分布函数低阶Hermite展开式, 创新性地提出了包含黏性热耗散和压缩功的耦合的双分布函数热晶格玻尔兹曼数值方法, 将能量场内温度的变化以动量源的形式引入晶格波尔兹曼动量演化方程, 实现了能量场与动量场之间的耦合. 研究了考虑黏性热耗散和压缩功的和不考虑的两种热自然对流模型, 重点分析了不同瑞利数和普朗特数下流场内的流动情况以及温度、速度和平均努赛尔数的变化趋势. 本文实验结果与文献结果一致, 验证了本文数值方法的可行性和准确性. 研究结果表明: 随着瑞利数和普朗特数的增大, 方腔内对流传热作用逐渐增强, 边界处形成明显的边界层; 考虑黏性热耗散和压缩功的模型对流作用相对增强, 黏性热耗散和压缩功对自然对流的影响在微尺度流动过程中不能忽略.     

Large-eddy simulation of equilibrium plasma-assisted combustion in supersonic flow

A large-eddy simulation (LES) model with a new localized dynamic subgrid closure for the magnetohydrodynamics (MHD) equations is used to investigate plasma-assisted combustion in supersonic flow. A 16-species and 74-reactions kinetics model is used to simulate hydrogen-air combustion and high-temperature air dissociation. The numerical model is validated with experimental data for non-reacting and reacting supersonic flow over a rearward-facing step. The creation of a plasma source near the step corner is shown to have a strong localized effect with the high temperature region resulting in an increase of the radical species concentration in the mixing region. This has the potential for enhancing combustion. In addition, downstream fuel–air mixing is improved, primarily by the creation of a strong baroclinic torque effect in the near field of the plasma source. Furthermore, by adding an uniform external magnetic field, the Lorentz force effect helps to further enhance mixing by lifting the shear layer and increasing fuel penetration by approximately 20%.     

Semi-analytical validation of a dynamic large-eddy simulation procedure for turbulent premixed flames via the G-equation

The performance of a dynamic subgrid model for the turbulent burning speed of a premixed flame is investigated for a series of idealized test cases where the flame front is wrinkled by a multiple-scale shear flow; a rigorous asymptotic subgrid model is also implemented for comparison. Explicit formulae for the flame wrinkled shape and turbulent speed are available to generate a reference database. The role of the subgrid wrinkling models is to achieve the same overall flame shape and propagation speed in a simulation where only the largest scales of the flow are explicitly accounted for. Very good results are obtained when the subgrid burning speed enhancement is estimated using the asymptotic subgrid model. On the other hand, the dynamic model attempts to exploit the scaling observable in the simulation to extrapolate the turbulent burning speed enhancement in the original system. The performance of this strategy is adequate for some regimes but poor for others; the source of the problem is traced back to the existence of a scaling transition that occurs as the flame propagating speed is adjusted during the large-eddy simulation. A modification to the scaling of the enhanced burning is implemented to account for the existence of the two distinct scaling ranges; it improves significantly the predictions of the dynamic model away from the transition, but results in the near-critical range remain predictably very poor compared with the rigorous asymptotic model results. These conclusions based on a priori performance for the reference steady data are confirmed by comparing unsteady large-eddy and direct simulations. Results based on rigorous mathematical tools are possible here because of the separation of length scales in the special class of idealized flow fields used in this study: their relevance to more realistic flows is also discussed.     

Large-eddy simulation of a bluff-body stabilised turbulent premixed flame using the transported flame surface density approach

A premixed propane–air flame stabilised on a triangular bluff body in a model jet-engine afterburner configuration is investigated using large-eddy simulation (LES). The reaction rate source term for turbulent premixed combustion is closed using the transported flame surface density (TFSD) model. In this approach, there is no need to assume local equilibrium between the generation and destruction of subgrid FSD, as commonly done in simple algebraic closure models. Instead, the key processes that create and destroy FSD are accounted for explicitly. This allows the model to capture large-scale unsteady flame propagation in the presence of combustion instabilities, or in situations where the flame encounters progressive wrinkling with time. In this study, comprehensive validation of the numerical method is carried out. For the non-reacting flow, good agreement for both the time-averaged and root-mean-square velocity fields are obtained, and the Karman type vortex shedding behaviour seen in the experiment is well represented. For the reacting flow, two mesh configurations are used to investigate the sensitivity of the LES results to the numerical resolution. Profiles for the velocity and temperature fields exhibit good agreement with the experimental data for both the coarse and dense mesh. This demonstrates the capability of LES coupled with the TFSD approach in representing the highly unsteady premixed combustion observed in this configuration. The instantaneous flow pattern and turbulent flame behaviour are discussed, and the differences between the non-reacting and reacting flow are described through visualisation of vortical structures and their interaction with the flame. Lastly, the generation and destruction of FSD are evaluated by examining the individual terms in the FSD transport equation. Localised regions where straining, curvature and propagation are each dominant are observed, highlighting the importance of non-equilibrium effects of FSD generation and destruction in the model afterburner.     

Numerical simulation of instantaneous flow structure of swirling and non-swirling coaxial-jet particle-laden turbulence flows

A subgrid scale two-phase second-order-moment (SGS-SOM) model based on the two-fluid continuum approach is presented for the analysis of the instantaneous flow structures of swirling and non-swirling coaxial-jet particle-laden turbulence flows. Since the interaction between the two-phase subgrid scale stresses and the anisotropy of two-phase subgrid scale stresses is fully considered, it is superior to the conventional subgrid scale model on the basis of single gas phase or together with their similar forms for the particle phase for not taken these characters thoroughly into account. The swirling numbers s=0.47 and s=0 of coaxial-jet particle-laden turbulence flows (measured by M. Sommerfeld, H.H. Qiu, Detailed measurements in a swirling particulate two-phase flow by a phase Doppler anemometer, Int. J. Heat Fluid Flow 12 (1991) 20-28) are numerically simulated by large eddy simulation using this model, together with a Reynolds-averaged Navier-Stokes model using the unified second-order-moment two-phase turbulence model (RANS-USM). The instantaneous results show that the multiple recirculating gas flow structure is similar to that of single-phase swirling flows; but the particle flow structure contains less vortices. Both SGS-SOM and RANS-USM predicted that the two-phase time-averaged velocities and the root-mean-square fluctuation velocities are validated and are in good agreement with the experimental results. It is seen that for the two-phase time-averaged velocities both the models give almost the same results, hence the RANS-USM modeling is validated by large eddy simulation. For the two-phase root-mean-square fluctuation velocities the SGS-SOM results are obviously better than the RANS-USM results.     

大雷诺数Navier-Stokes方程的两水平亚格子模型稳定化方法

基于两重网格离散方法,提出三种求解大雷诺数定常Navier-Stokes方程的两水平亚格子模型稳定化有限元算法.其基本思想是首先在一粗网格上求解带有亚格子模型稳定项的Navier-Stokes方程,然后在细网格上分别用三种不同的校正格式求解一个亚格子模型稳定化的线性问题,以校正粗网格解.通过适当的稳定化参数和粗细网格尺寸的选取,这些算法能取得最优渐近收敛阶的有限元解.最后,用数值模拟验证三种算法的有效性.     

Formation of zones of an inverse energy cascade in an axisymmetric jet

The results of investigation into the mechanism of formation of an inverse turbulence energy cascade in an axially symmetric submerged jet based on numerical modeling by the large eddy simulation. The flow structure is calculated using two models of subgrid turbulence in a broad range of Reynolds numbers with the imposed external harmonic low-amplitude perturbation and swirling. Calculations showed that upon imposing swirling, the effect of formation of the inverse cascade is suppressed, while upon imposing the external low-amplitude harmonic perturbations corresponding to the frequency of highest perceptivity, it is enhanced. We can assume that the regions with the inverse turbulence energy flow are formed there where the dynamics of large quasi-2D structures is determined by the mechanisms of combining eddies and involvement processes, while the tension mechanism of vortex tubes is suppressed. It is shown that the balance of these mechanisms can be controlled by means of imposition of low-amplitude harmonic perturbations.