- 1、本文档共6页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
实际流体的动量传递一纳维-斯托克斯方程.doc.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)
式中v=%/p—
文档评论(0)