第四章：对流环境中的非等温液滴蒸发模型与数值求解.docVIP

第4章 对流环境中的非等温液滴蒸发模型与数值求解 目前常见的液滴蒸发研究多集中于等温液滴蒸发[18,36,37]，即假设液滴内部导热率无限大。事实上，燃料液滴的导热率只比气体稍大，液滴从气液界面吸收的热量不会立即传到液滴内部，特别是当液滴在强对流环境中蒸发时，毕渥数（）相对增大，液滴表面温度比其内部升高更快，液滴内部的非等温性将对液滴蒸发过程产生影响。近年来，国外一些学者先后开展了非等温液滴蒸发模型的研究[9,10,30]，但对非等温液滴蒸发过程中的环境压力效应与膨胀效应缺乏较深入的研究。 等温液滴蒸发模型因为不考虑液滴内部的温度梯度，又称为零维蒸发模型；常见的非等温液滴蒸发模型因为考虑了液滴内部沿径向的一维导热，又称为一维蒸发模型。 液滴蒸发过程中客观存在着瞬态加热和平衡蒸发两个阶段[32]，在瞬态加热过程中，液滴从环境获得的热量主要用于液滴加热，液滴升温快，蒸发慢，因此，在液滴蒸发初期必然存在着由于受热而产生的膨胀现象，而目前许多关于非等温液滴蒸发的研究忽略了液滴蒸发过程中的膨胀现象[9,18,30]。 本章将在前文等温液滴蒸发研究的基础上，通过采用沿液滴内部径向瞬态导热的非等温液滴蒸发模型，分析非等温液滴蒸发数值求解方法与网格处理，并对计算方法与程序的可靠性给予验证。 4.1 物理模型与控制方程 单组分球对称液滴蒸发气体在液体中的扩散速率与液滴蒸发速率相比很小，认为气体在气液界面上的径向量为零；气相准稳态，满足Clausius-Clapeyron方程 (4-1) 考虑液滴内部温度梯度时，液滴表面能量守衡方程： (4-2) 液滴表面的传质关系式[18,47]： (4-3) 液滴内部一维瞬态导热微分方程： (4-4) 初始条件： ，， 边界条件： ， ， 4.2 离散与求解 非等温液滴蒸发是传热、传质、气体流动及液滴内部的导热通过液滴表面的质量守恒和能量守恒的耦合实现的，因此，非等温液滴蒸发求解时需要对有关方程进行离散，进而耦合求解。 4.2.1 方程离散 (1) 液滴表面质量守恒方程的离散 由式(4-1)可得 (4-5) 整理得 (4-6) 式(4-6)右边第一项为单位时间因液滴蒸发引起的半径减小，第二项为单位时间因液滴密度减小引起的半径变化，因此，液滴半径的变化是由质量的减小与密度的减小共同引起的。 将式(4-6)离散为如下形式： (4-7) (2) 液滴的区域离散 如图4.1所示，采用内结点法、均匀网格对液滴区域进行离散。图中结点所在的控制容积的两个边界分别为m、n，液滴中心到这两个边界的距离分别为rm和rn，节点i+1与i-1为结点的两个相邻结点。 (3) 液滴内部导热方程的离散 采用控制容积法和隐式格式离散液滴内部非稳态导热微分方程。在结点的控制容积内对式(4-4)从到积分： (4-8) 整理式(4-8)得： (4-9) 令  则式(4-9)可表示为： (4-10) (4) 液滴边界离散 液滴加热蒸发时，为获得边界温度，可以采用直接离散法或辅助源项法。 (a) 直接离散法 对式(4-2)隐式离散得： (4-11) 式中： —— 次外层节点温度，； —— 最外层节点与次外层结点之间距离，。 整理(4-11)得： (4-12) (b) 辅助源项法[49]： 液滴边界区域离散示意图如图4.2。s、p、e为液滴外层三个连续结点，rs、re分别为液滴中心到最外层控制容积两个边界的距离。将式(4-10)运用于边界处： (4-13) 式中： 为了在的代数方程中不出现未知的边界条件，需要利用已知的边界条件把去掉，为此，须对式(4-13)作如下变换： (4-14) 根据边界条件和Fourier定律： (4-15) 式中： —— 进入边界控制容积的热流密度，，规定进入为正，对于液滴中心点，。 将式(4-15)代入式(4-14)可得 (4-16) 式中： 求出p点温度后，可根据Fourier定律反求s点温度。本文采用辅助源项法。 4.2.2 蒸发求解 液滴内部的一维导热与外部的传热、传质通过液滴表面的质量、能量守恒方程实现耦合。因此，在完成液滴内部一维导热离散与液滴表面质量、能量守恒方程离散后，还需要分析液滴外部的传热、传质，并在此基础上完成整个液滴蒸发求解。 (1) 液滴外部的传