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

第七章不可压缩流体动力学基础在前面的章节中，我们学习了理想流体和粘性流体的流动分析，按照水力学的观点，求得平均量。但是，很多问题需要求得更加详细的信息，如流速、压强等流动参数在二个或三个坐标轴方向上的分布情况。 本章的内容介绍流体运动的基本规律、基本方程、定解条件和解决流体问题的基本方法。第一节流体微团的运动分析运动方式：①移动或单纯的位移（平移）②旋转③线性变形④角变形。位移和旋转可以完全比拟于刚体运动，至于线性变形和脚变形有时统称为变形运动则是基于液体的易流动性而特有的运动形式，在刚体是没有的。在直角坐标系中取微小立方体进行研究。X(b)AH(d)X(b)AH(d)、平移：如果图（a所示的基体各角点的质点速度向量完全相同时，贝u构成22了液体基体的单纯位移，其移动速度为Ux、Uy、Uz。基体在运动中可能沿直线也可能沿曲线运动，但其方位与形状都和原来一样（立方基体各边的长度保持不变）。二、线变形：从图（b）中可以看出，由于沿y轴的速度分量，B点和C点都比.：UyA.：UyA点和D点大了7；Uydy，而一上就代表dy二1时液体基体运动时，在单位时间内.:y沿y轴方向的伸长率。CUxCUx dUy-：Uz.z.x.z、角变形（角变形速度）Dd：'d:■.:ux一dydt：:ydy：=y角变形：亠d:■.:ux一dydt：:ydy：=y角变形：亠0：：Uydxdt:Uydt=d：_r-e6矽丿四、旋转（旋转角速度）■y：：U■y：：Ux•，z = -二=1即,当满足无涡流条件时,L、；Ujz_：X，满足柯西条件，就有:二一存在：=y即,当满足无涡流条件时,L、；Ujz_：X，满足柯西条件，就有:二一存在：=y「u.:x；：yUz；：ux ；：uz—. — Jfz ;：X1i二.k2czUxUyUz那么，代入欧拉加速度表达式，得：C（X=叫+Ux竺+UzTy+Uyg+Uz国y—UyBTOC\o"1-5"\h\zdt ct cX zCUy和y CUyGy= = +Uy 十Uy日Z+Uz^X+UxC^—Uz^x>dt ct cyCUZ 和z CUZctz=——=——+Uz——+Uy6lx+Ux6+Uy国x-UxCOydt ct cz各项含义：（1） 平移速度（2） 线变形运动所引起的速度增量（3） （4）角变形运动所引起的速度增量5）（6）微团的旋转运动所产生的速度增量流体微团的运动可分解为平移运动，旋转运动，线变形运动和角变形运动之和。――亥姆霍兹速度分解定理第二节有旋运动1、无涡流（势流）如在液体运动中，各涡流分量均等于零，即 「x=•；y=J：z=0，则称这种运动为无涡流。?U:即流速势。满足此条件的流动（无涡流）就叫势流。 （下一章作详细介绍）2、有涡流：如在液体运动中，涡流分量 「y及亠中间的任一个或全部不等于零，则这样的液体运动就叫做旋流或有涡流。自然界中的实际液体几乎都是这

种有涡的流动涡线：流场中一些假想的线，在所讨论的瞬时，涡线上各个质点的涡旋向量都与此线在该点处相切与流线同样的分析方法，得到涡线方程:dxdydz涡量：设流体微团的旋转角速度为■x,y,z,t，则门=2yYxi•iJj"zk称为涡量，是与空间坐标和时间有关的矢量函数其中门涡量，是与空间坐标和时间有关的矢量函数其中门x、门y和门z是涡量在x、y、z坐标上的投影。根据旋转角速度的定义，有::：:：Ux：：Uz：「Uy；=Ux:z.:z.:xjy哈米尔顿算子是矢量算子,.:x；:yz可知，Iu:哈米尔顿算子是矢量算子,.:x；:yz可知，Iu:x jyUxUy;:zUz那么，I-Iu=0就自然满足。或者写成,"x:'x ::y:z即涡量的定义使之自然满足涡量连续性微分方程。12；：Uz——Ijx汶:y12；：Uz——Ijx汶:y4.1例：已知某圆管（半径r。）中液体流动的流速分布为:J2 2 2Ux r。-yz4A试判断该流动是有涡流还是无涡流？并求涡线微分方程。J -z4」所以，该流动是有涡流将上三式代入涡线微分方程,dxdy'ydz将上三式代入涡线微分方程,dxdy'ydz，得:--zdyJ

z

4」ydyzdz=0dyJ

z

4」ydyzdz=0积分后，得到： z2亍涡线是和管轴同轴的同心圆。涡管：在涡量场中任意画一封闭曲线,通过这条曲线上的每一点所做出的涡线构成一管状的曲面，称为涡管成一管状的曲面，称为涡管涡通量：设A为涡量场中一开口曲面，微元面dA的外法线单位向量为n，涡量在n方向上的投影为，则面积积分J=a'」dA二a'」ndA=a'"Xdyd 1ydzd 1zdxdy称为涡通量。有旋运动的一个重要的运动学性质：在同一瞬间，通过同一涡管的各截面的涡通量相等。证明：我们知道，根据涡量的定义，可以很容易知道，涡量自然满足涡量连续性微分方程，即：二•丄•二=0，对这个微分方程在任意封闭体积上作积分， 也是满:x :y :z00足的，若任意体积取为，一段涡管和两个截面A1和A2，就有:二•丄•二dV；:xjy ；:z可以将体积分化成圭寸闭曲面积分:A1.A2A^dydz—dxdz'dxdy=ail'xdydz…::-ydxdz…::-zdxdy亠iC.xdydz亠Oydxdz亠Ozdxdy亠i::xdydz•】ydxdz•JdxdyTOC\o"1-5"\h\zA3 y-dA亠|OdA亠|O.dA其中门dA=0^A1 =A2 =A3 =A3- ndA亠Ii：：ndAA1 -A2一AIS鳥儿吩。所以，JndA「A2.dA 得证对于微元涡管，近似认为截面上各点的涡量为常数,\ ='.12A2性质：涡管不可能在流体内部开始或终止， 而只能在流体中自行封闭成涡环， 或者终止于和开始于边界面。龙卷风开始于地面，终止于云层。速度环量：在流场中任取一封闭曲线 s，则流速沿曲线s的积分：-二：uds二uxdxUydyuzdz称为曲线s上的速度环量，并规定积分沿 s逆时针方向绕行为s的正方向。（一） 斯托克斯定理根据斯托克斯公式，丨s='uds='uxdxuydyuzdz=f『匹一邑；ydz+包-邑»dJ邑-空〕dxdyaR®&丿 （往&丿 ex cyJ-门xdAxi」ydAy"zdAz-[dA二JA性质：沿任意封闭曲线s的速度环量等于通过以该曲线为边界的曲面 A的涡通量。一一斯托克斯定理。（二） 汤姆逊定理汤姆逊定理：在理想流体的涡量场中，如果质量力具有单值的势函数， 那么沿由流体质点所组成的封闭曲线的速度环量不随时间而变，即：d：dt则在dt时段内沿x轴从左侧面abed流入六面的液体质量为uxx1(dx)则在dt时段内沿x轴从左侧面abed流入六面的液体质量为uxx1(dx)2cPdx

P+——2&解释：速度环量=涡通量，所以，流体的涡旋具有不生、不灭的性质。第三节不可压缩流体连续性微分方程1、微分形式的连续性方程在推导这个方程式时，我们认为运动着的液体系连续地充满它所占据的空间 ，流动时不形成空隙，并且表征液体运动的各物理量也都是时间和空间的连续函数。在时间t,于流场中取一具有边长为 dx、dy、dz的微分六面体,在随后的一无限小段dt内,流进和流出该微分六面体的质量。流出 -流入=质量增量。Uy及Uz微分六面体形心A点的坐标为（x、y、z）,密度为Uy及UzTOC\o"1-5"\h\z'ux1 ■dx[ux一(dx)] dydzdt£x 2P ex2流出的液体质量为： [ux•二』dx）]「•二岂dydzdt'x2 x2质量的变化：(P、P+—dtdxdydz-Fdxdydz◎P=——dxdydzdtCt联立，得到：贸仏)- dxdydzdt讽My)...负仏） - dxdydzdt— dxdydzdtcyczcPdxdydzdtcPdPu ePuycPU」=0（—般形式的液体连续性方程）适合可压缩和不,:t 汶 ;:y:z可压缩液体。或，写成:+——.:y.:u=0土•丄•二=0（适合不可压缩液体，恒定流和非恒定流）.X 鋼:z,若保持其它是质量守恒定律在水力学中的表现形式。它表征着不可压缩液体在运动时连续性，,若保持其2、积分形式的液体连续性方程连续性方程写成矢量形式： 二-\■:-u=0其中为微分算子。体积积分:根据高斯公式， 'd u-ndA=0…v泾 “A对于恒定流，「IILiundA=0对于不可压缩，n对于不可压缩，n是液体边界的外法线方向考虑到速度和面积的方向，就可知:-考虑到速度和面积的方向，就可知:-u1dA1u2dA2=0，即，u1dA1=u2dA2 （微小流束的流量平衡）积分后，可以得到，wA=v2A2其中v、v2为各自断面上的断面平均流速。解:ux.:ux2xy& (x2+解:ux.:ux2xy& (x2+y22ux二-yx一22，x yuy二2 2，x y例：判断，流速为:uyX-2xyx2y2uz=0的流动是否满足连续性方程。uz=0，那么亠0；z满足且•凹•竺=0，所以，满足连续性方程。:x :y :z第四节以应力表示的粘性流体运动微分方程式、粘性流体的内应力表面力，9个分量:丨丨2Pxx -xy -xz-yxpyy yz-zx 、zy pzzX方向:、Fx=ma:Xdxdydz[(pxxxxdx;:xBxxdxxx;:x)]dydz[(yxcTyxdy)-(■yxdy)]dzdx亠[(.：•-乂zxdzzx)—(•CTzxdzzx;:z2)]dxdy或,dxdydzduxdt-Pxx：：zx?x;z1::pxx 1Jyxyx2二；-X?Xr;:x:y同理:yy.■xdudt：：zxduCTxydtdu:-y;:zyz：\z;:zdtdudt其中，密度「对于不可压缩流体是已知常量，单位质量力X,Y,Z也是已知常量；未知量为9个应力和三个速度分量。不容易求解。第五节应力和变形速度的关系、切应力和角应变速度的关系一元流动的牛顿内摩擦定律为:.=」巴 或可写为—二dr 切应力与流速梯度或直角变形速度的关系。dy dtd.■'是直角变形速度，它是角变形速度的 2倍，在xoy平面上，dt那么，对于三元流动的牛顿内摩擦定律，可以写成如下形式:zxyz二-zyzxyz二-zyjz ；:x旬y'+Z&丿六个切应力均可用粘性系数和直角变形速度的乘积来表示。法向应力和线变形速度的关系在理想流体中，同一点各方向的法向应力相等,P在理想流体中，同一点各方向的法向应力相等,Pxx=Pyy=Pzzi-P，P-°（代表是压应力）。在粘性流体中，粘性不仅产生与切应力有关的角变形速度，而且使线变形速度x方向上的伸长变形。x方向上的伸长变形。取边长为dx二dy的方形流体微团进行研究，先考虑方形微团在微团在x方向上作伸长变形时， BC伸长为BC，而对角线AC旋转至AC，使二产u生角变形，这样在AC面上产生切应力n，线变形-所产生的n就要有力来平衡，这样&就在AB面上产生了附加法向切应力 ,xx。根据这样的分析，就有 .n和.xx力的平衡,■xx0dxcos45xxdrABAC■xx0dxcos45xxdrABACCCsin450,2dxu其中，CC•是微团在x方向上伸长的长度， CC立dxdtex那么，化简后，兰1:：那么，化简后，兰1:：Uxdtdv?-n二2Jdt2；：Ux其中，dv是45度角的角变形速度,ex那么，附加法向应力和线变形速度的关系:xxx二2■—cux=2卩一-exCUy=2卩一->cux=2卩一-exCUy=2卩一->cy=2卩竺程J线变形运动使法向应力随伸长变形而减小，就有,综合一下，得到,Pyy二一PtXXyyPzz+2卩亜&

和y!+2卩——cy+2卩生czJ这就是粘性流体法向应力和线变形速度的关系。其中,Pt为理想液体的压强，它的大小与作用面方位无关。在粘性流体中，法向应力应该与方向有关了， 所以，定义任意一点上三个相互垂直平面上的法向应力的平均值的负值为粘性流体在该点的压强。—PxxPyy-Pzz1=Pt—PxxPyy-Pzz1=Pt3 ；:y■:z对于满足连续性方程的不可压缩流体,.:x:y•且=0，所以,.:zP=Pt。而对于可压缩流体，土•且•且代表的是质点的体积膨胀率， 与坐标的选择无关，而.:x :y :z压强P是空间坐标的函数，与方向无关。yy二-P2「显+邑3 ex cyzzyy二-P2「显+邑3 ex cyzz3 ：：x;:yjz第六节N-S方程I%I%将pxx的表达式代入以应力表示的粘性流体运动微分方程，得到:-p^-：Ux_2ex3+12hp矽＜若流体粘滞性是常数，那么,:y ；：xdUx1—dtX亠P ；:x「22」Ux匚 2"3J：：x222 Uy2.二疗uz"：：x：：z_2x■y_2UyP次£yJdxJr2Uz2；=zdUxdt==>；：P.:x2::Ux2；=x2::Ux2；=x+ 13 Wx ；:xUx(x对于不可压缩流体，有连续性方程,XPex-2x

-2

:X¥「UxC2-y同理，可得,Ux+v(—2x「X :x-2:'Uy'(2.xC2U■A1;:p1B.-；-斜12P-2:-Uxr2-2:'Uy-2■yrd口2-y+zz::Uy-)

f2 ):z它与连续性方程.:u_x:y :z-y:-Ux;xdUxdtdUx

dtdUydt

dUzdt=0联立,上式中,匹是流体质点的加速度,dt；：2Ux.:y22::Ux2工dUx

dt邑).z2x-2:y2J:UxPcz2dUx

dtdUxdt5-:yJUz;:z所以,＞不可压缩粘性流体的运动微分方程求得速度的三个分量和压强。对于欧拉的描述方法，竺F竺"竺F:x :y .'z22(IUx+v( -P；x：：x21CP 已UyY-A2「-y-x21；：P/:UzZ (r! .z ：：X2:UxT-=y2rUy2:y2:Uz2""

-=y2:z2；：Uy+2)=

jz2字)二之；Uyjt-：Uz当流体为理想液体时,第七节运动粘滞系数将式（2-13）代入，则上式化为UxuQjUy.x；x；：Ux；z.：Uy .：UyUy Uz■y:z£Uz 丄 CUzUy Uz-:y：z理想流体运动微分方程及其积分=0，n-s方程就简化为:1dp dUxX—!-：..xdti :.xjt丫亠—fti■'二y_ 1 cp cuzZ—再引用二和••的表达式，则又可化为1 ；：p；x;:p:y：：PdUydtduz4Uy土U-：Uxjx：:y(2-23)/UyUy;uz-Uyjx:Uyy-U;uz--uUz^y)-(Uy•-Uzy)(2-24)/UyMy(UzU■y：：Uz ：：Uz /-Uz- (U；-jz；:t ■1卫；-jy1 ：：pUyUx^)—(Uz'■一Uxz)(2-25).'zUyd)—(Uxy—Uyx)上式为欧拉方程式用整体运动性质（位移、变形及旋转）所表现的形式,说明液体在各力作用下而运动时，表现有位移、变形及旋转作用的可能。对于不可压缩的理想液体， 欧拉方程及连续方程提供出解决运动问题的四个独立条件。而论，可用以解决其中四个未知数，例如 Ux、Uy、Uz及PO、欧拉-哥罗米柯方程式F面我们再把欧拉方程变成具有涡旋分量的另一形式。由流速u22二Ux Uyuz2可求得2

：：u(—):X 2二Ux且uCX：x今将此式与欧拉方程式2—24）中的第一式相减则得21沖：：U()I i.X ::X同理rt因得X一丄兰PcX&U();:XUz(巴cz：：同理rt因得X一丄兰PcX&U();:XUz(巴cz：：UzjX2Uz，y-2Uy，z2(u^-y—Uy'z)CUy)-Uy匸Y1Up-：Uy ::_(2u)2u(X；Tz-UP.t :y21;:p：：uz ；：2uZ---(—)2比•X-uXP:Z.t :z2(2-26)今设体积力（X、Y、Z）为有力势的力，亦即体积力分量可由其势能二-二（X、确定如下:cn.：cn.：X-:y;:z将之代入式（2—26）,稍加整理则得£Ux;:t（二■ U_-U；'y£Ux;:t（二■ U_-U；'y):XP2（二2』.U_J=2(Uz，'x-U/'z)P2(■:2■-Pu_)=2(U/'y■uy'x):zP2；=tUz;:t2(2-27)式（2—27）即是在有势力作用下用涡旋分量表示的理想液体的运动方程式，又称欧拉-哥罗米柯方程式。三、理想液体伯努力积分形式的能量方程（ 25分钟）①由于水流是在重力场中运动，重力为有势力所以，存在一势能函数 所以，存在一势能函数 兀X:x .:ycJi:：z②水为不可压缩的液体,1P—(P)i■':.x ::Xj：?故上式可写为2.uX pu(二 )=2(Uy・・z—Uz「y)TOC\o"1-5"\h\z:：t T2/Uy-/Uy-—(-PJ=2(匕x-UxJ.:y? 22「Uz ：- p u( )=2(Ux，y-UyX)；:t :z ? 2能量方程的推导：基于上述柯罗米柯方程，通过观测:①第一种情况无旋流•无旋流•'x二-y则存在一流速势函数 =(x,y,z,t)，且汐Ux ,Uy ,UzTOC\o"1-5"\h\zex cy cz代入柯罗米柯方程,三(二)一二RX)=o.tx:x■- 2'(「)kp『)=0.t:yy: 2■r) 'r-—^o.t ：Z:Z : 2推出22:xft.yrt2pu、):2:xft.yrt2pu、):2=0=0“•二；z t=0砂\ +nrtP—=C(t) 能量方程:2单纯重力作用时，二-gz-gz・L=c(t).:t T2②第二种情况恒定流■U恒定流■Ux：：Uy:Uz = = =0;:t;:t ;:t(■::X■2pU于2(Uyz:y（二(■::X■2pU于2(Uyz:y（二■2pU)=2(Uzx2一Uxz).:z(■:)=2(Uxy积分:dxdydzdxdydzd(二UxUyUz如果行列式等于零，则「 2'数C，与t无关单纯重力作用下理想液体恒定流的能量方程2, pU■■:二gz时，gz CP2满足条件式1-4其中的一条，就可以满足行列式为零的要求2.1.•，x=片=入=02.—=――=—,这是一个流线的微分方程式，所以适用于同一个流线之上。UXUyUzdxdydz3.'x，y ，z，这是一条涡线的方程式，所以适用于同一个涡线上各点。Ux_巴_Uz4<x ，y ，z，恒定流中以流线与涡线相重合为特征的螺旋流，所以方程式也适用于恒定螺旋流中。根据以上论述，对于由上述1、2、3、根据以上论述，对于由上述1、2、3、4各种情况所限定的两个点上(点1及点2),则2y2g=Z22P2U22-y2g(2—34)沿流线「②沿涡线各项含义：位能、压能和动能。总的意义:液体在运动过程中，尽管其机械能可以相互转化,各项含义：位能、压能和动能。总的意义:液体在运动过程中，尽管其机械能可以相互转化,但总的机械能是守恒的，上式是普遍能量守恒原理在理想液体中的表现形式。对于实际液体:1 .:px-P;ufta2CU丄()2(U/ y:x2_Uy，z)/Uy.：t2cu丄()-2(Uxz:y2_Uz‘X)2DU丄()2(U/-x2■UX■y)积分形式势流222x2x(二.7-(.pu+——2■) Uz迄p2推出CJIp2::u.:ux—()-()-八Ux=2(Uzy—Uy・・z):X_=xP:-Xjt2CJI■y+EUy(£)2.7 2=2(Ux・z—Uz，x)2p；U ：UzL2—-—(.)一—( ) Uz=2(Uyz-Ux・y):z:z :z2 ：：t-d■- )P2-（竺dx；t和y OUz ―2 ，—2 _2+ dy+ dz)+十(可uxdx+可uydy uzdz)=2Ct CtdxUxdydz(l2UxdxF2uUzdz)--js沿流线:2pu-d( )P2AUx=ucos(u,x),A.dx=dscos(u,x)dx■Ux=u-dsdxdsUxUy'yUz：Ux.:tdx■Uy亍J竺dz.i-Uxds;tUxU;tds;tdsUzuds[

u'(u.i2222.：udsu2)ds(R为摩擦所作的功)COTOC\o"1-5"\h\zp u 刃 cR„(R为摩擦所作的功)CO，•”一d(兀+二十——)一——ds——ds=0p 2 acs-2.-dpu’ 1cu—(z h)ds ds=0js 2g'gft任意两点:2pu 1 ：-U乙 h ds2g1g7Zi2■u12g=Z22P2.u22gh'i2hi(为惯性消耗，hi为惯性水头)第八节流体流动的初始条件和边界条件粘性流体的基本方程是二阶偏微分方程， 现在的任务就由原本结果具体的流体问题,变成了解决这个粘性流体基本方程的解的问题。方程的解不仅需要给出描述流体运动的方程，而且非常重要的是，需要给出流动的初始条件(时间变量)和边界条件(空间变量)一、初始条件：方程组的解在初始时刻应满足的条件。在初始时刻 t二t0，给出：Ux(x,y,z;t。)=Uxo(x,y,zfUy(x,y,z;t°)=Uyo(x,y,z)Uzx,y,z;to二Uzox,y,zpx,y,z;to=pox,y,z:■、边界条件：在流场的边界，方程组的解应满足的条件。在固体接触面上(ux,uy,uz)=(ux,uy,uzW 粘性，无滑移边界条件(ux,Uy,uz)=0 固壁静止Ux,Uy,Uzfnf〕.Ux,Uy,Uzwn 理想流体，有滑移边界条件不同液体的分界面：两侧液体的速度和压强保持连续。Vf1—vf2，Pf1二Pf2