第5章粘性流体动力学基本方程组5.1粘性流体动力学基本方程流体运动所遵循的规律是由物理学三大守恒定律规定的,即质量守恒定律,动量守恒定律和能量守恒定律。这三大定律对流体运动的数学描述就是动力学基本方程组。但这个方程组是不封闭的,要使其封闭还需加上辅助的物性关系等。一般情况下,现在还求不出这个方程组的解析解,但研究这个方程组的性质却具有极其重要的意义,因为所有的流动现象都是由这个方程组所规定的。粘性流动的一个基本特征是流动的有旋性。因此研究涡的产生、输运和扩散就是很重要的了。这些性质也都是由流体动力学基本方程组所规定的。对流体运动的描述有两种方法,即拉格朗日法和欧拉法;对基本定理的数学表述也有两种方法,即积分形式和微分形式。本章将采用欧拉法和微分形式来表述基本方程。5.1.1质量守恒定律——连续方程连续方程是质量守恒定律在运动流体中的数学表达式。由于不涉及力的问题,因此粘性流体力学与非粘性流体方程完全相同,在非粘性流体中所做的推导和讨论在这里全部有效。考察流体通过一微元体的界面所引起的微元体内质量的变化问题。根据质量守恒定律,单位体积上通过微元体界面流出的质量流量即矢量u的散度u,它应等于微元体内单位时间单位体积所减少的质量:0tu(5.1.1)展开后得:0uvwtxyz(5.1.2)连续方程表示单位时间内流人流出微元体的质量必与密度变化相平衡。对于定常流,此式可变为:0uvwxyz(5.1.3)0u(5.1.4)对于不可压缩流,(5.1.2)式变为:0uvwxyz即iiux=0(5.1.5)由张量分析的知识可知,iiux是应变量张量的主对角线上三元素之和,恒为常数,表示微元体的体积变化率。式(5.1.5)表示总的体积变化率为零,与流体的不可压缩一致。-62-5.1.2动量守恒定律——运动方程粘性流体的运动方程是动量守恒定律对于粘性流体运动规律的数学表述,可由牛顿第二定律推出。以微元体为分析对象则可表述为:在惯性系中,流体微元体的质量和加速度的乘积等于该微元体所受外力的合力。对于流体运动应考虑两类外力:一为彻体力(用F来表示),它是作用在微元体上所有质量上的力,如重力;另一类为表面力(用P来表示),它是作用在微元体界面上的力,如压力、摩擦力等。运动方程可写成如下向量形式:DDtuFP(5.1.6)其中微分符号D()Diiutttxu(5.1.7)称为物质导数或随体导数,它所代表的是微团的某性质对时间的变化率。例如,DDtu是该微团的速度u随时间的变化率,即加速度,亦即D()Dttuuuu(5.1.8)从欧拉法的观点看,此式右端第一项由流动的非定常性引起,称为当地加速度;右端第二项由流场中速度分布的不均匀性引起,表示经Dt时间后由于微团空间位置的变化而引起的速度的变化,称为迁移加速度。式(5.1.6)中的彻体力F可表示成:xyzFFFFijk(5.1.9)在这里彻体力可以看成是已知的外力,而表面力则和流体速度场的变形情况有关。它决定了流体的应力状态。所以我们分别研究流体的应力和应变变化率后,将建立它们之间的关系。为了写出表面力的式子,我们从流体中取出正六面微元体(图5-1)。它的左下方的点的坐标为(x,y,z)。对于垂直于x轴的两个微元面上分别作用了如下的合应力(应力即单位面积上的作用力):xP和dxxxxPP图5-1微元体的应力张量这里的注足x表示x方向上的应力向量,则作用在垂直于x轴的微元面上的应力的合力为:(d)dddddddxxxxxyzyzxyzxxPPPP(5.1.10)同样可得作用在垂直于y轴和z轴的微元面上的应力的合力分别为:dddyxyzyP,dddzxyzzP于是可得作用于单位容积的表面力的合力向量为-63-(ddddddddd)/dddyxzyxzxyzxyzxyzxyzxyzxyzPPPPPPP(5.1.11)式中xP,yP和zP都是向量,还可以把它们沿三个坐标方向分解,即分解为正应力和平行于各微元面的切应力。例如,作用于与x轴垂直的微元面上的应力xP可分解为(图5-1):xxxyxzPijk(5.1.12a)同理有yy...