流体力学 第七章 不可压缩流体动力学基础

第七章 不可压缩流体动力学基础 在前面的章节中，我们学习了理想流体和粘性流体的流动分析，按照水力学的观点，求得平均量。但是，很多问题需要求得更加详细的信息，如流速、压强等流动参数在二个或三个坐标轴方向上的分布情况。本章的内容介绍流体运动的基本规律、基本方程、定解条件和解决流体问题的基本方法。第一节 流体微团的运动分析运动方式：移动或单纯的位移（平移）旋转线性变形角变形。位移和旋转可以完全比拟于刚体运动，至于线性变形和脚变形有时统称为变形运动则是基于液体的易流动性而特有的运动形式，在刚体是没有的。在直角坐标系中取微小立方体进行研究。一、平移：如果图（a）所示的基体各角点的质点速度向量完全相同时，则构成了液体基体的单纯位移，其移动速度为。基体在运动中可能沿直线也可能沿曲线运动，但其方位与形状都和原来一样（立方基体各边的长度保持不变）。二、线变形：从图（b）中可以看出，由于沿y轴的速度分量，B点和C点都比A点和D点大了，而就代表时液体基体运动时，在单位时间内沿y轴方向的伸长率。 ，三、角变形（角变形速度） 角变形： 四、旋转（旋转角速度） 即，那么，代入欧拉加速度表达式，得：各项含义：（1） 平移速度（2）线变形运动所引起的速度增量（3）（4）角变形运动所引起的速度增量（5）（6）微团的旋转运动所产生的速度增量流体微团的运动可分解为平移运动，旋转运动，线变形运动和角变形运动之和。 亥姆霍兹速度分解定理第二节 有旋运动1、无涡流（势流） 如在液体运动中，各涡流分量均等于零，即，则称这种运动为无涡流。 当满足无涡流条件时，满足柯西条件，就有：存在。 即流速势。满足此条件的流动（无涡流）就叫势流。（下一章作详细介绍）2、有涡流：如在液体运动中，涡流分量、及中间的任一个或全部不等于零，则这样的液体运动就叫做旋流或有涡流。自然界中的实际液体几乎都是这种有涡的流动。 涡线：流场中一些假想的线，在所讨论的瞬时，涡线上各个质点的涡旋向量都与此线在该点处相切。与流线同样的分析方法，得到涡线方程：涡量：设流体微团的旋转角速度为，则称为涡量，是与空间坐标和时间有关的矢量函数。其中、和是涡量在、坐标上的投影。根据旋转角速度的定义，有： 哈米尔顿算子是一矢量算子，可知， 那么， 就自然满足。或者写成， 即涡量的定义使之自然满足涡量连续性微分方程。例：已知某圆管（半径）中液体流动的流速分布为： 试判断该流动是有涡流还是无涡流？并求涡线微分方程。 所以，该流动是有涡流。将上三式代入涡线微分方程，得： 积分后，得到： 涡线是和管轴同轴的同心圆。涡管：在涡量场中任意画一封闭曲线，通过这条曲线上的每一点所做出的涡线构成一管状的曲面，称为涡管。涡通量：设A为涡量场中一开口曲面，微元面dA的外法线单位向量为，涡量在方向上的投影为，则面积积分 称为涡通量。有旋运动的一个重要的运动学性质：在同一瞬间，通过同一涡管的各截面的涡通量相等。证明：我们知道，根据涡量的定义，可以很容易知道，涡量自然满足涡量连续性微分方程，即： ，对这个微分方程在任意封闭体积上作积分，也是满足的，若任意体积取为，一段涡管和两个截面A1和A2，就有： 可以将体积分化成封闭曲面积分： 其中 所以， 得证对于微元涡管，近似认为截面上各点的涡量为常数， 性质：涡管不可能在流体内部开始或终止，而只能在流体中自行封闭成涡环，或者终止于和开始于边界面。龙卷风开始于地面，终止于云层。速度环量：在流场中任取一封闭曲线s，则流速沿曲线s的积分： 称为曲线s上的速度环量，并规定积分沿s逆时针方向绕行为的正方向。（一）斯托克斯定理根据斯托克斯公式， 性质：沿任意封闭曲线s的速度环量等于通过以该曲线为边界的曲面A的涡通量。斯托克斯定理。（二）汤姆逊定理汤姆逊定理：在理想流体的涡量场中，如果质量力具有单值的势函数，那么沿由流体质点所组成的封闭曲线的速度环量不随时间而变，即： 解释：速度环量=涡通量，所以，流体的涡旋具有不生、不灭的性质。第三节 不可压缩流体连续性微分方程1、微分形式的连续性方程在推导这个方程式时,我们认为运动着的液体系连续地充满它所占据的空间,流动时不形成空隙,并且表征液体运动的各物理量也都是时间和空间的连续函数。在时间,于流场中取一具有边长为、的微分六面体,在随后的一无限小段内,流进和流出该微分六面体的质量。流出-流入=质量增量。微分六面体形心A点的坐标为（、）,密度为,质点的速度分量为、及,则在时段内沿轴从左侧面流入六面的液体质量为 流出的液体质量为：质量的变化：联立，得到：（一般形式的液体连续性方程）适合可压缩和不可压缩液体。或，写成： （适合不可压缩液体，恒定流和非恒定流）它是质量守恒定律在水力学中的表现形式。它表征着不可压缩液体在运动时,若保持其连续性,则线性变形必系伸长现象与缩短现象同时发生。2、积分形式的液体连续性方程 连续性方程写成矢量形式： 其中为微分算子。体积积分：根据高斯公式，对于恒定流，对于不可压缩， n是液体边界的外法线方向考虑到速度和面积的方向，就可知：，即， （微小流束的流量平衡）积分后，可以得到， 其中、为各自断面上的断面平均流速。例：判断，流速为：，的流动是否满足连续性方程。解： ，那么 满足 ，所以，满足连续性方程。第四节 以应力表示的粘性流体运动微分方程式一、 粘性流体的内应力表面力，9个分量：X方向： 或，同理：其中，密度对于不可压缩流体是已知常量，单位质量力也是已知常量；未知量为9个应力和三个速度分量。不容易求解。第五节 应力和变形速度的关系一、 切应力和角应变速度的关系一元流动的牛顿内摩擦定律为： 或可写为， 切应力与流速梯度或直角变形速度的关系。是直角变形速度，它是角变形速度的2倍，在xoy平面上，那么，对于三元流动的牛顿内摩擦定律，可以写成如下形式： 六个切应力均可用粘性系数和直角变形速度的乘积来表示。二、法向应力和线变形速度的关系在理想流体中，同一点各方向的法向应力相等，（代表是压应力）。在粘性流体中，粘性不仅产生与切应力有关的角变形速度，而且使线变形速度也产生附加法向应力。使一点的法向应力与作用面方位有关。取边长为的方形流体微团进行研究，先考虑方形微团在x方向上的伸长变形。微团在x方向上作伸长变形时，伸长为，而对角线旋转至，使产生角变形，这样在面上产生切应力，线变形所产生的就要有力来平衡，这样就在面上产生了附加法向切应力。根据这样的分析，就有和力的平衡， 其中，是微团在x方向上伸长的长度，那么，化简后， 其中，是45度角的角变形速度，那么，附加法向应力和线变形速度的关系：综合一下，得到，线变形运动使法向应力随伸长变形而减小，就有，这就是粘性流体法向应力和线变形速度的关系。其中，为理想液体的压强，它的大小与作用面方位无关。在粘性流体中，法向应力应该与方向有关了，所以，定义任意一点上三个相互垂直平面上的法向应力的平均值的负值为粘性流体在该点的压强。对于满足连续性方程的不可压缩流体，所以，。而对于可压缩流体，代表的是质点的体积膨胀率，与坐标的选择无关，而压强是空间坐标的函数，与方向无关。 第六节 N-S方程将的表达式代入以应力表示的粘性流体运动微分方程，得到：若流体粘滞性是常数，那么， = 对于不可压缩流体，有连续性方程，所以，同理，可得，不可压缩粘性流体的运动微分方程它与连续性方程联立，求得速度的三个分量和压强。上式中，是流体质点的加速度，对于欧拉的描述方法，第七节 理想流体运动微分方程及其积分当流体为理想液体时，运动粘滞系数，n-s方程就简化为： （2-23）将式（213）代入，则上式化为 （2-24）再引用和的表达式，则又可化为 （2-25）上式为欧拉方程式用整体运动性质（位移、变形及旋转）所表现的形式，说明液体在各力作用下而运动时，表现有位移、变形及旋转作用的可能。对于不可压缩的理想液体，欧拉方程及连续方程提供出解决运动问题的四个独立条件。一般而论，可用以解决其中四个未知数，例如、及。二、欧拉哥罗米柯方程式下面我们再把欧拉方程变成具有涡旋分量的另一形式。由流速可求得今将此式与欧拉方程式（224）中的第一式相减则得因得同理 （2-26）今设体积力（、）为有力势的力，亦即体积力分量可由其势能（、）确定如下：将之代入式（226），稍加整理则得 （2-27）式（227）即是在有势力作用下用涡旋分量表示的理想液体的运动方程式，又称欧拉哥罗米柯方程式。三、理想液体伯努力积分形式的能量方程（25分钟）由于水流是在重力场中运动，重力为有势力所以，存在一势能函数，水为不可压缩的液体，为常数故上式可写为能量方程的推导：基于上述柯罗米柯方程，通过观测：第一种情况无旋流 则存在一流速势函数，且代入柯罗米柯方程，推出 能量方程单纯重力作用时，第二种情况恒定流 积分：如果行列式等于零，则，与无关单纯重力作用下理想液体恒定流的能量方程时，满足条件式1-4其中的一条，就可以满足行列式为零的要求1，无旋流或势流，并可用于流场空间内所有的空间点上。2这是一个流线的微分方程式，所以适用于同一个流线之上。3，这是一条涡线的方程式，所以适用于同一个涡线上各点。4，恒定流中以流线与涡线相重合为特征的螺旋流，所以方程式也适用于恒定螺旋流中。根据以上论述，对于由上述1、2、3、4各种情况所限定的两个点上（点1及点2），则得 （234）各项含义：位能、压能和动能。总的意义：液体在运动过程中，尽管其机械能可以相互转化，但总的机械能是守恒的，上式是普遍能量守恒原理在理想液体中的表现形式。对于实际液体：积分形式势流推出积分沿流线：； （为摩擦所作的功）任意两点： （为惯性消耗，为惯性水头）第八节 流体流动的初始条件和边界条件 粘性流体的基本方程是二阶偏微分方程，现在的任务就由原本结果具体的流体问题，转变成了解决这个粘性流体基本方程的解的问题。方程的解不仅需要给出描述流体运动的方程，而且非常重要的是，需要给出流动的初始条件（时间变量）和边界条件（空间变量）。一、 初始条件：方程组的解在初始时刻应满足的条件。在初始时刻，给出： 二、边界条件：在流场的边界，方程组的解应满足的条件。（1）在固体接触面上 粘性，无滑移边界条件 固壁静止 理想流体，有滑移边界条件（2） 不同液体的分界面：两侧液体的速度和压强保持连续。 ，（3） 液体和蒸汽的界面，若不考虑液面上饱和蒸汽中的动量、热量和质量交换时。 液体在平均液面垂直方向上的速度等于液面在垂直与平均液面方向上的高度随时间的变化率的相反数。（4） 自由液面，若忽略表面张力的影响，（5） 流道的入口和出口的边界条件入口和出口断面上的流速和压强的分布以及温度的变化，流量的输入速度等的边界条件，具体问题具体分析。第九节 不可压缩粘性流体紊流运动的基本方