用求解概率密度函数输运方程的方法模拟湍流自由射流

用求解概率密度函数输运方程的方法模拟湍流自由射流周向阳，郑楚光，马毓义（华中理工大学煤燃烧国家重点实验室武汉４３００７４）关键词概率密度函数；湍流；自由射流１前言和通常的湍流统计矩模型（如ｋ－ε模型或Ｒｅｙｎｏｌｄｓ应力模型等）相比，用求解速度和标量...     

DSM-LPDF两相湍流模型及旋流两相流动的模拟

本文由流体－颗粒速度的拉氏联合概率密度函数（ＰＤＦ）输运方程出发,用Ｓｉｍｏｎｉｎ建议的Ｌａｎｇｅｖｉｎ模型封闭颗粒所遇到流体瞬时速度的条件期望项,并用Ｍｏｎｔｅ Ｃａｒｌｏ方法直接求解 ＰＤＦ输运方程,将其和求解流体雷诺应力方程模型的有限差分方法结合,建立了雷诺应力－拉氏ＰＤＦ（ＤＳＭ－ＬＰＤＦ,简称ＤＬ）两相湍流模型．用此模型模拟了旋流数为０．４７的突扩旋流气粒两相流动,并与文献中ＰＤＰＡ实验和用类似于单相流动湍流模型封闭方法的时平均统一二阶矩（ＵＳＭ）模型的预报进行了对比．     

湍流两相流动有燃烧颗粒相概率密度函数输运方程理论

由有燃烧的湍流气粒两相流动的瞬态方程和统计力学概率密度函数概念出发,推导了有燃烧颗粒相的质量－动量－能量联合概率密度函数（ＰＤＦ）输运方程,并对方程中条件期望项用梯度模拟概念进行了模拟封闭。封闭后的ＰＤＦ方程可作为建立颗粒拟流体模型方程和封闭二阶矩模型的基础,也可以通过Ｍｏｎｔｅ－Ｃａｒｌｏ 法求解用以直接计算颗粒雷诺应力和湍流动能,以便和二阶 矩模型的结果相对照,改善二阶矩模型。     

k－ε－PDF两相湍流模型和台阶后方气粒两相流动的模拟

本文提出了两相湍流的ｋ－ε－ＰＤＦ模型。ＰＤＦ模型所得的湍流两相流动的统计平均方程，与雷诺时均方程有相似的形式，但ＰＤＦ模型可精确计算出颗粒相各方程的脉动关联项。本文将该模型用于预报台阶后方湍流两相流动，与ｋ－ε－ｋｐ模型相比，它可以更合理地预报出颗粒湍流的各向异性。     

用直接模拟数据检验RANS二阶矩燃烧模型

用谱方法对三维槽道不可压湍流反应流动进行了直接模拟(DNS),并对关联量精确输运方程中各项的贡献进行统计,发现产生项和耗散项起重要作用.进一步用DNS数据库对BANS二阶矩(SOM)湍流燃烧模型进行检验,发现扩散项和产生项的模拟值在大部分区域都和精确统计值接近,但在近壁区需要修正;耗散项的模拟有待改进.     

U-RANS／PDF方法计算钝体后漩涡脱落

本研究发展了U-RANS/PDF混合算法研究湍流和化学反应相互作用对燃烧稳定性的影响,采用有限体积和Monte Carlo相结合的方法在非结构网格中求解相容的U-RANS方程和脉动速度-湍流频率-标量的联合PDF方程.本文对钝体火焰驻定器后冷态流场进行了计算,结果表明此混合算法能够捕捉流场中非稳态的漩涡脱落现象.着重研究了湍流频率模型系数的改变对漩涡脱落频率以及拟序结构在动量输运中的作用的影响.     

钝体后湍流预混燃烧的PDF模拟

本文采用PDF方法对矩形燃烧室内钝体后的湍流预混火焰进行了数值模拟。脉动速度-频率-标量联合的PDF输运方程用Monte Carlo方法求解,质量、动量和能量的平均值由基于无结构网格的有限体积法求解,压力通过状态方程获得。PDF方程中所需的平均密度、平均速度和压力由有限体积法提供,并将用Monte Carlo方法求出的雷诺应力、化学反应源项和比热比传递给有限体积法。本文对丙烷和空气燃烧的不同简化化学反应机理进行了研究,并与实验结果进行比较,获得满意的结果。     

近壁区湍流脉动压力相关项的理论模型

提出在湍流边界层近壁区采用三维波的理论模型描述湍流相干结构,根据理论模型计算了Ｒｅｙｎｏｌｄｓ应力输运方程中的脉动速度与脉动压力梯度的相关项,理论计算结果与直接数值模拟（ＤＮＳ）符合很好。表明该理论方法不仅有益于对湍流机制的了解,而且可能为湍流的近壁模型化开辟一条新的途径。     

湍流射流扩散火焰中NOx排放量的数值模拟

本文采用k-ε模型和求解标量的联合概率密度函数输运方程相结合的方法，利用文献[1]提出的简化NOx生成机理，对甲烷湍流射流扩散火焰中NOx的排放量进行了计算。计算结果表明，采用这一方法对污染物排放进行数值模拟是行之有效的，很有发展潜力。     

湍流扩散火焰的非稳态火焰面模拟

采用稳态的和非稳态的火焰面模型同时对一个湍流甲烷射流扩散火焰进行了数值模拟,比较了两者对湍流平均火焰结构、活性自由基和污染物(氮氧化物)排放的模拟效果。速度场采用κ-ε模型计算,守恒标量混合物分数的分布通过其概率密度函数(PDF)输运方程的求解得到。稳态的火焰面结构由查询火焰面数据库得到,而非稳态的火焰面结构由火焰面方程和流场方程耦合求解来计算。采用详细的GRI—Mech 3.0机理描述甲烷的氧化和氮氧化物的形成。数值模拟结果和实验数据作了广泛的对比,验证了火焰面模型对湍流扩散燃烧的定量模拟能力。     

Radiatively generated maximal mixing scenario for the Higgs boson mass and the least fine-tuned minimal supersymmetric standard model

We argue that given the experimental constraints on the Higgs boson mass the least fine-tuned parameter space of the minimal supersymmetric standard model is with negative top-squark masses squared at the grand unification scale. While the top-squark mass squared is typically driven to positive values at the weak scale, the contribution to the Higgs boson mass squared parameter from the running can be arbitrarily small, which reduces the fine-tuning of electroweak symmetry breaking. At the same time the top-squark mixing is necessarily enhanced and the maximal mixing scenario for the Higgs boson mass can be generated radiatively even when starting with negligible mixing at the unification scale. This highly alleviates constraints on possible models for supersymmetry breaking in which fine-tuning is absent.     

A priori assessment of the subgrid scale viscous/scalar dissipation closures in compressible turbulence

An appraisal is made of several subgrid scale (SGS) viscous/scalar dissipation closures via a priori analysis of direct numerical simulation data in a temporally evolving compressible mixing layer. The effects of the filter width, the compressibility level and the Schmidt number are studied for several models. Based on the scaling of SGS kinetic energy, a new formulation for SGS viscous dissipation is proposed. This yields the best overall prediction of the SGS viscous dissipation within the inertial subrange. An SGS scalar dissipation model based on the proportionality of the turbulent time scale with the scalar mixing time scale also performs the best for the filter widths in the inertial subrange. Two dynamic methods are implemented for the determination of the model coefficients. The one based on the global equilibrium of dissipation and production is shown to be more satisfactory than the conventional dynamic model.     

On the role of mixing models in the simulation of MILD combustion using finite-rate chemistry combustion models

The present work shows an in-depth analysis about the role of mixing models on the simulation of MILD combustion using a finite-rate combustion model, the Partially Stirred Reactor approach (PaSR). Different approaches of increasing complexity are compared: a simple model based on a fraction of the integral time scale, a fractal-based mixing model and a dynamic mixing model based on the resolution of transport equations for scalar variance and dissipation rate. The approach is validated using detailed experimental data from flames stabilized on the Adelaide Jet-in-Hot Co-flow (JHC) burner at different fuel-jet Reynolds numbers (5k, 10k and 20k) and different co-flow oxygen dilution levels (3%, 6% and 9%). The results indicate the major role of mixing models to correctly handle turbulence/chemistry interactions and clearly indicate the superior performances of the dynamic mixing model over the other tested approaches.     

An a priori assessment of the Partially Stirred Reactor (PaSR) model for MILD combustion

Moderate or Intense Low-oxygen Dilution (MILD) combustion has drawn increasing attention as it allows to avoid the thermo-chemical conditions prone to the formation of pollutant species while ensuring high energy efficiency and fuel flexibility. MILD combustion is characterized by a strong competition between turbulent mixing and chemical kinetics so that turbulence-chemistry interactions are naturally strengthened and unsteady phenomena such as local extinction and re-ignition may occur. The underlying physical mechanisms are not fully understood yet and the validation of combustion models featuring enhanced predictive capabilities is required. Within this context, high-fidelity data from Direct Numerical Simulation (DNS) represent a great opportunity for the assessment and the validation of combustion closure formulations. In this study, the performance of the Partially Stirred Reactor (PaSR) combustion model in MILD conditions is a priori assessed on Direct Numerical Simulations (DNS) of turbulent combustion of MILD mixtures in a cubical domain. Modeled quantities of interest, such as heat release rate and reaction rates of major and minor species, are compared to the corresponding filtered quantities extracted from the DNS. Different submodels for the key model parameters, i.e., the chemical time scale τ_c and the mixing time scale τ_mix, are considered and their influence on the results is evaluated. The results show that the mixing time scale is the leading scale in the investigated cases. The best agreement with the DNS data regarding the prediction of heat release rate and chemical source terms is achieved by the PaSR model that employs a local dynamic approach for the estimation of the mixing time scale. An overestimation of the OH species source terms occurs in limited zones of the computational domain, characterized by low heat release rates.     

A turbulence characteristic length scale for compressible flows

The current RANS models are generally established and calibrated under incompressible condition and these kinds of models could succeed in predicting many features of incompressible flows. However, these models extended to the high-speed, compressible flows are always less accurate. In the paper, a compressible von Kármán length scale is proposed for compressible flows considering the variable densities. It contains no empirical coefficients and is based on phenomenological theory. In the turbulent kinetic equation, the extra unclosed terms induced by non-constant densities are treated as dissipation terms and the equation is closed algebraically via the introduction of the von Kármán length scale. The original and the proposed von Kármán length scale lead to two different kinds of SAS (scale adaption simulation) models, KDO (turbulence kinetic energy dependent only) and CKDO (compressible KDO), respectively. Compressible mixing layer with significant compressibility is studied within standard k–?, k–ω, KDO turbulence models and their compressible versions. The compressibility effects such as the reduced mixing layer thickness, growth rate and turbulence intensity can be reproduced by CKDO. The new length scale can improve the performances of the model in predicting the mixing layer thickness, stream-wise velocity and Reynolds shear stresses when the convective Mach number is 0.8. Besides, the new length scale also leads to accurate computed growth rate when the convective Mach number ranges from 0.1 to 1.0.     

Diesel engine emissions and combustion predictions using advanced mixing models applicable to fuel sprays

An advanced mixing model was applied to study engine emissions and combustion with different injection strategies ranging from multiple injections, early injection and grouped-hole nozzle injection in light and heavy duty diesel engines. The model was implemented in the KIVA-CHEMKIN engine combustion code and simulations were conducted at different mesh resolutions. The model was compared with the standard KIVA spray model that uses the Lagrangian-Drop and Eulerian-Fluid (LDEF) approach, and a Gas Jet spray model that improves predictions of liquid sprays. A Vapor Particle Method (VPM) is introduced that accounts for sub-grid scale mixing of fuel vapor and more accurately and predicts the mixing of fuel-vapor over a range of mesh resolutions. The fuel vapor is transported as particles until a certain distance from nozzle is reached where the local jet half-width is adequately resolved by the local mesh scale. Within this distance the vapor particle is transported while releasing fuel vapor locally, as determined by a weighting factor. The VPM model more accurately predicts fuel-vapor penetrations for early cycle injections and flame lift-off lengths for late cycle injections. Engine combustion computations show that as compared to the standard KIVA and Gas Jet spray models, the VPM spray model improves predictions of in-cylinder pressure, heat released rate and engine emissions of NOx, CO and soot with coarse mesh resolutions. The VPM spray model is thus a good tool for efficiently investigating diesel engine combustion with practical mesh resolutions, thereby saving computer time.     

Numerical consistency check between two approaches to radiative corrections for neutrino masses and mixings

We briefly outline the two popular approaches on radiative corrections to neutrino masses and mixing angles, and then carry out a detailed numerical analysis for a consistency check between them in MSSM. We find that the two approaches are nearly consistent with a discrepancy factor of 4.2% with running vacuum expectation value (VEV) (13% for scale-independent VEV) in mass eigenvalues at low-energy scale but the predictions on mixing angles are almost consistent. We check the stability of the three types of neutrino models, i.e., hierarchical, inverted hierarchical and degenerate models, under radiative corrections, using both approaches, and find consistent conclusions. The neutrino mass models which are found to be stable under radiative corrections in MSSM are the normal hierarchical model and the inverted hierarchical model with opposite CP parity. We also carry out numerical analysis on some important conjectures related to radiative corrections in the MSSM, viz., radiative magnification of solar and atmospheric mixings in the case of nearly degenerate model having same CP parity (MPR conjecture) and radiative generation of solar mass scale in exactly two-fold degenerate model with opposite CP parity and non-zero U_e3 (JM conjecture). We observe certain exceptions to these conjectures. We find a new result that both solar mass scale and U_e3 can be generated through radiative corrections at low energy scale. Finally the effect of scaledependent vacuum expectation value in neutrino mass renormalisation is discussed     

煤粉炉喷钙脱硫过程的总体模型与模拟

1前言我国燃煤中约有三分之一用于火力发电,减少火电厂SO。排放是治理酸雨和SO。污染的重中之重。炉内喷钙脱硫结合尾部增湿活化技术（LIFAC）因流程简单、投资省、运行费用低而且效率较高,特别适用于燃烧中、低硫煤锅炉的改造,是一种具有广泛应用前景的技术［1］本文?..     

Viscous fluid cosmological models in LRS Bianchi type V universe with varyingΛ

Some LRS Bianchi type V viscous-fluid cosmological models are investigated, in which the coefficient of shear viscosity is considered as proportional to the scale of expansion in the model. This leads toA=Bn, whereA andB are metric potentials,n being a constant. The coefficient of bulk viscosity is also assumed to be a power function of mass density. The cosmological constant is found to be a decreasing function of time, which is supported by results from recent type Ia supernovae observations. Some physical aspects of the models are also discussed.     

Renormalisation group improved leptogenesis in family symmetry models

We study renormalisation group (RG) corrections relevant for leptogenesis in the case of family symmetry models such as the Altarelli–Feruglio A₄$A_4$ model of tri-bimaximal lepton mixing or its extension to tri-maximal mixing. Such corrections are particularly relevant since in large classes of family symmetry models, to leading order, the CP violating parameters of leptogenesis would be identically zero at the family symmetry breaking scale, due to the form dominance property. We find that RG corrections violate form dominance and enable such models to yield viable leptogenesis at the scale of right-handed neutrino masses. More generally, the results of this paper show that RG corrections to leptogenesis cannot be ignored for any family symmetry model involving sizeable neutrino and τ Yukawa couplings.