《高等流体力学》习题集学习资料
《高等流体力学》习
题集
《高等流体力学》复习题
一、 基本概念
1. 什么是理想流体?正压流体,不可压缩流体? [答]:教材P57
当流体物质的粘度较小,同时其内部运动的相对速度也不大,所产生的粘性应力比起其它类型的力来说可以忽略不计时,可把流体近似地看为是无粘性的,这样无粘性的流体称为理想流体。
内部任一点的压力只是密度的函数的流体,称为正压流体。
流体的体积或密度的相对变化量很小时,一般可以看成是不可压缩的,这种流体就被称为不可压缩流体。
2. 什么是定常场;均匀场;并用数学形式表达。
[答]:如果一个场不随时间的变化而变化,则这个场就被称为定常场。其数学表达式为:
)(r ??=
如果一个场不随空间的变化而变化,即场中不显含空间坐标变量r ,则这个场就被称
为均匀场。其数学表达式为:)(t ??=
3. 理想流体运动时有无切应力?粘性流体静止时有无切应力?静止时无切应力是否无粘
性?为什么?
[答]:理想流体运动时无切应力。
粘性流体静止时无切应力。但是,静止时无切应力,而有粘性。因为,粘性是流体的固有特性。
4.流体有势运动指的是什么?什么是速度势函数?无旋运动与有势运动有何关系?[答]:教材P119-123
如果流体运动是无旋的,则称此流体运动为有势运动。
φ的速度梯度来表对于无旋流动来说,其速度场V总可以由某个速度标量函数(场)),(t r
φ称为速度场V的速度势函数。
示,即φ
r
V,则这个标量函数(场)),(t
=
?
无旋运动与有势运动的关系:
势流运动与无旋运动是等价的,即有势运动是无旋的,无旋运动的速度场等同于某个势函数的梯度场。
5.什么是流函数?存在流函数的流体具有什么特性?(什么样的流体具有流函数?)[答]:
6.平面流动中用复变位势描述的流体具有哪些条件(性质)?
[答]:教材P126-127
理想不可压缩流体的平面无旋运动,可用复变位势描述。
7. 什么是第一粘性系数和第二粘性系数?在什么条件下可以不考虑第二粘性系数?Stokes
假设的基本事实依据是什么? [答]:教材P89
第一粘性系数μ:反映了剪切变形对应力张量的贡献,因此称为剪切变形粘性系数; 第二粘性系数μ’:反映了体变形对应力张量的贡献,因而称为体变形粘性系数。 对于不可压缩流体,可不考虑第二粘性系数。
Stokes 假设的基本事实依据:平均法向正应力ε就是压力函数的负值,即体变形粘性系数03
2
=+='λμμ。
8. 从运动学观点看流体与固体比较有什么不同? [答]:教材P55
若物质分子的平均动能远小于其结合能,即E mv ?<<221
,这时物质分子间所形成的对偶
结构十分稳定,分子间的运动被严格地限定在很小的范围内,物质的分子只能在自己的平衡位置周围振动。这时物质表现为固态。
若物质分子的平均动能与其结合能大致相等,即E mv ?≈221
,其分子间的对偶结构不断
地遭到破坏,又不断地形成新的对偶结构。这时,物质分子间不能形成固定的稳定对偶结构,而表现出没有固定明确形状的液态。
若物质分子的平均动能远大于其结合能,即E mv ?>>221
,物质几乎不能形成任何对偶结
构。这时,物质表现为气态。
9. 试述流体运动的Helmholts 速度分解定律。 [答]:教材P65
可变形流体微团的速度分解:流体微团一点的速度可分解为平动速度分量与转动运动分量和变形运动分量之和,这称为流体微团的Helmholts 速度分解定理
r S r V V δδ?+?+=0
10. 流体微团有哪些运动形式?它们的数学表达式是什么? [答]:V δδω?+?+=0 1)平动运动:0V = 2)转动运动:δω? rot 2
1
=ω 3)变形运动:δ?
11. 描述流体运动的基本方法有哪两种?分别写出其描述流体运动的速度、加速度的表达
式。 [答]:教材P58-60
描述流体运动的基本方法:
1)拉格朗日方法:对流体介质的每一质点进行跟踪,着眼于流体介质中的每个质点,需要
对流体介质中的每个质点进行区别。 各质点速度表达式:t
t c b a t c b a ??=
)
,,,(),,,( 各质点加速度表达式:2
2)
,,,(),,,(t t c b a t c b a V
??=?
2)欧拉方法:定点观察描述流场的运动,着眼于空间的定点,而不是流体质点。
速度表达式:332132321213211321),,,(),,,(),,,(),,,(),(e t x x x u e t x x x u e t x x x u t x x x V t r V V ++===
加速度表达式:V V t V V t x u u t u V t t t dt d j i j i )(??+??
=??+??=???+??=+??=???+??=
12. 什么是随体导数(加速度)、局部导数(加速度)及位变导数(加速度)?分别说明
0=dt v d ,0=??t v 及()0=??v v
的物理意义? [答]:教材P60
随体导数:流体质点在其运动过程中的加速度所对应的微商,叫做随体导数; 局部导数:流体位置不变时的加速度所对应的微商,叫做局部导数;
位变导数:质点位移所造成的加速度所对应的微商,叫做位变导数。
物理意义:0=dt v
d :随体导数为0,流体质点在其运动过程中的加速度为0;
0=??t v
:局部导数为0,流体位置不变时的加速度为0,流体是定常流动; ()0=??v v :位变导数为0,流体质点位移所造成的加速度为0,流体速度分布均匀。
13. 什么是流体的速度梯度张量?试述其对称和反对称张量的物理意义。 [答]:教材P65-67
对流体微团M ,其中o r 处的速度为0V ,那么处的速度可以表示为 j j
x x V δ??+
=0,或者j j i i i x x u u u δ??+
=0, 即)(0V ??+=δ。这里,x u
j
i ?=??为二阶张量,是速度的梯度,因此称之为速度梯度张量。
速度梯度张量分解为对称和反对称部分:S A x u
V i
j +=??=?
反对称张量的物理意义:
反对称张量表征了流体微团旋转运动,所对应的矢量ω为流体微团的角速度矢量。
k ijk z v y w z u x w z v y w y u x v z u x w y u x v A ωεωωωωωω=???
?
?
??---=???????
?
?
?
???-??-??-??-??-????-??-??-????-??=0000) (21) (21) (210
) (21) (21)
(2101
2
1323
V rot e e e z y x 2
1
321=++=ωωωω
对称张量的物理意义:
对称张量表征了流体微团的变形运动。其中,对角线上的元素()321 , , εεε表示了流体单元
微团在3个坐标轴上的体变形分量,而三角元素??? ??32121 ,21 ,2
1
θθθ表示了流体单元微团在3个
坐标平面上的角变形分量的一半。
???
????
? ??=????????
?
?
?????+????+????+??????+????+????+????=31212
32312
121212
1212
1
) (21) (21) (21) (21) (21) (21 εθθθεθθθεz w z
v y w z u x w z v y w y v y u x v z u x w y
u x v x u
A y
u
??x
w ??-z
v ??-反对称部分
Z z
?
14. 流体应力张量的物理意义是什么?它有什么性质? [答]:教材P71
流体应力张量的物理意义:
应力张量表示了坐标面的三个面力密度矢量z y x p p p
, ,的九个分量}{ij p 组成的一二阶张
量,即为面力密度张量。
应力张量的性质:应力张量是对称张量,具有对称性 应力张量具有二阶对称张量的性质
(1) 应力张量的几何表示为应力椭球面,即二次型
1222)(222=+++++=??zx p yz p xy p z p y p x p r P r zx yz xy zz yy xx
(2) 应力张量有三个互相垂直的主轴方向,即是应力椭球的三个对称的直径的方向。在主轴坐标系下,应力张量具有标准形式:
????
? ?
?='0
00'0
00'332211p p p P (3) 应力张量的三个不变量为:
反对称部分
15. 某平面上的应力与应力张量有什么关系?nm mn p p =的物理含义是什么? [答]:教材P71
应力n p 与应力张量P 的关系:P n p n p ij n ?=?= ,即:空间某点处任意平面上的应力等于这点处的应力张量与该平面法向单位矢量的左向内积。
nm mn p p =的物理意义:
i ji j j ji i j ij i n nm n p m m p n m p n m p m P n p ===?=??=
)(
mn m p n p n P m =?=??=
)(
应力张量的对称性,使得在以n 为法线的平面上的应力n
p
在 m
方向上的投影等于(=)在以m 为法线的平面上的应力m p 在 n
方向
上的投影。
16. 流体微团上受力形式有哪两种?它们各自用什么形式的物理量来表达? [答]:教材P68-71
(1)质量力,也称体力,这种力作用在物质中每个质点上,其大小与每个质点的质量成正比。作用于某物质体上质量力的合力将通过该物质体的质心。
δτρδ)(f = , ?=τ
δτρ)(r F f )(为质量力密度,与位置有关。
(2)面力,作用于流体微团表面S 上的力。
S p n δδ= , ?=S
n S p δ n p 为面力分布密度,P p p ij n ?=?=
17. 什么是广义的牛顿流体和非牛顿流体?
?????+--++=---++=++=22311212332312232211331231233221132
12231223221111333322233
22111p p p p p p p p p p p p p p p I p p p p p p p p p I p p p I
[答]:教材P86-87
牛顿内摩擦定律:流体微团的运动变形的的大小与其上所受的应力存在线性关系。 遵从或近似遵从牛顿内摩擦定律的一类流体称为牛顿流体。不遵从牛顿内摩擦定律的流体称为非牛顿流体。
广义牛顿内摩擦定律:偏应力张量的各分量与速度梯度张量的各分量间存在线性关系。 遵从或近似遵从广义牛顿内摩擦定律的一类流体称为广义牛顿流体。
18. 试述广义牛顿内摩擦定律的物理意义及相应的数学表达式? [答]:教材P87
广义牛顿内摩擦定律的物理意义:偏应力张量的各分量与速度梯度张量的各分量间存在线性关系。
数学表达式:lk ijkl lk ijkl k
l
ijkl
ij a c s c x v c +=??=τ,其中,二阶张量lk s 和lk a 市速度梯度张量的对称和反对称部分,而四阶张量ijkl c 称为动力粘性系数张量。
19. 什么是层流运动、紊流(湍流)运动和临界雷诺数?圆管中层流和紊流运动的速度分
布规律是什么? [答]:
层流流动是平稳有规律的流动状态,流体介质各部分之间分层流动,互不掺混,流体内部的微团具有连续而平滑的迹线,流场中各种有关物理量(参数)的变化较为缓慢,表现出明显的连续性和平稳性。
湍流流动是极不规则的流动形态,流体介质各部分之间,各层之间有着剧烈的掺混,其流体内部微团的运动迹线很不规则,杂乱无章,表征流体运动状态的各种物理量也表现出不同程度的跃变和随机性。
雷诺数:流体运动中,惯性力与粘性力的无量纲比值 μ
ρν
vd
vd
=
=Re 下临界雷诺数:从湍流状态到层流状态的转折点; 上临界雷诺数:从层流状态到湍流状态的转折点。 圆管中层流和紊流运动的速度分布规律: 层流:)(422
0r R l
p p u l --=
μ (1) 定常流动的速度沿径向的分布规律,由式(1)可以看出,流动截面上的速度分布是一抛物回转面。 湍流:光滑圆管中的速度分布:
394.5)lg(756.5**+=ν
yU U u
粗糙圆管中的速度分布与光滑圆管中的速度分布相同,只是改变方程的常数。
20. 流体的阻力可分为哪几种?管路中的阻力通常分为哪几种?
[答]:粘性时产生阻力的根本原因,依据阻力产生的不同机理,可分为:摩擦阻力和压差阻力。
管路中的阻力通常分为:沿程阻力(即摩擦阻力)和局部阻力。
21. 试说明粘性流体流动的三个基本性质。 [答]:教材P170-174 (1)粘性运动的有旋性
粘性流体运动时,有旋是绝对的,粘性流体的无旋运动是不存在的。
(2)运动过程中有能量的损耗性
在粘性流动中永远伴随着机械能的损耗。这部分能量转换成热能形式传递给流体介质及相邻的固壁,使其温度升高而耗散。 (3)粘性涡旋运动的扩散性
在粘性流体中,涡旋强的地方要向涡旋弱的地方传送涡量,直至涡量相等为止。
22. 使流体涡量产生变化的因素有哪些?其中哪些是流体运动的内在因素,哪些是外在因
素?
[答]:流体涡量产生变化的因素有:(1)质量力无势;(2)流体不正压;(3)粘性剪切应力;(4)流体微团的体积变化;(5)流体涡线微元的变形(涡线的拉伸、压缩、扭曲)。
其中,流体运动的外在因素为:(1)质量力无势;(2)流体不正压;(3)粘性剪切应力。
内在因素为:(4)流体微团的体积变化;(5)流体涡线微元的变形(涡线的拉
伸、压缩、扭曲)。
23. 试说明层流边界层和湍流边界层的速度分布特征。 [答]:层流边界层:层流边界层内的速度分布呈线性分布规律;
湍流边界层:分为层流底层和湍流核心区。层流底层内的速度分布呈线性分布,湍流核心区速度分布呈对数分布规律。
24. 试述雷诺应力j i u u ''-ρ的物理意义及其与分子粘性应力的异同。 [答]:教材P230
雷诺应力j i u u ''-ρ的物理意义:在湍流运动中,由脉动速度引起的应力,称之为雷诺应力。 雷诺应力与分子粘性应力的异同:
相同:都是由于分子动量传递产生的应力,都是剪切应力。
不同:(1)引起动量传递的原因不同(雷诺应力:分子脉动;分子粘性应力:分子热运动);
(2)分子粘性应力与粘性这一物质固有属性有关,而雷诺应力取决于流体的流动特性,与流场性质有关,与所处位置和时均速度有关。
25. 试述平板湍流边界层的结构及其速度分布特征。 [答]:教材P241-242
结构:沿壁面法向,在板面附近有层流子层流区,其速度呈线性分布(01.0=<δ
y
),而后
为很小的过渡区,接着为湍流核心区。
结构:层流子层流区 过渡区 湍流核心区; 内层:粘性底层 过渡区 湍流核心区; 外层:粘性顶层及边界层其余部分。 速度分布特征:
层流子层流区(80*
≤<
ν
yU ):ν
*
*yU U u x =,速度呈线性分布; 过渡区:308*
≤<
ν
yU 湍流核心区(ν
*
30yU <):
9.4)lg(6.5**+=ν
yU U u x ,速度呈对数分布。
二、 推导及证明
1. 根据质量守恒定律推导连续性方程。
[证明]:教材P78-79
根据物理学中的质量守恒定律,由某封闭的物质面S 所围成的体积τ中的物质在运动过程中不消灭也不创生,即使说,在运动过程中由物质面S 所围成的体积τ中的流体介质的质量保持不变,是守恒的。
在体元素δτ中,若流体介质的密度为ρ,那么其质量就为ρδτδ=m ,于是有限体积τ中的质量m 为
?=τ
ρδτm (1)
根据质量守恒定律的物理含义:体积τ中的质量m 在其运动过程中保持不变,这意味着,
质量m 的随体导数为零,即
0)(==?τ
ρδτdt d
dt dm (2) 由物质体元素的随体导数表达式 V dt d ??=δτδτ 知 δτρρ
ρδτττ
????+=)()(V dt d dt d (3)
于是由式(2)有 0)(=??+?δτρρ
τV dt
d (3)
即 0])([=??+???δτρρ
τV t (4)
考虑到奥-高公式(??=S
V
n V div S u δ)有
0????=+??=?+??S v t S t n S δρδτρδρδτρτ (5) 式(3)到式(5)都可称之为积分形式的连续方程。
由式(3)和式(4)的被积函数为零可直接得到微分形式的连续方程:
0=??+dt d ρρ
(6) 0(=??+??V t ρρ
(7) 2.
根据动量定律推导出微分形式的运动方程。
[证明]:教材P80-81
封闭曲面S 所围成的体积τ中流体物质体的动量为体积分:?τ
δτρV
,
其变化率就是体积分的随体导数:
?τ
δτρ)(V dt d
, 而该物质体τ上所受外力为其上的质量力:?τ
δτρF 和面力:S P n S p S
S
n δδ
???= , 由动量定理得:??????+=+=s
s n S P n F S p F V dt d
δδτρδδτρδτρτττ ,
因为
δτρρ
δτρδτ
ρδτρρδτρδτρδτρτττττττ)()()()()(V dt
d V dt V d V V dt
d V dt V d dt d V V dt d V dt d
??++=??++=+=??????? 由连续性方程知,0)(=??+?δτρρτV dt
d V
,
所以??????+=+=s
s n S P n F S p F dt V
d δδτρδδτρδτρτττ
,
又????+??=??+??=s
n S v t V t V dt V d δρδτρδτρδτρτττ)( 得到S p F S V v t V
S n s
n δδτρδρδτρ
ττ
????+=+??
由奥-高公式????=?τ
δτδP S P n S
所以,δτδτρδτρτττP F dt
V
d ??+=???
于是得到微分形式的动量方程P F dt
V
d ??+=
ρρ 3.
根据能量守恒定律推导出微分形式的能量方程。
[证明]:教材P83-85
4.
试推导出运动方程的Bernoulli 积分和lagrange 积分。
[证明]:教材P108-109 5.
在不可压缩流体中,若流线是11c f =和22c f =两曲面的交线。试证明:
()()2121,f f f f F V ???=
,其中F 是1f 和2f 所决定的函数。 [证明]:设 构成曲线坐标系,于是 满足:
由题设:流线是 两曲面的交线,那么速度场 的方向将同时垂直于
的梯度方向。因此: 于是速度场可以表示为: 即要证明F 不显含 f 3 :
321,,f f f 321,,f f f 0
)
,,()
,,(321321≠??x x x f f f 11c f =22c f =V
21 ,f f 2
1//f f V ???
))(,,(21321f f f f f F V ???= 03
=??f F
)])(,,([21321=?????=??f f f f f F V
而
又
所以, 也就是说F 中不显含 f 3。于是,有:()()2121,f f f f F V ???=
。 此题得证。
6.
证明不可压缩理想流体作二维定常流动时,忽略质量力,其流函数ψ和涡旋Ω满足
()()
0,,=?Ω?y x ψ,若Ω为常数,则压力方程为=Ω++ψρ22v p 常数。 [证明]:由: P F V V t V ??+=?Ω+?+??ρ1)2( 2
理想、定常、忽略质量力 p V V ?-=?Ω+?ρ1
)2(2
0)2(2=?Ω++?V p
V ρ
两边取旋度
0)])2([2=?Ω++???V p
V ρ
0)()()()()(=??Ω-Ω??+Ω??-??Ω=?Ω??V V V V V
不可压 0)(=??ΩV 蜗旋场无源 0)(=Ω?? V 二维流动 0)(=??ΩV 0)(=Ω??
V 0=Ω??V 0=?Ω
??+?Ω??
y
v x u )( ),,()(),,(2132121321f f f f f F f f f f f F V ?????+?????=??
)()()(211221=?????-?????=?????f f f f f f )(),,(21321f f f f f F V ?????=??
)()(
21332211f f f f F f f F f f F ???????+???+???=0)(2133
=???????=f f f f F 3
3231
3322212
3
121
1
1
x f x f x f x f x f x f x f x f x f ??????????????????=?????)(213f f f 0
)
,,(),,(321321≠??=x x x f f f 03
=??f F
0=?Ω
????-?Ω????y x x y ψψ 0)
,()
,(=?Ω?y x ψ
由: 0)2(2=?Ω++?V p
V ρ
(1)
)()()(x y x y y x z e x e y e v e u e v e u e V ??+??Ω=-Ω=+?Ω=?Ωψψ
由于,Ω为常数:)()()(ψψψΩ?=?Ω?+?Ω?=
?Ωx y e x e y V
代入(1),得: 0)2(2=Ω++?ψρp
V
两端积分,得:
C p
V =Ω++ψρ
22 7.
进行圆管中流体摩擦试验时,发现圆管中沿轴向的压降p ?是流速u 、密度ρ、粘性
系数μ、管长l 、管内径d 及管壁粗糙度d
h
k ?=的函数,而且p ?与l 成正比。试用因次分析方法证明2
2
1u d l p ρλ
=?,其中()Re ,k λλ=为无因次系数。 [证明]:由题意可假设存在关系 γβαρλu d l k p 1 Re),(=? (1)
相应各量的量纲(因次)为:[]2]][[][T L M p =? αα][][L d = β
βρ??????=3][L M γ
γ
??
????=T L u ][ 式(1)对应量纲的协调条件为:γγβαβT L M T L M -+-+=][][][][][][31-2-11 于是,对于M 量纲,有: 1=β
T 量纲,有: 2-=-γ 2=γ
L 量纲,有: 131-=+-+γβα 1-=α
将:1-=α 1=β 2=γ 带入(1)式,得:2
2
1u d l p ρλ
=?
此题得证。 8.
试从运动方程:P F dt V d ??+=
ρρ 和本构关系)3
1
(2I V S pI P ??-+-=μ 推导出:粘性不可压缩流体的运动方程为: V p F dt V d
?+?-=νρ
1
如果体力有势即G F -?= 则有: Ω?=??Ω-Ω
νV dt
d )( [证明]:(1)将本构关系带入运动方程: ) 31
(2I V S p F dt V d
??-??+?-=μρρ 考虑到不可压缩流体 0=??V 上式为:S p F dt V
d ??+?-=νρ
21
V x u x x u x u x x u x x u x u x S i i j i j j i i i j i j i i j i
?=????+??=????+????=??+????=??2
1)(21)(21)(2122
V p F dt V d
?+?-=νρ
1
(2)考虑到体力有势:V p G dt V d
?+?--?=νρ
1
)()2()(V V V
V t V V V t V dt V d
???-??+??=??+??= V V p
G V t V
?=Ω?-++?+??νρ
)2(2 两边取旋度:
V V t
V
???=Ω???-????ν)(
Ω??-??Ω-??Ω+Ω??=Ω??? )()()()()(V V V V V 旋度无源:0)(=Ω?? V , 不可压0)(=??ΩV
所以,Ω??-??Ω=Ω???
)()()(V V V
Ω?=Ω??+??Ω-?Ω
?
ν)()(V V t
由于dt
d V t =Ω??+?Ω
?
)( 所以,可得:Ω?=??Ω-Ω
νV dt
d )( 证毕。 9.
证明对粘性不可压缩流体定常运动,若外力有势ξ,则有:
222)2
1
)(1(ζξρν=++??-?V p s V
其中s 为沿流线的弧元素,ζ为涡量,ν为运动粘性系数,V 为流体速度,p 为压力函数,
ρ 为密度。
[证明]:Ω??-+-?=Ω??-?-=?Ω+?+??
νρ
ξνρ)(1)2(2p
p F V V t V
Ω??-=++?+?Ω νρ
ξ)2(2p
V V (1)
在(1) 两边同以流线切线方向的单位向量s e
作左向内积:
)()2()]2([22Ω???-=++??=++?+?Ω?
V V V p s p V V e s νξρρξ
)()2
(2
Ω???-=++?? V V p s V ξρν (2) (1) 两边同求散度:
0)]2([2=++?+?Ω??ρ
ξp V V
0)2()(22=++?+?Ω??ρ
ξp V V