流体主要计算公式.docx
《流体主要计算公式.docx》由会员分享,可在线阅读,更多相关《流体主要计算公式.docx(30页珍藏版)》请在冰点文库上搜索。
![流体主要计算公式.docx](https://file1.bingdoc.com/fileroot1/2023-6/22/4cff4dcb-7802-4daf-af46-fdafacb3b2a9/4cff4dcb-7802-4daf-af46-fdafacb3b2a91.gif)
流体主要计算公式
主要流体力学事件有:
∙1738年瑞士数学家:
伯努利在名著《流体动力学》中提出了伯努利方程。
∙1755年欧拉在名著《流体运动一般原理》中提出理想流体概念,并建立了理想流体基本方程和连续方程,从而提出了流体运动解析方法,同时提出了速度势概念。
∙1781年拉格朗日首先引进了流函数概念。
∙1826年法国工程师纳维,1845年英国数学家、物理学家斯托克思提出了著名N-S方程。
∙1876年雷诺发现了流体流动两种流态:
层流和紊流。
∙1858年亥姆霍兹指出了理想流体中旋涡许多基本性质及旋涡运动理论,并于1887年提出了脱体绕流理论。
∙19世纪末,相似理论提出,实验和理论分析相结合。
∙1904年普朗特提出了边界层理论。
∙20世纪60年代以后,计算流体力学得到了迅速发展。
流体力学内涵不断地得到了充实及提高。
理想势流伯努利方程
(3-14)
或
(3-15)
物理意义:
在同一恒定不可压缩流体重力势流中,理想流体各点总比能相等即在整个势流场中,伯努利常数C均相等。
(应用条件:
“ ”所示)
符号说明
物理意义
几何意义
单位重流体位能(比位能)
位置水头
单位重流体压能(比压能)
压强水头
单位重流体动能(比动能)
流速水头
单位重流体总势能(比势能)
测压管水头
总比能
总水头
二、沿流线积分
1.只有重力作用不可压缩恒定流,有
2.恒定流中流线及迹线重合:
沿流线(或元流)能量方程:
(3-16)
注意:
积分常数C,在非粘性、不可压缩恒定流流动中,沿同一流线保持不变。
一般不同流线各不相同(有旋流)。
(应用条件:
“ ”所示,可以是有旋流)
流速势函数(势函数)
观看录像>>
• 存在条件:
不可压缩无旋流,即
或
必要条件
存在全微分dϕ
直角坐标
(3-19)
式中:
ϕ——无旋运动流速势函数,简称势函数。
• 势函数拉普拉斯方程形式
对于不可压缩平面流体流动中,将(3-19)式代入连续性微分方程(3-18),有:
或
(3-20)
适用条件:
不可压缩流体有势流动。
点击这里练习一下
极坐标
(3-21)
流函数
1.流函数
存在条件:
不可压缩流体平面流动。
直角坐标
连续性微分方程:
必要条件
存在全微分dy
(3-22)
式中:
y——不可压缩流体平面流动流函数。
适用范围:
无旋流、有旋流、实际流体、理想流体不可压缩流体平面流动。
流函数拉普拉斯方程形式
对平面势流,有
,则
或
(3-23)
适用条件:
不可压缩流体平面有势流动。
极坐标
(3-24)
2.流函数物理意义
(1)流函数等值线
就是流线。
得平面流线方程(3-1):
,得证。
(2)不可压缩流体平面流动中,任意两条流线流函数之差dy等于这两条流线间所通过单位宽度流量dq。
AB断面所通过流量:
图3-26
粘性流体运动微分方程
1.粘性流体特点
(1)实际流体面积力包括:
压应力和粘性引起切应力。
切应力由广义牛顿内摩擦定律确定:
(2)实际流动流体任一点动压强,由于粘性切应力存在,各向大小不等,即pxx≠pyy≠pzz。
任一点动压强由式(2-5)为:
(3-11)
第三节 流体动力学基本方程式
一、连续性微分方程
在流场内取一微元六面体(如图3-23),边长为dx,dy,dz,中心点O流速为(ux,uy,uz)
以x轴方向为例:
图3-23
左表面流速
右表面流速
所以 单位时间内x方向流出流进质量流量差:
x方向:
同理可得:
y方向:
z方向:
质量守恒定律:
单位时间内流出及流入六面体流体质量差之总和应等于六面体内因密度变化而减少质量,即:
(3-6)
(1)流体连续性微分方程一般形式
由(3-6)式可得
(3-7)
适用范围:
理想流体或实际流体;恒定流或非恒定流;可压缩流体或不可压缩流体。
(2)可压缩流体恒定流动连续性微分方程
当为恒定流时,有,则(3-7)式为
(3-8)
适用范围:
理想、实际、可压缩、不可压缩恒定流。
(3)不可压缩流体连续性微分方程
当为不可压缩流时,有
,则(3-7)式为
(3-9)
物理意义:
不可压缩流体单位时间内流入单位空间流体体积(质量),及流出流体体积(质量)之差等于零。
适用范围:
理想、实际、恒定流或非恒定流不可压缩流体流动。
二、理想流体运动微分方程
理想流体动水压强特性及静水压强特性相同:
从理想流体中任取一(x,y,z)为中心微元六面体为控制体,边长为dx,dy,dz,中心点压强为p(x,y,z),如图3-24。
图3-24
受力分析(x方向为例):
1.表面力
因为理想流体,所以t=0
左表面
右表面
2.质量力
单位质量力在各坐标轴上分量为X,Y,Z,所以x方向质量力为Xρdxdydz
由牛顿第二运动定律
,x方向有:
理想流体运动微分方程(欧拉运动微分方程)
(3-10)
适用范围:
恒定流或非恒定流,可压缩流或不可压缩流体。
若加速度
等于0,则上式就可转化为欧拉平衡微分方程(2-6)式
三、粘性流体运动微分方程
1.粘性流体特点
(1)实际流体面积力包括:
压应力和粘性引起切应力。
切应力由广义牛顿内摩擦定律确定:
(2)实际流动流体任一点动压强,由于粘性切应力存在,各向大小不等,即pxx≠pyy≠pzz。
任一点动压强由式(2-5)为:
(3-11)
2.实际流体运动微分方程式
图3-25
同样取一微元六面体作为控制体,如图3-25。
x向受力
左右向压力、上下向切力、前后面切力、质量力
x方向(牛顿第二运动定律
)
考虑条件:
1)不可压缩流体连续性微分方程(3-9):
2)切应力及主应力关系表达式(3-11)。
可得不可压缩粘性流体运动微分方程:
纳维-斯托克斯方程(Navier-Stokes,N-S)方程
(3-12)
拉普拉斯算符
,例:
想一想:
N-S方程及欧拉运动微分方程有何联系?
N-S方程是不可压缩粘性流体运动微分方程,而欧拉运动微分方程则是理想流体运动微分方程。
当流动流体运动粘度等于0,即为理想流体时,N-S方程即为欧拉运动微分方程。
第四节 欧拉运动微分方程积分
由于欧拉运动微分方程是一个一阶非线性偏微分方程组(迁移加速度三项中包含了未知数及其偏导数乘积),因而至今还无法在一般情况下积分,只能在一定条件下积分。
欧拉运动微分方程组(3-10)各式分别乘以dx,dy,dz(流场任意相邻两点间距ds坐标分量),然而相加得:
(3-13)
一、在势流条件下积分
考虑条件
1.恒定流:
;
2.均匀不可压缩流体,即ρ=const,;
3.质量力只有重力,即X=Y=0,Z=-g;
4.有势流动,满足式(3-5):
;
因此,(3-13)式中各项为:
(考虑欧拉加速度表达式(3-3))
(引入有势流动条件4)
由以上得:
积分得:
第一节 流态判别
一、两种流态运动特征
1883年英国物理学家雷诺(ReynoldsO.)通过试验观察到液体中存在层流和紊流两种流态。
1.层流 观看录像>>
层流(laminarflow),亦称片流:
是指流体质点不相互混杂,流体作有序成层流动。
特点:
(1)有序性。
水流呈层状流动,各层质点互不混掺,质点作有序直线运动。
(2)粘性占主要作用,遵循牛顿内摩擦定律。
(3)能量损失及流速一次方成正比。
(4)在流速较小且雷诺数Re较小时发生。
2.紊流 观看录像>>
紊流(turbulentflow),亦称湍流:
是指局部速度、压力等力学量在时间和空间中发生不规则脉动流体运动。
特点:
(1)无序性、随机性、有旋性、混掺性。
流体质点不再成层流动,而是呈现不规则紊动,流层间质点相互混掺,为无序随机运动。
(2)紊流受粘性和紊动共同作用。
(3)水头损失及流速1.75~2次方成正比。
(4)在流速较大且雷诺数较大时发生。
二、雷诺实验
如图6-1所示,实验曲线分为三部分:
(1)ab段:
当υ<υc时,流动为稳定层流。
(2)ef段:
当υ>υ''时,流动只能是紊流。
(3)be段:
当υc<υ<υ''时,流动可能是层流(bc段),也可能是紊流(bde段),取决于水流原来状态。
图6-1
图6-2
观看录像一>> 观看录像二>> 观看录像三>>
实验结果(图6-2)数学表达式
层流:
m1=1.0,hf=k1v,即沿程水头损失及流线一次方成正比。
紊流:
m2=1.75~2.0,hf=k2v1.75~2.0 ,即沿程水头损失hf及流速1.75~2.0次方成正比。
层流:
紊流:
流态判别
一、两种流态运动特征
1883年英国物理学家雷诺(ReynoldsO.)通过试验观察到液体中存在层流和紊流两种流态。
1.层流
层流(laminarflow),亦称片流:
是指流体质点不相互混杂,流体作有序成层流动。
特点:
(1)有序性。
水流呈层状流动,各层质点互不混掺,质点作有序直线运动。
(2)粘性占主要作用,遵循牛顿内摩擦定律
牛顿内摩擦定律
a.牛顿内摩擦定律:
液体运动时,相邻液层间所产生切应力及剪切变形速率成正比。
即
(N/m2,Pa) (1-6)
τ—粘性切应力,是单位面积上内摩擦力。
说明:
1)流体切应力及剪切变形速率,或角变形率成正比。
——区别于固体重要特性:
固体切应力及角变形大小成正比。
2)流体切应力及动力粘度μ成正比。
3)对于平衡流体du/dy=0,对于理想流体μ=0,所以均不产生切应力,即t=0。
b.牛顿平板实验及内摩擦定律
图1-1 流体绝对粘度
设板间y向流速呈直线分布,即:
则:
实验表明,对于大多数流体满足:
引入动力粘度μ,则得牛顿内摩擦定律
(1-7)
式中:
流速梯度
代表液体微团剪切变形速率。
线性变化时,即;非线性变化时,
即是u对y求导。
证明:
在两平板间取一方形质点,高度为dy,dt时间后,质点微团从abcd运动到a′b′c′d′。
由图1-2得:
图1-2
则:
说明:
流体切应力及剪切变形速率,或角变形率成正比。
(3)能量损失及流速一次方成正比。
(4)在流速较小且雷诺数Re较小时发生。
2.紊流
紊流(turbulentflow),亦称湍流:
是指局部速度、压力等力学量在时间和空间中发生不规则脉动流体运动。
特点:
(1)无序性、随机性、有旋性、混掺性。
流体质点不再成层流动,而是呈现不规则紊动,流层间质点相互混掺,为无序随机运动。
(2)紊流受粘性和紊动共同作用。
(3)水头损失及流速1.75~2次方成正比。
(4)在流速较大且雷诺数较大时发生。
三、层流、紊流判别标准——临界雷诺数
临界雷诺数
上临界雷诺数:
层流→紊流时临界雷诺数,它易受外界干扰,数值不稳定。
下临界雷诺数:
紊流→层流时临界雷诺数,是流态判别标准,它只取决于水流边界形状,即水流过水断面形状。
变直径管流中,细断面直径d1,粗断面直径d2=2d1,则粗细断面雷诺数关系是 。
圆管流
(5-1)
层流
紊流
明渠流
(5-2)
式中:
R——水力半径,R=A/P;
A——过水断面面积;
P——湿周,即断面中固体边界及流体相接触部分周长。
问题:
雷诺数及哪些因数有关?
其物理意义是什么?
当管道流量一定时,随管径加大,雷诺数是增大还是减小?
答案:
雷诺数及流体粘度、流速及水流边界形状有关。
Re=惯性力/粘滞力,
随d增大,Re减小。
.为什么用下临界雷诺数,而不用上临界雷诺数作为层流及紊流判别准则?
答:
上临界雷诺数不稳定,而下临界雷诺数较稳定,只及水流过水断面形状有关。
3.当管流直径由小变大时,其下临界雷诺数如何变化?
答:
不变,临界雷诺数只取决于水流边界形状,即水流过水断面形状。
三、紊流基本方程
对N-S方程(3-12)和连续性方程(3-9)进行时间平均即可得出紊流时均流动方程。
连续性方程
(6-20)
N-S方程(x方向)
(6-21)
式中:
——由于脉动产生附加法应力
统称为雷诺应力
——由于脉动产生附加切应力
它们是紊流传输项,也是造成紊流动量交换及质点混掺主要原因。
在紊流边界层外侧或紊流扩散中,雷诺应力远远超过粘性切应力。
边界层概念
一、边界层提出
1.边界层(boundarylayer):
图6-17
亦称附面层,雷诺数很大时,粘性小流体(如空气或水)沿固体壁面流动(或固体在流体中运动)时壁面附近受粘性影响显著薄流层,如图6-17。
判断:
边界层内流体流动及粘性底层流体流动都属于层流。
你回答:
对
错
2.流场求解可分为两个区进行
根据边界层概念,可将流场求解可分为两个区进行:
边界层内流动必须计入流体粘性影响可利用动量方程求得近似解。
边界层外流动视为理想流体流动,可按势流求解。
二、层流边界层和紊流边界层
1.边界层描述
普兰特把贴近于平板边界存在较大切应力,粘性影响不能忽略薄层称为边界层,图6-18。
边界中水流同样存在两种流态:
层流和紊流。
图6-18
2.边界层厚度
边界层厚度δ(boundarylayerthickness):
自固体边界表面沿其外法线到纵向流速ux达到主流速U099%处,这段距离称为边界层厚度。
边界层厚度顺流增大,即δ是x函数。
3.转捩点,临界雷诺数
转捩点:
在x=xcr处边界层由层流转变为紊流过渡点。
临界雷诺数:
(6-45)
特点:
临界雷诺数大小及来流脉动程度有关,脉动强,
小。
层流边界层及紊流边界层(图6-19)
层流边界层(laminarboundarylayer):
当边界层厚度d较小时,边界层内流速梯度很大,粘滞应力作用也很大,这时边界层内流动属于层流,这种边界层称为层流边界层。
紊流边界层(turbulenceboundarylayer):
当雷诺数达到一定数值时,边界层中层流经过一个过渡区后转变为紊流,就成为紊流边界层。
在紊流边界层内,最紧靠平板地方,dux/dy仍很大,粘滞力仍起主要作用,其流态仍为层流,所以紊流边界层内有一粘性底层。
图6-19
光滑平板边界层
临界雷诺数范围:
临界雷诺数并非常量,而是及来流扰动程度有关,如果来流受到扰动,脉动强,流态改变在较低雷诺数就会发生。
边界层厚度
层流边界层
紊流边界层
(6-46)
(6-47)
5.边界层特点
(1)边界层厚度为一有限值(当ux→0.99u时)
(2)边界层厚度沿程增加(δ=δ(x))
(3)边界层内:
;边界层外:
按理想流体或有势流动计算。
(4)边界层分层流边界层和紊流边界层。
边界层分离
1.边界层分离(separationofboundarylayer):
因压强沿流动方向增高,边界层内流体从壁面离开现象称边界层分离。
观看录像>>
平板绕流边界层分离,如图6-20。
压强梯度保持为零,即dp/dx=0
无论板有多长,都不会发生分离,这时边界层只会沿流向连续增厚。
压强沿程增大,即p2>p1或梯度dp/dx>0
边界层迅速地增厚,压强增大(流速减小)和阻力增大使边界层内动量减小,如两者共同作用在一足够长距离,致使边界层内流体流动停滞下来,分离便由此而生,自分离点B起,边界流线必脱离边界,其下游近壁处形成回流(或涡旋),在分离点:
(6-48)
(6-49)
图6-20
点击这里练习一下!
2.尾流
尾流:
分离流线及物体边界所围下游区域,如图6-21。
减小尾流主要途径:
使绕流体型尽可能流线型化。
观看录像>>
图6-21
1.流体流动两种形态(层流和紊流)特点。
(质点是否掺混,运动是否有序,水头损失及流速间关系)
2.层流、紊流判别标准——下临界雷诺数Rec Rec只取决于边界形状(过水断面形状)。
对圆管流Rec<2300时为层流。
3.均匀流基本方程:
τ0=ρgRJ τ=ρgR'J
4.不可压缩恒定均匀圆管层流
圆管层流流速呈旋转抛物面分布:
。
圆管层流最大流速:
圆管层流断面平均流速:
断面平均流速是最大流速为2倍。
圆管层流水头损失:
,即水头损失及流速一次方成正比,沿程阻力系数λ=64/Re。
5.紊流特点:
无序性、耗能性、扩散性。
时均化处理紊流。
瞬时流速=时均流速+脉动流速
6.紊流切应力:
7.紊流流速分布
a.近壁处:
,线性分布
b.紊流核心区:
,对数分布
粘性底层厚度:
,随Re增大而减小
8.能量损失,
恒定紊流能量方程
一、水流阻力及水头损失
产生流动阻力和能量损失根源:
流体粘性和紊动。
1.水头损失两种形式
(1)沿程阻力和沿程水头损失
沿程阻力(frictionaldrag):
当限制流动固体边界使流体作均匀流动时,流动阻力只有沿程不变切应力,该阻力称为沿程阻力。
沿程水头损失(frictionalheadloss):
由沿程阻力作功而引起水头损失称为沿程水头损失。
(2)局部阻力和局部水头损失 观看录像>>
局部阻力(localresistance):
液流因固体边界急剧改变而引起速度分布变化,从而产生阻力称为局部阻力。
局部水头损失(localheadloss):
由局部阻力作功而引起水头损失称为局部水头损失。
(3)特点
沿程水头损失hf:
主要由于“摩擦阻力”所引起,随流程增加而增加。
在较长直管道和明渠中是以hf为主流动。
局部阻力水头损失hj :
主要是因为固体边界形状突然改变,从而引起水流内部结构遭受破坏,产生漩涡,以及在局部阻力之后,水流还要重新调整结构以适应新均匀流条件所造成。
例“弯头”,“闸门”,“突然扩大”等。
(4)水头损失叠加原理
水头损失叠加原理:
流段两截面间水头损失为两截面间所有沿程损失和所有局部损失总和。
即:
(6-26)
式中:
n——等截面段数;
m——局部阻力个数。
不同固体边界下水头损失如图6-11:
图6-11
2.沿程水头损失公式
(1)魏斯巴赫(Weisbach)公式
实验表明:
(6-27)
式中:
λ——沿程阻力系数。
R——水力半径,R=A/P。
适用范围:
适用于任意形状等截面流道恒定均匀流。
(2)圆管流达西-魏斯巴赫公式(简称为D-W公式)
圆管R=d/4,则
(6-28)
适用范围:
适用于圆管紊流或层流,为恒定均匀管流通用