第3章守恒方程
全面理解物理守恒定律对成功分类、求解并应用其结论是很重要的。因此本章我们从质量、动量、热能和化学物质组分守恒的控制方程推导开始。推导过程基于流体力学的观点,读者会发现本章这些内容会与黏性流体力学教材中的内容有相当大的重叠。然而,我们会在两个方面与传统的表述有所不同。首先,因为主要考虑化学反应流,会保留很多在单**动中可以忽略的特性。其次,因为主要考虑轴对称流动,大部分数学表述是基于圆柱坐标系,而非笛卡儿坐标系。尽管这种选择增添了一些复杂性,但这能凸显那些会在笛卡儿坐标系中被忽视的重要问题。
我们的总���目标是处理欧拉坐标系下用空间坐标和时间作为自变量表述的偏微分方程形式的守恒方程。这种方法将用于系统的守恒定律的概念与固定在空间中有流体流过的控制体的行为联系起来。对于一个系统,即具有一定质量的流体,在其上能应用守恒定律,如质量守恒、动量守恒(F=ma)和能量(热力学**定律)守恒。但在实际情况下,追踪整个系统的流动和无数流体微团的相互作用是不可能的。幸运的是,如2.3节所讨论的,用随体导数的概念是可行的,它能定量地将系统的守恒定律与固定的控制体联系起来。
守恒定律通常基于很简单、很直观的概念。其表述为: 系统广延变量的积累速率等于通过系统边界面的净传递速率(流入减去流出)加上内部的净生成率(增加减去减少),即
dNdtsystem=(N·in-N·out)+N·gen(3.1)
对于一个含有固定质量的系统,必然不会有流体流进或流出系统。如果有的话,系统的质量就会改变,违反其定义。因此,N不能通过系统边界被“对流”。当然,我们相当清楚,流体流动能够传递质量、动量和能量,流体的这些行为必须被重视。
对于某一瞬时占据固定于空间的控制体上的系统,随体导数提供了系统和控制体之间必需的联系:
dNdtsystem=ρDηDtδV=(N·in-N·out)+N·gen(3.2)
其中,δV代表控制体的体积。应该强调一下,N·代表通过控制面的“非流动”的传递,例如分子扩散过程。通过控制体表面的对流传递明显包含在随体导数中。
一般来说,我们要找寻的微分方程呈以下形式:
ρDηDt=∑n·i(η)(3.3)
其中,n·=N·/δV代表每单位体积的N的输运率。我们尽可能地将n·(η)以η的显函数与空间导数形式来表达,即找到从η作因变量的偏微分方程。
3.1质量连续性
一定体积流体的质量守恒方程在任何流动情况都是必需的。当N是质量m时,相关的强度变量(每单位质量的广延变量)是η=1。那就是说,η是单位质量1。在这种情况下,控制体内没有质量的产生和消耗。例如,化学反应中,可能产生或消耗个别的化学组分,但整体上,没有质量产生或消耗。**能通过控制体表面传递净质量的方法是通过对流。尽管个别物质会因分子运动通过控制体表面扩散,但在此过程中也不会有净质量传递。这部分内容会在后面的质量传递章节中深入讨论。
我们注意到,两相流中对于某一给定相会有源项和汇项。例如,考虑水滴在潮湿的空气中的蒸发过程。对于液相和蒸汽相都能写出质量守恒方程。液体转变为蒸汽(反之亦然)会形成源项,即在某一相中有物质的产生(或消耗)。当然,在相与相之间,质量一定是守恒的,没有质量的净产生和净消耗。
由于质量连续方程已经在式(3.2)的随体导数形式的导数中使用了,它就对直接推导连续方程本身不再有直接作用。质量守恒方程的应用只能简单地得到一个不重要的等式。作为替代,我们从已经得出的积分形式的式(2.30)出发得到
dmdtsystem=∫CVρt+·ρVdV=0(3.4)
假设一个无限小的微元控制体,假设被积函数在控制体内是相同的,那么,连续性方程就能被写成微分方程形式,如:
ρt+·ρV=0(3.5)
连续性方程更常见的形式是引入密度的随体导数。使用随体导数的定义,表示为
DρDt=ρt+V·(ρ)=ρt+(V·)ρ(3.6)
以及乘积求微分的链式法则
·(ρV)=ρ·V+V·(ρ)=ρ·V+(V·)ρ(3.7)
连续性方程可被写成包含随体导数的形式,如
DρDt+ρ·V=0(3.8)
在上边的等式中,要注意到
V·(ρ)=def(V·)ρ(3.9)
是可交换使用的。等号右边的形式在流体力学文献中是*常见的,但左边的形式可能更容易理解。
对不可压缩流体,ρ是常数,根据式(3.8),很明显,速度场的散度为零,即
·V=0(3.10)
这在2.7节中讨论过,并且不可压缩流体也不会有净体积缩胀:
1VDVDt=·V(3.11)
尽管不可压缩流体很容易变形,但它的体积并不发生改变。
在圆柱坐标系中,用随体导数的定义,连续性方程可给定为
ρt+ρuz+1rrρvr+1rρwθ=0(3.12)
这是一个具有双曲线特性的一阶微分方程。
3.1.1流函数
对于稳态(不随时间变化)的二维流动,流函数的概念有很大的实用性。流函数的推导恰好满足连续性方程。在圆柱坐标中,有两种二维情况值得研究: rz平面,称为轴对称坐标,以及rθ平面,称为极坐标。
流函数的物理意义是流体沿着流线流动,流线即等流函数线。由于根据定义,流动不能穿过流线,因此任意两条流线之间的质量流量为常数。此外,两条流线之间流量的大小是由这两条流线上的流函数的数值不同决定的。