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

1、第七章不可压缩流体动力学基础在前面的章节中,我们学习了理想流体和粘性流体的流动分析,按照水力学的观点,求得平均量。但是,很多问题需要求得更加详细的信息,如流速、压强等流动参数在二个或三个坐标轴方向上的分布情况。本章的内容介绍流体运动的基本规律、基本方程、定解条件和解决流体问题的基本方法。第一节流体微团的运动分析运动方式:移动或单纯的位移(平移旋转线性变形角变形。位移和旋转可以完全比拟于刚体运动,至于线性变形和脚变形有时统称为变形运动则是基于液体的易流动性而特有的运动形式,在刚体是没有的。在直角坐标系中取微小立方体进行研究。 一、平移:如果图(a所示的基体各角点的质点速度向量完全相同时,则构成了

2、液体基体的单纯位移,其移动速度为z y x u u u 、。基体在运动中可能沿直线也可能沿曲线运动,但其方位与形状都和原来一样(立方基体各边的长度保持不变。二、线变形:从图(b 中可以看出,由于沿y 轴的速度分量,B 点和C 点都比A 点和D 点大了dyyu y ,而yu y 就代表1=dy 时液体基体运动时,在单位时间内沿y 轴方向的伸长率。xu x ,yu y ,zu z 三、角变形(角变形速度d d d d DCABCDBAdtyu dydt dy yu d x x= dtxudxdt dx xu d yy=+=-d d 2d d -=角变形: +=+=-=x u yu d d d y

3、xz 212 +=x u z u z x y 21+=y u z u z y x 21 四、旋转(旋转角速度-=-=y u xu x y z 21-=z u y u y zx21 即, -=x u z u z x y 21zyxu u u z y x k j i =21那么,代入欧拉加速度表达式,得:z x x x xx xz y y z z y y y y y y y x z z x x z z z z zz z y x x y y x x y du u u u u u u u dt t xu u u u u u u u dt t y u u u u u u u u dttz =+-=+-=

4、+-各项含义:(1 平移速度(2线变形运动所引起的速度增量 (3(4角变形运动所引起的速度增量 (5(6微团的旋转运动所产生的速度增量流体微团的运动可分解为平移运动,旋转运动,线变形运动和角变形运动之和。 亥姆霍兹速度分解定理第二节 有旋运动1、无涡流(势流如在液体运动中,各涡流分量均等于零,即0=z y x ,则称这种运动为无涡流。当满足无涡流条件时,y zx z y x u u y zu u z x u u x y =,满足柯西条件,就有:x y z u x u y u z =存在。即流速势。满足此条件的流动(无涡流就叫势流。(下一章作详细介绍2、有涡流:如在液体运动中,涡流分量x 、y

5、及z 中间的任一个或全部不等于零,则这样的液体运动就叫做旋流或有涡流。自然界中的实际液体几乎都是这种有涡的流动。涡线:流场中一些假想的线,在所讨论的瞬时,涡线上各个质点的涡旋向量都与此线在该点处相切。 y与流线同样的分析方法,得到涡线方程:zyxdzdydx=涡量:设流体微团的旋转角速度为(t z y x ,则k j i z y x+=2称为涡量,是与空间坐标和时间有关的矢量函数。其中x 、y 和z 是涡量在x 、y 、z坐标上的投影。根据旋转角速度的定义,有: zu yu y z x -= xu zu z x y -= yu xu x y z -=哈米尔顿算子是一矢量算子,k zj y i

6、x +=,可知,-+ -+ -=yu x u j x u z u i z u y u u u u z y x k j i u x y z x y z zyxu=那么,(0=u就自然满足。 或者写成,0=+zyxz y x即涡量的定义使之自然满足涡量连续性微分方程。例:已知某圆管(半径0r 中液体流动的流速分布为:(22204zy rJ u x +-= 0=yu0=z u试判断该流动是有涡流还是无涡流?并求涡线微分方程。021=-=z u y u y zx z J x u z u z x y -= -=421y Jy u x u x yz = -=421 所以,该流动是有涡流。将上三式代入涡线微

7、分方程,zy x dzdydx=,得:yJ d z zJ d y 44=-0=+z d z y d y 积分后,得到: C y z =+22涡线是和管轴同轴的同心圆。涡管:在涡量场中任意画一封闭曲线,通过这条曲线上的每一点所做出的涡线构成一管状的曲面,称为涡管。涡通量:设A 为涡量场中一开口曲面,微元面dA 的外法线单位向量为n,涡量在n方向上的投影为n ,则面积积分+=Az y x An Ad x d y d z d x d y d z dA A d J称为涡通量。有旋运动的一个重要的运动学性质:在同一瞬间,通过同一涡管的各截面的涡通量相等。 证明:我们知道,根据涡量的定义,可以很容易知道,

8、涡量自然满足涡量连续性微分方程,即:0=+zyxz y x ,对这个微分方程在任意封闭体积上作积分,也是满足的,若任意体积取为,一段涡管和两个截面A1和A2,就有:0=+dV zyxvz y x可以将体积分化成封闭曲面积分:+321A A A z y x d x d yd x d z d y d z +=321A z y x A z y x A z y x dxdydxdz dydz dxdydxdz dydz dxdy dxdz dydz+=321A A A A d A d A d其中03=A A d(+=21A A dA n dA n021=+-=A n A n dA dA所以,=21A

9、n A n dA dA 得证对于微元涡管,近似认为截面上各点的涡量为常数, 2211A A =性质:涡管不可能在流体内部开始或终止,而只能在流体中自行封闭成涡环,或者终止于和开始于边界面。龙卷风开始于地面,终止于云层。速度环量:在流场中任取一封闭曲线s ,则流速沿曲线s 的积分:+=sz y x sdz u dy u dx u s d u称为曲线s 上的速度环量,并规定积分沿s 逆时针方向绕行为s的正方向。 (一斯托克斯定理 根据斯托克斯公式,+=sz y x ss dz u dy u dx u s d u -+-+ -=A xy z x y z d x d y y u x u d x d z

10、 x u zu d y d z z u y u A AAz z y y x x J A d dA dA dA =+=性质:沿任意封闭曲线s 的速度环量等于通过以该曲线为边界的曲面A 的涡通量。斯托克斯定理。(二汤姆逊定理汤姆逊定理:在理想流体的涡量场中,如果质量力具有单值的势函数,那么沿由流体质点所组成的封闭曲线的速度环量不随时间而变,即:0=dtd y解释:速度环量=涡通量,所以,流体的涡旋具有不生、不灭的性质。第三节 不可压缩流体连续性微分方程1、微分形式的连续性方程在推导这个方程式时,我们认为运动着的液体系连续地充满它所占据的空间,流动时不形成空隙,并且表征液体运动的各物理量也都是时间和

11、空间的连续函数。在时间t ,于流场中取一具有边长为dx 、dy 、dz 的微分六面体,在随后的一无限小段dt 内,流进和流出该微分六面体的质量。流出-流入=质量增量。微分六面体形心A 点的坐标为(x 、y 、z ,密度为,质点的速度分量为x u 、y u 及z u ,则在dt 时段内沿x 轴从左侧面abcd 流入六面的液体质量为xdx -2 12dx xdx +2dydzdt dx x dx x u u x x -221( 流出的液体质量为:dydzdt dx x dx x u u x x +221( 质量的变化:dxdydzdt tdxdydz dxdydz dt t =-+联立,得到:(d

12、xdydzdt tdxdydzdt zu dxdydzdt yu dxdydzdt x u z y x =-0=+zu yu xu tz y x (一般形式的液体连续性方程适合可压缩和不可压缩液体。或,写成: 0=+z u y u xu dt d zy x 0=+zu yu xu z y x (适合不可压缩液体,恒定流和非恒定流它是质量守恒定律在水力学中的表现形式。它表征着不可压缩液体在运动时,若保持其连续性,则线性变形必系伸长现象与缩短现象同时发生。 2、积分形式的液体连续性方程 连续性方程写成矢量形式:(0=+u t 其中为微分算子。体积积分:(0=+d ut v 根据高斯公式,(0=+A

13、 vdAn u d t对于恒定流,(0=AdA n u对于不可压缩,(0=AdA n un 是液体边界的外法线方向考虑到速度和面积的方向,就可知:02211=+-dA u dA u ,即,2211dA u dA u = (微小流束的流量平衡积分后,可以得到,2211A v A v = 其中v 、2v 为各自断面上的断面平均流速。 例:判断,流速为:22yx y u x +-=,22yx x u y +=,0=z u 的流动是否满足连续性方程。解: 22yx y u x +-=,22yx x u y +=,0=z u ,那么2222yxxyxu x += 2222yxxyyu y +-=0=zu

14、 z满足0=+zu yu xu z y x ,所以,满足连续性方程。第四节 以应力表示的粘性流体运动微分方程式一、粘性流体的内应力表面力,9个分量: zzzyzx yz yyyx xz xy xxp p p X 方向:x xma F=dydzdx xp p dx xp p dxdydz X xx xx xx xx 2(2(-+d x d y dz zdz zdzdx dy ydy yzx zx zx zx yx yx yx yx 2(2(2(2(-+-+dtdu dxdydz X x =dtdu zyxp X x zx yx xx =+或,dt du z y xp X xzx yx xx=+(

15、11 同理:dt du z y yp Y y zyxy yy=+(11 dtdu z y xp Z zxz yz zz=+(11 其中,密度对于不可压缩流体是已知常量,单位质量力Z Y X ,也是已知常量;未知量为9个应力和三个速度分量。不容易求解。第五节 应力和变形速度的关系一、切应力和角应变速度的关系一元流动的牛顿内摩擦定律为:dydu = 或可写为,dtd = 切应力与流速梯度或直角变形速度的关系。dtd 是直角变形速度,它是角变形速度的2倍,在xoy 平面上,xu yu dtd y x z +=2那么,对于三元流动的牛顿内摩擦定律,可以写成如下形式: += += +=z u y u x

16、 u zu x u y uy zzyyz z xxz zx y xyx xy 六个切应力均可用粘性系数和直角变形速度的乘积来表示。二、法向应力和线变形速度的关系在理想流体中,同一点各方向的法向应力相等,0,-=p p p p p zz yy xx (代表是压应力。在粘性流体中,粘性不仅产生与切应力有关的角变形速度,而且使线变形速度zu y u x u zy x ,也产生附加法向应力。使一点的法向应力与作用面方位有关。 取边长为dy dx =的方形流体微团进行研究,先考虑方形微团在x 方向上的伸长变形。 微团在x 方向上作伸长变形时,BC 伸长为C B ',而对角线AC 旋转至C A &

17、#39;,使产生角变形,这样在AC 面上产生切应力n ,线变形xu x 所产生的n 就要有力来平衡,这样就在AB 面上产生了附加法向切应力xx 。根据这样的分析,就有n 和xx 力的平衡,dx dx n xx 245cos 0= n xx 2=dxC C ACAB d d 245sin sin 0'=其中,C C '是微团在x 方向上伸长的长度,dxdt xu C C x ='那么,化简后,x u dtd x =21 其中,d 是45度角的角变形速度,xu dtd xn =2那么,附加法向应力和线变形速度的关系: xu x xx =2综合一下,得到,=z u y u x

18、 u z zz y yy x xx 222线变形运动使法向应力随伸长变形而减小,就有,+-=+-=+-=z u p p y u p p x u p p z t zzy t yy x t xx 222这就是粘性流体法向应力和线变形速度的关系。其中,t p 为理想液体的压强,它的大小与作用面方位无关。在粘性流体中,法向应力应该与方向有关了,所以,定义任意一点上三个相互垂直平面上的法向应力的平均值的负值为粘性流体在该点的压强。(+-=+-=z u y u xu p p p pp zx t zz xxy yy3231 对于满足连续性方程的不可压缩流体,0=+zu yu xu z y x ,所以,t p

19、 p =。而对于可压缩流体,zu yu xu z y x +代表的是质点的体积膨胀率,与坐标的选择无关,而压强p 是空间坐标的函数,与方向无关。 +-+-=z u y u x u x u p p z y x x xx 322 +-+-=z u y u xu y u p p z y x y yy 322 +-+-=z u yu xu zu p p z y x z zz 322 第六节 N-S 方程将xx p 的表达式代入以应力表示的粘性流体运动微分方程,得到:dt du x u z u z x u yu y z u y u x u x u p x X x z xy xzx x = + + +-+

20、-+113221y若流体粘滞性是常数,那么,zx u yx u xu xu xpX z y x x -+-22222232323221dtdu zx u zu yx u yu x z x y x =+222222 = >zx u yx u xu xu x pX z y x x +-2222223131311dtdu zu yu x x x =+2222(31122zu yu xu xxu xpX z y x x +-dtdu zu yu x x x =+2222对于不可压缩流体,有连续性方程,0=+zu yu xu z y x ,所以,dtdu zu yu xu xpX x x x x

21、=+-(1222222同理,可得,=+-=+-=+-dt du z u y u x u z p Z dt du z u y u x u y p Y dtdu z u y u x u x pX z z z z y y y y x xxx(1(1(1222222222222222222不可压缩粘性流体的运动微分方程 它与连续性方程0y =+zu yu xu z x 联立,求得速度的三个分量和压强。上式中,dtdu x 是流体质点的加速度,对于欧拉的描述方法, zu u yu u xu u tu dtdu x zx y x x x x +=+=+-+=+-+=+-z u u yu u xu u tu

22、 zu yu xu zpZ z u u y u u x u u t u zu yu xu y pY zu u y u u x u u t u zu yu xu x pX z zz yz xz z z z y z y y y x y y y y x zx y x x x x x x (1(1(1222222222222222222第七节 理想流体运动微分方程及其积分当流体为理想液体时,运动粘滞系数0=,n-s 方程就简化为:111x y z du pX x dt du pY y dt du pZ z dt -=-=-= (2-23将式(2-13代入,则上式化为111xxxx xyzy y y y

23、 x y z z z z z x y z u u u u pX u u u x t x y z u u u u p Y u u u y t x y z u u u u p Z u u u z t x y z -=+-=+-=+(2-24 再引用和的表达式,则又可化为1(1(1(x x xy z z y y z z y y y y z x x z z x x z z z z x y y x x y y x u u pX u u u u u xtxu u pY u u u u u y t y u u pZ u u u u u z t z -=+-=+-=+-(2-25上式为欧拉方程式用整体运动性质

24、(位移、变形及旋转所表现的形式,说明液体在各力作用下而运动时,表现有位移、变形及旋转作用的可能。对于不可压缩的理想液体,欧拉方程及连续方程提供出解决运动问题的四个独立条件。一般而论,可用以解决其中四个未知数,例如x u 、y u 、z u 及p 。 二、欧拉-哥罗米柯方程式下面我们再把欧拉方程变成具有涡旋分量的另一形式。由流速2222x y z u u u u =+可求得 2(2y x z xyzu u u uu u u xxxx=+今将此式与欧拉方程式(2-24中的第一式相减则得21(2(22y x x x z z y x z y y zpuX x x u u u u u u u t z x

25、xyu u u t-=+-=+-因得21(2(2x z y y z u puX u u x tx -=+-同理 221(2(21(2(2yx z zx z y x x yu puY u u y t y u p u Z u u z t z -=+-=+-(2-26 今设体积力(X 、Y 、Z 为有力势的力,亦即体积力分量可由其势能=(x 、y 、z 确定如下:,X x=-,Y y=-Z z=-将之代入式(2-26,稍加整理则得222(2(2(2(2(2(2xy z z y y z x x z z x y y x u p uu u t x u p u u u t y u p uu u t z +=-+=-+=-(2-27 式(2-27即是在有势力作用下用涡旋分量表示的理想液体的运动方程式,又称欧拉-哥罗米柯方程式。三、理想液体伯努力积分形式的能量方程(25分钟 由于水流是在重力场中运动,重力为有势力所以,存在一势能函数xX -=,yY -=,zZ -=水为不可压缩的液体,为常数(1p x xpi =故上式可写为 (22(2y z z y x u u up x t u -=+(22(2z x x z y u u upy t u