高等流体力学
第一章 流体力学的基本概念
连续介质:流体是由一个紧挨着一个的连续的质点所组成的,没有任何空隙的连续体,即所谓的连续介质。
流体质点:是指微小体积内所有流体分子的总和。 欧拉法质点加速度:a x =
du x ∂u x ∂u ∂u ∂u
=+u x x +u y x +u z x dt ∂t ∂x ∂y ∂z
a y =
du y dt
=
∂u y ∂t
+u x
∂u y ∂x
+u y
∂u y ∂y
+u z
∂u y ∂z
a z =
du z ∂u z ∂u ∂u ∂u =+u x z +u y z +u z z dt ∂t ∂x ∂y ∂z
d
表示。在dt
质点的随体导数:质点携带的物理量随时间的变化率称为质点的随体导数,用欧拉法描述中的任意物理量Q 的质点随体导数表述如下:
dQ ∂Q ∂Q
=+u k dt ∂t ∂k
式中Q 可以是标量、矢量、张量。质点的随体导数公式对任意物理量都成立,故将质点的
随体导数的运算符号表示如下:
d ∂∂=+u k dt ∂t ∂k
其中
∂∂
称为局部随体导数,u k 称为对流随体导数,即在欧拉法描述的流动中,物理量∂t ∂k
的质点随体导数等于局部随体导数与对流随体导数之和。 体积分的随体导数:
d ∂Φ⎡∂Φ⎤⎡d Φ⎤
()ΦdV =dV +Φu dS =+div Φv dV =+Φdivv dV v v s n v v ⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰⎢⎥⎢⎥dt ∂t ⎣∂t ⎦⎣dt ⎦
d ∂a ⎡∂a ⎤⎡da ⎤ ()adV =dV +au dS =+div av dV =+adivv v ⎰⎰⎰v ∂t s n ⎰⎰⎰v ⎢⎰⎰⎰v ⎢⎥⎥dV dt ⎰⎰⎰⎣∂t ⎦⎣dt ⎦
变形率张量: D ij =
11 ε12ε13
ε21 ε22 ε23
ε31 ε32 ε33
其中εii 表示所在方向的线性变形率,其余εij (i ≠j )为角变形率。D ij 为变形张量。
i ⎪ εij = + 2⎝∂x j ∂x i ⎪⎭
1⎛∂u ∂u j ⎫
旋转角速度: 0 -ωz R ij =
y ωz 0 -ωx -ωy ωx 0
1⎛∂u y ∂u x ⎫1⎛∂u z ∂u y ⎫1⎛∂u x ∂u z ⎫
⎪-ωz = ⎪ωy =2 ∂z -∂x ⎪ωx =2 ∂y -∂z ⎪⎪ 2 ∂x ∂y ⎝⎭⎝⎭⎝⎭
1⎛∂u y ∂u x ⎫1⎛∂u x ∂u z ⎫
⎪-判断有旋流和无旋流:ωx =ωy =ωz =0,ωz = =0,=-ω ⎪=0 y ⎪2 ∂x ∂y 2∂z ∂x ⎝⎭⎝⎭
ωx =
∂u y ∂u x ∂u x ∂u z 1⎛∂u z ∂u y ⎫∂u z ∂u y
⎪-==0 , ==
2⎝∂y ∂z ⎪∂x ∂x ∂y ∂z ∂y ∂z ⎭
涡量与速度环量的关系:涡量,流体力学中多用涡量来表示流体微团的旋转。定义旋转角速度的两倍为涡量,即Ωk =2ωk ;速度环量,速度沿封闭曲线的积分称为速度环量,通常用
Γ来表示,Γ=l v ⋅dl 。在笛卡尔坐标系下为Γ=l u x dx +u y dy +u z dy 。涡量与速度环
量的关系,数学表示如下:
⎰⎰
s
Ω⋅ndS =
l
v ⋅dl 。说明通过面的涡通量等于沿边界的速度
环量。 应力张量:
1、切应力的特性:切应力互等定律,即作用在两相互垂直平面且与该平面的交线相垂直的切应力大小都是相等的。表述如下:τxy =τyx ,τyz =τzy ,τzx =τxz
2、压应力的特性:压应力的大小与其作用面的方位有关,三个相互垂直方向的压应力一般是不相等的,即p xx ≠p yy ≠p zz 。但在几何关系上可以证明,同一点上,三个相互垂直面的压应力之和,与该组垂直面的方位无关,即p xx +p yy +p zz 值总保持不变。在实际流体中,任何三个相互垂直面上的压应力的平均值定义为动水压强,以p 表示,则
1
p =(p xx +p yy +p zz )
3
牛顿流体的本构方程:将应力张量σij 与变形张量εij 联系起来的方程称为本构方程 用张量的形式表示:σij =-p δij +2μεij σij =-p δij +μ 这就是不可压缩牛顿流体的本构方程。写成分量形式
⎛∂u i ∂u j
+ ∂x
⎝j ∂x i ⎫⎪ ⎪⎭
σ11
∂u y ∂u x ∂u z
=-p +2μ σ22=-p +2μ σ33=-p +2μ
∂x ∂z ∂y
⎛∂u z ∂u y ⎫⎛∂u x ∂u y ⎫⎛∂u x ∂u z ⎫ ⎪=μ ∂y +∂x ⎪⎪ σ13=σ31=μ ∂z +∂x ⎪ σ23=σ32=μ ∂y +∂z ⎪
⎝⎭⎝⎭⎝⎭
σ12=σ21
第二章 流体运动的基本方程
微分形式的连续性方程的表达式:
∂ρ∂(ρu i )∂ρ∂(ρu x )∂(ρu y )∂(ρu z )+=0 +++=0; ∂t ∂x i ∂t ∂x ∂y ∂z
不可压缩流体的确切定义,理解其含义:
d ρ∂ρ∂ρ =+u k
dt ∂t ∂x k
∂ρ
=0只是指密度是恒定不变的,但流体质点的密度换可以随流动中位置发生变化。只有∂t
满足上式,密度质点才能保持不变。 即
∂ρ
表明质点密度在时间上恒定不变。 ∂t
u k
∂ρ
表明质点的密度不随流动中位置的变化而变化 ∂x k
N-S 方程的各种表示形式:
(1)
du i 1∂p
=f i -+υ∇2u i dt ρ∂x i
dv 1
=f -∇p +υ∇2v dt ρ
(2)
(3)
dv 1
=f -gradp +υ∇2v dt ρ
(4)
∂u i ∂u 1∂p +u j i =f i -+υ∇2u i dt ∂x j ρ∂x i
du x 1∂p
=f x -+υ∇2u x dt ρ∂x
(5)
du y dt
=f y -
1∂p
+υ∇2u y ρ∂y
du z 1∂p
=f z -+υ∇2u z dt ρ∂z
∂u x ∂u x ∂u x ∂u x ∂2u x ∂2u x ∂2u x 1∂p
(6)+u x +u y +u z =f x -+υ(2+2+2)
∂t ∂x ∂y ∂z ρ∂x ∂x ∂y ∂z
∂u y
∂2u y ∂2u y ∂2u y 1∂p +u x +u y +u z =f y -+υ(2++2) ∂t ∂x ∂y ∂z ρ∂y ∂x ∂y 2∂z
∂u y
∂u y
∂u y
∂u z ∂u z ∂u z ∂u z ∂2u z ∂2u z ∂2u z 1∂p +u x +u y +u z =f z -+υ(2++2) 2∂t ∂x ∂y ∂z ρ∂z ∂x ∂y ∂z
流体的能量包括哪几种形式,并对各种形式解释,写出单位质量流体能量的表达式
运动流体的能量包括内能、动能和势能三种形式 内能是指分子运动的动能和分子间结合的能量,它随温度而变化。单位质量流体所含有的内能用e 1表示。
1v 22
若质量为∆m 的流体,其速度为v ,则动能为∆mv ,因此单位质量的动能e k =
22
势能来源于保守力场。一般情况下,作用于流场的保守力是重力场,因此流体的势能取决于位置的高度。设z 为某一个基准面以上的高程,则单位质量的势能可表示为e p =gz 则单位质量流体的能量方程可写为:
v 2
e= e1++gz
2
流体运动微分形式的基本方程组由哪些方程组成,通常有几个未知量,方程组是否封闭。 连续性方程、N-S 运动方程和能量方程。共12未知量,而方程组只有5个,因而不封闭的。对于不可压缩流体,如何求解速度场、压强场以及温度场,说明其求解步骤。
对于不可压缩均质流体,ρ为常数,则有连续性方程和运动方程即可求解v 和p, 然后再由能量方程求解温度场。
第三章 势流运动
势流运动控制方程及求解步骤:
∇2ϕ=0
∂ϕv 2p +++gz =f (t ) ∂t 2ρ
势流求解的常用方法
求解势流最常用的方法有流网法、势流叠加法、复变函数法以及数值计算法等。 速度势函数与流函数 势函数和速度的关系
u x =
∂ϕ ∂x
u y =
∂ϕ ∂y
流函数和速度的关系
u x =
∂ψ ∂y ∂ψ
∂x
u y =-
复势与复速度 f (z )=ϕ(x , y )+i ψ(x , y ) 它的实数部分是速度势函数,虚数部分是流函数 复势的导数为
df ∂ϕ∂ψ
=+i =u x -iu y dz ∂x ∂x
称复势的导数为复速度,其实数部分是x 向的分速度,其虚数为y 向的分速度的负值。
恒定平面势流的解析方法
1、以速度势函数ϕ为未知函数 2、以流函数ψ为未知函数
3、以复势f (z )为未知函数,其本身又含有三种方法:奇点法、镜像法和保角变化法 保角变化法的思路 设在z =x +iy 平面上一个复杂的流动边界,借助于某一解析变换函数一般为复势已知的典型流动(如圆柱绕流),ξ=g (z )变换到ζ=ξ+i η平面上另外的流动,
因为对于这些简单形状的物体在ξ平面上的解是已知的,则通过这种变换可以得到复杂图形的复势。然后再通过z =g
-1
(ζ),将ζ平面变换为z 平面。
第四章黏性流体运动
基本方程及求解途径
∂u x ∂u y ∂u z
连续性方程:++=0
∂x ∂y ∂z
运动方程:
ρ
du x ∂p =-+μ∇2u x dt ∂x
ρ
du y dt
=-
∂p
+μ∇2u y ∂y
ρ
du z ∂p
=-+μ∇2u z dt ∂z
求解途径:1、解析解2、近似解3、数值解 黏性流体运动的基本性质 1、黏性流体运动的有旋性 2、机械能量的损耗性 3、涡量的扩散性
黏性流体运动的解析解
∂2u x ∂2u x 1∂p
+==c
μ∂x ∂y 2∂z 2
两平行平板间的层流
u x =-
1dp 21y 1h -y 2+U +U
2μdx 2h 2
()
泊肃叶流
u x =-
1dp 2
h -y 2
2μdx
()
哈根-泊肃叶流
u x =-
1∂p 2
R -r 2
4μ∂x
()
小雷诺数流动近似解的思路
雷诺数小意味着黏性力对流动起主导作用,而惯性力则是次要因素,作为零级近似可以将惯性力全部舍去。如果是一级近似可以保留非线性惯性力项中的主要部分而将次要部分略去。 大雷诺数流动近似解的思路 求解大雷诺数问题的基本思路,把整个流场分为外部流体流动和边界层内流体流动。前面一部分流动属于理想流体的范围,运动方程为欧拉方程。第二部分流动属于黏性流体范围,运动方程为N-S 方程,由于边界层的厚度δ比特征长度小的多,而且x 方向速度分量沿法向的变化比切向大得多,所以N-S 方程在边界层内可以得到相当大的简化。简化后的方程称为普朗特边界层方程。 边界层的概念
边界附近必须考虑黏性作用的很薄流层称为边界层。 1、靠近边界 2、有涡
3、速度梯度大 边界层的厚度
名义厚度,通常定义当地流速u x (x , y )等于0.99U (x )时的y 值为边界层厚度δ,这样定义的边界层厚度也称名义厚度。
位移厚度,在固体壁面附近边界层中,由于流速受到壁面的阻止而降低,使得在这个区域所通过的流量较之理想流体的所通过的流量减少,相当于边界层的固体壁面向流动内部移动了一个距离δ1后理想流体流动所通过的流量。这个距离δ1称为边界层位移厚度。 边界层方程的相似性解的概念
当边界层方程具有相似性解时,其流速u x (x , y )分布具有以下性质:如果把任意x 断面的流速分布图形u x ~y 的坐标用相应的尺度化为无量纲坐标,则任意x 断面的流速分布图形均相同。
边界层分离现象
第五章 紊流运动
紊流的特征及其分类 特征:
1、不规则性 2、扩散性 3、连续性 4、耗能性 5、三维有涡性 分类:
1、各向同性均匀紊流
2、剪切紊流 又分为自由剪切紊流和壁面剪切紊流 壁面剪切紊流的发生过程及紊流结构 在靠近边壁很近的黏性底层中,平面上具有顺流向的低速带和高速带相间的带状结构。低速带随时均流动向下游移动时,其下游头部缓慢上举,与壁面间的距离逐渐增大,常形成横向漩涡。旋涡上下产生压差使漩涡顶着低速带上升,涡旋本身则变形成为马蹄形涡,头部上举后进入流速较大的流层,马蹄形涡发生拉伸变形。马蹄形涡的头部上举最终形成低速流体突然向上层的高速水流“喷射”,同时,高速水流乘机俯冲而入,这个过程称为“清扫”,清扫过后流速分布恢复正常,拐点消失。低速流体的喷射和高速流体的清扫是猝发现象的两个主要组成部分,清扫过后,又是新的低速带的出现而重复以上各个阶段。在发生猝发现象的地点,其下游将出现局部的紊流斑。紊流斑随主流向下游扩展,最后紊流部分占据了全部板宽,发展为充分发展紊流。
时间平均法和系综平均法的概念
时间平均法是将随机变量的瞬时值在一定的时段内进行平均。
系综平均法是是在同样条件下重复进行多次试验,任取其中足够多次的量测值做算术平均,所得的函数值具有确定性。
紊流运动方程—雷诺方程,雷诺方程的形式与N-S 方程的区别,雷诺应力项的意义
⎫∂u i ∂u 1∂p 1∂⎛ μ∂u i -ρu i ' u ' j ⎪+f i +u j i =-+
⎪∂t ∂x j ρ∂x i ρ∂x j ⎝∂x j ⎭
区别:1、雷诺方程中增加了一项-ρu i u j ,它代表了紊流对时均流动产生的影响,称为雷诺应力
2、物理量为平均值
紊流模型的用途,紊流模型通常有哪几类(零方程模型、一方程模型,二方程模型,其他模型)
应用时均紊流的连续性方程和雷诺方程解决紊流问题时,未知数共10个,远超过方程的数目。这就造成了时均紊流方程的不封闭。因此紊流模型是根据紊流的运动规律以寻求附加的
' '
条件和关系式,从而使方程组封闭可解。
零方程模型(涡黏性模型、混合长度模型、涡量传递模型) 一方程模型—k 方程模型 二方程模型—k-ε方程模型 紊动动能k
k =
1' ' u u 2i i
能量耗散率ε
∂u i ' ∂u i '
ε=ν
∂x j ∂x j
第六章 涡旋运动
涡旋的运动学性质
1、涡管中任一横截面上的涡通量保持为一常数,或者给定瞬间流入涡管的涡通量等于流出涡管的涡通量。
2、涡管不能在流体中产生或消失。 涡旋的形成
既然涡旋的不生不灭是在理想、正压、外力有势三个条件具备的时候才成立,那么涡旋运动的产生或消失的必然来源于这三个条件没有得到完全满足
黏性流体中涡旋的产生1、流场中存在固体边界2、两股流速不同的流体汇合所形成的间断面3、自由表面
非正压流体中涡旋的形成 外力无势所产生的涡旋运动
高等流体力学
第一章 流体力学的基本概念
连续介质:流体是由一个紧挨着一个的连续的质点所组成的,没有任何空隙的连续体,即所谓的连续介质。
流体质点:是指微小体积内所有流体分子的总和。 欧拉法质点加速度:a x =
du x ∂u x ∂u ∂u ∂u
=+u x x +u y x +u z x dt ∂t ∂x ∂y ∂z
a y =
du y dt
=
∂u y ∂t
+u x
∂u y ∂x
+u y
∂u y ∂y
+u z
∂u y ∂z
a z =
du z ∂u z ∂u ∂u ∂u =+u x z +u y z +u z z dt ∂t ∂x ∂y ∂z
d
表示。在dt
质点的随体导数:质点携带的物理量随时间的变化率称为质点的随体导数,用欧拉法描述中的任意物理量Q 的质点随体导数表述如下:
dQ ∂Q ∂Q
=+u k dt ∂t ∂k
式中Q 可以是标量、矢量、张量。质点的随体导数公式对任意物理量都成立,故将质点的
随体导数的运算符号表示如下:
d ∂∂=+u k dt ∂t ∂k
其中
∂∂
称为局部随体导数,u k 称为对流随体导数,即在欧拉法描述的流动中,物理量∂t ∂k
的质点随体导数等于局部随体导数与对流随体导数之和。 体积分的随体导数:
d ∂Φ⎡∂Φ⎤⎡d Φ⎤
()ΦdV =dV +Φu dS =+div Φv dV =+Φdivv dV v v s n v v ⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰⎢⎥⎢⎥dt ∂t ⎣∂t ⎦⎣dt ⎦
d ∂a ⎡∂a ⎤⎡da ⎤ ()adV =dV +au dS =+div av dV =+adivv v ⎰⎰⎰v ∂t s n ⎰⎰⎰v ⎢⎰⎰⎰v ⎢⎥⎥dV dt ⎰⎰⎰⎣∂t ⎦⎣dt ⎦
变形率张量: D ij =
11 ε12ε13
ε21 ε22 ε23
ε31 ε32 ε33
其中εii 表示所在方向的线性变形率,其余εij (i ≠j )为角变形率。D ij 为变形张量。
i ⎪ εij = + 2⎝∂x j ∂x i ⎪⎭
1⎛∂u ∂u j ⎫
旋转角速度: 0 -ωz R ij =
y ωz 0 -ωx -ωy ωx 0
1⎛∂u y ∂u x ⎫1⎛∂u z ∂u y ⎫1⎛∂u x ∂u z ⎫
⎪-ωz = ⎪ωy =2 ∂z -∂x ⎪ωx =2 ∂y -∂z ⎪⎪ 2 ∂x ∂y ⎝⎭⎝⎭⎝⎭
1⎛∂u y ∂u x ⎫1⎛∂u x ∂u z ⎫
⎪-判断有旋流和无旋流:ωx =ωy =ωz =0,ωz = =0,=-ω ⎪=0 y ⎪2 ∂x ∂y 2∂z ∂x ⎝⎭⎝⎭
ωx =
∂u y ∂u x ∂u x ∂u z 1⎛∂u z ∂u y ⎫∂u z ∂u y
⎪-==0 , ==
2⎝∂y ∂z ⎪∂x ∂x ∂y ∂z ∂y ∂z ⎭
涡量与速度环量的关系:涡量,流体力学中多用涡量来表示流体微团的旋转。定义旋转角速度的两倍为涡量,即Ωk =2ωk ;速度环量,速度沿封闭曲线的积分称为速度环量,通常用
Γ来表示,Γ=l v ⋅dl 。在笛卡尔坐标系下为Γ=l u x dx +u y dy +u z dy 。涡量与速度环
量的关系,数学表示如下:
⎰⎰
s
Ω⋅ndS =
l
v ⋅dl 。说明通过面的涡通量等于沿边界的速度
环量。 应力张量:
1、切应力的特性:切应力互等定律,即作用在两相互垂直平面且与该平面的交线相垂直的切应力大小都是相等的。表述如下:τxy =τyx ,τyz =τzy ,τzx =τxz
2、压应力的特性:压应力的大小与其作用面的方位有关,三个相互垂直方向的压应力一般是不相等的,即p xx ≠p yy ≠p zz 。但在几何关系上可以证明,同一点上,三个相互垂直面的压应力之和,与该组垂直面的方位无关,即p xx +p yy +p zz 值总保持不变。在实际流体中,任何三个相互垂直面上的压应力的平均值定义为动水压强,以p 表示,则
1
p =(p xx +p yy +p zz )
3
牛顿流体的本构方程:将应力张量σij 与变形张量εij 联系起来的方程称为本构方程 用张量的形式表示:σij =-p δij +2μεij σij =-p δij +μ 这就是不可压缩牛顿流体的本构方程。写成分量形式
⎛∂u i ∂u j
+ ∂x
⎝j ∂x i ⎫⎪ ⎪⎭
σ11
∂u y ∂u x ∂u z
=-p +2μ σ22=-p +2μ σ33=-p +2μ
∂x ∂z ∂y
⎛∂u z ∂u y ⎫⎛∂u x ∂u y ⎫⎛∂u x ∂u z ⎫ ⎪=μ ∂y +∂x ⎪⎪ σ13=σ31=μ ∂z +∂x ⎪ σ23=σ32=μ ∂y +∂z ⎪
⎝⎭⎝⎭⎝⎭
σ12=σ21
第二章 流体运动的基本方程
微分形式的连续性方程的表达式:
∂ρ∂(ρu i )∂ρ∂(ρu x )∂(ρu y )∂(ρu z )+=0 +++=0; ∂t ∂x i ∂t ∂x ∂y ∂z
不可压缩流体的确切定义,理解其含义:
d ρ∂ρ∂ρ =+u k
dt ∂t ∂x k
∂ρ
=0只是指密度是恒定不变的,但流体质点的密度换可以随流动中位置发生变化。只有∂t
满足上式,密度质点才能保持不变。 即
∂ρ
表明质点密度在时间上恒定不变。 ∂t
u k
∂ρ
表明质点的密度不随流动中位置的变化而变化 ∂x k
N-S 方程的各种表示形式:
(1)
du i 1∂p
=f i -+υ∇2u i dt ρ∂x i
dv 1
=f -∇p +υ∇2v dt ρ
(2)
(3)
dv 1
=f -gradp +υ∇2v dt ρ
(4)
∂u i ∂u 1∂p +u j i =f i -+υ∇2u i dt ∂x j ρ∂x i
du x 1∂p
=f x -+υ∇2u x dt ρ∂x
(5)
du y dt
=f y -
1∂p
+υ∇2u y ρ∂y
du z 1∂p
=f z -+υ∇2u z dt ρ∂z
∂u x ∂u x ∂u x ∂u x ∂2u x ∂2u x ∂2u x 1∂p
(6)+u x +u y +u z =f x -+υ(2+2+2)
∂t ∂x ∂y ∂z ρ∂x ∂x ∂y ∂z
∂u y
∂2u y ∂2u y ∂2u y 1∂p +u x +u y +u z =f y -+υ(2++2) ∂t ∂x ∂y ∂z ρ∂y ∂x ∂y 2∂z
∂u y
∂u y
∂u y
∂u z ∂u z ∂u z ∂u z ∂2u z ∂2u z ∂2u z 1∂p +u x +u y +u z =f z -+υ(2++2) 2∂t ∂x ∂y ∂z ρ∂z ∂x ∂y ∂z
流体的能量包括哪几种形式,并对各种形式解释,写出单位质量流体能量的表达式
运动流体的能量包括内能、动能和势能三种形式 内能是指分子运动的动能和分子间结合的能量,它随温度而变化。单位质量流体所含有的内能用e 1表示。
1v 22
若质量为∆m 的流体,其速度为v ,则动能为∆mv ,因此单位质量的动能e k =
22
势能来源于保守力场。一般情况下,作用于流场的保守力是重力场,因此流体的势能取决于位置的高度。设z 为某一个基准面以上的高程,则单位质量的势能可表示为e p =gz 则单位质量流体的能量方程可写为:
v 2
e= e1++gz
2
流体运动微分形式的基本方程组由哪些方程组成,通常有几个未知量,方程组是否封闭。 连续性方程、N-S 运动方程和能量方程。共12未知量,而方程组只有5个,因而不封闭的。对于不可压缩流体,如何求解速度场、压强场以及温度场,说明其求解步骤。
对于不可压缩均质流体,ρ为常数,则有连续性方程和运动方程即可求解v 和p, 然后再由能量方程求解温度场。
第三章 势流运动
势流运动控制方程及求解步骤:
∇2ϕ=0
∂ϕv 2p +++gz =f (t ) ∂t 2ρ
势流求解的常用方法
求解势流最常用的方法有流网法、势流叠加法、复变函数法以及数值计算法等。 速度势函数与流函数 势函数和速度的关系
u x =
∂ϕ ∂x
u y =
∂ϕ ∂y
流函数和速度的关系
u x =
∂ψ ∂y ∂ψ
∂x
u y =-
复势与复速度 f (z )=ϕ(x , y )+i ψ(x , y ) 它的实数部分是速度势函数,虚数部分是流函数 复势的导数为
df ∂ϕ∂ψ
=+i =u x -iu y dz ∂x ∂x
称复势的导数为复速度,其实数部分是x 向的分速度,其虚数为y 向的分速度的负值。
恒定平面势流的解析方法
1、以速度势函数ϕ为未知函数 2、以流函数ψ为未知函数
3、以复势f (z )为未知函数,其本身又含有三种方法:奇点法、镜像法和保角变化法 保角变化法的思路 设在z =x +iy 平面上一个复杂的流动边界,借助于某一解析变换函数一般为复势已知的典型流动(如圆柱绕流),ξ=g (z )变换到ζ=ξ+i η平面上另外的流动,
因为对于这些简单形状的物体在ξ平面上的解是已知的,则通过这种变换可以得到复杂图形的复势。然后再通过z =g
-1
(ζ),将ζ平面变换为z 平面。
第四章黏性流体运动
基本方程及求解途径
∂u x ∂u y ∂u z
连续性方程:++=0
∂x ∂y ∂z
运动方程:
ρ
du x ∂p =-+μ∇2u x dt ∂x
ρ
du y dt
=-
∂p
+μ∇2u y ∂y
ρ
du z ∂p
=-+μ∇2u z dt ∂z
求解途径:1、解析解2、近似解3、数值解 黏性流体运动的基本性质 1、黏性流体运动的有旋性 2、机械能量的损耗性 3、涡量的扩散性
黏性流体运动的解析解
∂2u x ∂2u x 1∂p
+==c
μ∂x ∂y 2∂z 2
两平行平板间的层流
u x =-
1dp 21y 1h -y 2+U +U
2μdx 2h 2
()
泊肃叶流
u x =-
1dp 2
h -y 2
2μdx
()
哈根-泊肃叶流
u x =-
1∂p 2
R -r 2
4μ∂x
()
小雷诺数流动近似解的思路
雷诺数小意味着黏性力对流动起主导作用,而惯性力则是次要因素,作为零级近似可以将惯性力全部舍去。如果是一级近似可以保留非线性惯性力项中的主要部分而将次要部分略去。 大雷诺数流动近似解的思路 求解大雷诺数问题的基本思路,把整个流场分为外部流体流动和边界层内流体流动。前面一部分流动属于理想流体的范围,运动方程为欧拉方程。第二部分流动属于黏性流体范围,运动方程为N-S 方程,由于边界层的厚度δ比特征长度小的多,而且x 方向速度分量沿法向的变化比切向大得多,所以N-S 方程在边界层内可以得到相当大的简化。简化后的方程称为普朗特边界层方程。 边界层的概念
边界附近必须考虑黏性作用的很薄流层称为边界层。 1、靠近边界 2、有涡
3、速度梯度大 边界层的厚度
名义厚度,通常定义当地流速u x (x , y )等于0.99U (x )时的y 值为边界层厚度δ,这样定义的边界层厚度也称名义厚度。
位移厚度,在固体壁面附近边界层中,由于流速受到壁面的阻止而降低,使得在这个区域所通过的流量较之理想流体的所通过的流量减少,相当于边界层的固体壁面向流动内部移动了一个距离δ1后理想流体流动所通过的流量。这个距离δ1称为边界层位移厚度。 边界层方程的相似性解的概念
当边界层方程具有相似性解时,其流速u x (x , y )分布具有以下性质:如果把任意x 断面的流速分布图形u x ~y 的坐标用相应的尺度化为无量纲坐标,则任意x 断面的流速分布图形均相同。
边界层分离现象
第五章 紊流运动
紊流的特征及其分类 特征:
1、不规则性 2、扩散性 3、连续性 4、耗能性 5、三维有涡性 分类:
1、各向同性均匀紊流
2、剪切紊流 又分为自由剪切紊流和壁面剪切紊流 壁面剪切紊流的发生过程及紊流结构 在靠近边壁很近的黏性底层中,平面上具有顺流向的低速带和高速带相间的带状结构。低速带随时均流动向下游移动时,其下游头部缓慢上举,与壁面间的距离逐渐增大,常形成横向漩涡。旋涡上下产生压差使漩涡顶着低速带上升,涡旋本身则变形成为马蹄形涡,头部上举后进入流速较大的流层,马蹄形涡发生拉伸变形。马蹄形涡的头部上举最终形成低速流体突然向上层的高速水流“喷射”,同时,高速水流乘机俯冲而入,这个过程称为“清扫”,清扫过后流速分布恢复正常,拐点消失。低速流体的喷射和高速流体的清扫是猝发现象的两个主要组成部分,清扫过后,又是新的低速带的出现而重复以上各个阶段。在发生猝发现象的地点,其下游将出现局部的紊流斑。紊流斑随主流向下游扩展,最后紊流部分占据了全部板宽,发展为充分发展紊流。
时间平均法和系综平均法的概念
时间平均法是将随机变量的瞬时值在一定的时段内进行平均。
系综平均法是是在同样条件下重复进行多次试验,任取其中足够多次的量测值做算术平均,所得的函数值具有确定性。
紊流运动方程—雷诺方程,雷诺方程的形式与N-S 方程的区别,雷诺应力项的意义
⎫∂u i ∂u 1∂p 1∂⎛ μ∂u i -ρu i ' u ' j ⎪+f i +u j i =-+
⎪∂t ∂x j ρ∂x i ρ∂x j ⎝∂x j ⎭
区别:1、雷诺方程中增加了一项-ρu i u j ,它代表了紊流对时均流动产生的影响,称为雷诺应力
2、物理量为平均值
紊流模型的用途,紊流模型通常有哪几类(零方程模型、一方程模型,二方程模型,其他模型)
应用时均紊流的连续性方程和雷诺方程解决紊流问题时,未知数共10个,远超过方程的数目。这就造成了时均紊流方程的不封闭。因此紊流模型是根据紊流的运动规律以寻求附加的
' '
条件和关系式,从而使方程组封闭可解。
零方程模型(涡黏性模型、混合长度模型、涡量传递模型) 一方程模型—k 方程模型 二方程模型—k-ε方程模型 紊动动能k
k =
1' ' u u 2i i
能量耗散率ε
∂u i ' ∂u i '
ε=ν
∂x j ∂x j
第六章 涡旋运动
涡旋的运动学性质
1、涡管中任一横截面上的涡通量保持为一常数,或者给定瞬间流入涡管的涡通量等于流出涡管的涡通量。
2、涡管不能在流体中产生或消失。 涡旋的形成
既然涡旋的不生不灭是在理想、正压、外力有势三个条件具备的时候才成立,那么涡旋运动的产生或消失的必然来源于这三个条件没有得到完全满足
黏性流体中涡旋的产生1、流场中存在固体边界2、两股流速不同的流体汇合所形成的间断面3、自由表面
非正压流体中涡旋的形成 外力无势所产生的涡旋运动