流体动力学基本方程Word文档下载推荐.docx
- 文档编号:18372761
- 上传时间:2022-12-15
- 格式:DOCX
- 页数:30
- 大小:96.17KB
流体动力学基本方程Word文档下载推荐.docx
《流体动力学基本方程Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《流体动力学基本方程Word文档下载推荐.docx(30页珍藏版)》请在冰豆网上搜索。
2)对于不可压缩球对称流动(如三维空间中的点源产生的流动)
则有4r2V(r,t)m(t),
2
即V(r)r2,其中m(t)代表点源强度(单位时间发出的流体
体积)。
例2、均质不可压缩流体(密度为)从圆管(半径为R)入口端以
速度V0流入管内,经过一定距离后,圆管内流体的速度发展为抛物型剖面,即
。
通
常称这种流动为圆管的入口流。
试求当管内流动发展为抛物型剖面时的最大速度Vm。
解:
如图,将整个入口段取为控制体,对不可压缩流体有:
界面VvdSv0,由于管壁无渗透故上式
2R
可写为:
V0R2V2rdr,可得Vm2V0。
II动量方程流体团所受合外力=该流体团的质量
II-1方程的导出
1直角坐标系下推导微分形式的动量定理
t时刻,考虑一个正六面体形状的流体微团,如图所示,该流体微团刻所占控制体CV,其边界CS。
受力分析:
体力合力=Fd
面力合力
x
px
y,z
sx
p
xx
x2
y
py
x,y
2y,z
yx,y
sy
z
pz
x,y,z
zx,y,z
sz
x,y,z
vdV
vF
vpx
vpy
vpz
CSpvndS
即dVv
分量形式:
dvx
Fx
dvy
dvz
pxy
pyy
pzy
pxx
pyx
pzx
pxxpyx
Fy
或写成
dvi
Fipji
xj
或
F
P。
Fz
P意义:
单位体积流体团所受面力的合力。
2积分形式的动量定理的导出
考虑体系,该流体团t时刻所占控制体
CV,
其边界CS。
由动量定理有
ddt
vVd
Fd
CSvpndS
利用输运定理可得dV
V
vv
VV
vS。
dttCV
于是得到积分形式动量定理:
VtCV
S
Fdò
CSpndS
该定理的应用:
经常应用于求流体与边界的相互作用力。
例题1求流体作用于闸门上的力。
(设渠宽w)
取控制体如图所示,根据假定只需讨论动量方程的
x方向分量方程。
x方向动量通量
V12wD1
V22wD2
x方向合外力
D1
wPa
g(D1
y)dy
D2
0wPag(D2y)dy(hD2)Pa
闸门受合力=R
Pa(D1
h)
代入动量方程方程得
22
w(V2D2V1D1)
2gw(D1
故
122Rgw(D1D2)
w(V2D1V1D2)
D
注:
求R时可直接设Pa
0。
注微分形式的动量定理也可由积分形式的动量定理导出,推导过程如下:
vv
dvdVdVvdVVdtdtdtdt
ddm
其中ddtddtm0,因而得到
dvdVdV
dtdtCVdt
上式表明:
流体团总动量的变化率
=组成该流体团的流体质点的动量变化率之和。
另外,乙CSpvndS
nvPdS
综上可得
dV
Pd,
0,再考虑到系统大小形状的任意性可得
vdVdt
尽管得到了流动的动量方程,但是不像经典力学有了动量定理就可以求解质点运动一样,流体运动的动量方程中应力张量等于什么我们还不知道,并且速度的随体导数同时包含空间导数和时间导数,使得我们不仅需要初始条件,还需要边界条件才能确定一个具体流动。
3兰姆—葛罗米柯形式的动量方程
VvV2vvv
rotVVFP
t2
II-2地转参照系下的动量方程就很多空间和时间尺度都较小的流动而言,地球参照系通常课近似看作惯性系。
但是对于大尺度的流体运动问题,必须考虑地球自转的影响。
在海洋和大气的大尺度运动问题中,通常把地心看成惯性参照系,地球相对于地心有自转运动。
我们在此介绍地转参照系下的动量方程,为将来学习物理海洋学、地球流体动力学等打基础。
地球上运动质点的绝对速度VvaVvrVve,其中Vvr代表质点相对于地球表面的运动速度,牵连速度
vvvv
Vevrv(牵连速度=地球表面上该质点所在位置绕地心的自转速度),为地球自转角速度。
绝对加速度:
wvawvrwvewvc,
vvdv
其中wvr代表相对加速度,牵连加速度wveddt
vrv,科氏加速度wvc2vVr。
动量方程:
dVr
wvewvc
其中dVr
i
xi
因为真实力与参照系无关,
般情况下可以忽略地球自转角速度的变化,
认为
0,于是有
VrVvrVvr
vvrv2vVvr。
III.能量方程
III-1能量方程的推导:
t时刻流体团所占控制体CV,其边界CS,能量平衡关系式:
(UV22)
FV
V2
U
=
U+
乙CSpvnVv
vs=n
P
vVs
乙CSnvP
乙CSkT
s
k
T
所以得到能量方程微分形式
:
UV
因为
V2ddV2
U=U
2dtdt2
VsPVs?
PV
CSCS
vvv
FV(PV)(kT)q,
v其中(PV)
pjiVipji
(pjiVi)VipjiiVipjisji
xxxx
jjjj
pjiaji。
由于旋转运动张量
A是反对称张量,而应力张量
P是对称张量,故有pjiaji0(因pij是对称张量)
记pjisjipijsji
pjiv
P:
S。
另外VijiV(
P),于是有如下形式的能量方程:
d(U
V22)
FVV(P)P:
(kT)q。
t时刻系统能量增加率
1
外力的功率
其中
单位时间内通过边界流入的热量
3
单位时间内从外界吸收的其他能量
4
(1)
(U
),U代表单位质量流体的内能(分子热运动动能+分子间相互作用势能)
(2)
CSpnVs
(3)
vò
CSf
vs
kTsv,fv为热流强度,根据付利叶热传导定律对各向同性流体
vf
(4)设单位时间内单位质量流体从外界吸收的辐射能为q,则(4)qd
故能量方程积分形式为:
乙CSpvnVvsCSkTsvCVq
CSnCSCV
方程中各项意义分析:
代表单位体积流体能量变化率;
FvVv代表作用在单位体积流体微团上的体力的功率;
Vv(P)代表作用在单位体积流体微团表面的面力的合力的功率;
(kT)q代表单位时间内单位体积流体微团通过热传导和辐射吸收从外界获得的能量。
III-2动能方程
将动量方程
dVFv
P两边同时点积Vv得:
vdVvvvVdVFVV(P)。
dtF
其中VvdVv
1d(VV)
1dV,故有动能定理
2dt
dV2vv
V(P)。
dt2
单位体积流体微团动能变化率=作用于该微团上的体力的功率+作用于该微团上的合面力的功率。
III-3热流量方程:
dUP:
S(kT)q
面力的功率包含两项VvP(P:
S),其中合面力的功率Vv(P)转化为系统的宏观运动动能,
另一部分P:
S转化为系统的内能。
尽管系统内部的应力是内力,但是粘性应力必然导致机械能的耗散。
如果系统要维持定常状态,必须有外力对系统做功,补充其机械能损耗。
参考本章后面的例题。
o-xyz到旋转后的坐标系o-xyz,二阶
IV.本构方程数学预备:
记VvE,根据二阶张量定义,将坐标系旋转,从原坐标系张量E的张量元满足变换:
Vj
im
jn
其中变换矩阵
vivjvvkvkvvivjv
逆变换:
Vn
minj
xm
本构方程的导出
1应力张量分解:
Pp
P——偏应力张量,代表运动流体的应力张量与各向同性应力张量(记为p)的差异。
记作PP是对称二阶张量。
2线性假设(Newton粘性定律的推广,对于剪切流动u1kx2,21u1)
x2偏应力产生于速度场的不均匀性。
线性假设:
假设偏应力张量各分量与速度梯度张量的各分量成线性关系:
uk。
ijcijkl。
xl
cijkl是四阶张量,满足变换关系cijklimjnkplqcmnpq。
cijkl是由81个系数组成的一组系数,这组系数确定了偏应力张量各张量元与速度梯度张量各张量元之间的关系,由于偏应力张量和速度梯度张量都满足二阶张量定义,于是有
cuqcul
ijimjnmnimjncmnpqimjncmnpqkplq
xpxk
可知cijklimjnkplqcmnkl。
数学上定义,由81个元素组成的量,若其元素满足该变换的则称之为四
阶张量。
3各向同性流体及其四阶张量cijkl的表达式
3-1各向同性流体:
若在原坐标系
o-xyz和旋转后的坐标系o-xyz中偏应力张量分别表示为
cijklul和ij
ulcijkl,若
ul
ul则应当有ij
ij,于是要求cijklcijkl。
xk
考虑一个特例来理解流体粘性的各向同性:
水池中插入并移动平板引起的两个纯剪切流动的粘性应力大小与平板放置方向无关。
只要加上一个速度梯度,就对应一个粘性应力,粘性系数与速度梯度的方向无关。
u1
21c21kl
c2121
3-2对于各向同性流体,可以证明(参见吴书p75)四阶张量cijkl可表示为
cijkl
ijklikjliljk,其中,,是标量,即
whenijklwhenijklwhenikjl
wheniljk
3-3偏应力张量是对称张量cijkl
cjikl,于是
于是
ijkl
ikjliljk。
另外,由上式还可知cijklcijlk。
4分解
akl,于是
ij
cijklsklcijklaklcijklskl
偏应力与变形运动相关联。
如果流体只有旋转运动而没有变形运动,那么偏应力张量=
5将cijkl的表达式带入上式,得
ijij
klskl
ikjliljkskl
skkij
2sij
最后得到:
pijpij
1s
2s
pij
3skkij
其中sij
代表无体积变化的纯剪切运动,
skki
j代表各向同性膨胀运动。
3kkij
6Stokes假设
对于不可压缩流体,
skkij=0。
对于可压缩流体
skk
ij表示流体发生膨胀或收缩时引起的法向应力,
被称为第二粘性系数或膨胀粘性系数。
Stokes假设:
系统处于准热力学平衡状态时,可近似认为0。
7的意义
考虑纯剪切运动,uky,粘性应力212s21
u,可知
为动力学粘性系数。
8p的意义设流体满足Stokes假设,可以证明
p作用于球形微团上的法应力的平均值。
So,it'
sameasureofthelocalintensityofthe“squeezing”ofthefluid.
roverall
证明:
Theaveragevalueofthenormalcomponentofthestressonasurfaceelementatpositiondirectionsofthenormalntotheelementis
pnnnPnd
nipijnjsindd.
Sincenv(sin
cos
sinsin
cos
),
pnn
ni
pijnjsin
dd
pij
ijpiip
或者在球坐标系下
vn
(1,0,0)
ver
nn
evrP
versin
prrp
Hence,pcharacterizingthefluidpressureinamovingfluidwhichisanalogoustothestaticfluidpressureinthesensethatit'
sameasureofthelocalintensityofthe‘squeezing'
ofthefluid.
(关于p与热力学压强的关系,建议学生查庄礼贤《流体力学》对应章节。
)
9关于偏应力张量P
Ageneralrelativemotionnearanypointmayberepresentedasthesuperpositionoftwosimpleshearingmotion,eachofwhichgivesrisetoatangentialstressdeterminedbyandthecorrespondingvelocitygradient,togetherwitharigidrotationandanisotropicexpansion,neitherofwhichhasaneffect(inafluidofisotropic
structure)onthenon-isotropicpartofthestress'
tensoranddij
ij)mayofcauseberegarded
3ij
sijandasymmetrical
1)平板上的切应力212s21
u0,平板所受总阻力h
221bl2blu0。
21h
astheonlypossiblelineartensorialrelation,involvingonescalarparameter,betweentensordijwhosediagonalelementshavezerosum.(以上8和9)引自Batchlor,1994)本构方程(广义牛顿公式)的适用范围:
1)大多数液体;
2)非高温、非高频振动的气体;
非牛顿流体:
油漆、橡胶、蜂蜜、血液、沥青等。
例1写出纯剪切流动偏应力张量各分量
例2吴书p203,23
2)y3h/2处流体内摩擦力为0。
例3吴书p203,22柱坐标系下应力张量的表达式见p190。
除przpzru外,应力张量其他非对角元均为零。
ru2管壁处的切应力przu2cr0,单位长圆管对流体的阻力prz2r04cr02。
与圆管共轴的半径为r
r0/2的单位长流体柱表面的总摩擦力
cr2。
V流体力学基本方程组
V-1完备的微分形式流体力学基本方程组
(V)
0,
P,
dU
(k
T)q,
PpI
IdivVIdivV,(
pp
T.
内能
UUT,V,
具体函数形式由热力学理论给出。
对于完全气体
UCvT。
V-2N-S方程
将PpI2(S
s3kkI)代入动量方程即得:
Fvp
[2(S
ΔI)],其中Δ
skk。
当流场温度变化不大时,
近似为常数,
故有
[2
(S
I)]
V,
uj
ui
V。
xixi
最后得到
2Vv
V)。
又,若流体不可压缩,方程化为
N—S方程:
2Vv。
又,若流体粘性可略,方程化为理想流体Euler方程:
10
V-3耗散函数
耗散函数——单位时间内粘性导致的单位体积流体机械能转化成的内能。
P:
Spijsji
sij
vdivV
jsji
3i
psii
2sijsji
(si
i)2p
v2v2v
V2S:
S(V)2pV
v其中pV为压缩功,
而2
:
S23(
v2
V)2为粘性力的功,它将导致机械能转化成内能。
定义耗散函数
2S
(V)2,
它等于单位时间内由于粘性应力做功导致的机械能转化成的
内能。
它可以化成如下形式:
4(s12s13
3[(s11s22)
(s11
s33)(s22
s33)]。
可见,恒大于或等于零。
这说明粘性力做功总是使机械能转化成内能,这个过程不可逆。
例题:
拖动上板引起的剪切运动uky,k
①写出该流动的耗散函数。
常数。
设平板面积间距h,忽略边缘效应,
Sk
0,
Vv0,
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 流体动力学 基本 方程