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

第七章不可压缩流体动力学基础在前面的章节中，我们学习了理想流体和粘性流体的流动分析，按照水力学的观点，求得平均量。但是，很多问题需要求得更加详细的信息，如流速、压强等流动参数在二个或三个坐标轴方向上的分布情况。本章的内容介绍流体运动的基本规律、基本方程、定解条件和解决流体问题的基本方法。第一节流体微团的运动分析运动方式：①移动或单纯的位移（平移）②旋转③线性变形④角变形。位移和旋转可以完全比拟于刚体运动，至于线性变形和脚变形有时统称为变形运动则是基于液体的易流动性而特有的运动形式，在刚体是没有的。在直角坐标系中取微小立方体进行研究。一、平移：如果图（a）所示的基体各角点的质点速度向量完全相同时，则构成了液体基体的单纯位移，其移动速度为zyxuuu、、。基体在运动中可能沿直线也可能沿曲线运动，但其方位与形状都和原来一样（立方基体各边的长度保持不变）。二、线变形：从图（b）中可以看出，由于沿y轴的速度分量，B点和C点都比A点和D点大了dyyuy，而yuy就代表1dy时液体基体运动时，在单位时间内沿y轴方向的伸长率。xux，yuy，zuz三、角变形（角变形速度）ddddDCABCDBAdtyudydtdyyudxxdtxudxdtdxxudyydd2dd角变形：xuyudddyxz212xuzuzxy21yuzuzyx21四、旋转（旋转角速度）yuxuxyz21zuyuyzx21即，xuzuzxy21zyxuuuzyxkji21那么，代入欧拉加速度表达式，得：zxxxxxxzyyzzyyyyyyyxzzxxzzzzzzzyxxyyxxyduuuuuuuudttxuuuuuuuudttyuuuuuuuudttz各项含义：（1）平移速度（2）线变形运动所引起的速度增量（3）（4）角变形运动所引起的速度增量（5）（6）微团的旋转运动所产生的速度增量流体微团的运动可分解为平移运动，旋转运动，线变形运动和角变形运动之和。——亥姆霍兹速度分解定理第二节有旋运动1、无涡流（势流）如在液体运动中，各涡流分量均等于零，即0zyx，则称这种运动为无涡流。当满足无涡流条件时，yzxzyxuuyzuuzxuuxy，满足柯西条件，就有：xyzuxuyuz存在。即流速势。满足此条件的流动（无涡流）就叫势流。（下一章作详细介绍）2、有涡流：如在液体运动中，涡流分量x、y及z中间的任一个或全部不等于零，则这样的液体运动就叫做旋流或有涡流。自然界中的实际液体几乎都是这种有涡的流动。涡线：流场中一些假想的线，在所讨论的瞬时，涡线上各个质点的涡旋向量都与此线在该点处相切。321xoy与流线同样的分析方法，得到涡线方程：zyxdzdydx涡量：设流体微团的旋转角速度为tzyx,,,，则kjizyx2称为涡量，是与空间坐标和时间有关的矢量函数。其中x、y和z是涡量在x、y、z坐标上的投影。根据旋转角速度的定义，有：zuyuyzxxuzuzxyyuxuxyz哈米尔顿算子是一矢量算子，kzjyix，可知，yuxujxuzuizuyuuuuzyxkjiuxyzxyzzyxu那么，0u就自然满足。或者写成，0zyxzyx即涡量的定义使之自然满足涡量连续性微分方程。例：已知某圆管（半径0r）中液体流动的流速分布为：22204zyrJux0yu0zu试判断该流动是有涡流还是无涡流？并求涡线微分方程。021zuyuyzxzJxuzuzxy421yJyuxuxyz421所以，该流动是有涡流。将上三式代入涡线微分方程，zyxdzdydx，得：yJdzzJdy440zdzydy积分后，得到：Cyz22涡线是和管轴同轴的同心圆。涡管：在涡量场中任意画一封闭曲线，通过这条曲线上的每一点所做出的涡线构成一管状的曲面，称为涡管。涡通量：设A为涡量场中一开口曲面，微元面dA的外法线单位向量为n，涡量在n方向上的投影为n，则面积积分AzyxAnAdxdydzdxdydzdAAdJ称为涡通量。有旋运动的一个重要的运动学性质：在同一瞬间，通过同一涡管的各截面的涡通量相等。证明：我们知道，根据涡量的定义，可以很容易知道，涡量自然满足涡量连续性微分方程，即：0zyxzyx，对这个微分方程在任意封闭体积上作积分，也是满足的，若任意体积取为，一段涡管和两个截面A1和A2，就有：0dVzyxvzyx可以将体积分化成封闭曲面积分：321AAAzyxdxdydxdzdydz321AzyxAzyxAzyxdxdydxdzdydzdxdydxdzdydzdxdydxdzdydz321AAAAdAdAd其中03AAd21AAdAndAn021AnAndAdA所以，21AnAndAdA得证对于微元涡管，近似认为截面上各点的涡量为常数，2211AA性质：涡管不可能在流体内部开始或终止，而只能在流体中自行封闭成涡环，或者终止于和开始于边界面。龙卷风开始于地面，终止于云层。速度环量：在流场中任取一封闭曲线s，则流速沿曲线s的积分：szyxsdzudyudxusdu称为曲线s上的速度环量，并规定积分沿s逆时针方向绕行为s的正方向。（一）斯托克斯定理根据斯托克斯公式，szyxssdzudyudxusduAxyzxyzdxdyyuxudxdzxuzudydzzuyuAAAzzyyxxJAddAdAdA性质：沿任意封闭曲线s的速度环量等于通过以该曲线为边界的曲面A的涡通量。——斯托克斯定理。（二）汤姆逊定理汤姆逊定理：在理想流体的涡量场中，如果质量力具有单值的势函数，那么沿由流体质点所组成的封闭曲线的速度环量不随时间而变，即：0dtdyxdydzd'c'a'b'Azdxdcabxozy解释：速度环量=涡通量，所以，流体的涡旋具有不生、不灭的性质。第三节不可压缩流体连续性微分方程1、微分形式的连续性方程在推导这个方程式时,我们认为运动着的液体系连续地充满它所占据的空间,流动时不形成空隙,并且表征液体运动的各物理量也都是时间和空间的连续函数。在时间t,于流场中取一具有边长为dx、dy、dz的微分六面体,在随后的一无限小段dt内,流进和流出该微分六面体的质量。流出-流入=质量增量。微分六面体形心A点的坐标为（x、y、z）,密度为,质点的速度分量为xu、yu及zu,则在dt时段内沿x轴从左侧面abcd流入六面的液体质量为xdx21()2xxuudxx1()2xxuudxxxdx2dydzdtdxxdxxuuxx2)]21([流出的液体质量为：dydzdtdxxdxxuuxx2)]21([质量的变化：dxdydzdttdxdydzdxdydzdtt联立，得到：dxdydzdttdxdydzdtzudxdydzdtyudxdydzdtxuzyx0zuyuxutzyx（一般形式的液体连续性方程）适合可压缩和不可压缩液体。或，写成：0zuyuxudtdzyx0zuyuxuzyx（适合不可压缩液体，恒定流和非恒定流）它是质量守恒定律在水力学中的表现形式。它表征着不可压缩液体在运动时,若保持其连续性,则线性变形必系伸长现象与缩短现象同时发生。2、积分形式的液体连续性方程连续性方程写成矢量形式：0ut其中为微分算子。体积积分：0dutv根据高斯公式，0AvdAnudt对于恒定流，0AdAnu对于不可压缩，0AdAnun是液体边界的外法线方向考虑到速度和面积的方向，就可知：02211dAudAu，即，2211dAudAu（微小流束的流量平衡）积分后，可以得到，2211AvAv其中1v、2v为各自断面上的断面平均流速。例：判断，流速为：22yxyux，22yxxuy，0zu的流动是否满足连续性方程。解：22yxyux，22yxxuy，0zu，那么2222yxxyxux2222yxxyyuy0zuz满足0zuyuxuzyx，所以，满足连续性方程。第四节以应力表示的粘性流体运动微分方程式一、粘性流体的内应力表面力，9个分量：zzzyzxyzyyyxxzxyxxpppX方向：xxmaFdydzdxxppdxxppdxdydzXxxxxxxxx)]2()2[(dxdydzzdzzdzdxdyydyyzxzxzxzxyxyxyxyx)]2()2[()]2()2[(dtdudxdydzXxdtduzyxpXxzxyxxx或，dtduzyxpXxzxyxxx)(11同理：dtduzyypYyzyxyyy)(11dtduzyxpZzxzyzzz)(11其中，密度对于不可压缩流体是已知常量，单位质量力ZYX,,也是已知常量；未知量为9个应力和三个速度分量。不容易求解。第五节应力和变形速度的关系一、切应力和角应变速度的关系一元流动的牛顿内摩擦定律为：dydu或可写为，dtd切应力与流速梯度或直角变形速度的关系。dtd是直角变形速度，它是角变形速度的2倍，在xoy平面上，xuyudtdyxz2那么，对于三元流动的牛顿内摩擦定律，可以写成如下形式：zuyuxuzuxuyuyzzyyzzxxzzxyxyxxy六个切应力均可用粘性系数和直角变形速度的乘积来表示。二、法向应力和线变形速度的关系在理想流体中，同一点各方向的法向应力相等，0,pppppzzyyxx（代表是压应力）。在粘性流体中，粘性不仅产生与切应力有关的角变形速度，而且使线变形速度zuyuxuzyx,,也产生附加法向应力。使一点的法向应力与作用面方位有关。取边长为dy