第4章流体动力学基础

1、第四章 流体动力学基础本章在流体运动学的基础上，加进动力学因素，对运动流体的应力状态作进一步分析，定义应力张量，并给出应力张量和变形率张量之间的联系。建立不可压缩流体运动微分方程 N-S方程。对理想流体运动微分方程 欧拉方程在恒定条件下沿流线积分得到恒定元流的能量方程 伯努利方程，进而推广到总流，得到恒定总流的能量方程。将动量守恒定律用于恒定总流得到恒定总流的动量方程。41运动流体的应力状态l 在静止流体里，无论是理想还是粘性流体，流体质点只能承受压应力，即静水压强。任一点上的静水压强与作用方向无关，只是位置的函数。这说明静止流体的应力状态可由一个静压强（数量场）来描述。l 在运动的流体中，既

2、可能有压应力又可能有切应力。把流体在运动状态下的压应力叫做动水压强，以示与静水压强的区别。l 在运动的理想流体里，由于没有粘滞性的作用，虽有质点的相对运动，也不会有切应力，因此理想流体中只有动水压强，而且可用分析静水压强特性的同样方法推证：任一点的动水压强在各方向上的大小都相等，和静水压强有同样的特性。l 在运动的实际流体中，由于粘滞性作用，既有压应力又有切应力。任意一点处的应力是矢量，而且还与作用面方向有关。所以把法向为的作用面上的应力矢量表示为，这里我们定义法线的正方向为受力面的外法向，即法向应力为正表示流体受拉。应力矢量的分量形式为，其中每一个分量的两个脚标的含义是：前一个表示作用面方向

3、；后一个表示应力分量之投影方向。由此，也可知等的含义。l 由如下九个量组成的二阶张量，称为应力张量，记为 主对角线上的三个元素是法应力分量，其它是切应力分量。可以证明这个张量是对称的，所以它只有六个独立的分量。l 有了应力张量P，任意方位作用面上的应力都可知道，为：，如法向为的作用面上应力的y方向的分量为 l 运动流体中的每一点都对应一个应力张量，有了这个应力张量，即可知道该点处任意方位作用面上的应力，可见运动流体的应力状态可由应力张量来描述。l 应力张量主对角线上三个元素之和 是坐标变换中的不变量，即其值不随坐标轴的转动而改变，任意三个相互垂直的作用面上的法应力之和都是相同的。于是可定义 为

4、流体的动压强。它由场点唯一对应，而与作用面的方位无关。所以运动流体中存在一动压强场，它是数量场。要注意p并非任意方位作用面上真正的压应力.l 各向同性的不可压缩牛顿流体的应力和变形速率之间存在线性关系： 其中 它是牛顿内摩擦定律在三维情况下的推广，称为广义牛顿内摩擦定律。42流体运动微分方程一. 以应力表示的流体运动微分方程l 在流场中取出一个空间六面体微元, 按照欧拉观点表述动量守恒原理：单位时间微元内动量的增加必等于单位时间净流入微元的动量加上微元内流体所受合力。l 在单位时间里，净流入微元左右一对表面的y方向的动量为，净流入前后和上下两对表面的y方向的动量分别为 和 .l 作用于六面体微

5、元表面沿y方向的表面力有：左右一对面元法向力 ；前后一对面元切向力 ；上下一对面元切向力 .相加得沿y方向的总表面力 l 作用于六面体微元的沿y方向的质量力为 l 单位时间六面体微元内y方向动量的增加为 l 根据动量守恒原理，得到 将上式的左侧展开 易知上式右侧第一项为，根据连续方程第二项为零。所以 同理可得x,z 两个方向上的方程： 这就是以应力表示的流体的运动方程。X,Y,Z表示质量力的三个分量。二. 不可压缩粘性流体的运动微分方程 N-S方程l 将广义牛顿内摩擦定律代入运动微分方程,即有 最后一个等号是由于用了不可压流体的连续方程。同理可得 x,z方向上的方程，并可合并成如下的矢量式：

6、. 这就是不可压粘性流体的N-S方程，式中 是拉普拉斯算子，表示对变量求调和量。l 不可压粘性流体的N-S方程表明了时变惯性力、位变惯性力、质量力、压差力和粘性力之间的平衡。三理想流体的运动微分方程欧拉方程l 理想流体忽略粘性作用，流体中没有切应力，运动微分方程简化为： 称为欧拉方程。它表明了时变惯性力、位变惯性力、质量力、压差力之间的平衡。l 流体静止时，只受质量力、压差力的作用，运动方程退化为欧拉平衡方程 四. 流体动力学定解问题和解法概述l 基本微分方程组前面导出的微分形式流体运动方程连同连续方程，形成对流体运动的基本控制方程组，是求解流速场和压力场的理论基础。四个方程可求四个未知量：p

7、和，方程组是封闭的。但由于运动方程是二阶偏微分方程，其中的位变惯性力（常称为对流项）是非线性的，解析求解非常困难。l 解法概述 只有在极少数简单流动的情况下，N-S方程才有解析解。而绝大部分流动都不能直接对N-S方程解析求解，我们只能抓住问题的主要方面，作相应的简化，才能进行进一步的解析处理。忽略粘性，作理想流体假设，从流动的维数上作简化，都是常见的手段。如果流动是有势流动，解析处理就有更多的便利条件。后面我们就将分门别类地对各种流动进行求解方法的讨论。应该强调，各种简化都是在基本方程的基础上进行的，所以深入理解方程中各项的物理意义是非常重要的。l 初始条件和边界条件 流体运动基本方程还要加上

8、初始条件和边界条件才能形成流体动力学的定解问题。流体运动所遵循的动力学方程是普遍的，因此流动的个性就体现在初始条件和边界条件上面。初始条件是对不恒定流动指定初始时刻流场的速度和压强分布，边界条件是指运动方程的解在流场的边界上必须满足的运动学和动力学条件。常见的边界条件有：固壁条件和液体的自由表面条件。理想流体的固壁条件称为可滑移条件，即流体不能穿越固壁，但可有切向相对运动，所以 . 实际（粘性）流体的固壁条件称为不可滑移条件，即附着在固壁上的流体质点与固壁不能有相对运动，所以 .这里和分别表示附着在固壁上的流体质点与固壁上相应点的速度。液体的自由表面动力学条件为自由表面上压强为常数（大气压）。

9、五. 理想流体的运动微分方程的积分l 伯努利积分 讨论在恒定条件下理想流体运动方程沿流线的积分。将理想流体的运动方程式写成： ， ， 在流线上沿流动方向取一段弧长，因为恒定流动的流线不随时间变化，迹线与流线重合， 将上述两组各三个等式的左右边分别相乘，然后相加得等式： 上式左边可改写为： 如质量力是有势的，有势函数W，即 ，则右边前三项是力势函数W的全微分： 右边后三项为： ，对于密度为常数的不可压缩流体 经上述变换后，得到 或 . 上式的积分 称为伯努利积分。 伯努利积分表明：在理想流体的恒定流动中，同一流线上各点的值是一个常数。其中W是力势函数，是不可压缩流体的密度。 从推导过程看，积分是

10、在流线上进行的，所以不同的流线可以有各自的积分常数，将它记作，称为流线常数。l 欧拉积分 以上讨论了伯努利积分, 其成立的条件是：理想，恒定，不可压，质量力有势。现在再加无旋（有势）条件，导出理想流体运动方程的欧拉积分。 因为有理想、恒定的条件,所以x方向运动方程为: 再由无旋条件得 ，即 . 同理， . 将三式分别乘上，然后相加，得到 . 上式的积分 称为欧拉积分。 欧拉积分表明：在理想流体的恒定无旋流动中，流场中各点的值是一个常数。其中W是力势函数，是不可压缩流体的密度。 从推导过程看，可在流场中任取，所以积分常数C称为通用常数。 表面上看，伯努利积分和欧拉积分很相似，但两者的适用条件和使

11、用范围是不同的。43恒定总流的能量方程一. 恒定元流的能量方程l 伯努利积分在流线上成立，也在元流上成立。重力场中，理想、不可压流体恒定元流的1-1、2-2两个断面上，有： .l 伯努利积分是欧拉方程的各项单位质量流体受力（加速度）取了势函数而来的。力势函数是能量量纲（力对位移作积分），所以各项都是单位重量流体所具有的能量。是单位重量流体所具有的位置势能（简称单位位置势能），是单位重量流体所具有的压强势能（简称单位压强势能）。是单位重量流体所具有的总势能（简称单位总势能）。为单位重量流体所具有的动能（简称单位动能）。三项之和为单位重量流体的总机械能（简称单位总机械能）。所以伯努利积分表示沿程机

12、械能的守恒，可称为元流的能量方程。l 伯努利积分中的各项也都是长度量纲，可称为水头。 是位置水头，是压强水头，是测压管水头，为速度水头，三项之和 为总水头。l 毕托管利用测压管和测速管得到总水头和测压管水头之差 速度水头，可用来测量流场中某点的流速。要注意理解为什么测压管开口方向应与流速垂直，而测速管开口方向则应迎着流速方向。二. 恒定总流的能量方程l 将恒定元流能量方程累加可得到恒定总流能量方程，即： .l 若过水断面A取在渐变流段中，则其上的测管水头 可视为常数，这样 .l 速度水头 在过水断面A上是不均匀的，设 . 称为动能修正系数，易知 .l 总流通过渐变流段中过水断面的能量通量为 .

13、l 对理想不可压流体恒定总流，不考虑S侧上的能量损失，控制体内能量不变，所以通过各过水断面的能量流量相同，由连续方程决定了重量流量 沿程不变，所以 总流能量方程的上述一维表达要求1-1、2-2两个过水断面处于渐变流段中。l 考虑实际（粘性）流体流动时的能量损失， 断面1-1是上游断面，断面2-2是下游断面，h1-2为总流在断面1-1和2-2之间平均每单位重量流体所损耗的机械能，称为水头损失。三. 恒定总流能量方程的几何表示 水头线l 将恒定总流能量方程的各项水头沿程变化的情况几何表示出来，称为水头线。将三项水头之和称为总水头，记作H .可分别画出测管水头线和总水头线。l 实际（粘性）流体的总水

14、头线肯定是沿程下降的，将水头线的斜率 称为水力坡度，表示单位重量流体在单位长度流程上损失的平均水头。l 水头线不一定是直线，也不一定连续（有局部损失时会有间断）。四. 恒定总流能量方程的应用l 恒定总流能量方程表明三种机械能相互转化和总机械能守恒的规律，由此可根据具体流动的边界条件求解实际总流问题。其实，在实际生活中我们经常会遇到三种机械能相互转化的例子。l 利用恒定总流能量方程时，上下游过水断面必须在渐变流段中选取，并且选在尽可能多地集中已知量和未知量的位置。一般可求解某个水头，或者可求解两断面上测管水头或速度水头之差。l 有能量输入或输出的恒定总流能量方程为 . 注意Hm 应为单位重量流体

15、增加或减少的能量。l 对于四周通大气的流柱，过水断面的测管水头就取为断面形心处的位置水头。44 恒定总流动量方程l 对总流管中上游过水断面1-1和下游过水断面2-2之间的一段有限体积的控制体表述动量守恒原理如下：单位时间控制体内动量的增加必等于单位时间净流入控制体的动量加上控制体内流体所受合力。l 限于讨论恒定流动，因此控制体内动量不会改变，则单位时间净流入控制体的动量加上控制体内流体所受合力为零，即单位时间净流出控制体的动量等于控制体内流体所受合力。对恒定总流的情况，动量守恒原理为：单位时间通过两个过水断面净流出控制体的动量等于该段总流内的流体所受合力。l 单位时间里通过总流过水断面的动量

16、，根据恒定总流的动量守恒原理，应有 l 为把总流动量方程的表达一维化，假设断面A是在渐变流段中的一个过水断面，则断面上各点的方向一致。用断面平均流速v代替u，定义的大小为v，方向为的方向，用代替，那么总流过水断面的动量通量成为 ，它和精确的动量通量方向一致，大小并不相等。为使两者相等，需要补上一个修正系数，从而使动量通量的大小写成 . 可见，这个修正系数应为 ，称为动量修正系数。l 一维化的恒定总流动量方程为 或 其中，进口断面A1，出口断面A2位于渐变流段中。l 讨论 ： = 质量力（重力）侧壁表面力合力两过水断面上表面力合力。侧壁表面力合力往往是欲求的。过水断面取在渐变流段中，其上的表面力没有切向分量，只有法向力，可用断面形心处的相对压强作为平均压强计算。l 恒定总流动量方程中以流出控制体即下游出口断面的动量流量为正，反之为负，这里有一个取正负号的问题。另外，动量方程本身是矢量式，实际求解时须分解成分量式，又有一个向坐标轴投影的问题，也涉及到正负号的选取。应予特别注意。l 应用恒定总流动量方程，只要计算边界上的动量流量就可求取流体所受外力的合力，避开了这段流动内部的细节，在一些关心受力的总效果，而不关心应力分布的问题中显示出便利之处。l 恒定总流动量方程中的 是