基于人工神经网络的湍流大涡模拟方法

大涡模拟方法(LES)是研究复杂湍流问题的重要工具,在航空航天、湍流燃烧、气动声学、大气边界层等众多工程领域中具有广泛的应用前景.大涡模拟方法采用粗网格计算大尺度上的湍流结构,并用亚格子(SGS)模型近似表达滤波尺度以下的流动结构对大尺度流场的作用.传统的亚格子模型由于只利用了单点流场信息和简单的函数关系,在先验验证中相对误差较大, 在后验验证中耗散过强. 近几年来,机器学习方法在湍流建模问题中得到了越来越多的应用.本文介绍了基于人工神经网络(ANN)的湍流亚格子模型的最新进展.详细地讨论了人工神经网络混合模型、空间人工神经网络模型和反卷积人工神经网络模型的构造方法.借助于人工神经网络强大的数据插值能力,新的亚格子模型的先验精度和后验精度均有显著提升. 在先验验证中,新模型所预测的亚格子应力的相关系数超过了0.99,在预测精度上远高于传统的大涡模拟模型. 在后验验证中,新模型对各类湍流统计量和瞬态流动结构的预测都优于隐式大涡模拟方法、动态Smagorinsky模型、动态混合模型等传统模型.因此, 人工神经网络方法在发展复杂湍流的先进大涡模拟模型中具有很大的潜力.      

流对波边界层湍流运动的影响

采用大涡模拟方法和Smagorinsky亚格子模型，求解三维Navier-Stokes方程，研究了波流边界层中的湍流特性．将大涡模拟结果与相应的直接数值模拟结果和实验数据进行比较，吻合较好．获得了不同波雷诺数，不同波流比情况下的大涡模拟数据库，并由此分析了波流边界层中各种湍流统计量，如速度廓线、剪应力、湍流强度等的变化规律．     

速度估计模型在大涡模拟中的应用

介绍了速度估计模型的基本思想及其在物理空间的实现。速度估计模型依靠湍流大尺度的非线性作用估计出小尺度,从而可直接求解亚格子应力项而不需要额外的模型。本文采用该模型对不同雷诺数下的各向同性衰减湍流进行了模拟,并与直接数值模拟、理论分析和其他亚格子模型的大涡模拟结果进行了比较。也初步考察了网格分辨率、不同精度的紧致格式对大涡模拟结果的影响。     

基于人工神经网络的湍流大涡模拟方法

大涡模拟方法 (LES)是研究复杂湍流问题的重要工具,在航空航天、湍流燃烧、气动声学、大气边界层等众多工程领域中具有广泛的应用前景.大涡模拟方法采用粗网格计算大尺度上的湍流结构,并用亚格子(SGS)模型近似表达滤波尺度以下的流动结构对大尺度流场的作用.传统的亚格子模型由于只利用了单点流场信息和简单的函数关系,在先验验证中相对误差较大,在后验验证中耗散过强.近几年来,机器学习方法在湍流建模问题中得到了越来越多的应用.本文介绍了基于人工神经网络(ANN)的湍流亚格子模型的最新进展.详细地讨论了人工神经网络混合模型、空间人工神经网络模型和反卷积人工神经网络模型的构造方法.借助于人工神经网络强大的数据插值能力,新的亚格子模型的先验精度和后验精度均有显著提升.在先验验证中,新模型所预测的亚格子应力的相关系数超过了0.99,在预测精度上远高于传统的大涡模拟模型.在后验验证中,新模型对各类湍流统计量和瞬态流动结构的预测都优于隐式大涡模拟方法、动态Smagorinsky模型、动态混合模型等传统模型.因此,人工神经网络方法在发展复杂湍流的先进大涡模拟模型中具有很大的潜力.     

大涡模拟在内燃机中应用的研究进展

大涡模拟 能够以比较合理的计算成本, 提供更多详细的湍流信息, 故近年来已经广泛地应用于科学和工程领域, 并也成为内燃机缸内湍流流动与燃烧过程模拟计算的最有潜力的数值方法.综合现有研究成果, 对内燃机中大涡模拟的研究进展和模拟方法进行了比较全面的评述.介绍了大涡模拟的基本概念、方法、亚网格模型,着重讨论了内燃机缸内冷态流场、燃油喷雾过程以及两相液雾湍流燃烧大涡模拟的国内外研究进展.最后论述了大涡模拟在内燃机应用中当前需要解决的问题及其发展趋势.     

低密度流体界面不稳定性大涡模拟

采用拉伸涡亚格子尺度应力模型对湍流输运中的亚格子作用项进行模式化处理,发展了适用于可压多介质黏性流动和湍流的大涡模拟方法和代码MVFT（multi-viscous flow and turbulence）。利用MVFT代码对低密度流体界面不稳定性及其诱发的湍流混合问题进行了数值模拟。详细分析了扰动界面的发展,流场中冲击波的传播、相互作用、湍流混合区边界的演化规律,以及流场瞬时密度和湍动能的分布和发展。数值模拟获得的界面演化图像和流场中波系结构与实验结果吻合较好。三维和二维模拟结果的比较显示,两者得到的扰动界面位置、波系及湍流混合区边界基本一致,只是后期的界面构型有所不同,这也正说明湍流具有强三维效应。     

基于全局动态Vreman模型的复杂非定常湍流的大涡模拟

在湍流数值模拟方法中,大涡模拟方法可以提供丰富的大涡旋信息,已逐渐成为复杂湍流问题数值研究的重要方法。而大涡模拟中,最重要的一环是尽量准确地构建能反映流场物理本质特征的亚格子应力模型。基于该思想,将一种新型的大涡模拟亚格子应力模型-Vreman亚格子应力模型用于高雷诺数三维后台阶流动的求解,计算结果与实验结果进行对比分析结果较吻合,验证了该模型的可靠性。这是对该模型用于无任何均匀流动方向的高雷诺数复杂湍流非定常流动的首次检验,计算结果优于基于传统的Smagorinsky涡粘性的动态亚格子模型。     

槽道湍流近壁区多尺度输运特性研究

利用槽道湍流直接数值模拟的数据库和离散正交子波，对近壁湍流的多尺度输运特性进 行了研究. 通过在流向和展向分别进行子波多尺度分解，得到了近壁区湍动能在流向和展向 多尺度传输的不同性质，发现流向传输以能量的反传为主，而在展向能量存在明显的正传， 并且当过滤尺度较大时以正传为主. 近壁湍流能量传输的各向异性为进一步构造各向异 性大涡模拟亚格子模式提供了必要的参考.     

燃烧室两相流场亚网格燃烧模型的研究

在三维任意曲线坐标系下采用不同的亚网格燃烧模型对环形燃烧室火焰筒气液两相湍流瞬态反应流进行大涡模拟．计算中所采用的数学度模型有：k方程亚网格尺度模型估算亚网格湍流黏性；热通量辐射模型估算辐射换热，分别采用亚网格EBU燃烧模型（E-A model）、亚网格二阶矩输运方程模型（SOM）和亚网格二阶矩代数模型（SOM-A）估算化学反应速率．并在非交错网格系统下气相采用SIMPLE算法和混合差分格式求解，液相采用Lagrange处理，并用PSIC算法对其进行求解．通过实验结果和计算结果的比较，表明在三维任意曲线坐标系下对燃烧室火焰简两相湍流油雾燃烧流场进行大涡模拟，3种不同的亚网格燃烧模型都能真实反映两相湍流化学反应流流动及实际燃烧过程，而采用亚网格二阶矩输运方程模型稍优于其他两种亚网格燃烧模型．     

基于人工神经网络的亚格子应力建模

亚格子(SGS)应力建模在湍流大涡模拟(LES)中有着极为重要的作用. 传统亚格子应力模型存在相对误差较大、耗散过强等问题. 近年来, 计算机技术的发展使得人工神经网络(ANN)等机器学习方法逐渐成为亚格子应力建模型的新研究范式. 本文着重考虑滤波宽度及雷诺数影响, 在不可压缩槽道湍流中建立了亚格子应力的ANN模型. 该模型以滤波后的直接数值模拟(fDNS)流场物理量及滤波尺度为输入信息, 相应滤波尺度下的亚格子应力为输出量. 通过对不同滤波尺度及不同雷诺数数据的训练, ANN模型能够给出与直接数值模拟(DNS)高度吻合的亚格子应力. 此外, 模型在亚格子耗散等非ANN建模量上也有着优异的预测性能, 与基于DNS获得的对应物理量的相关系数大都在0.9以上, 较梯度模型及Smagorinsky模型有明显提升. 在后验测试中, ANN模型对流向平均速度剖面的预测同样优于梯度模型、Smagorinsky模型及隐式大涡模拟(ILES)等传统LES模型. 在脉动速度均方根预测方面, 除了某些法向位置外ANN模型的性能整体上相对其他3个模型有所提升. 然而, 随着网格尺度的增大ANN模型预测的结果与fDNS结果的偏差逐渐增大. 总之, ANN方法在发展高精度亚格子应力模型上具有很大的潜力.      

模型燃烧室内瞬态紊流的大涡模拟

本文采用三种不同亚网格尺度模型对带有V型稳定器的模型燃烧室二维瞬态紊流流动进行了大涡模拟。并在交错网格系下用SIMPLE算法和混合差分格式求解离散方程。数值研究拟不同型式入口速度分布和不同亚网格尺度模型下模型燃烧室二维瞬态紊流流场。计算结果表明不同入口速度分布和不同亚网格尺度模型对瞬态流场和出口速度分布有一定的影响。本文通过数值模拟,揭示了V型稳定器后旋涡的产生和脱落过程。通过计算结果及实验数据的比较可知,本文采用的亚网格尺度模型可以用来模拟模型燃烧室紊流流场及稳定器后面回流区的流动情况。     

Comparative study of implicit and subgrid‐scale model large‐eddy simulation techniques for low‐Reynolds number airfoil applications

A computational study of a high‐fidelity, implicit large‐eddy simulation (ILES) technique with and without the use of the dynamic Smagorinsky subgrid‐scale (SGS) model is conducted to examine the contributions of the SGS model on solutions of transitional flow over the SD7003 airfoil section. ILES without an SGS model has been shown in the past to produce comparable and sometimes favorable results to traditional SGS‐based large‐eddy simulation (LES) when applied to canonical turbulent flows. This paper evaluates the necessity of the SGS model for low‐Reynolds number airfoil applications to affirm the use of ILES without SGS‐modeling for a broader class of problems such as those pertaining to micro air vehicles and low‐pressure turbines. It is determined that the addition of the dynamic Smagorinsky model does not significantly affect the time‐mean flow or statistical quantities measured around the airfoil section for the spatial resolutions and Reynolds numbers examined in this study. Additionally, the robustness and reduced computational cost of ILES without the SGS model demonstrates the attractiveness of ILES as an alternative to traditional LES. Published 2012. This article is a US Government work and is in the public domain in the USA.     

Study of a Filtered Flamelet Formulation for Large Eddy Simulation of Premixed Turbulent Flames

Despite significant advances in the understanding and modelling of turbulent combustion, no general model has been proposed for simulating flames in industrial combustion devices. Recently, the increase in computational possibilities has raised the hope of directly solving the large turbulent scales using large eddy simulation (LES) and capturing the important time-dependant phenomena. However, the chemical reactions involved in combustion occur at very small scales and the modelling of turbulent combustion processes is still required within the LES framework. In the present paper, a recently presented model for the LES of turbulent premixed flames is presented, analysed and discussed. The flamelet hypothesis is used to derive a filtered source term for the filtered progress variable equation. The model ensures proper flame propagation. The effect of subgrid scale (SGS) turbulence on the flame is modelled through the flame-wrinkling factor. The present modelling of the source term is successfully tested against filtered direct numerical simulation (DNS) data of a V-shape flame. Further, a premixed turbulent flame, stabilised behind an expansion, is simulated. The predictions agree well with the available experimental data, showing the capabilities of the model for performing accurate simulations of unsteady premixed flames.     

Joint RANS/LES Approach to Premixed Flame Modelling in the Context of the TFC Combustion Model

We present an original timesaving joint RANS/LES approach to simulate turbulent premixed combustion. It is intended mainly for industrial applications where LES may not be practical. It is based on successive RANS/LES numerical modelling, where turbulent characteristics determined from RANS simulations are used in LES equations for estimation of the subgrid chemical source and viscosity. This approach has been developed using our TFC premixed combustion model, which is based on a generalization of the Kolmogorov’s ideas. We assume existence of small-scale statistically equilibrium structures not only of turbulence but also of the reaction zones. At the same time, non-equilibrium large-scale structures of reaction sheets and turbulent eddies are described statistically by model combustion and turbulence equations in RANS simulations or follow directly without modelling in LES. Assumption of small-scale equilibrium gives an opportunity to express the mean combustion rate (controlled by small-scale coupling of turbulence and chemistry) in the RANS and LES sub-problems in terms of integral or subgrid parameters of turbulence and the chemical time, i.e. the definition of the reaction rate is similar to that of the mean dissipation rate in turbulence models where it is expressed in terms of integral or subgrid turbulent parameters. Our approach therefore renders compatible the combustion and turbulent parts of the RANS and LES sub-problems and yields reasonable agreement between the RANS and averaged LES results. Combining RANS simulations of averaged fields with LES method (and especially coupled and acoustic codes) for simulation of corresponding nonstationary process (and unsteady combustion regimes) is a promising strategy for industrial applications. In this work we present results of simulations carried out employing the joint RANS/LES approach for three examples: High velocity premixed combustion in a channel, combustion in the shear flow behind an obstacle and the impinging flame (a premixed flame attached to an obstacle).     

Subgrid Linear Eddy Mixing and Combustion Modelling of a Turbulent Nonpremixed Piloted Jet Flame

A linear eddy model for subgrid mixing and combustion has been coupled to a large eddy simulation of the turbulent nonpremixed piloted jet flame (Sandia Flame D). For the combustion reaction, simplified, single-step, irreversible, Arrhenius kinetics are used. The large scale and the subgrid structure of the flow are compared with experimental observations and, where appropriate, with a flamelet model of the flame. The main objective of this work is to demonstrate the feasibility of the LES-LEM approach for determining the structure of the subgrid scalar dissipation rate and the turbulence-chemistry interactions. The results for the large- and subgrid-scale structure of the flow show a reasonable agreement with the experimental observations.     

Shear‐improved Smagorinsky modeling of turbulent channel flow using generalized Lattice Boltzmann equation

Generalized Lattice Boltzmann equation (GLBE) was used for computation of turbulent channel flow for which large eddy simulation (LES) was employed as a turbulence model. The subgrid‐scale turbulence effects were simulated through a shear‐improved Smagorinsky model (SISM), which is capable of predicting turbulent near wall region accurately without any wall function. Computations were done for a relatively coarse grid with shear Reynolds number of 180 in a parallelized code. Good numerical stability was observed for this computational framework. The results of mean velocity distribution across the channel showed good correspondence with direct numerical simulation (DNS) data. Negligible discrepancies were observed between the present computations and those reported from DNS for the computed turbulent statistics. Three‐dimensional instantaneous vorticity contours showed complex vortical structures that appeared in such flow geometries. It was concluded that such a framework is capable of predicting accurate results for turbulent channel flow without adding significant complications and the computational cost to the standard Smagorinsky model. As this modeling was entirely local in space it was therefore adapted for parallelization. Copyright © 2010 John Wiley & Sons, Ltd.     

Subgrid Modeling of Turbulent Premixed Flames in the Flamelet Regime

Large eddy simulation (LES) models for flamelet combustion are analyzed by simulating premixed flames in turbulent stagnation zones. ALES approach based on subgrid implementation of the linear eddy model(LEM) is compared with a more conventional approach based on the estimation of the turbulent burning rate. The effects of subgrid turbulence are modeled within the subgrid domain in the LEM-LES approach and the advection (transport between LES cells) of scalars is modeled using a volume-of-fluid (VOF) Lagrangian front tracking scheme. The ability of the VOF scheme to track the flame as a thin front on the LES grid is demonstrated. The combined LEM-LES methodology is shown to be well suited for modeling premixed flamelet combustion. The geometric characteristics of the flame surfaces, their effects on resolved fluid motion and flame-turbulence interactions are well predicted by the LEM-LES approach. It is established here that local laminar propagation of the flamelets needs to be resolved in addition to the accurate estimation of the turbulent reaction rate. Some key differences between LEM-LES and the conventional approach(es) are also discussed. This revised version was published online in July 2006 with corrections to the Cover Date.