Improvement of the pressure algorithm of the stand-alone PDF method to treat unsteady three-dimensional turbulent reacting flows

This paper deals with the particle-mesh probability density function (PDF) method. It shows how an existing but less precise pressure algorithm for the stand-alone method can be improved. The present algorithm is able to handle the general case of an unsteady three-dimensional turbulent reacting flow. The transport equation of the joint PDF of velocity and composition is solved with a particle method. Open boundary conditions are realized and for statistical reasons a simple but effective particle splitting procedure is applied.

Based on a simple configuration, the properties of the presented improved pressure algorithm are analysed. It is shown which numerical condition must be taken care of so that the algorithm is able to correct the particle positions such that the normalization condition is fulfilled as accurately as specified.

To verify the algorithm the combustion of a methane–air mixture enclosed in an open simulation volume is calculated. It is shown that the simple particle splitting algorithm works very effectively in the studied case. The behaviour of the improved pressure algorithm is examined by different calculations. To analyse the convergence of the algorithm, the particle number per cell and the grid spacing are varied. To demonstrate the accuracy, a statistically stationary inflow/outflow configuration is used and the numerical solution is compared to an analytical one. For a less symmetric test case, the previous unsteady combustion problem is simulated, including an additional mean velocity in one direction.

The presented improved pressure algorithm provides the opportunity to calculate unsteady three-dimensional turbulent reacting flows with a stand-alone method, and offers an alternative to the complex hybrid finite-volume/particle PDF method.     

A LES/PDF simulator on block-structured meshes

A block-structured mesh large-eddy simulation (LES)/probability density function (PDF) simulator is developed within the OpenFOAM framework for computational modelling of complex turbulent reacting flows. The LES/PDF solver is a hybrid solution methodology consisting of (i) a finite-volume (FV) method for solving the filtered mass and momentum equations (LES solver), and (ii) a Lagrangian particle-based Monte Carlo algorithm (PDF solver) for solving the modelled transport equation of the filtered joint PDF of compositions. Both the LES and the PDF methods are developed and combined to form a hybrid LES/PDF simulator entirely within the OpenFOAM framework. The in situ adaptive tabulation method [S.B. Pope, Computationally efficient implementation of combustion chemistry using in situ adaptive tabulation, Combust. Theory Model. 1 (1997), pp. 41–63; L. Lu, S.R. Lantz, Z. Ren, and B.S. Pope, Computationally efficient implementation of combustion chemistry in parallel PDF calculations, J. Comput. Phys. 228 (2009), pp. 5490–5525] is incorporated into the new LES/PDF solver for efficient computations of combustion chemistry with detailed reaction kinetics. The method is designed to utilise a block-structured mesh and can readily be extended to unstructured grids. The three-stage velocity interpolation method of Zhang and Haworth [A general mass consistency algorithm for hybrid particle/finite-volume PDF methods, J. Comput. Phys. 194 (2004), pp. 156–193] is adapted to interpolate the LES velocity field onto particle locations accurately and to enforce the consistency between LES and PDF fields at the numerical solution level. The hybrid algorithm is fully parallelised using the conventional domain decomposition approach. A detailed examination of the effects of each stage and the overall performance of the velocity interpolation algorithm is performed. Accurate coupling of the LES and PDF solvers is demonstrated using the one-way coupling methodology. Then the fully two-way coupled LES/PDF solver is successfully applied to simulate the Sandia Flame-D, and a turbulent non-swirling premixed flame and a turbulent swirling stratified flame from the Cambridge turbulent stratified flame series [M.S. Sweeney, S. Hochgreb, M.J. Dunn, and R.S. Barlow, The structure of turbulent stratified and premixed methane/air flames I: Non-swirling flows, Combust. Flame 159 (2012), pp. 2896–2911; M.S. Sweeney, S. Hochgreb, M.J. Dunn, and R.S. Barlow, The structure of turbulent stratified and premixed methane/air flames II: Swirling flows, Combust. Flame 159 (2012), pp. 2912–2929]. It is found that the LES/PDF method is very robust and the results are in good agreement with the experimental data for both flames.     

Comparisons of Lagrangian and Eulerian PDF methods in simulations of non-premixed turbulent jet flames with moderate-to-strong turbulence-chemistry interactions

Transported probability density function (PDF) methods have been applied widely and effectively for modelling turbulent reacting flows. In most applications of PDF methods to date, Lagrangian particle Monte Carlo algorithms have been used to solve a modelled PDF transport equation. However, Lagrangian particle PDF methods are computationally intensive and are not readily integrated into conventional Eulerian computational fluid dynamics (CFD) codes. Eulerian field PDF methods have been proposed as an alternative. Here a systematic comparison is performed among three methods for solving the same underlying modelled composition PDF transport equation: a consistent hybrid Lagrangian particle/Eulerian mesh (LPEM) method, a stochastic Eulerian field (SEF) method and a deterministic Eulerian field method with a direct-quadrature-method-of-moments closure (a multi-environment PDF-MEPDF method). The comparisons have been made in simulations of a series of three non-premixed, piloted methane–air turbulent jet flames that exhibit progressively increasing levels of local extinction and turbulence-chemistry interactions: Sandia/TUD flames D, E and F. The three PDF methods have been implemented using the same underlying CFD solver, and results obtained using the three methods have been compared using (to the extent possible) equivalent physical models and numerical parameters. Reasonably converged mean and rms scalar profiles are obtained using 40 particles per cell for the LPEM method or 40 Eulerian fields for the SEF method. Results from these stochastic methods are compared with results obtained using two- and three-environment MEPDF methods. The relative advantages and disadvantages of each method in terms of accuracy and computational requirements are explored and identified. In general, the results obtained from the two stochastic methods (LPEM and SEF) are very similar, and are in closer agreement with experimental measurements than those obtained using the MEPDF method, while MEPDF is the most computationally efficient of the three methods. These and other findings are discussed in detail.     

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

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

A coupled pressure-based computational method for incompressible/compressible flows

Pressure-based flow solvers couple continuity and linearized truncated momentum equations to derive a Poisson type pressure correction equation and use the well known SIMPLE algorithm. Momentum equations and the pressure correction equation are typically solved sequentially. In many cases this method results in slow and often difficult convergence. The current paper proposes a novel computational algorithm, solving for pressure and velocity simultaneously within a pressure-correction coupled solution approach using finite volume method on structured and unstructured meshes. The method can be applied to both incompressible and subsonic compressible flows. For subsonic compressible flows, the energy equation is also coupled with flow field and the density of fluid is obtained by equation of state. The procedure eliminates the pressure correction step, the most expensive component of the SIMPLE-like algorithms. The proposed coupled continuity-momentum-energy equation method can be used to simulate steady state or transient flow problems. The method has been tested on several CFD benchmark cases with excellent results showing dramatically improved numerical convergence and significant reduction in computational time.     

两相壁面射流PDF方程有限分析方法及数值模拟

提出求解位置-速度相空间中高维两相流PDF(probability density function)方程的有限分析方法,将位置-速度相空间颗粒PDF方程约化到速度空间,并解析求解,颗粒的位置PDF用轨道方法求解.对壁面射流两相流动进行数值模拟,并与颗粒雷诺应力轨道方法进行比较计算,结果优于颗粒雷诺应力轨道方法.     

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

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

改进的多重网格法重建含遮拦物的干涉波前

阐述了由单张干涉图求取相位主值图像及最小二乘相位恢复。推导出基于梯度拟合的泊松方程及其离散化形成后,详细讨论了多重网格法对此方程求解的原理和过程。结合路径无关与路径相关算法的优点,设计了自适应最优路径法作为多重网格法的预处理,改进迭代的初始条件,提高了重建精度,同时大大加快了收敛速度。实验结果证明改进的多重网格法对含遮挡物的干涉图能获得很好的波前重建效果。     

Probabilistic solution of nonlinear oscillators excited by combined Gaussian and Poisson white noises

The stationary probability density function (PDF) solution of the stochastic response of nonlinear oscillators is investigated in this paper. The external excitation is assumed to be a combination of Gaussian and Poisson white noises. The PDF solution is governed by the generalized Kolmogorov equation which is solved by the exponential-polynomial closure (EPC) method. In order to evaluate the effectiveness of the EPC method, different nonlinear oscillators are considered in numerical analysis. Nonlinearity exists either in displacement or in velocity for these nonlinear oscillators. The impulse arrival rate, mono-modal PDF and bi-modal PDF are also considered in this study. Compared to the PDF given by Monte Carlo simulation, the EPC method presents good agreement with the simulated result, which can also be observed in the tail region of the PDF solution.     

The probabilistic solution of stochastic oscillators with even nonlinearity under poisson excitation

The probabilistic solutions of nonlinear stochastic oscillators with even nonlinearity driven by Poisson white noise are investigated in this paper. The stationary probability density function (PDF) of the oscillator responses governed by the reduced Fokker-Planck-Kolmogorov equation is obtained with exponentialpolynomial closure (EPC) method. Different types of nonlinear oscillators are considered. Monte Carlo simulation is conducted to examine the effectiveness and accuracy of the EPC method in this case. It is found that the PDF solutions obtained with EPC agree well with those obtained with Monte Carlo simulation, especially in the tail regions of the PDFs of oscillator responses. Numerical analysis shows that the mean of displacement is nonzero and the PDF of displacement is nonsymmetric about its mean when there is even nonlinearity in displacement in the oscillator. Numerical analysis further shows that the mean of velocity always equals zero and the PDF of velocity is symmetrically distributed about its mean.     

Multi-environment probability density function method for modelling turbulent combustion using realistic chemical kinetics

A computational fluid dynamics (CFD) tool for performing turbulent combustion simulations that require finite-rate chemistry is developed and tested by modelling a series of bluff-body stabilized flames that exhibit different levels of finite-rate chemistry effects ranging from near equilibrium to near global extinction. The new modelling tool is based on the multi-environment probability density function (MEPDF) methodology and combines the following: the direct quadrature method of moments (DQMOM); the interaction-by-exchange-with-the-mean (IEM) mixing model; and realistic combustion chemistry. Using DQMOM, the MEPDF model can be derived from the transport PDF equation by depicting the joint composition PDF as a weighted summation of a finite number of multi-dimensional Dirac delta functions in the composition space. The MEPDF method with multiple reactive scalars retains the unique property of the joint PDF method of treating chemical reactions exactly. However, unlike the joint PDF methods that typically must resort to particle-based Monte-Carlo solution schemes, the MEPDF equations (i.e. the transport equations of the weighted delta-peaks) can be solved by traditional Eulerian grid-based techniques. In the current study, a pseudo time-splitting scheme is adopted to solve the MEPDF equations; the reaction source terms are computed with a highly efficient and accurate in-situ adaptive tabulation (ISAT) algorithm. A 19-species reduced mechanism based on quasi-steady state assumptions is used in the simulations of the bluff-body flames. The modelling results are compared with the experimental data, including mixing, temperature, major species and important minor species such as CO and NO. Compared with simulations using a Monte-Carlo joint PDF method, the new approach shows comparable accuracy.     

Time-averaging strategies in the finite-volume/particle hybrid algorithm for the joint PDF equation of turbulent reactive flows

The influence of time-averaging on bias is investigated in the finite-volume/particle hybrid algorithm for the joint PDF equation for statistically-stationary turbulent reactive flows. It is found that the time-averaging of the mean fluctuating velocity (TAu) leads to the same variances of the fluctuating velocity before and after the velocity correction, whereas without TAu the estimates are different, and an additional numerical dissipation rate is introduced for the turbulent kinetic energy (TKE). When 100 particles per cell are used without TAu, a large bias error is found to be involved in the unconditional statistics of the statistically-stationary solutions of two tested turbulent flames, the Cabra H₂/N₂ lifted flame and the Sandia piloted flame E. The use of TAu reduces this bias dramatically for the same number of particles per cell. The conditional statistics in these flames, however, are hardly affected by TAu. To a large extent, the effect of the bias error on the unconditional statistics is similar to the effect of increasing the model constant C _{ω 1} in the stochastic turbulence frequency model.     

A quadrature-based LES/transported probability density function approach for modeling supersonic combustion

The joint-scalar probability density function (PDF) approach provides a comprehensive framework for large eddy simulation (LES) based combustion modeling. However, currently available stochastic approaches for solving the high-dimensional PDF transport equation can be error prone and numerically unstable in highly compressible shock-containing flows. In this work, a novel Eulerian approach called the direct quadrature method of moments (DQMOM) is developed for evolving the PDF-based supersonic combustion model. The DQMOM technique uses a set of scalar transport equations with specific source terms to recover the PDF. The new technique is coupled to a compressible LES solver through the energy equation. The DQMOM approach is then used to simulate two practical flow configurations: a supersonic reacting jet and a cavity-stabilized supersonic combustor. Comparisons with experimental data demonstrate the predictive accuracy of the method.     

湍流扩散燃烧的数值研究-PDF方法和火焰面模型的性能比较

采用标量概率密度函数(PDF)方法、稳态和非稳态火焰面模型三种方法对一个值班湍流CH_4/O_2/N_2射流扩散火焰(Sandia Flame D)进行数值计算,以比较不同燃烧模型的性能。PDF方法通过计算反应标量的PDF输运方程来得到标量分布,而火焰面模型只求解单标量混合物分数的PDF方程,组分和温度分布通过火焰面方程的求解或者火焰面数据库的插值得到。计算结果和实验数据对比表明PDF方法计算结果最好但计算量相当大,稳态火焰面模型则反之。综合而言,非稳态火焰面模型的预测结果相对稳态模型有了非常大的改进,而计算量仍然容易接受,非常适合工程应用。     

各向异性颗粒PDF输运方程及其颗粒湍流扩散

根据Minier的模型,获得颗粒所见流体脉动速度的朗之万模型,并推导颗粒所见流体湍流脉动速度自关联函数,该函数具有各向异性的特征,进一步由流体脉动速度朗之万模型和颗粒运动方程获得相应的颗粒PDF输运方程.为描述在网格生成湍流中的颗粒扩散运动,将颗粒PDF输运方程作相应的简化,获得一个具有颗粒湍流各向异性的解析解,将方程的解与Wellsand Stock的颗粒湍流扩散实验结果进行比较,获得了很好的结果.揭示了具有各向异性特征的颗粒PDF输运方程的优越性.     

高精度有限差分在湍流直接数值模拟中的应用

本文采用高精度有限差分法对矩形管道内充分发展湍流换热进行了直接数值模拟,湍流雷诺数Rer和普朗特数分别为400(Rem=6200)和0.71,压力泊松方程分别采用两种不同的离散格式(二阶和四阶中心差分离散).结果表明,同二阶中心差分格式相比,高精度有限差分可以在较少的网格下得到较好的结果;压力泊松方程采用四阶中心差分或二阶中心差分对计算结果的影响甚小,但采用二阶中心差分时,可以节省大量的计算时间.     

两相湍流色噪声扩维法PDF模型及数值模拟

由颗粒运动的朗之万方程出发,对流体脉动速度采用扩维方法,得到两个不同层次的PDF输运方程.通过对颗粒运动方程求解和高斯分布假设,解决PDF方程的封闭问题,获得颗粒二阶矩模型,然后将颗粒应力方程简化成代数方程,建立代数应力模型.将对流扩散方程的有限分析法运用到求解两相流模型中,对壁面两相射流进行数值模拟,对比分析数值结果与实验结果.     

轴对称钝体后湍流扩散燃烧的PDF模拟

在无结构网格中，对轴对称钝体驻定的湍流扩散火焰进行数值模拟．采用有限容积(FV)／Monte Carlo(MC)混合算法求解湍流燃烧问题的混合算法．Monte Carlo法求解脉动速度一标量．频率的联合概率密度函数方程，有限容积法求解平均质量、动量和能量方程．求解的两组方程是相容的，合理的耦合方式可以减少统计偏差，计算精度和效率显著优于单独的颗粒方法．文中对化学反应采用层流火焰面模型，并将数值计算结果与实验结果作了比较和分析．     

湍流旋流燃烧瞬时速度PDF的实验测量

本文应用三维激光粒子动态分析仪(PDA)对分级进风旋流燃烧室内的湍流燃烧过程进行了实验测量。通过瞬时速度场的测量结果,获得了燃烧室内各测点上气体轴向与切向瞬时速度的概率密度函数(PDF),揭示了热态速度场的湍流脉动特性。     

A Eulerian PDF scheme for LES of nonpremixed turbulent combustion with second-order accurate mixture fraction

Monte Carlo simulations of joint probability density function (PDF) approaches have been developed in the past largely with Reynolds averaged Navier Stokes (RANS) applications. Current interests are in the extension of PDF approaches to large eddy simulation (LES). As LES resolves accurately the large scales of turbulence in time, the Monte Carlo simulation and the flow field need to be tightly coupled. A tight coupling can be achieved if the consistency between the scalar field solution obtained via finite-volume (FV) methods and that from the stochastic solution of the PDF is ensured. For nonpremixed turbulent flames with two distinct streams, the local reactive mixture is described by the mixture fraction. A Eulerian Monte Carlo method is developed to achieve a second-order accuracy in the instantaneous filtered mixture fraction that is consistent with the corresponding FV. The performances of the proposed scheme are extensively evaluated using a one-dimensional model. Then, the scheme is applied to two cases with LES. The first one is a non-reacting mixing flow of two different fluids. The second case is the Sandia piloted turbulent flame D with a steady state flamelet model. Both results confirm the consistency of the proposed method to the level of filtered mixture fraction.