粘性流体力学精确解.docVIP

粘性流体力学（朱克勤许春晓） 第三章粘性流体运动方程的精确解

PAGE1

PAGE53

第三章粘性流体运动方程的精确解

自从建立了以方程为核心的粘性流体运动方程组，人们就开始致力于寻求在各种定解条件下粘性流体运动方程的精确解。精确解对于深入认识和分析流体运动规律具有重要意义，它也为检验各类数值方法的可靠性和精确度提供重要的依据，作为基本流场的精确解也是流动稳定性分析的出发点。

方程是一个非线性的二阶偏微分方程，在许多情况下，求粘性流体运动的精确解是一件非常困难的工作。迄今为止，可以得到精确解的流动例子非常有限关于方程精确解的

关于方程精确解的详细述评可见参考文献[2]。

本章首先介绍平行剪切流，由于这时方程的非线性项为零，寻求这一类流动的精确解，数学上处理起来比较容易。然后将讨论包含有非线性项的方程精确解的几个著名例子。

图3.1平行平板间的剪切流

3.1平行平板间的定常剪切流

考虑相距的两块无限大平行平板，下板静止，上板沿水平方向以匀速运动。平板间的不可压缩粘性流体在恒定的压力梯度和运动平板的作用下沿方向作定常运动。

在图3.1所示的直角坐标系中，流体运动速度只有沿方向的分量，并且只是坐标的函数。这里流线相互平行，称为平行剪切流。连续方程自动满足，由于非线性项为零，方程(2.5.11)在方向的投影为

(3.1.1)

边界条件

， (3.1.2)

将常微分方程(3.1.1)积分两次，并利用边界条件确定积分常数可得到速度分布

(3.1.3)

其中第一部分是由于压力梯度引起的，为抛物线分布，与流体粘性系数有关；第二部分是由于平板运动引起的，为线性分布，与流体粘性系数无关。

根据速度分布可以求出通过单位宽度平板间的体积流率

(3.1.4)

平均流速是

(3.1.5)

相应的切应力

(3.1.6)

当时，切应力的最小值和最大值分别在上板面和下板面处出现。

在粘性流体力学中，将单纯由压力梯度引起的流动称为Poiseuille流动，比如管道中的流动；将由运动固壁引起的流动称为Couette流动，比如本节讨论的平行壁间的剪切流。Couette流动允许压力梯度存在，没有压力梯度的Couette流动称为简单Couette流动。

图3.2平行板间的分层剪切流

在以上问题中，如果没有压力梯度，平板间是两层互不掺混的粘性流体，厚度分别为和，粘性系数分别为和，且。这时动量方程仍是方程(3.1.1)，求解时以流体的界面为界将流场分成两个区域，积分后共有四个积分常数需要确定，除了原有的两个固壁边界条件外，在界面处增加速度和切应力连续的两个条件，便可解得速度分布

(3.1.7)

粘性大的流体层速度梯度较小。应该特别指出的是，在流体的界面处，两侧流场的涡量是不连续的。

3.2同轴圆筒间的定常流

3.2.1同轴旋转圆筒间的Couette流

图3.3同轴圆筒间的流动

考虑两个同轴圆筒间的粘性流体，内筒的外径为，外筒的内径为。圆筒作匀速旋转，角速度分别为和。旋转式圆筒型粘度计内的流动可归入这一类流动（见图3.3）。

假定圆筒足够长，可以忽略圆筒底部壁面的影响。在筒壁带动下流体作轴对称定常运动，引入柱坐标系。本问题的流场中，流体运动速度只有沿方向的一个分量，速度只与坐标有关。这时，连续方程自动满足，方程的非线性项为零，在方向上的投影式(2.5.15b)可简化为

(3.2.1)

满足筒壁边界条件

(3.2.2)

(3.2.3)

的速度分布为

(3.2.4)

由(2.5.18)得到涡量场

(3.2.5)

它的特点是与无关。在两种情况下涡量场为零，其一是内外筒半径和旋转角速度满足关系时，这是涡量场变号的临界情况，这时的速度分布为。其二是在的同时让，这时速度分布退化为，这相当于一个涡核半径为的位势涡；根据1.2.3节的分析，其涡量为零是由于流场中曲率涡量和切变涡量正好相互抵消的结果。

进一步计算切应力

(3.2.6)

这表明内筒壁面上受到的摩擦切应力较大，外筒壁面上受到的摩擦切应力较小。但是在单位长度的圆筒面上，内外筒面受到的摩擦力矩是相同的，等于

(3.2