紊流理论紊流模型 ppt课件
合集下载
河流动力学第二章 水流的紊动(1)
![河流动力学第二章 水流的紊动(1)](https://img.taocdn.com/s3/m/e6d246bb87c24028915fc3dc.png)
u*
(y 0 )
说明:在紊流核心区(y>0),紊流流速呈对数规律分布。
在水面处 y H , u umax 代入得
umax u 1 ln H
u*
y
最大流速umax 往往未知,且该式无法反应边界粗糙情
况
紊流存在不同的壁面类型
一、光壁的流速分布
y ,u u
u 1 ln yu* A
向着波谷方向的压力
这两个垂直于流层的压力将促使这 个流层的波幅更加增大
波幅增大到一定程度后,动水压力形成的力 偶和切力产生的力偶,将促使涡体形成。
在涡体上部,旋转方向和上部流速方向一致,流速加 大而压强变小,下部则流速减小而压强增大,这样就 产生了一个压差即升力,迫使涡体从一个流层进入另 一个流层而混掺。
u
y
5.75 lg(30.2 )
u*
ks
式中 χ 为修正系数,它是 ks / 的函数,见下图。
χ与 k s / 的关系
沿水深积分可得紊流的断面平均流速 光壁条件下的断面平均流速公式
U u*
3.25 5.75 lg Ru*
糙壁条件下的断面平均流速
U u*
6.25
5.75 lg
R Ks
引进校正值 以上两式合并为
(1)忽略粘滞切应力,近壁处切应力为一常量,且 其值等于边壁处的切应力 0,即= 0 ;
(2)混合长度l随着离边壁的距离y呈线性变化,即l=y。
0
l 2 ( dux
dy
)2
2y2 ( dux
dy
)2
u*
0 y dux
dy
dux 1 dy
u* y
摩阻流速如何确定?
ux 1 ln y C
(y 0 )
说明:在紊流核心区(y>0),紊流流速呈对数规律分布。
在水面处 y H , u umax 代入得
umax u 1 ln H
u*
y
最大流速umax 往往未知,且该式无法反应边界粗糙情
况
紊流存在不同的壁面类型
一、光壁的流速分布
y ,u u
u 1 ln yu* A
向着波谷方向的压力
这两个垂直于流层的压力将促使这 个流层的波幅更加增大
波幅增大到一定程度后,动水压力形成的力 偶和切力产生的力偶,将促使涡体形成。
在涡体上部,旋转方向和上部流速方向一致,流速加 大而压强变小,下部则流速减小而压强增大,这样就 产生了一个压差即升力,迫使涡体从一个流层进入另 一个流层而混掺。
u
y
5.75 lg(30.2 )
u*
ks
式中 χ 为修正系数,它是 ks / 的函数,见下图。
χ与 k s / 的关系
沿水深积分可得紊流的断面平均流速 光壁条件下的断面平均流速公式
U u*
3.25 5.75 lg Ru*
糙壁条件下的断面平均流速
U u*
6.25
5.75 lg
R Ks
引进校正值 以上两式合并为
(1)忽略粘滞切应力,近壁处切应力为一常量,且 其值等于边壁处的切应力 0,即= 0 ;
(2)混合长度l随着离边壁的距离y呈线性变化,即l=y。
0
l 2 ( dux
dy
)2
2y2 ( dux
dy
)2
u*
0 y dux
dy
dux 1 dy
u* y
摩阻流速如何确定?
ux 1 ln y C
紊流理论(紊流模型)
![紊流理论(紊流模型)](https://img.taocdn.com/s3/m/c64338a4770bf78a652954c3.png)
切细节,但在实际工程中具有重要意义的并不是紊 流的一切细节,而是紊流对于时间的平均效应。 – 雷诺(Osborne Reynolds)建议用统计方法将N -S方程取时间平均。但取平均的过程产生了新问 题:方程增加了新未知项,时均方程组不再封闭, 因此各种类型的紊流模型应运而生。 – 紊流模型可定义为一组方程,这组方程确定时均流 方程中的紊动输运项,从而封闭时均流动方程组( 零方程、单方程、双方程)。
雷诺方程数值模拟(RANS)
• 单方程模型
– 将上式i方向乘以u’j,j方向乘以u’i,然后相加, 进行时间平均得
雷诺方程数值模拟(RANS)
• 单方程模型
雷诺方程数值模拟(RANS)
• 单方程模型
– 令i=j,将k带入,可得
雷诺方程数值模拟(RANS)
• 单方程模型
– 可见二者相同。
雷诺方程数值模拟(RANS)
紊流理论(紊流模型)
概述
• 紊流运动的瞬时连续方程和运动方程 • 连续方程 • 运动方程
概述
• 直接求解连续方程和运动方程,即 DNS(Direct Numerical Simulation)。
• 原则上讲,DNS并无理论上的困难。
– 一方面,描述紊流运动的精确的微分方程已经得出 ,即N-S方程;从数学观点看来,紊流就是N-S 方程的通解。
概述
• 紊流模型分类
– 直接数值模拟(DNS)。直接求解N-S方程,必 须采用很小的时间步长与空间步长。
– 大涡模拟(LES)。 DNS模拟大尺度运动,亚网格 上模拟小尺度运动。
– Reynolds时均方程法(RANS)。将N-S方程对 时间平均,并通过一些假定建立模型,是目前工程 中所采用的基本方法。根据对雷诺应力的处理方式 不同,分为基于涡粘性假设的模型和应力输运模型 两类,根据引入方程的数量,前者又分为零方程、 单方程和双方程模型。
雷诺方程数值模拟(RANS)
• 单方程模型
– 将上式i方向乘以u’j,j方向乘以u’i,然后相加, 进行时间平均得
雷诺方程数值模拟(RANS)
• 单方程模型
雷诺方程数值模拟(RANS)
• 单方程模型
– 令i=j,将k带入,可得
雷诺方程数值模拟(RANS)
• 单方程模型
– 可见二者相同。
雷诺方程数值模拟(RANS)
紊流理论(紊流模型)
概述
• 紊流运动的瞬时连续方程和运动方程 • 连续方程 • 运动方程
概述
• 直接求解连续方程和运动方程,即 DNS(Direct Numerical Simulation)。
• 原则上讲,DNS并无理论上的困难。
– 一方面,描述紊流运动的精确的微分方程已经得出 ,即N-S方程;从数学观点看来,紊流就是N-S 方程的通解。
概述
• 紊流模型分类
– 直接数值模拟(DNS)。直接求解N-S方程,必 须采用很小的时间步长与空间步长。
– 大涡模拟(LES)。 DNS模拟大尺度运动,亚网格 上模拟小尺度运动。
– Reynolds时均方程法(RANS)。将N-S方程对 时间平均,并通过一些假定建立模型,是目前工程 中所采用的基本方法。根据对雷诺应力的处理方式 不同,分为基于涡粘性假设的模型和应力输运模型 两类,根据引入方程的数量,前者又分为零方程、 单方程和双方程模型。
圆管中的层流、紊流运动ppt课件
![圆管中的层流、紊流运动ppt课件](https://img.taocdn.com/s3/m/f30093e66137ee06eff9182a.png)
圆管中的层流、 紊流运动
第五章 流动阻力与水头损失
• • • • §5–1 概述 §5–2 粘性流体的流动型态 §5–3 均匀流基本方程 §5–4 圆管中的层流运动
•
•
§5–5 圆管中的紊流运动
§5–6 局部水头损失
§5–4 圆管中的层流运动
一、流速分布
r y r r0 u
r
r0
umax
v
1、圆管层流的流速分布
l
~
粘性底层
~ ~ ~ ~ ~ ~
~ ~
紊流核心
~ ~ ~ ~ ~
l
粘性底层
紊流核心
2、粘性底层的流速分布
dux o dy
uo
u du x 0 0 dy 0
u / 0 0 0 / 0
令Hale Waihona Puke u 0 / *与流速量纲相同,称剪切流速
du 1 dy x u* k y
u 1 x ln y C ( y ) 0 u * k
说明:在紊流核心区(y>0),紊流流速呈对数规律分布。
流量qV=2.5×10-4m3/s,求hf
[例2] =850kg/m3的油在管径100mm, =0.18×10-4m2 /s的管中以v =0.0635m /s的速 度运动,求:(1)管中心处的最大流速; (2)在离管中心r=20mm 处的流速; (3)沿程阻力因数 ; (4)管壁切应力0及每km管长的水头损失。
[例3] 应用细管式粘度计测定油的粘度,已知细管直径d=6mm,测量段长l=2m 如图。实 测油的流量qV =77×10-6m3 /s ,水银压差计的读值hp=0.3m,油的密度=900kg/m3 。试
第五章 流动阻力与水头损失
• • • • §5–1 概述 §5–2 粘性流体的流动型态 §5–3 均匀流基本方程 §5–4 圆管中的层流运动
•
•
§5–5 圆管中的紊流运动
§5–6 局部水头损失
§5–4 圆管中的层流运动
一、流速分布
r y r r0 u
r
r0
umax
v
1、圆管层流的流速分布
l
~
粘性底层
~ ~ ~ ~ ~ ~
~ ~
紊流核心
~ ~ ~ ~ ~
l
粘性底层
紊流核心
2、粘性底层的流速分布
dux o dy
uo
u du x 0 0 dy 0
u / 0 0 0 / 0
令Hale Waihona Puke u 0 / *与流速量纲相同,称剪切流速
du 1 dy x u* k y
u 1 x ln y C ( y ) 0 u * k
说明:在紊流核心区(y>0),紊流流速呈对数规律分布。
流量qV=2.5×10-4m3/s,求hf
[例2] =850kg/m3的油在管径100mm, =0.18×10-4m2 /s的管中以v =0.0635m /s的速 度运动,求:(1)管中心处的最大流速; (2)在离管中心r=20mm 处的流速; (3)沿程阻力因数 ; (4)管壁切应力0及每km管长的水头损失。
[例3] 应用细管式粘度计测定油的粘度,已知细管直径d=6mm,测量段长l=2m 如图。实 测油的流量qV =77×10-6m3 /s ,水银压差计的读值hp=0.3m,油的密度=900kg/m3 。试
第十章紊流
![第十章紊流](https://img.taocdn.com/s3/m/a241674b77232f60ddcca146.png)
2017/9/30
西安交通大学流体力学课程组
8
时均法运算性质
设 f,
f ,g g 为紊流时均参数,
脉动参数, f , g 为瞬时参数
(1)
f f
1 f T
T 2 T t 2 t
1 fd f T
T 2 T t 2 t
d f
时均值的时均值仍为原时均值
(2) f g f g
ui ui ui ui 2ui 1 p uj uj uj uiuj t x j x j x j x j xi x j x j
2017/9/30 西安交通大学流体力学课程组 21
脉动速度方程
瞬时量方程
ui ui 2 ui 1 p uj t x j xi x j x j
ui u 1 p uj i t x j xi x j ui u u i j x j
fg f g f g
f f x x
u j x j 0
ui ui ui 2 ui 1 p uj uj t x j x j xi x j x j
ui uj ( uiuj ) ui (u iuj ) x j x j x j x j
f f f
如
2017/9/30
u u u
p p p
7
西安交通大学流体力学课程组
紊流的时间平均法
严格来说,时均平均法只适用定常紊流,实际上已推广 用于非定常紊流 定常紊流是指时均特性不随时间变化的紊流流动 将紊流流动分为两部分,即:时均流动和脉动运动
时均流动代表主流,关注的重点是时均流动特性 脉动流动反映紊流的实质,对时均流动一切特性都产生 影响
紊流理论(基本理论)-60页精选文档
![紊流理论(基本理论)-60页精选文档](https://img.taocdn.com/s3/m/00b0d75d52ea551811a68707.png)
– 来流的紊动度和壁面的粗糙程度都影响转捩的发生,粗糙表 面的边界层更容易发展为紊流,因此会导致阻力危机的提前 发生。
– 因此,绕流阻力主要 取决于物体背流面尾 流中的负压,与背流 面的形状关系密切。 这一问题直到边界层 理论提出后才得到解930年,普朗特建立了层流稳定性理论。 – 层流稳定性的基本点是:层流流动总会受到一些扰动,可能
层流稳定性理论
• 层流到紊流的转捩-壁面边界层流动
– 以经典的圆柱绕流为例。
– 粘性流动壁面无滑移,产生边界层,在背流面发生分离,形 成一个由漩涡组成的尾流区。
层流稳定性理论
• 层流到紊流的转捩-壁面边界层流动
– 临界Re附近,边界层流动由层流向紊流转捩,分离点下移, 尾流区缩小,形状阻力降低,产生了阻力危机。
第一项
第二项
第三项
=fx1144P x2414px3' 1(42xU 42 442yU 24)24(42xu42' 442yu32' )
第四项
第五项
U u' UU Uu' u' U u' u' V U V u' v' Uv' u'
14t 2 43t 1 4x4 4 44x2 4 4x4 4 43x 1 4y4 4 44y2 4 4y4 4 43y
uUu',vVv',wWw',pPp'
层流稳定性理论
问题:对于这样一个主流流动,主流满足N-S方程, 叠加后的流动也满足,那么扰动将随时间放大还是衰 减。
uti uj xuij fi 1 xpi x2 j uxi j
u t
u u x
v u y
fx
– 因此,绕流阻力主要 取决于物体背流面尾 流中的负压,与背流 面的形状关系密切。 这一问题直到边界层 理论提出后才得到解930年,普朗特建立了层流稳定性理论。 – 层流稳定性的基本点是:层流流动总会受到一些扰动,可能
层流稳定性理论
• 层流到紊流的转捩-壁面边界层流动
– 以经典的圆柱绕流为例。
– 粘性流动壁面无滑移,产生边界层,在背流面发生分离,形 成一个由漩涡组成的尾流区。
层流稳定性理论
• 层流到紊流的转捩-壁面边界层流动
– 临界Re附近,边界层流动由层流向紊流转捩,分离点下移, 尾流区缩小,形状阻力降低,产生了阻力危机。
第一项
第二项
第三项
=fx1144P x2414px3' 1(42xU 42 442yU 24)24(42xu42' 442yu32' )
第四项
第五项
U u' UU Uu' u' U u' u' V U V u' v' Uv' u'
14t 2 43t 1 4x4 4 44x2 4 4x4 4 43x 1 4y4 4 44y2 4 4y4 4 43y
uUu',vVv',wWw',pPp'
层流稳定性理论
问题:对于这样一个主流流动,主流满足N-S方程, 叠加后的流动也满足,那么扰动将随时间放大还是衰 减。
uti uj xuij fi 1 xpi x2 j uxi j
u t
u u x
v u y
fx
紊流理论 ppt课件
![紊流理论 ppt课件](https://img.taocdn.com/s3/m/ed84a876b90d6c85ec3ac6d8.png)
hy
y h
普适流速分布律(卡门)
• 进一步改写为
umauxu1ln1
hy
y h
• 上式即为卡门流速分布亏损律。
• 管道中心,y=0处,梯度为0,卡门相似 理论不满足;y=h处,速度趋于无穷, 不合理。因为公式中未考虑粘性切应力。
• 公式应扣除管道中心和边壁的区域。
普适流速分布律(普朗特)
• 从普朗特混掺长度理论出发
u5.75logy5.5
明渠紊流流速分布与分区结构
• 明渠紊流的流速分布可用对数律来表示, 但是:
– 公式的两个常数确存在差异。
– 自由水面附近存在速度下降的现象。
– 70年代以来激光测速技术的发展,发现存在 分区结构:贴近边壁存在粘性底层,向上依 次存在过渡区,对数区,在明渠紊流上部还 存在尾流律区。
– 物理量的分析:流速与压强沿x方向无变化,即对x 的偏导数为0。垂向流速V很小,相对之下可以忽略。
明渠紊流流速分布与分区结构
• 二维恒定明渠紊流雷诺方程
UU xV U ygsinx Pxu2 +yuv+2U( 2 yU 2)
• 即 U V xV V ygcosy Pxuv+yv2 +2V
u 1 yu*
ln
u*
C1
普适流速分布律(普朗特)
• 无量纲数 u u , y yu*
u*
•
则无量纲公式变为: u
1
ln
y
C1
• 式中系数为卡门常数,常数是与壁面有
关的量,需要通过试验确定。
• 尼古拉兹对光滑圆管进行试验得到卡门 常数为0.4,C1为5.5,则光滑圆管中紊 流时均流速分布为:
2
l2
流体力学 5.2紊流
![流体力学 5.2紊流](https://img.taocdn.com/s3/m/08e777eedd88d0d233d46a83.png)
IV:紊动粘性应力对紊动位移做功,引起能量传 递(扩散)。此项小,是次要的。
V:紊动粘性应力对紊动变形做功,紊动能转化为 热能,亦称粘性耗损。(此项总为负)。是最后 消耗紊动能量的途径,是不能忽略的重要项。
第三章 紊流的基本方程与紊流模型
§3.4 紊流的涡量方程 前面讲紊流必是
三维的,即三个方向 瞬时流速都有梯度, 因此紊流必然是有 涡流动,其涡量
第三章 紊流的基本方程与紊流模型
注意:
(1)混合长假定与实际不符,与流体连续 性矛盾。
(2)紊流由大小不同的涡旋组成,L代表 哪一种涡不明确。
( 附3近)。假设u’1 u’2也不全符合实际。如壁面
(4)在圆管中心及明渠水面无速度梯度 (则 u’1=0) ,而实测之u’1并不很小 。
第三章 紊流的基本方程与紊流模型
但此项总为负(正应力与切应力作功综合 结果),故亦称紊动能产生项。
第三章 紊流的基本方程与紊流模型
※VI:时均粘性应力对时均变形作功 能量损失项(转变为热能)
ui X j
u j X i
ui X j
Sij Sij
总为负
第三章 紊流的基本方程与紊流模型
讨论:雷诺应力和粘性力对位移和变形 都要作功,那么两种力作功何者为大?
(
ui
)做
功
,引
起
的
IV:粘性应力对变形做功,使机械能转化 为热能损失
第三章 紊流的基本方程与紊流模型
紊流时均化的能量方程 将式(3.20)时均化 ui (N S)
式(3.21)包括时均能、脉动能(脉动值 与脉动值相关)、关联能 (时均值与脉动 值相关)
第三章 紊流的基本方程与紊流模型
V:紊动粘性应力对紊动变形做功,紊动能转化为 热能,亦称粘性耗损。(此项总为负)。是最后 消耗紊动能量的途径,是不能忽略的重要项。
第三章 紊流的基本方程与紊流模型
§3.4 紊流的涡量方程 前面讲紊流必是
三维的,即三个方向 瞬时流速都有梯度, 因此紊流必然是有 涡流动,其涡量
第三章 紊流的基本方程与紊流模型
注意:
(1)混合长假定与实际不符,与流体连续 性矛盾。
(2)紊流由大小不同的涡旋组成,L代表 哪一种涡不明确。
( 附3近)。假设u’1 u’2也不全符合实际。如壁面
(4)在圆管中心及明渠水面无速度梯度 (则 u’1=0) ,而实测之u’1并不很小 。
第三章 紊流的基本方程与紊流模型
但此项总为负(正应力与切应力作功综合 结果),故亦称紊动能产生项。
第三章 紊流的基本方程与紊流模型
※VI:时均粘性应力对时均变形作功 能量损失项(转变为热能)
ui X j
u j X i
ui X j
Sij Sij
总为负
第三章 紊流的基本方程与紊流模型
讨论:雷诺应力和粘性力对位移和变形 都要作功,那么两种力作功何者为大?
(
ui
)做
功
,引
起
的
IV:粘性应力对变形做功,使机械能转化 为热能损失
第三章 紊流的基本方程与紊流模型
紊流时均化的能量方程 将式(3.20)时均化 ui (N S)
式(3.21)包括时均能、脉动能(脉动值 与脉动值相关)、关联能 (时均值与脉动 值相关)
第三章 紊流的基本方程与紊流模型
水力学课件 第4章层流和紊流、液流阻力和水头损失
![水力学课件 第4章层流和紊流、液流阻力和水头损失](https://img.taocdn.com/s3/m/17c3d4c56e1aff00bed5b9f3f90f76c661374cbe.png)
13
实验结果——关于流态
1. vc΄> vc 2. v< vc 为层流
v > vc΄ 为紊流 3. vc <v< vc΄ 为过渡区
14
实验结果——关于hf与v的关系 lg hf lg k m lg v
取反对数得:hf kvm
AB段 (层流):
m 1(1 45 ) ; hf ~ v1
DE段 (紊流):
(2)紊流过渡粗糙区 ( , Re)
d
结论:
① 沿程水头损失系数既和Re有关也 和相对粗糙度有关
4.9.1人工粗糙管的试验研究— 尼古拉兹试验
3紊流区 lg Re 3.6
(3)紊流粗糙区
()
d
结论:
① λ和Re无关,只和相对粗糙度有关; ② hf是v的2次方
讨论
紊流分区与壁面分类关系:
Re vd
——雷诺数
Rec
vc d
为下临界雷诺数;
Rec
vcd
为上临界雷诺数。
G
对于圆管,临界雷诺数相对稳定:
Rec 2300
17
雷诺数的物理意义:惯性力与粘性力的比
F
V
dv dt
L3 U T
L2U 2
UL
T A du L2 U LU
dy
L
对于非圆管:
Re vR
过 水 断 面 上 , 水 流 与 固 体 边 界 接 触 的 长 度 , 称 为 湿 周 , 用 表 示 。
l
( z1
p1 g
)
(z2
p2 g
)
'
l
gA' gR'
( z1
p1 g
)
实验结果——关于流态
1. vc΄> vc 2. v< vc 为层流
v > vc΄ 为紊流 3. vc <v< vc΄ 为过渡区
14
实验结果——关于hf与v的关系 lg hf lg k m lg v
取反对数得:hf kvm
AB段 (层流):
m 1(1 45 ) ; hf ~ v1
DE段 (紊流):
(2)紊流过渡粗糙区 ( , Re)
d
结论:
① 沿程水头损失系数既和Re有关也 和相对粗糙度有关
4.9.1人工粗糙管的试验研究— 尼古拉兹试验
3紊流区 lg Re 3.6
(3)紊流粗糙区
()
d
结论:
① λ和Re无关,只和相对粗糙度有关; ② hf是v的2次方
讨论
紊流分区与壁面分类关系:
Re vd
——雷诺数
Rec
vc d
为下临界雷诺数;
Rec
vcd
为上临界雷诺数。
G
对于圆管,临界雷诺数相对稳定:
Rec 2300
17
雷诺数的物理意义:惯性力与粘性力的比
F
V
dv dt
L3 U T
L2U 2
UL
T A du L2 U LU
dy
L
对于非圆管:
Re vR
过 水 断 面 上 , 水 流 与 固 体 边 界 接 触 的 长 度 , 称 为 湿 周 , 用 表 示 。
l
( z1
p1 g
)
(z2
p2 g
)
'
l
gA' gR'
( z1
p1 g
)
高二物理竞赛课件:流体力学的紊流流动
![高二物理竞赛课件:流体力学的紊流流动](https://img.taocdn.com/s3/m/bad4778559f5f61fb7360b4c2e3f5727a5e924d0.png)
0.0002
0.0001 0.000,05
2 3 4 5 6 810 7 2 0.000001
d
0.000,01
3
456
8
10 8
0.000005
d
[例] 沿程损失:已知管道和流量求沿程损失
已知: d=200mm , l=3000m 的旧无缝钢管, ρ=900 kg/m3, Q=90T/h,
lg
Re
d
2
1.42
lg
1.273
qV
2
V.湍流平方阻力区 λ=f(ε/d )
1
1/2
2lg
d
2
1.74
莫迪图
(用于计算新的工业管道 ) (600 Re 108)
0.1 0.09 0.08 0.07
0.06
0.05
0.04
0.03 0.025
0.02
0.015
0.01 0.009 0.008
夏天
hf2
2
l d
V2 2g
0.0385
3000 0.2
0.884 2
2 9.81
23.0m
(油柱)
[例] 沿程损失:已知管道和压降求流量
已知: d=10cm , l=400m 的旧无缝钢管比重为0.9,
=10 -5 m2/s 的油
p 800 KPa
求: 管内流量qv
解:
p 800 103
或
32.8d
Re 1 2
d ——管径 ——沿程损失系数
(4) 水力光滑与水力粗糙
管壁粗糙凸出部分的平均高度叫做管壁的绝对粗糙度(ε) ε /d 称为相对粗糙度
(a)
水力光滑 δ>ε 光滑管
0.0001 0.000,05
2 3 4 5 6 810 7 2 0.000001
d
0.000,01
3
456
8
10 8
0.000005
d
[例] 沿程损失:已知管道和流量求沿程损失
已知: d=200mm , l=3000m 的旧无缝钢管, ρ=900 kg/m3, Q=90T/h,
lg
Re
d
2
1.42
lg
1.273
qV
2
V.湍流平方阻力区 λ=f(ε/d )
1
1/2
2lg
d
2
1.74
莫迪图
(用于计算新的工业管道 ) (600 Re 108)
0.1 0.09 0.08 0.07
0.06
0.05
0.04
0.03 0.025
0.02
0.015
0.01 0.009 0.008
夏天
hf2
2
l d
V2 2g
0.0385
3000 0.2
0.884 2
2 9.81
23.0m
(油柱)
[例] 沿程损失:已知管道和压降求流量
已知: d=10cm , l=400m 的旧无缝钢管比重为0.9,
=10 -5 m2/s 的油
p 800 KPa
求: 管内流量qv
解:
p 800 103
或
32.8d
Re 1 2
d ——管径 ——沿程损失系数
(4) 水力光滑与水力粗糙
管壁粗糙凸出部分的平均高度叫做管壁的绝对粗糙度(ε) ε /d 称为相对粗糙度
(a)
水力光滑 δ>ε 光滑管
流体在管道中的流动PPT幻灯片课件
![流体在管道中的流动PPT幻灯片课件](https://img.taocdn.com/s3/m/d9a1e012f242336c1eb95eee.png)
•
de=4R=4×0.173=0.693m
•
V=Q/A=0.2/0.48=0.417m/s
10
材料工程基础
• 矩形 a·3a=3a2=0.48m2 a=0.4m b=1.2m
R ab 0.48 0.15 m 2(a b) 2 1.6
• de=4×0.15=0.6m V=Q/A=0.417m/s
则
hf
p1 p2
g
测压管中的水柱高差△P即为有效截面1-1和2-2 间的压头损失。
14
材料工程基础
图4-3 水平等 直管道中水头损失
15
材料工程基础
伯努利(能量)方程实验
16
材料工程基础
【例4-1】 管道直径 d 100mm,输送水的流量 qV 0.01
m3/s,水的运动黏度 1106 m2/s,求水在管中的流动状 态?若输送 1.14104 m2/s的石油,保持前一种情况下的流 速不变,流动又是什么状态?
6
材料工程基础
4.1.2 雷诺数
流体的流动状态与流速、管径和流体的黏性等物理性质有关。
uc d
引入比例系数 Rec
uc Rec d Rec d 或
Rec
ucd
Rec 称为临界雷诺数,是一个无量纲数。
7
材料工程基础
流体在任意形状截面的管道中流动时,雷诺数的形式是
Re ude
材料工程基础
第四节 流体在管道中的流动 一维定常流动
4.1 流体的两种流动状态 4.2 圆管中流体的层流流动 4.3 圆管中流体的紊流流动 4.4 流动阻力损失 4.5 管路计算
1
流体力学5.1紊流ppt课件
![流体力学5.1紊流ppt课件](https://img.taocdn.com/s3/m/3054a933cc7931b765ce154a.png)
北京工业大学市政学科部——马长明
高等流体(水)力学讲稿
4
第四讲
紊(湍)流运动基础
(4)充分发展的紊流研究 紊流的发生:剪切层的存在—产生涡。(剪切层的类型图示)
对充分发展的紊流研究分类: 1)自由剪切紊流:剪切层由流速间断面引起,紊动发展 不受边壁的限制; 2)边壁剪切紊流:剪切层由边壁附着引起,紊动发展受 边壁限制。 3)均匀各向同性紊流:作为理论研究的假想模型。流动 中无速度梯度,也无剪切应力。
稳定流动性图示
(3)紊流是否有规律—紊流结构的实验研究 1)壁面剪切紊流的拟序结构(Quasi-order Structure) 猝发(Burst)过程:低速带,上升马蹄涡,喷射,清扫。
拟序结构图片
拟序图示1 轴对称射流
拟序图示2 尾流
2)自由紊流的相干结构(Coherent Structure)
二维差混合层
北京工业大学市政学科部——马长明 高等流体(水)力学讲稿
6
第四讲
紊(湍)流运动基础
二、基于统计理论的紊流运动方程
1、基本统计量 对不可压缩流动问题,基本未知量ui,p可认为是具有一定统计规律 的随机变量,即:可表示为
u u u i i i
p p p
其中上划线表示平均值,上ui’ 、p’表示脉动值,称为“涨落” (Fluctuation)。 (1)系综平均与时间平均 1)系综平均(ensemble average)
北京工业大学市政学科部——马长明 高等流体(水)力学讲稿
3
第四讲
紊(湍)流运动基础
3、紊流问题研究分类
理论上要求回答:紊流的发生、发展物理机理和结构。 应用上要求回答:紊流的平均流场、阻力、能耗与扩散的定量确定。 (1)N-S方程是否能用于描述紊流运动? (2)紊流的发生——流动稳定性研究(仅有几个线性流动解) 在给定边界条件的小扰动值下,求解线性化后的N-S方程,对不同的 雷诺数,由扰动是否衰减,来确定临界雷诺数。
流体力学层流紊流ppt课件
![流体力学层流紊流ppt课件](https://img.taocdn.com/s3/m/da6b5025a55177232f60ddccda38376baf1fe0dc.png)
和
R0
r0 2
圆管均匀流断面上的切应力
' 0
在半径
r
是线性分布
的,在管轴中心 r=0 处为零。在边壁 r=r0 处最大。 18
§5-2层流紊流及能量损失:均匀流沿程损失
二、沿程损失的通用公式
根据均匀流基本方程,总流的 hf 取决于边壁上的平均摩擦切应力 τ0。若能确定τ0 的大小,则易得到 hf 的变化规律。
R A
A
dH 4R 4
湿周
11
§5-1层流紊流及能量损失:层流及紊流
5、非圆管雷诺数ReH
Re H
VdH
12
§5-1层流紊流及能量损失:层流及紊流
三、紊流的成因
发生条件:
1、有初始扰动
层流:粘性力起主导作用
2、粘性对质点运动的束缚降低 紊流:惯性力起主导作用
13
§5-1层流紊流及能量损失:层流及紊流
3、圆管流态判断
Re Vd
Re Rec 2000 层流 Re Rec 2000 过渡流或紊流
4、非圆管流态判断
(1)特征量——描述运动的某一属性、相对独立的量。 例如管流的平均流速和管径、缝隙流里的平均流速和 缝隙高度。
(2)湿周——断面上固体边界与流体接触的周长,
(3)水力半径 (4)当量直径
根据实验,圆管均匀流边壁上的摩擦切应力τ0 与五个因素有关: 断面平均流速 V 水力直径 d
流体密度ρ 流体的动力粘度μ
壁面的粗糙高度 ks
19
§5-2层流紊流及能量损失:均匀流沿程损失
达西公式
由量纲和谐原理得
0
8
V
2
f Re, ks
d
hf
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
– 他们的共同贡献是指出了封闭 Reynolds 方程或 Reynolds 应力的封闭表达式应从湍流脉动场的性 质去寻找。称为一阶封闭格式。
– 周培源(1945)和Rotta(1951),绕过 Boussinesq 涡粘性假定,提出了一个描述紊流切 应力张量演化的微分方程,即雷诺应力张量,得到 了应力输运模型,也称为二阶封闭或者二阶矩封闭 模型。70年代得到应用。
– 塞弗曼(P G Saffmann)在1977年曾经预言,大 约在本世纪末,高度发展的计算机将有足够的能力 用数值计算方法求解紊流的精确方程,但在较近的 未来,精确地求解紊流问题无疑是不可能的。
– 到现在,DNS算法已有很大发展。
概述
• DNS之外的其他方法
– 求解时均N-S方程(即雷诺方程)。 – 作为紊流的通解, N-S方程描述了流体运动的一
雷诺方程数值模拟(RANS)
• 雷诺时均方程的封闭问题
– 将瞬时值写成时均值+脉动值uu u,代入连
续方程和运动方程,并对方程两边取平均:
ux uy uz 0 x y z
u ti uj x uij fi 1 x pi xj x uij uiuj
– 引入数学方程或者代数公式,确定所产生的新未知
– 1945年,Prandtl 提出一个涡粘性依赖于紊动能k 的模型,建议使用一个偏微分方程的模型对精确的 k方程进行近似,得到了单方程模型,即k方程模型。
概述
• 雷诺时均方程模型发展史
– 1942年,Kolmogorov提出了第一个完整的紊流模 型,除了k方程,还引入了另外一个参数ε,能量耗 散率,得到了双方程模型,即k-ε模型。70年代得 到应用。
项雷诺应力 u iu j ,而不引入新的变量。
雷诺方程数值模拟(RANS)
• 涡粘性假设(紊动粘性)
– 对应层流中切应力与流速梯度关系的公式: d u dy
– 引入一个涡粘度 t ,将紊流中的雷诺应力与流场 中的时均流速梯度建立下述关系:
t uvt dduyt dduy
t
-uiuj
t
xuij
xi
t /t
– 对于热输运问题,称为紊动普朗特数;对于质量输 运问题,称为紊动施密特(Schmidt)数。实验表明, 在流场中各点,甚至在不同型式的水流中,紊动普 朗特数几乎不变。
概述
• DNS之外的其他方法
– 大涡模拟LES。 – DNS可以获得尺度大于网格尺度的紊流结构,但却
无法模拟小于该网格尺度的紊动结构。 – 大涡模拟的思路是: DNS模拟大尺度紊流运动,而
利用亚网格尺度模型模拟小尺度紊流运动对大尺度 紊流运动的影响。 – LES较DNS拟占计算机的内存小,模拟需要的时间 也短,并且能够得到较雷诺平均模型更多的信息。
uj xi
雷诺方程数值模拟(RANS)
• 涡粘性假设(紊动粘性)
– 当i=j时,假设不合理,引入湍流脉动所产生的压力
t
பைடு நூலகம்
-uiuj
t xuji
uj xi
ptij
– 紊流涡粘度与层流中的粘度相对应,也可称为表观
粘度。粘度是流体本身的物理特性,与流动情况无
关。但是涡粘度则不是流体的物理性质,而是紊流
二、紊流模型
• 1、概述 • 2、雷诺方程数值模拟(RANS) • 3、大涡模拟(LES) • 4、直接数值模拟(DNS)
概述
• 紊流运动的瞬时连续方程和运动方程
• 连续方程 • 运动方程
u 0,ui 0 xi
du f 1 p 2u
dt
ui t
uj
ui xj
fi
1 p
xi
2ui
xjxj
切细节,但在实际工程中具有重要意义的并不是紊 流的一切细节,而是紊流对于时间的平均效应。 – 雷诺(Osborne Reynolds)建议用统计方法将N -S方程取时间平均。但取平均的过程产生了新问 题:方程增加了新未知项,时均方程组不再封闭, 因此各种类型的紊流模型应运而生。 – 紊流模型可定义为一组方程,这组方程确定时均流 方程中的紊动输运项,从而封闭时均流动方程组 (零方程、单方程、双方程)。
概述
• 紊流模型分类
– 直接数值模拟(DNS)。直接求解N-S方程,必 须采用很小的时间步长与空间步长。
– 大涡模拟(LES)。 DNS模拟大尺度运动,亚网格 上模拟小尺度运动。
– Reynolds时均方程法(RANS)。将N-S方程对 时间平均,并通过一些假定建立模型,是目前工程 中所采用的基本方法。根据对雷诺应力的处理方式 不同,分为基于涡粘性假设的模型和应力输运模型 两类,根据引入方程的数量,前者又分为零方程、 单方程和双方程模型。
概述
• 雷诺时均方程模型发展史
– 1895年,Reynolds发表对紊流研究结果的文章。
– 1897年,Boussinesq((布辛涅司克))涡粘性假定。
– 他们都没有尝试对雷诺方程进行系统的求解,许多 粘性流动的机理尚未清楚。
– 1925年,Prandtl 提出了混合长度理论来计算涡粘 性,为早期的研究奠定了基础。早期贡献最为显著 的还有von Karman,1930 年提出相似性假定。此 类模型并没有引入微分方程,称为零方程或者代数 模型。
的一种流动特性,决定于流动的时均流速场和边界
条件。
– 引入涡粘性假设并未构成紊流模型,只是提供了构 造紊流模型的基础,但使模拟紊动应力问题转化为 确定的分布。
雷诺方程数值模拟(RANS)
• 紊动扩散概念
– 将紊动热(或质量)输运与紊动动量输运直接类比, 假设热(或质量)输运与被输运的量有关:
ui'
概述
• 直接求解连续方程和运动方程,即 DNS(Direct Numerical Simulation)。
• 原则上讲,DNS并无理论上的困难。
– 一方面,描述紊流运动的精确的微分方程已经得出, 即N-S方程;从数学观点看来,紊流就是N-S方 程的通解。
– 另一方面,数值计算方法的发展,已足以求解N-S 方程。
– 但是,现代计算机的储存的能力和运算速度尚不足 以求解任何一个实际的紊流问题。
概述
• DNS实际应用中的困难
– 紊流运动所包含的单元,比流动区域的尺度要小得 多。为了用数值计算方法求解紊动单元的运动要素, 数值计算的网格必须比紊动单元的尺度更小。如此 之多的网格点贮存各种变量,远远超过了现代计算 机的内存容量,而且,随着网格点的增加,算术运 算的次数显著增多,所需计算的时间令人望而却步。
– 周培源(1945)和Rotta(1951),绕过 Boussinesq 涡粘性假定,提出了一个描述紊流切 应力张量演化的微分方程,即雷诺应力张量,得到 了应力输运模型,也称为二阶封闭或者二阶矩封闭 模型。70年代得到应用。
– 塞弗曼(P G Saffmann)在1977年曾经预言,大 约在本世纪末,高度发展的计算机将有足够的能力 用数值计算方法求解紊流的精确方程,但在较近的 未来,精确地求解紊流问题无疑是不可能的。
– 到现在,DNS算法已有很大发展。
概述
• DNS之外的其他方法
– 求解时均N-S方程(即雷诺方程)。 – 作为紊流的通解, N-S方程描述了流体运动的一
雷诺方程数值模拟(RANS)
• 雷诺时均方程的封闭问题
– 将瞬时值写成时均值+脉动值uu u,代入连
续方程和运动方程,并对方程两边取平均:
ux uy uz 0 x y z
u ti uj x uij fi 1 x pi xj x uij uiuj
– 引入数学方程或者代数公式,确定所产生的新未知
– 1945年,Prandtl 提出一个涡粘性依赖于紊动能k 的模型,建议使用一个偏微分方程的模型对精确的 k方程进行近似,得到了单方程模型,即k方程模型。
概述
• 雷诺时均方程模型发展史
– 1942年,Kolmogorov提出了第一个完整的紊流模 型,除了k方程,还引入了另外一个参数ε,能量耗 散率,得到了双方程模型,即k-ε模型。70年代得 到应用。
项雷诺应力 u iu j ,而不引入新的变量。
雷诺方程数值模拟(RANS)
• 涡粘性假设(紊动粘性)
– 对应层流中切应力与流速梯度关系的公式: d u dy
– 引入一个涡粘度 t ,将紊流中的雷诺应力与流场 中的时均流速梯度建立下述关系:
t uvt dduyt dduy
t
-uiuj
t
xuij
xi
t /t
– 对于热输运问题,称为紊动普朗特数;对于质量输 运问题,称为紊动施密特(Schmidt)数。实验表明, 在流场中各点,甚至在不同型式的水流中,紊动普 朗特数几乎不变。
概述
• DNS之外的其他方法
– 大涡模拟LES。 – DNS可以获得尺度大于网格尺度的紊流结构,但却
无法模拟小于该网格尺度的紊动结构。 – 大涡模拟的思路是: DNS模拟大尺度紊流运动,而
利用亚网格尺度模型模拟小尺度紊流运动对大尺度 紊流运动的影响。 – LES较DNS拟占计算机的内存小,模拟需要的时间 也短,并且能够得到较雷诺平均模型更多的信息。
uj xi
雷诺方程数值模拟(RANS)
• 涡粘性假设(紊动粘性)
– 当i=j时,假设不合理,引入湍流脉动所产生的压力
t
பைடு நூலகம்
-uiuj
t xuji
uj xi
ptij
– 紊流涡粘度与层流中的粘度相对应,也可称为表观
粘度。粘度是流体本身的物理特性,与流动情况无
关。但是涡粘度则不是流体的物理性质,而是紊流
二、紊流模型
• 1、概述 • 2、雷诺方程数值模拟(RANS) • 3、大涡模拟(LES) • 4、直接数值模拟(DNS)
概述
• 紊流运动的瞬时连续方程和运动方程
• 连续方程 • 运动方程
u 0,ui 0 xi
du f 1 p 2u
dt
ui t
uj
ui xj
fi
1 p
xi
2ui
xjxj
切细节,但在实际工程中具有重要意义的并不是紊 流的一切细节,而是紊流对于时间的平均效应。 – 雷诺(Osborne Reynolds)建议用统计方法将N -S方程取时间平均。但取平均的过程产生了新问 题:方程增加了新未知项,时均方程组不再封闭, 因此各种类型的紊流模型应运而生。 – 紊流模型可定义为一组方程,这组方程确定时均流 方程中的紊动输运项,从而封闭时均流动方程组 (零方程、单方程、双方程)。
概述
• 紊流模型分类
– 直接数值模拟(DNS)。直接求解N-S方程,必 须采用很小的时间步长与空间步长。
– 大涡模拟(LES)。 DNS模拟大尺度运动,亚网格 上模拟小尺度运动。
– Reynolds时均方程法(RANS)。将N-S方程对 时间平均,并通过一些假定建立模型,是目前工程 中所采用的基本方法。根据对雷诺应力的处理方式 不同,分为基于涡粘性假设的模型和应力输运模型 两类,根据引入方程的数量,前者又分为零方程、 单方程和双方程模型。
概述
• 雷诺时均方程模型发展史
– 1895年,Reynolds发表对紊流研究结果的文章。
– 1897年,Boussinesq((布辛涅司克))涡粘性假定。
– 他们都没有尝试对雷诺方程进行系统的求解,许多 粘性流动的机理尚未清楚。
– 1925年,Prandtl 提出了混合长度理论来计算涡粘 性,为早期的研究奠定了基础。早期贡献最为显著 的还有von Karman,1930 年提出相似性假定。此 类模型并没有引入微分方程,称为零方程或者代数 模型。
的一种流动特性,决定于流动的时均流速场和边界
条件。
– 引入涡粘性假设并未构成紊流模型,只是提供了构 造紊流模型的基础,但使模拟紊动应力问题转化为 确定的分布。
雷诺方程数值模拟(RANS)
• 紊动扩散概念
– 将紊动热(或质量)输运与紊动动量输运直接类比, 假设热(或质量)输运与被输运的量有关:
ui'
概述
• 直接求解连续方程和运动方程,即 DNS(Direct Numerical Simulation)。
• 原则上讲,DNS并无理论上的困难。
– 一方面,描述紊流运动的精确的微分方程已经得出, 即N-S方程;从数学观点看来,紊流就是N-S方 程的通解。
– 另一方面,数值计算方法的发展,已足以求解N-S 方程。
– 但是,现代计算机的储存的能力和运算速度尚不足 以求解任何一个实际的紊流问题。
概述
• DNS实际应用中的困难
– 紊流运动所包含的单元,比流动区域的尺度要小得 多。为了用数值计算方法求解紊动单元的运动要素, 数值计算的网格必须比紊动单元的尺度更小。如此 之多的网格点贮存各种变量,远远超过了现代计算 机的内存容量,而且,随着网格点的增加,算术运 算的次数显著增多,所需计算的时间令人望而却步。