燃烧温度计算程序_液体推进剂的燃烧（液滴的蒸发）

液体火箭发动机使用的液体推进剂包括液体燃料和液体氧化剂。对于一般的液体燃料，其燃烧形式包括：

1.液面燃烧。首先液面蒸发（扩散），然后在液面形成燃烧火焰。由于液态蒸发时间较长，故为扩散燃烧。

2.预蒸发燃烧。液体先蒸发为气体，在喷射到燃烧室中进行燃烧。如果没有预先混合氧化剂，则为扩散燃烧，如打火机气体的燃烧；如果预先与氧化剂混合，则为预混燃烧，如家用液化气的燃烧。

3.雾化燃烧。液面通过雾化器形成小液滴，然后再进行燃烧。由于雾化液滴增大了燃烧面积，从而强化了燃烧。

液体燃料还可以通过灯芯的方式实现燃烧。利用灯芯的毛细作用使得燃油再灯芯表面形成燃料蒸气，然后和空气混合燃烧。这种燃烧方式功率小，只适合小功率的燃烧器或家庭生活所用，如蜡烛、煤油灯、酒精灯等。

液体火箭发动机中的燃烧方式是雾化燃烧，故主要讨论这种燃烧方式。雾化燃烧的基本单元是液滴，把单个液滴的燃烧规律按统计规律叠加在一起，便形成总的雾化燃烧。因此，液滴的燃烧是基础。

液滴的蒸发

液滴的燃烧首先要经过液滴的蒸发过程。液滴的蒸发分为低温蒸发和高温蒸发，低温蒸发主要是由浓度梯度引起的，类似于斯蒂芬流，而高温蒸发的驱动力不仅仅是浓度，还有温度梯度。这里以单个小液滴在静止条件下的蒸发为例阐述液滴高温蒸发的特点和规律

液滴的蒸发模型与控制方程

研究液滴蒸发的主要目的是确定液滴的蒸发速率，定义为液滴单位表面积上向外蒸发的液体流量，用

表示。由质量流率公式

，蒸发速率

为

^[1]（蒸发速率定义式）

液滴蒸发的基本模型为：

1.液滴为规律球形，如下图所示，这时常采用球坐标系统。

液滴蒸发模型

2.液体蒸气的流动假设为一维、准定常、无黏、不可压。

3.液滴的蒸发动力只考虑热传导，而忽略其他因素如热辐射等，并忽略所有外力。

4.各物性参数均为常数，且路易斯数

。

严格上讲，液滴由于不断蒸发其大小在不断变化，应为非定常，但蒸发速率较慢，可假设为准定常；同样由于蒸发流动速率较慢，常假设为不可压，即

，同时可忽略静温和总温的差异。一维球坐标控制方程可从积分形式的变截面一维定常流动控制方程中推导，这里直接给出结果。由于忽略所有外力，因此没有动量方程。

混合物连续方程由

或

^[2]可得为

（连续方程式1）

式中：球表面积

。对于液滴表面，有

，R为液滴表面半径，代入上式，得

由蒸发速率定义式，结合上式可知

（连续方程式2）

组分连续方程由

,忽略化学反应，可得为

代入连续方程式1，并考虑球表面积

和连续方程式2，得

（组分连续方程式）

能量方程只考虑热传导对蒸发的作用，由能量方程式，忽略静温和总压的差异，可得

考虑球表面积

和连续方程式，则

^[3]（能量方程式）

比较液滴蒸发控制方程中的组分连续方程式和能量方程式，可见它们在形式上相似。

液滴蒸发速率的解

由于液滴蒸发控制方程具有相似性，只需求其中一个方程，另一个方程则可采用类比方法得到，引入无量纲守恒标量

（守恒标量定义式）

式中：

为环境温度；q为单位质量液滴从初始温度

蒸发为蒸汽所需的能量；

为环境液滴蒸气质量百分比。各参数的变化特征如下图所示。

液滴蒸发参数变化

用守恒标量代入组分连续方程式和能量方程式，得

（式1）

（式2）

由路易斯数定义式

,以及假设

，可知

可见式1和式2具有相同的形式，其解的形式也应相似，只需求解其中的一个方程即可。以能量方程为例，对式2积分一次，得

(a)

式中：

为积分常数，可利用下述液滴表面处的边界层条件确定，即

时

式中：

可由液滴表面的热平衡条件求得——该条件认为，液滴表面蒸发需吸收的能量

应等于周围介质传给液滴的热量

,即

考虑连续方程式和球表面公式

，则

(b)

边界条件式（b）代入积分式（a），可得积分常数

。代入式（a），得一次积分结果为

分离变量

积分得

其中

为积分常数，利用下述液滴环境边界条件确定，即

可得

，则

由于蒸发速率是在液滴表面，故取

，代入上式，整理得

令

（传质数定义式）

则

（求解的结果式）

类似地，可对组分连续方程式进行积分。这时需要考虑扩散物质质量守恒边界条件，即

，即

按照上述同样的积分过程，可得蒸发速率为

（式3）

其中

（式4）

比较求解的结果式和式3，由于两个形式相同，可知

，故通用B代替，则液滴的蒸发速率为

（液滴的蒸发速率式）

式中：

称为传质数，或传质驱动力，又称斯柏尔丁（D.B.Spalding）数。

分析液滴的蒸发速率式可知，影响液滴蒸发速率的因素主要由三方面：一是液滴本身的物性参数，如

；二为液滴的大小即R,液滴越小，蒸发速率越快，反之，则越慢；三即为传质数B。由于传质数在对数项中，因此，传质数对蒸发速率的影响相对较弱。

在液滴蒸发与燃烧 计算中，除了需要确定传质数B外，还需要确定蒸气的物性参数如

等（假设为常数）所对应的温度。罗（C.K.Law）和威廉姆斯（F.A.Williams）提出了简单的近似方法——取燃料沸点温度

和环境为温度

的平均值为特性温度即

则对应的物性参数为

其中导热系数考虑了液滴燃料导热系数

和燃料环境物质导热系数

的综合影响

传质数的计算

由传质数定义式和式4，并考虑守恒标量的定义式，由于

故

液滴蒸发所需的热量q包括两部分：一是将液滴从初始温度

升高到

所需的能量；二是液滴在

下蒸发为蒸气所需的汽化潜热（或称汽化焓）

，即

^[4]

由于

，传质数B的计算公式可统一表示为

（传质数式）

由于液滴的表面参数如

和

未知，故传质数的计算还需要考虑更多的模型，一般根据环境温度

相对于沸点

的相对大小来选择不同的计算方法。

1.当环境温度远高于液滴沸点，即

时，

。由于

，则

，即传质速率无限快，液滴瞬间雾化，相对传质数可由

直接计算。燃烧室正常燃烧即为这种情况。

2.在液滴的环境温度

相对于沸点

不是很高的情况下，可假设

，从液滴蒸发参数变化图温度变化中可以看出，这种假设带来的误差不大。由传质数式，将

代入，可得传质数为

3.当环境温度远低于液滴沸点，即

时，

，则

，相应传质数可由

计算，但表面质量分数

还未知，需要采用液滴蒸发表面满足饱和蒸气热力学平衡条件。燃烧室低温点火燃烧即为这种情况。

4.一般情况下，可利用液滴蒸发表面满足饱和蒸气热力学平衡的条件来求解液滴表面参数。

假设液滴蒸发表面满足饱和蒸气热力学平衡条件，根据克劳修斯-柯拉柏隆方程，液滴的饱和蒸气压

为

^[5]（克劳修斯-柯拉柏隆方程式）

一些液体燃料的饱和蒸气压于温度的对应关系见下表，利用表中数据可以拟合得到克劳修斯-柯拉柏隆方程中的系数

。

一些液体燃料的饱和蒸气压于温度的对应关系

设液滴的饱和蒸气混合物（包括饱和蒸气于其他混合物）的总密度为

，则

式中：

为饱和蒸气的密度，

为除开饱和蒸气意外的混合物密度。设混合蒸气为完全气体，且其温度均匀，则理想气体热状态方程代入上式可得

或

pM=

式中：p为饱和蒸气混合物总的压强，M为摩尔质量

考虑到总压于分压的关系：

，则

在表面w处，克劳修斯-柯拉柏隆方程式代入上式得

把上式代入传质数式，得到关于

的方程，可以解出相应的表面温度

，再代回传质数式从而得到该温度下的传质数。

参考

1. ^下标W代表液滴表面
2. ^这是一维定常流动的连续方程，即在一每一个横截面上的质量流率dot{m}都是相等的。
3. ^c_p为液体燃料蒸气的比定压热容
4. ^cd为液滴的比热容
5. ^c1、c2为液滴的物性常数