10液滴蒸发与燃烧.pptVIP

* * 气相能量 且 * * 液滴动量(速度) 且 * * 除了上面的方程外，还需要几个纯代数方程和C-H-O-N系统复杂平衡关系才能使方程组封闭。 要求解这个问题，上面的控制方程很容易使用一些子程序包[43,44]来进行积分。 * * 要计算气相的平衡属性，可以使用第2章提到的方法，也可以用一些已有的程序或子程序，例如参考文献[41]就可以达到这一目的。Olikara和Borman代码是使用起来很方便，它在计算平衡组分时可以计算不同的偏微分（例如(?h/?T),?h/??等。） * * 数学求解 给定一套气相和液相的初始条件，我们希望求解下面的函数： 气相： * * 液相： 接下来，我们会使用简单守恒定律（质量，动量和能量守恒）来求得上面列举函数。 * * 分析 对于这个问题，我们对稳态，稳定流动的控制体进行分析。所选择的控制体为沿燃烧器轴向长?x的一段（如图10.16a）。 质量守恒 图10.10b显示了进出控制体的液相和气相的质量流动。因为控制体中的质量没有变化，流入和流出控制体的总质量应该相等： * * 其中 和 分别是液体和气体的质量流率（kg/s）， 是单位长度进入控制体的二次风的质量流率（kg/s-m）。我们假设 是一个已知的关于x的函数。 整理上面的内容，可得： * * 取极限?x?0 ，结合微分的定义，我们可以得到下面的总质量守恒的方程： * * 方程10.85积分后可得 的表达式： 接下来我们只考虑液相；这样的话，如图10.17a所示的控制体就只包括液体。 是单位长度从液相进入气相的质量流率，即燃料蒸发速率。 * * * * 则液体燃料在x+?x处离开控制体的速率比在x处进入控制体的速率少 ?x ，即： 整理，并取?x ?0 ，得： * * 液体流过燃烧器的流率与单位时间进入其中的液滴数 及单个液滴的质量 md有关，则： 或 其中?l和D分别是液滴密度和直径。方程10.88微分后可得： * * D2对x的微分与对时间的微分可以通过液滴速度vd联系在一起，即dx=vddt ；则： * * 在前面我们曾得到过dD2/dt的关系式。将方程10.17代入方程10.90，可得： 其中 K 是蒸发系数(方程10.18). * * 单位时间进入燃烧器的液滴的数目很容易和初始燃料流率和假定的初始液滴大小D0联系在一起，即, 或 * * 上面的守恒分析得到了一个一元常微分方程，方程10.91，可以用来求解D2(x)或D(x)。知道了液滴的直径后， 可以通过方程10.86及一些其他条件来计算。 要得到液滴速度与坐标的关系vd(x)，需要应用液滴的动量守恒方程，下面将会讨论这些。气相的速度可以用下式表示： * * 其中?g是气相的密度， A是燃烧器的截面积。密度与气体温度和压力的关系可以通过理想气体方程得到： * * 其中 求解Tg需要气相能量守恒，而且，由于气相的成份在向下游流动时持续变化，它的分子量MWg也随之而变化。代入方程10.94和10.95到方程10.93之中，可得： * * 气相能量守恒 参见图10.17b，按照与分析质量守恒相同的步骤，控制体的能量守恒可以表达为： 其中我们假设没有热交换。展开上述微分，并认识到液体焓为常数，则方程10.97重新整理后可得： * * 由于我们想要得到温度分布T（x），所以将焓与温度和其他热力学变量之间建立关系： * * 其中各变量之间的关系可以由平衡来得到。应用链式定则，注意到为Pg常数， dhg/dx可以表达为： 其中下标g已经从T和?中去掉。 * * 由方程10.98和10.100的右半部分相等，来求解dT/dx ，可得： 如果知道了T ， P ， ?的平衡状态，气体焓hg ，和它的偏微分?hg/?T及?hg/?? ，都可以计算。Olikara和Borman[4]平衡模式（参见符录F）可以用来求解这些计算。 * * 方程10.101可以用方程10.85来消去 作进一步的简化，即 为了认清楚我们做了什么和要做什么，在这儿做一个简单的总结是很必要的。迄今为止，我们已经得到了dD2/dx和dTg/dx的一次常微分方程。 * * 给定初始条件，这些方程可以通过积分（数学上的）来得到D(x)和Tg(x) 。剩下的就是要找到求解? 和 d?/dx的关系式，以及再得到一个不同的可以积分求解液滴速度vd(x)的方程。 * * 气相成分 现在我们的目标是求当量比的轴向分布?(x) 。这一计算只是质量守恒的一个子集，因为所要算的只是任一位置的燃料和氧化剂的质量。 假设在入口处x=0有一