内燃机喷雾欧拉-拉格朗日数值研究.docxVIP

基金项目：国家自然科学基金资助项目. 中国工程热物理学会 燃烧学 学术会议论文 编号：184083 内燃机喷雾欧拉-拉格朗日数值研究 李亮1,解茂昭1,贾明1 （大连理工大学能源与动力工程学院,大连116024） （TelEmail:liliang2009223@163.com） 摘要：目前，关于内燃机喷雾的研究通常采用拉格朗日模型；然而，欧拉模型则相对较少。本文则在构建的欧拉-拉格朗日模型下，对ECN喷雾实验进行研究。结果显示：在150Mpa喷射压力下，喷嘴附近存在明显的连续欧拉液相。然而，在液核的下游，离散项液滴的破碎、碰撞、融合及聚集等均对界面面积密度、粒子直径的分布产生明显影响。定量分析知，在LES湍流模型下，除轴距6mm附近的狭窄区间燃料体积分数存在较大偏差外，其对整个喷雾的轴向和径向分布均做到了很好的预测。 关键字：欧拉-拉格朗日；喷雾；LES；界面面积密度；多相流 0 引言 众所周知，喷雾的初次雾化发生在液核周边。在这一区域，流体以连续相的状态射入，气液两相之间存在强烈的相互作用，而并非离散方法（DDM）中所描述的以众多液体包的形式射入。而且，实验中在该区域很少发现球形的液滴，而只能观测到一个气液界面为高度褶皱的连续相流体。因此，严格的说，DDM框架下的多种亚模型[1-3]与真实喷雾存在较大差别，且在数值建立过程中其对入射粒子量、喷射锥角等参数的设定具有很强的人为经验性。比较可行的方法是采用连续欧拉法对喷雾近场进行研究，同时把气液界面的面积密度定义为一个函数，并导出精确的输运方程。通过输运方程的求解来描述气液界面和喷雾下游离散液滴的行为。Vallet[4]正是基于这种思想，首次提出了液体射流的ELSA（Eulerian-Lagrangian spray and atomization）雾化模型；此模型主要针对的是高雷诺数和高韦伯数的流体计算，并逐步演化出多种应用模型。 Blokkel et al.[5]采用Vallet提出的方法对整个喷雾过程进行了数值研究，模拟了由近场连续相到远场离散项的转变，即假定两个液滴间的平均自由间距达到平均液滴直径的2倍时发生模型的转变，最终得出了非常满意的结果，而且当燃烧室气体密度增加时模型同样可以很精确的预测喷雾锥角的变化。Beau[6]同样采用了均匀模型，并将其命名为ELSA（Eulerian-Lagrangian Spray Atomization）模型，作者假设在液相体积分数小于0.1时发生欧拉到拉格朗日的转换。Lebas[7, 8]对比ELSA模型和DNS模型的结果来确定模型的常数和参数值，并加入了蒸发模型，考虑了切应力、破碎、碰撞、汽化等因素对界面密度输运方程的影响,同时研究了不同环境气体密度、温度、射流压力等条件下ELSA模型适用性。Ning[9]在Vallet基础上进一步增加了平衡蒸发模型，并同时修正了标准κ-ε湍流模型，以考虑可压缩性和涡管拉伸的影响。 Belhadef[10]所构建的方法中，通过修改界面密度输运方程中的生成和耗散项，利用欧拉均相模型对轴对称旋流喷射进行模拟，结果显示索特平均直径与试验值很好的匹配。García-Oliver[11, 12]研究和验证了蒸发和无蒸发条件下可压缩Σ?Yliq模型，通过修改标准κ-ε湍流输运方程中的耗散系数，可以精确地预测不同入射压力和环境气体密度下液体和喷雾的贯穿距、轴向速度和混合质量分数分布。而且，随着环境气体密度和入射压力的增加，预测精度也相应增加。 从先前的研究可见，欧拉—拉格朗日模型在喷雾中的应用相对较少。本文在先前基于OpenFOAM所构建的欧拉模型[13]的基础上，参照Vallet理论和前人的研究，构建欧拉-拉格朗日模型。在Vallet的最初模型中，压力是通过状态方程直接求解[4]；然而，在我们当前构建的模型中，我们则采用PISO算法来迭代求解压力泊松方程；由于该模型假设液相向气体环境的弥散主要是由于湍流扩散控制；即，湍流模型的选择将对计算产生影响；因此，本文对动量方程分别进行了RANS和LES封闭。最后，应用此模型对柴油喷射特性进行数值研究，模拟对象为ECN（Engine combustion Network） 喷雾A[14, 15]2106753#试验，此喷雾已广泛应用于发动机各种喷雾模型的验证。 1 理论模型 Vallet[16]在1998年提出了一种新的欧拉-拉格朗日计算模型，此模型的基本出发点是：高雷诺数和