FVMC混合方法求解RANSPDF方程研究.pdf

中国工程热物理学会 燃烧学 学术会议论文 编号：084205 FV/MC混合方法求解RANS/PDF方程 朱旻明，陈义良，韩省思，叶桃红 （中国科学技术大学热科学和能源工程系， 合肥 230027） （Tel: 0551 －3601650， Email: zhmm@ustc.edu.cn ） 摘要：本文在非结构网格中，采用有限体积法(FV)/Monte Carlo 法(MC)相结合的混合算法求解湍流燃 烧的 PDF 方程。在分析混合算法数值误差特性的基础上，提出了 FV 法和 MC 法采用不同网格划分的 双层网格方法。分析了双层网格技术减小数值误差的原理，讲述了双层网格系统中颗粒场和平均场的 数据交换。并以湍流射流扩散火焰为算例，说明在混合算法中采用双层网格划分，与单一网格划分相 比，能够极大的提高计算效率和精度。 关键词：概率密度函数；Monte Carlo ；有限体积法；非结构网格；双层网格 0 引言 概率密度函数(PDF)方法通过求解流场中单点联合概率密度函数方程模拟湍流燃 烧，它能够精确处理对流项和化学反应源项，是研究燃烧过程中污染物排放、燃烧装置 点火、熄火等问题的有力工具[[1]] 。概率密度函数方法的主要缺点是计算量大，首先， Monte Carlo方法求解概率密度函数方程需要大量的样本；其次，虽然原则上可以通过直 接积分详细的化学反应速率方程组得到化学热力学参数的概率密度函数，但是详细的化 学反应机理方程组刚性问题严重，求解困难且计算量大，目前用于PDF方法中的化学反 应计算方法主要有简化机理法，与火焰面模型相结合法和自适应建表方法等。Pope等人 提出将有限体积 (FV)和Monte Carlo(MC)相结合的混合算法求解湍流燃烧的概率密度函 数方程，并在结构网格中发展了此算法[[2]] 。在该算法中，用MC方法求解脉动速度-频率 -标量联合PDF方程，给出湍流输运项（如雷诺应力，热通量）及化学反应源项。将概率 密度函数方程积分得到的平均连续、速度和能量方程则用FV方法求解，给出平均压强、 速度和密度。与以往的混合算法相比，它满足了模型层次上的相容性，且合理的耦合方 [[2]] 式减少了偏差，使计算精度和效率都显著优于单独的颗粒方法 。 为了面向工程应用，适应工程实际中复杂的几何结构，我们将该算法用于非结构网 格，研制了适用于二维湍流燃烧问题的计算程序，并取得了一定的成功[[3]] 。混合算法的 复杂性导致了数值误差受多种因素的影响，研究其误差特性对于提高精度和减小计算量 有重要意义。本文根据混合算法的数值误差特性，提出了有限体积法和Monte Carlo法采 用不同网格划分的双层网格方法，并通过理论分析和数值实验说明了在混合算法中采用 双层网格能够进一步提高计算效率。 本文首先介绍了混合算法，接着分析混合算法的数值误差特性，第 3节提出双层网 格划分技术，第 4 节对湍流射流扩散火焰作了数值模拟，最后给出结论。 资助项目：国家自然科学基金 (No. [[2]] 1 求解概率密度函数方程的相容的混合算法 1.1 求解的方程及耦合关系 u ω Φ 的联合PDF 定义为 流场中单点密度加权平均的脉动速度 －频率 －标量