液雾燃烧的全尺度直接数值模拟方法.docVIP

*基金项目：国家自然科学基金（Nos. 中国工程热物理学会 燃烧学 学术会议论文 编号： 液雾燃烧的全尺度直接数值模拟方法 邵长孝，柴敏，王海鸥，罗坤，樊建人 （浙江大学能源清洁利用国家重点实验室，杭州，310027） （ Tel Email: shaocx@zju.edu.cn） 摘 要 过去的研究表明，雾化燃烧的模拟精度受雾化模型影响较大。本文发展了一种适用于液雾燃烧的全尺度直接数值模拟方法，直接求解雾化、蒸发和燃烧过程，其模拟结果的可靠性在不同的算例中都得到了检验。本文方法为下一步液雾湍流燃烧的高精度全尺度直接数值模拟研究提供了数值基础。 关键词 直接数值模拟；全尺度；雾化；蒸发；燃烧 0 前言 雾化燃烧广泛存在于航空发动机燃烧室等推进领域。雾化燃烧是一个多相、多尺度、多物理场耦合的复杂过程，其中涉及到液体燃料的雾化、蒸发、燃料蒸汽/空气混合、燃料蒸汽的燃烧等。雾化燃烧的复杂性对实验测量技术提出了严峻的挑战。随着计算方法的发展和计算能力的提升，数值模拟正在成为研究雾化燃烧过程的重要手段。 传统的雾化燃烧直接数值模拟方法假定液体燃料完全雾化（入口采用液滴粒径分布简化）或采用大量参数控制的雾化模型。例如，罗坤等[1]采用点源假设的液滴替代模型燃烧室中液体燃料，Irannejad等[2]和Esclapez等[3]采用随机雾化模型来模拟液体燃料的雾化燃烧过程。然而，雾化过程的近似处理会对雾化燃烧的模拟结果产生影响，例如Greenberg[4]发现不同的液滴粒径分布会对火焰锋面产生重要影响，Som和Aggarwal [5]研究发现不同的雾化模型对燃烧特性产生不同的影响。因此，为了减小粒径分布假设和雾化模型精度对燃烧特性的影响，有必要发展一种可以同时考虑雾化、蒸发和燃烧过程的全尺度直接数值模拟方法。 近些年来，高精度的界面追踪方法被用来直接模拟液体燃料的雾化过程，例如Level Set方法[6]，Volume of Fluid方法[7]，Front Tracking方法[8]等。但是，同时对雾化、蒸发和燃烧过程进行全尺度模拟的方法并未见报道。对于该全尺度方法的实现，有以下两个难点：（1）界面处传热传质阶跃条件对浓度方程、温度方程和动量方程离散的处理，（2）耦合燃烧过程的处理。 本文的目标是解决以上难点，提出一种雾化、蒸发和燃烧耦合的全尺度直接数值模拟方法，其中气液界面由最近发展的Level Set方法[9]追踪，界面处的阶跃条件由虚拟流体法（Ghost Fluid Method，GFM）实现，气相燃烧由正庚烷的单步化学反应机理来描述。本文的全尺度雾化燃烧直接数值模拟方法将为下一步雾化湍流燃烧的研究提供高可信度的基础。 1 数值模拟方法 1.1 控制方程 本文采用变密度、低马赫数近似的Navier-Stokes方程组： , (1) , (2) 其中，、、p、和S分别表示密度、速度、压力、动力粘度和应力张量。 本文还求解燃料（F），氧化剂（O）和产物（P）质量分数的标量输运方程： , (3) , (4) , (5) 其中，和分别为质量扩散系数和化学反应源项。 最后需要求解整个计算区域（包含气相和液相区域）的温度输运方程： , (6) 其中，为比热容，为导热率，为燃烧化学反应源项。 1.2 Level Set方法 本文使用Level Set方法处理界面。该方法可将界面隐式地表示为某连续的Level Set方程的零等值面。通常情况下被设为符号距离函数，即： , (7) 其中是界面上离x最近的位置。Level Set方程为正值时代表液相，为负值时代表气相。 然而，上述标准Level Set方法存在质量不守恒的问题。本文采用双曲正切函数来减小质量损失： ,