高等流体力学之三个特解

1、一、涡的衰减情况推导一、涡的衰减情况推导 由涡旋的传输方程知道，当流体具有粘性、由涡旋的传输方程知道，当流体具有粘性、非正压或者质量力无势时，均将破坏涡旋的守非正压或者质量力无势时，均将破坏涡旋的守恒。粘性、非正压与质量力无势这三者中，尤恒。粘性、非正压与质量力无势这三者中，尤以粘性流体为经常性起作用的因素。因为，对以粘性流体为经常性起作用的因素。因为，对于实际的流体，运动时总是呈现粘性。因此，于实际的流体，运动时总是呈现粘性。因此，粘性流体一般来讲是有旋的，而且其涡旋的大粘性流体一般来讲是有旋的，而且其涡旋的大小可以随时间产生、发展、衰减、消失。涡旋小可以随时间产生、发展、衰减、消失。涡旋还

2、会扩散，自涡旋强度大得地方向涡旋强度弱还会扩散，自涡旋强度大得地方向涡旋强度弱的地方扩散，直至涡旋强度均衡为止。涡旋强的地方扩散，直至涡旋强度均衡为止。涡旋强度的扩散性质决定了很多流体运动的物理现象；度的扩散性质决定了很多流体运动的物理现象；因此，研究涡旋在粘性流体中运动的规律具有因此，研究涡旋在粘性流体中运动的规律具有重要的实际意义。重要的实际意义。 下面以一空间孤立涡线为例，从涡旋传输方程下面以一空间孤立涡线为例，从涡旋传输方程具体分析涡旋扩散的规律。具体分析涡旋扩散的规律。 设在无边界的粘性流体中有一强度设在无边界的粘性流体中有一强度 为的无穷长为的无穷长直涡线。不难证实，此涡线引起的运

3、动与非粘性直涡线。不难证实，此涡线引起的运动与非粘性流体情形相似，运动是无旋的，涡线周围各处流体情形相似，运动是无旋的，涡线周围各处 的的 ，流体质点以，流体质点以 的速度作定长圆周运的速度作定长圆周运 动。差别只在于，在非粘性流体中，由于没有粘动。差别只在于，在非粘性流体中，由于没有粘性内摩擦阻力，因此，该直涡线的强度能够永远性内摩擦阻力，因此，该直涡线的强度能够永远保持不变，且不会向周围流体中扩散，不需要外保持不变，且不会向周围流体中扩散，不需要外加能量来维持质点的定常圆周运动；但在粘性流加能量来维持质点的定常圆周运动；但在粘性流体中，由于粘性的缘故，涡旋将衰减下去，要维体中，由于粘性的缘

4、故，涡旋将衰减下去，要维持这种运动，就必须有外加能量，例如用涡旋的持这种运动，就必须有外加能量，例如用涡旋的无穷长细柱体来供给涡源。无穷长细柱体来供给涡源。00 02Vr 现在讨论的问题是，假定在某现在讨论的问题是，假定在某 时刻，外加涡时刻，外加涡源突然中断，分析该直涡线的扩散（衰减）情况。源突然中断，分析该直涡线的扩散（衰减）情况。 根据质量力有势的不可压缩粘性流体的涡旋传输根据质量力有势的不可压缩粘性流体的涡旋传输方程方程 上述方程也可以写成下面的形式：上述方程也可以写成下面的形式： (1-1) 如图如图3-2所示，沿直涡线去所示，沿直涡线去oz轴，则有轴，则有()()VVt 0 xy

5、zk 0t DDtVv 由于运动的对称性和平面运动中速度由于运动的对称性和平面运动中速度V沿沿 方向方向的微商为零，故的微商为零，故 于是，（于是，（1-1）式可改写成）式可改写成 略去下标略去下标z，写成，写成 (1-2)z()0V ()0Vzzt t 采用极坐标，上式可写成采用极坐标，上式可写成 (1-3) 初始条件为：初始条件为： (1-4a) 边界条件为：边界条件为： (1-4b)()rtrrt0t 0r 0 0t r 0 解方程（解方程（1-3）可得）可得 (1-5) 式中，常数式中，常数A可用沿周线的速度环量等于该圆周所可用沿周线的速度环量等于该圆周所围得涡管强度这一条件来确定。设

6、在任意时刻沿围得涡管强度这一条件来确定。设在任意时刻沿半径为半径为r的圆周上的速度环量为的圆周上的速度环量为 ，同一时刻，同一时刻，半径为半径为r的涡管强度应该为的涡管强度应该为24rtAet 02rrdr 因此可得因此可得 (1-6) 将将 ， 代入上式，可得代入上式，可得224400224(1)rrrrttArdrerdrAet 0t 0 04A 于是得到涡量分布为于是得到涡量分布为 (1-7) 速度环量为速度环量为 (1-8) 由于由于 (1-9)2044rtet 240(1)rte 202CV dsV rdrV n 把（把（1-9）式代入上式，得到速度分布为）式代入上式，得到速度分布为

7、 (1-10) 由（由（1-7）式可以看到，在初始时刻）式可以看到，在初始时刻t=0，各处（，各处（r0）的运）的运 动都是无旋的，在动都是无旋的，在t0的任何时刻整个空间立即产生涡旋，的任何时刻整个空间立即产生涡旋， 分布情况可由（分布情况可由（1-7）式代表即涡旋随距离）式代表即涡旋随距离r的增加而单调地的增加而单调地 下降。在中心处（下降。在中心处（r=0）的涡旋随时间增长而单调地下降。）的涡旋随时间增长而单调地下降。 而在离中心某一距离处（而在离中心某一距离处（r=a）的点上，涡旋起初增加，达）的点上，涡旋起初增加，达 到一个极大值以后，开始降低，一直降到当到一个极大值以后，开始降低，

8、一直降到当 时的时的 值。在值。在r=a处，处， 随随t的变化规律如图的变化规律如图3-3所示。所示。204(1)4rtVet 同样根据（同样根据（1-10）式可以画出不通时刻速度随）式可以画出不通时刻速度随r的的变化规律。如图变化规律。如图3-4所示所示. 最后，还可以注意到，在任意最后，还可以注意到，在任意r处，当处，当时时 ， 。换句话说，初始时刻由于无。换句话说，初始时刻由于无 源涡线在粘性流体中引起的运动，随时间过程而衰源涡线在粘性流体中引起的运动，随时间过程而衰 减下去，直至运动停止，所有涡旋的动能都耗散变减下去，直至运动停止，所有涡旋的动能都耗散变 为热。相反，对于任意时刻为热。

9、相反，对于任意时刻t，当，当 时，时， ， ，也就是说，运动在直涡，也就是说，运动在直涡 线处逐渐消失。线处逐渐消失。t 00V t 0 0V 二、旋转圆盘附近的流动二、旋转圆盘附近的流动 Von-Karman研究了旋转圆盘附近的流动，研究了旋转圆盘附近的流动，并用并用N-S方程做了解析解。如图所示，假定方程做了解析解。如图所示，假定半径为无限大的平面圆盘在不可压缩流体半径为无限大的平面圆盘在不可压缩流体中以等加速度旋转，忽略质量力。如图中以等加速度旋转，忽略质量力。如图3-11所示，由于粘性，圆盘带动圆盘附近的所示，由于粘性，圆盘带动圆盘附近的流体旋转，离心力的作用使流体产生径向流体旋转，离

10、心力的作用使流体产生径向分速，压力下降。为补充径向分速流出的分速，压力下降。为补充径向分速流出的流体，自然出现轴向分速，最终形成轴对流体，自然出现轴向分速，最终形成轴对称螺旋形流动。称螺旋形流动。 在惯性圆柱坐标系中，速度分量为在惯性圆柱坐标系中，速度分量为uz，u2 ，u流动为轴对称，故流动各量不随周向角流动为轴对称，故流动各量不随周向角变化。基本微分方程为变化。基本微分方程为()()10rzruurrz（2-1a） 222221()rrrrrrzuuuuuupuurzrrrrrz （2-1b） 2222()rrzuu uuuuuuurrzrrrz（2-1c） 222211()zzzzzrz

11、uuuuupuurzzrrrz （2-1d） 边界条件边界条件 方程方程2-1中，中， ur 、 u和和uz是二阶微商，要有七个是二阶微商，要有七个边界条件才能确定。边界条件不够，因此在求解边界条件才能确定。边界条件不够，因此在求解之前必须对流动做进一步的分析，作出合理的假之前必须对流动做进一步的分析，作出合理的假定。定。 显然，决定流动速度和鸭梨分布的因素是圆盘旋显然，决定流动速度和鸭梨分布的因素是圆盘旋转速度转速度，流体粘性系数，流体粘性系数及空间点的坐标及空间点的坐标r和和z，故可得速度和压力为故可得速度和压力为0z 0zruuurz 0ruu（2-2） ( , , , )iiuur z

12、 （2-3） ( , , , )ppr z 下面估算圆盘旋转带动的粘性层厚度下面估算圆盘旋转带动的粘性层厚度。因为距。因为距旋转轴距离为旋转轴距离为r的粘性层流体单位体积所受的离心的粘性层流体单位体积所受的离心力为力为r2，所以在底面积为，所以在底面积为drds、高、高为粘性层为粘性层厚度的流体元上，所受的离心力为厚度的流体元上，所受的离心力为r2 drds 。同一流体元还受到圆盘面上切应力同一流体元还受到圆盘面上切应力w的作用，这的作用，这个力与流体滑动方向相反，且与周向速度成一角个力与流体滑动方向相反，且与周向速度成一角度，比如度，比如。该流体元的切应力的径向分量必与。该流体元的切应力的径

13、向分量必与离心力平衡。离心力平衡。 或或2sinwdrdsrdrds 2sinwr （2-4） 另一方面，切应力的周向分量必须正比于壁面上另一方面，切应力的周向分量必须正比于壁面上周向速度的轴向梯度，其数量级（以周向速度的轴向梯度，其数量级（以“”表示数表示数量级）关系为量级）关系为 从以上两个方程中消去从以上两个方程中消去 ，得，得 因为圆盘半径无限大，所以紧贴壁面流体滑动的因为圆盘半径无限大，所以紧贴壁面流体滑动的方向与半径方向与半径r无关，故圆盘带动的流体厚度为无关，故圆盘带动的流体厚度为coswr （2-5） w2tan（2-6） 将上式带入式（将上式带入式（2-4），得圆盘上得摩擦应

14、力为），得圆盘上得摩擦应力为 为了求解方程（为了求解方程（2-1），根据切向应力引起），根据切向应力引起u、离、离心力引起心力引起ur和圆盘无限大的流动特点，和圆盘无限大的流动特点， Von-Karman假设假设 wr （2-7） ( , , )rurfz ( , , )urgz ( , , )zuhz 0( , , )pppz （2-8） 对上式应用对上式应用定理，选定理，选、为量纲独立量，得为量纲独立量，得 可见可见 是无量自变量是无量自变量 ，用用表示为表示为 因此，式（因此，式（2-8）可写成）可写成1rur 2ur 3zu 04pp 5/z ，5z（2-9） ( )ruFr( )uG

15、r( )zuH0( )ppP， （2-10） 将式将式（2-9）和（）和（2-10）分别代入方程（）分别代入方程（3-1）各式，）各式，因为因为 则方程（则方程（3-1）各式变为）各式变为zz2222z220FF HFG20GHGFG220PFFH（2-11） 边界条件为边界条件为00FHP1G 0FGH； （2-12） 先从（先从（2-11）的前三式求）的前三式求F、G、H，然后代入第，然后代入第四式求四式求P。方程（。方程（2-11）是非线性常务分方程，不）是非线性常务分方程，不易求解，可用数值积分或幂级数法。易求解，可用数值积分或幂级数法。Karman用用数值积分法求的近似解，以后数值积

16、分法求的近似解，以后Cocbran得到更精得到更精确的数值结果，计算结果表示于图确的数值结果，计算结果表示于图3-12和表和表3-1中。中。 计算表明，当计算表明，当很小时，随很小时，随增大，增大，F和和G迅速趋于迅速趋于零，具有边界层特性。零，具有边界层特性。ur和和u只在圆盘附近明显，只在圆盘附近明显，在粘性影响范围内，压力变化的量级为在粘性影响范围内，压力变化的量级为。流。流体从远处吸向圆盘，并向周围抛出，旋转圆盘相体从远处吸向圆盘，并向周围抛出，旋转圆盘相当于离心泵的作用。各函数的边界值为当于离心泵的作用。各函数的边界值为00F 0.510F 1G 0.616G 50.01G ：，：

17、（相当于边界层外边界） 0.8838H ： 虽然上述结果是假设圆盘无限大的条件下导出的，虽然上述结果是假设圆盘无限大的条件下导出的，但只要圆盘半径比边界层厚度打得多，可忽略圆但只要圆盘半径比边界层厚度打得多，可忽略圆盘周边的影响，仍可用上述结果计算有限大圆盘盘周边的影响，仍可用上述结果计算有限大圆盘的摩擦力矩。壁面摩擦应力为的摩擦力矩。壁面摩擦应力为3122(0)()(0)0.616zwrGrz 3122(0)()0.51rzrwrz （2-13） 半径为半径为R的圆盘两面所受的力矩和力矩系数分别为的圆盘两面所受的力矩和力矩系数分别为312422024(0)1.938RzMr drR 2223

18、.87/ 2ReMMCR（2-14） 式中式中2Re/R 半径为半径为R圆盘上上不离心力甩出的体积流量为圆盘上上不离心力甩出的体积流量为202()0.8838rRqRudzR（2-15） 该流量等于轴向流入的流量。该流量等于轴向流入的流量。 图图3-13所示，当所示，当 时，理论结果与实验很时，理论结果与实验很吻合，说明在层流条件下，上述推论完全正确；吻合，说明在层流条件下，上述推论完全正确；当当 时，理论偏离实验，说明流动变为时，理论偏离实验，说明流动变为湍流，湍流摩擦阻力大。按湍流速度分布的七分湍流，湍流摩擦阻力大。按湍流速度分布的七分之一次方规律计算，旋转圆盘的力矩系数为之一次方规律计算

19、，旋转圆盘的力矩系数为5Re3 10 5Re3 10 150.146ReMC（2-16） 三、绕圆球流动情况三、绕圆球流动情况 图图3-16表示小雷诺数绕圆球流动。在球坐标表示小雷诺数绕圆球流动。在球坐标系中，由于流动的对称性，故横向速度为零，系中，由于流动的对称性，故横向速度为零，其余两个速度分量为其余两个速度分量为 和和 ，压力为，压力为 ( , )ru r( , )ur( , )p r 因而基本方程组为因而基本方程组为 （3-1）2211()(sin )0sinrr uurrr2222222cot211120()(sin)sinrrruuuuuprrrrrrrrr 2222211120(

20、)(sin)sinsinruuuuprrrrrrrr 边界条件边界条件 （3-2） 根据边界条件和流动对称性，可用分离变量法求解：根据边界条件和流动对称性，可用分离变量法求解： rR0ruur cosruVsinuV 设方程的解具有以下形式设方程的解具有以下形式 （3-3） 将式（将式（3-3）代入方程（）代入方程（3-1）各式，整理后）各式，整理后得得 （3-4）cos)(sin)(cos)(u321rrfpprfurf2211 13211)(420)(2rfffrffrfff 边界条件边界条件 （3-5） 方程（方程（3-4）是线性常微分方程组，可先由）是线性常微分方程组，可先由第一式和第

21、三式求出第一式和第三式求出 和和 与与 的关系，的关系，再将和代入第二式而得再将和代入第二式而得082121,:0:fVffrffRr08881 12 13 14frfrfrfr2f3f1f 解上式得解上式得 ，再求得，再求得 和和 ，即，即 （3-6） 根据边界条件，式（根据边界条件，式（3-6）中的各积分常数）中的各积分常数为为 1f2f3fDrrAfDrCrBrAfDrCrBrAf1022223232231VRA321RVB23VC0D 于是：按式（于是：按式（3-3），速度和压力分别为），速度和压力分别为 （3-7）根据应力关系式，球面上的法向应力和切向应力根据应力关系式，球面上的法向

22、应力和切向应力分别为分别为)2231 (cos33rRrRVur3331sin (1)44RRuVrr 23cos2V Rppr()(2)rrrr Rr Rupr 1()()rrr Rr Ruuurrr 将速度关系式（将速度关系式（3-7）代入，得）代入，得 （3-8） 法向应力和切向应力对球产生的阻力为法向应力和切向应力对球产生的阻力为 （3-9）31()cos2rrr Rr RVppR 31()sin2rr RVR 20(sincos )2sin6xrrrDRdRV 阻力系数为阻力系数为 （3-10） 式（式（3-10）称为）称为Stokes公式公式 从图从图3-17看出，在子午平面内看出

23、，在子午平面内Stokes解得流线前解得流线前后对称。但从图后对称。但从图3-18看出，正应力不对称，前大后看出，正应力不对称，前大后小，形成压差阻力；切向应力前后按正弦规律变化，小，形成压差阻力；切向应力前后按正弦规律变化，也产生了阻力。也产生了阻力。2224/ 2ReDdDCVR2RedV R 讨论：讨论： （1）Stokes全部忽略惯性项的近似不适用于整个全部忽略惯性项的近似不适用于整个流场，因为惯性项与粘性项之比的量级在不同径向流场，因为惯性项与粘性项之比的量级在不同径向位置处有显著差别。小位置处有显著差别。小Re数时，惯性项与粘性项数时，惯性项与粘性项之比，只在物面附近很小，在远离物面附近处却不之比，只在物面附近很小，在远离物面附近处却不是这样。故是这样。故Stokes近似只适用于小近似只适用于小Re数空间绕流数空间绕流的物面附近的流场，是一级近似；对于平面流动，的物面附近的流场，是一级近似；对于平面流动，Stokes近似根本得不到解。近似根本得不到解。 （2）和对速度分布没有影响，他们的影响只在阻）和对速度分布没有影响，他们的影响只在阻力系数中得