湍流燃烧中的概率密度函数方法

用求解速度和化学热力学参数联合概率密度函数(pdf)输运方程的方法,计算湍流燃烧问题时,湍流输运和化学反应等过程可以精确计算,无须模拟。它还可以提供比统计矩模型方法更多的信息,因此是一个很有潜力的方法。本文给出了湍流燃烧概率密度函数的输运方程,扼要地介绍了目前应用较多的随机过程模型,以及求解概率密度函数方程的Monte Carlo算法。最后引用两个例子说明概率密度函数方法的优越性。     

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

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

高阶矩湍流模型研究进展及挑战

高阶矩模型是湍流模式理论研究中的难点和前沿. 自周培源先生首次建立一般湍流的雷诺应力输运方程起, 为了更精确的预测复杂流动, 人们从未间断过对高阶矩模型的研究. 尤其进入新世纪以来, 随着计算机硬件水平的飞跃和高精度数值算法的突破, 湍流模拟方法正由RANS向LES转变. 而无论对于RANS框架、LES框架还是两者混合, 高阶矩模式都是其中先进封闭模式的代表. 基于此, 本文对高阶矩模型的发展情况进行了总结, 重点包括高阶矩模型中各项的建模方式、尺度提供方程的演化历程和数值求解技术的关键需求. 然后, 通过几类典型湍流问题展示了其相对于传统涡黏模型的优势, 并且给出了部分CFD软件对高阶矩模型的集成情况. 最后对高阶矩湍流模型未来面临的挑战和发展的方向进行了展望.      

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

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

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

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

旋转圆管湍流的大涡模拟数值研究

采用大涡模拟（LES）方法,并结合动力学亚格子尺度应力（SGS）模型,通过数值求解柱坐标系下的滤波Navier-Stokes方程,研究了绕管轴旋转圆管内的湍流流动特性.为验证计算的可靠性,以及动力学SGS模型对于旋转湍流的适用性,将大涡模拟计算所得的结果,与相应的直接模拟（DNS）结果和实验数据进行了对比验证,吻合良好.进一步对旋转圆管湍流的物理机理进行了探讨,研究了湍流特性随旋转速率的变化规律.当旋转速率增加时,湍流流动有层流化的发展趋势.基于湍动能变化的关系,分析了旋转效应对湍流脉动生成的抑制作用.     

大涡模拟及其在湍流燃烧中的应用

大涡模拟作为一种研究湍流流动和湍流燃烧的有效手段,在国 际上已经得到广泛应用。本文在回顾了大涡模拟(LES)的基本思想及其 实施方法的基础上着重介绍了前人在大涡模拟的亚格子湍流模式和亚 格子燃烧模式中的研究成果,同时给出了采用不同亚格子模式的大涡 模拟在湍流燃烧中的应用实例,指出了大涡模拟在湍流燃烧中的重要 作用,为大涡模拟的进一步发展和应用提供参考。     

NO_x生成湍流反应率数值模拟的进展

利用计算机数值模拟来优选燃烧器或炉内降低氮氧化物的各种措施,如空气分级,燃料分级（再燃烧）,烟气再循环,煤和生物质混烧等,预测其效果,受到国际上越来越大的重视．如何建立合理而经济的数学模型,是其中的关键问题之一．根据近年来国际上发表的不同研究者的研究结果,对湍流燃烧氮氧化物生成湍流反应率数值模拟的研究现状进行了评述,并且讨论了有待解决的问题和今后的发展方向．     

NO_x生成湍流反应率数值模拟的进展

利用计算机数值模拟来优选燃烧器或炉内降低氮氧化物的各种措施,如空气分级,燃料分级（再燃烧）,烟气再循环,煤和生物质混烧等,预测其效果,受到国际上越来越大的重视．如何建立合理而经济的数学模型,是其中的关键问题之一．根据近年来国际上发表的不同研究者的研究结果,对湍流燃烧氮氧化物生成湍流反应率数值模拟的研究现状进行了评述,并且讨论了有待解决的问题和今后的发展方向．     

湍流涡致振动频率特性

用数值模拟方法对固定圆柱湍流涡脱落频率与弹性圆柱湍流涡致振动频率特性进行了研究,湍流计算模型采用标准κ-ε模型,压力泊松方程提法基于非交错网格系统.研究结果表明:固定圆柱湍流绕流涡脱落频率基本不随雷诺数而变,对于同一固有频率弹性圆柱,涡振频率基本不随雷诺数而变;对于某一固定雷诺数流动涡振频率在一定范围内与系统固有频率有关.     

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.     

Modelling of a Premixed Swirl-stabilized Flame Using a Turbulent Flame Speed Closure Model in LES

This paper proposes a combustion model based on a turbulent flame speed closure (TFC) technique for large eddy simulation (LES) of premixed flames. The model was originally developed for the RANS (Reynolds Averaged Navier Stokes equations) approach and was extended here to LES. The turbulent quantities needed for calculation of the turbulent flame speed are obtained at the sub grid level. This model was at first experienced via an test case and then applied to a typical industrial combustor with a swirl stabilized flame. The paper shows that the model is easy to apply and that the results are promising. Even typical frequencies of arising combustion instabilities can be captured. But, the use of compressible LES may also lead to unphysical pressure waves which have their origin in the numerical treatment of the boundary conditions.     

不同构型和多点喷射对短环型燃烧室燃烧性能的影响

对短环形流燃烧室内有较强回流的湍流旋流流动进行了模拟,并从两个方面(燃烧室构型和多点喷射)对燃烧室性能的影响进行了分析。计算中采用Reynoldes应力湍流模型(RSM)、EBU-Arrheniue湍流燃烧模型和离散坐标辐射模型描述其燃烧流动,液相采用Lagrange法处理,气相采用SIMPLE法求解。研究表明:在燃料和空气总流量不变的情况下,燃烧室构型对燃烧室出口平均温度影响不大,对出口温度分布、燃烧室内空气流场有比较大的影响。喷嘴数目的改变对出口处的平均温度和平均速度影响不大,但是对出口截面处的温度分布影响比较大,在局部范围可能产生温度比较高的热斑。     

Large Eddy Simulation of Reactive Two-Phase Flow in an Aeronautical Multipoint Burner

Because of compressibility criteria, fuel used in aeronautical combustors is liquid. Their numerical simulation therefore requires the modeling of two-phase flames, involving key phenomena such as injection, atomization, polydispersion, drag, evaporation and turbulent combustion. In the present work, particular modeling efforts have been made on spray injection and evaporation, and their coupling to turbulent combustion models in the Large Eddy Simulation (LES) approach. The model developed for fuel injection is validated against measurements in a non-evaporating spray in a quiescent atmosphere, while the evaporation model accuracy is discussed from results obtained in the case of evaporating isolated droplets. These models are finally used in reacting LES of a multipoint burner in take-off conditions, showing the complex two-phase flame structure.     

铝粉湍流火焰诱导激波现象的实验研究

粉尘湍流火焰诱导激波问题是工业灾害研究中的重要课题．本文在自行设计的大型卧式燃烧管内，对铝粉火焰诱导激波现象进行了实验研究，测试了湍流火焰阵面前压缩波到激波的转捩过程，并将实验结果与数值模拟结果进行了比较．     

Large-Scale Turbulent Vortical Structures inside a Sudden Expansion Cylinder Chamber

A large eddy simulation (LES) is performed for turbulent flow around a bluff body inside a sudden expansion cylinder chamber, a configuration which resembles a premixed gas turbine combustor. To promote turbulent mixing and to accommodate flame stability, a flame holder is installed inside the combustion chamber. The Smagorinsky model and the Lagrangian dynamic subgrid-scale model are employed and tested. The calculated Reynolds number is 5,000 based on the bulk velocity and the diameter of inlet pipe. The simulation code is constructed by using a general coordinate system based on the physical contravariant velocity components. The predicted turbulent statistics are evaluated by comparing with the laser-doppler velocimetry (LDV) measurement data. The agreement of LES with the experimental data is shown to be satisfactory. Emphasis is placed on the time-dependent evolutions of turbulent vortical structures behind the flame holder. The numerical flow visualizations depict the behavior of large-scale vortices. The turbulent behavior behind the flame holder is analyzed by visualizing the sectional views of vortical structure. This revised version was published online in July 2006 with corrections to the Cover Date.     

LES of Low to High Turbulent Combustion in an Elevated Pressure Environment

A subgrid scale flame surface density combustion model for the Large Eddy Simulation (LES) of premixed combustion is derived and validated. The model is based on fractal characteristics of the flame surface, assuming a self similar wrinkling of the flame between smallest and largest wrinkling length scales. Experimental and direct numerical simulation databases as well as theoretical models are used to derive a model for the fractal parameters, namely the cut-off lengths and the fractal dimension suitable in the LES context. The combustion model is designed with the intent to simulate low as well as high Reynolds number premixed turbulent flame propagation and with a focus on correct scaling with pressure. The combustion model is validated by simulations of turbulent Bunsen flames with methane and propane fuel at pressure levels between 0.1 MPa and 2 MPa and at turbulence levels of $0 < u^{\prime }/s_{L}^{0} < 11$ , conditions typical for spark ignition engines. The predicted turbulent flame speed is in a very good agreement with the experimental data and a smooth transition from resolved flame wrinkling to fully modelled, nearly subgrid-only wrinkling is realized. Evaluating the influence of mesh resolution shows a predicted mean flame surface and turbulent flame speed independent of mesh resolution for cases with 9–86 % resolved flame surface. Additional simulations of a highly turbulent jet flame at 0.1 MPa and 0.5 MPa and the comparison with experimental data in terms of flame shape, velocity field and turbulent fluctuations validates the model also at conditions typical for gas turbines.     

采用详细化学反应机理模拟燃烧污染物生成的研究进展

结合燃烧污染物数值模拟研究的发展和现状,阐述了采用详细化学反应机理进行模拟的方法,重点介绍了详细化学反应机理模型以及机理简化方法,概括综述了湍流燃烧模拟方法,并结合具体实例,展示了这个领域取得的成果,最后通过分析,指出了发展前景和可以取得突破性进展的方向.     

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).