实际流体的动量传递一纳维-斯托克斯方程.doc

实际流体的动量传递一纳维-斯托克斯方程 实际流体具有粘性，故也称为粘性流体。与冶金机械欧拉方程的推导思路类似，实际流体的运动方程也可从牛顿第二运动定律导出，式(1-2-69)仍然适用。所不同的是，实际流体所受到的表面力除法向力外，还有切向力，对单位面积而言，则分别称为法应力和切应力。切应力即流体流动时的粘性应力。从量纲的角度考虑，应力与动量通量相当，因此上述表面应力均可表示为!，其下标i表示动量的传递方向，j表示流体的流动方向，即动量的方向。 下面进行流体的受力分析。 图1-2-23为直角坐标系下实际流体中一微元六面体，其边长分别为dx，dy和dz。此微元体的每个面上都受到与之毗邻的来自外部流体的表面力的作用，这种表面力表现为：由流体单位元体线变形而产生的法向力和由流体粘性引起的切向粘性力。图1-2-23示出了微元体在x轴方向的受力情况。 图1-2-23微元体受力示意图微元体在x轴方向所受质量力为： pXdxdydz (1-2-80) 微元体在x轴方向所受法向力为左、右两侧面所受法向力之代数和，即 [(!**+^fdx)-]dydz=^Mxdydz 在x方向所受切向粘性力为上、下和前、后四个面上的4个x方向之切向粘性力的代数和，即 [(!+x+管办)-!yx]dxdz+[(!zx+#fdz)-!zx]dxd+=(#^+)dxd+dz因此，微元体在x轴方向受到的总表面力为： (#f+賢+#f)dxdydz (1-2-81) 微元体沿x轴方向总的净动量率为： []dxdydz (1-2-82) #x #y #z 微元体内动量率的积累为： jOUx !t ■dxdydz (puxux)(puxu%)!x+!y 因此，根据牛顿第二运动定律[式1-2-69]，有将上式右端的各项展开，并运用连续性方 程，可得.+丄(+!#yx+!#zx)=!Ux O\3x!y!z/!t 同理，有 !(puxuz)!(pux)!z+!t (1-2-83) (1-2-84) !ux!ux!uxDux !T+UxaT+UyaT+Uz + Y+- z+- !x!y!z/ !xy!yy ) !x!y!z/ (1-2-84a) (l-2-85a) (l-2-85b) (l-2-85c) 办1!z/—!t!x,i!x?z!z—Dt DuxDt DUy ：Dt !yz )_Duz ,!x+!y+!z)Dt 式(1-2-85)即为以应力形式表示的实际流体的运动方程。 对不可压缩流体，上述三个运动微分方程中只有密度和质量力的三个分量X，Y，Z为己知，尚有9个应力分量和3个速度分量共12个未知数，加上连续性方程也只有四个方程，欲求解12个未知数是不可能的，所以还必须设法找出这些未知量之间的其他关系。 对一维流动的牛顿型流体，其切应力与法向速度梯度的关系可用牛顿粘性定律来描述。当流体作三维运动时，情况要复杂得多，每一切应力都与其作用面上两个方向的速度梯度有关，其关系为(推导过程从略）： #xy_#yz_$/9ux9u%\(+!x)1-)-86a）#yz_#zy_$/!Uy!uz、 \!z!y)1-)-86b）#zx_#zx_^(!z?+axz)(1-2-86c) 而法应力则由如下两部分组成:一部分由流体的静压力产生，其结果使流体微元承受压缩应力而发生体积形变;另一部分由粘性应力作用产生，其结果是使流体微元在法线方向上承受拉伸或压缩应力，发生线性变形。各法应力与静压力和速度梯度之间的关系如下：。!ux #xx_-p+2$ax-■2 (!ux3u2\ 3 $!x!y!z/1-)-87a）。!uy tyy_-9+2$!%-■2 /!ux!!3u2\ 3 $!x+!y+!z/1-)-87b）。！uz #zz_-p+2$az-■2 /!ux!Uy!uz、 3 $!x+!y+!z/1-)-87c）将方程式（1-2-86)和式(1-2-87a)代入方程式(1-2-85a)，得 一丄!92!2!x2(a2!!2Uy^Uz) Dto!x+%!x23%!x2+!x!y+!x!zj /!2uxc!u^/ c!u2\ +v(!%2+!x!y)+v(!z2+axaz) 整理后，可得 Du# Dt 1p /+ux 2x 丄1 u u, ( ) -!#+v(X2*)+孓(#+式+)(1-2-88a) 同理，有 11_ol- Y -- 丨lly+2Uy++)+1 X2Ty2T 5(1+,+)d-2-88b) y\xyz/ Du^v1p/2z2zuz、1(2uUyUz\ { 、 :z-—+v|^+^+^r)+yv(#+亥+) (1-2-88c) Dt 将以上三式写成向量形式为： Du%FDtFg .丄2.p+#22,u+\v$(2.u) (1-2-89