流体力学第八章 粘性不可压缩流体绕物体的流动

合集下载

流体力学第八章 粘性流体动力学基础

流体力学第八章  粘性流体动力学基础
对形心取矩,忽略了 质量力引起的力矩:
力矩方程为:
pz
pz z
dz
zy
zy
z
dz
z
py
yz
形心 dx
dy
zy
pz
yz
yz
y
dy
py
p y y
dy
y
yzdz
dy 2
( yz
yz
y
dy)dz
dy 2
zydy
dz 2
( zy
zy
z
dz)dy
dz 2
0
略去高阶小量后得: yz zy
本问题是N-S 方 程的精确解之一
在上述条件下,流动将是二元的,质量力可略 去不计,N-S方程和连续方程可简化为:
vx
vx x
vy
vx y
1
p x
(
2vx x2
2vx y 2
)
(a)
vx
2vy y2
2vy z2
)
(8--12)
vz t
vx
vz x
vy
vz y
vz
vz z
Z
1
p z
(
2vz x2
2vz y2
2vz z 2
)
N-S方程的矢量形式:
可压缩
vr
(vr
)vr
r F
1
p
2vr
( vr )
(8-13)
t
3
不可压缩
vr
(vr
)vr
r F
1
p
2vr
t
(8-14)
讨论 1.方程(8-12)的求解:
dz

流体力学第8章中文版课件

流体力学第8章中文版课件

Chapter 8: External flows
14
8.3 绕淹没体的流动
分离前的湍流边 界层 分离前的层流 边界层
2013-11-25
Chapter 8: External flows
15
8.3 绕淹没体的流动
2013-11-25
Chapter 8: External flows
16
8.3 绕淹没体的流动
W FD
sphere volume CD V 2 A
4 3 1 S water R CD V 2R 2 3 2
1 2
8RS water V 3C D
2013-11-25
1/ 2
8 0.15 1.02 9800 3 1.20 CD
Re
VD


129 0.3 2.42 10 6 1.6 10 5
V 129 m/s
2013-11-25 Chapter 8: External flows 20
8.3 绕淹没体的流动
求解:(b) 对于球在水中的下落情况,则必须考虑施加在球体上的与阻力FD 同方向的浮力 B 的作用:
如果物体形状上有一 个突然的变化,分离 点将出现在形状突然 变化点或其附近。 另外,分离后流 体在某一个位臵 上又会重新附着 在物体上。
2013-11-25
Chapter 8: External flows
10
8.2 分离
在分离点的上游,壁面附 在分离点的下游,壁面附 近的 x方向上的速度分量 近的 x方向上的速度分量在 负 x 方向,因此在正 x 方向,因此 壁面上 壁面上的 的 u/y一定是负的。 u/y是正的。

流体力学 第八章 绕流运动

流体力学 第八章 绕流运动

第八章绕流运动一、应用背景1、问题的广泛存在性:在自然界和工程实际中,存在着大量的流体绕物体的流动问题(绕流问题),如:飞机在空气中的飞行、河水流过桥墩、大型建筑物周围的空气流动、植物护岸(消浪,船行波),粉尘颗粒在空气中的飞扬和沉降,水处理中固体颗粒污染物在水中的运动。

(一种:流体运动;另外一种:物体运动),我们研究,将坐标系固结于物体上,将物体看成静止的,讨论流体相对于物体的运动。

2、问题的复杂性上一章的内容中可以看出,流体力学的问题可以归结为求解在一定边界条件和初始条件下偏微分方程组的求解。

但描述液体运动的方程式非常复杂的:一方面,是方程的非线性性质,造成方程求解的困难;另一方面,复杂的边界条件和初始条件都给求解流体力学造成了很多麻烦。

迄今为止,只有很少数的问题得到了解决。

平面泊萧叶流动,圆管coutte流动等等。

而我们所要解决的绕流问题正是有着非常复杂的边界条件。

3、问题的简化及其合理性流体力学对此的简化则是,简化原方程,建立研究理想液体的势流理论。

实际液体满足势流运动的条件:粘性不占主导地位,或者粘性还没有开始起作用。

正例:远离边界层的流体绕流运动、地下水运动、波浪运动、物体落入静止水体中,水的运动规律研究。

反例:研究阻力规律、能量损失、内能转换等等。

圆柱绕流(经典之一)半无限长平板绕流(经典之二)分成两个区域:一个区域是远离边界的地方,此区域剪切作用不明显,而且流体惯性力的影响远远大于粘性力的影响(理想液体)(引导n-s方程);另一个是靠近边界的地方(附面层,粘性底层),此区域有很强烈的剪切作用,粘性力的影响超强,据现代流体力学的研究表明,此区域是产生湍流的重要区域,有强烈的剪切涡结构,但此区域只有非常薄的厚度。

此区域对绕流物体的阻力、能量耗损、扩散、传热传质都产生重要影响。

4、本章的主要研究内容(1) 外部:理想液体,(简化方法,求解方式)、(2) 内部:附面层理论,(简化方法,求解方式,求解内容,现象描述) (3) 两者的衔接。

粘性流体的不可压缩流动

粘性流体的不可压缩流动

Chapter 9-1 粘性不可压缩流体流动§1概述一、粘性不可压缩流动模型1、关于粘性 粘性摩擦的存在必导致绕流阻力的存在,运动的衰减及涡量的扩散。

在大e R 数下,惯性力>>粘性力,采用理想流体模型,理想流体理论对不脱体绕流情况下的升力,压力分布和速度分布给出了符合实际的结果,但在阻力等与粘性效应相关的问题上却无能为力。

因而,在研究阻力等起源于粘性的现象时须抛弃理想流体假设。

在小e R 数和中e R 数情况下,粘性作用不可忽略。

2、关于不可压缩流动(流体的压缩性对流动的影响可略)液体压缩系数小,一般可认为不可压缩(极端情况如激波等除外)。

气体在低速运动(速度远小于声速)、非定常时速度变化缓慢,且重力方向上流场的尺度<10km 时,可略其压缩性。

(当研究对流层(~10km )内大气运动时,不能忽略重力场引起的压缩效应)。

3、基本方程组和边界条件均质不可压缩流体.const ρ=,且温度变化小,const μ=,故有20V dV pF V dt γρ⎫∇⋅=⎪⎬∇=-+∇⎪⎭求速度和压力场的完备方程组。

能量方程22:dUk T S S dtρμ=∇+ 用于求温度场 本构方程 2P p I S μ=-+ 用于求应力边界条件:在固壁表面上,流体的法向和切向速度分别等于固体表面的对应速度分量。

在自由表面上,0, 0nn n p p p τ=-=。

二、粘性流动分类,求解问题的几种途径层流:流体运动规则、稳定,各部分分层流动互不掺混,质点轨迹光滑。

脉线清晰 湍流:流体运动极不规则、极不稳定,伴有高频扰动,各部分激烈掺混,质点轨迹杂乱无章。

决定流动状态的参数是e R 数(Batchlor page255),e R <<2000 一定是层流,此时粘性力足以保持流动的稳定。

层流:极少有准确解(某些特殊的简单问题,非线性方程得以简化) 近似解法:大e R 数,边界层理论小e R 数,部分或全部忽略惯性力。

第八章 粘性不可压缩流体的层流运动

第八章 粘性不可压缩流体的层流运动

8.6 不可压缩粘性流体在无穷长直圆管内流。

由实验知,其璧面传热系数h 与圆管的直径D ,热传导系数k,流体的平均速度U ,密度ρ,粘度系数μ和流体比热c 有关,其中h 具有h/D 的量纲。

试由量纲分析证明 P r ).(R e ,f Nu = 式中khD Nu =叫做努塞尔特(Nusselt )数,μρUD =Re 是雷诺数,kc μ=Pr 是普朗特数。

解:由题意:,,,,,(][c U k D f h μρ=此式中有n=6个物理量,其中含4=r 个基本量纲,按π定理可简化为2=-r n 个无量纲间的函数关系。

记质量,长度,时间和温度的基本量纲分别为K T L M ,,,写出各量的量纲如下:[]L D =,[][]13)/(--==KMLTLK W k ,[]1-=LT U ,[]3-=ML ρ,11][--=TML μ,[]13--=⎥⎦⎤⎢⎣⎡=KMT D k h ,122][-=KT L c 。

现取D ,k ,U ,ρ为基本量,将其余各量与这些基本量组合成无量纲量。

例如,设]ξγβαρ][][][][U k D h =,列出此式两侧的量纲有:ξγβαβγβξβ3313-++---+--=LKTMKMT显然两侧的幂次应该分别相等:⎪⎪⎩⎪⎪⎨⎧=-++-=--=--=+031331ξγβαβγβξβ解得⎪⎪⎩⎪⎪⎨⎧===-=0011ξγβα,即[]][][1k D h -=,于是khD Nu =构成一个无量纲量。

同理: ),,,,,(][1c U k D h f μρ=,取μ,,,k U D 为基本量,将其余各量与这些基本量组合成无量纲量。

设[]ξγβαμρ][][][][k U D =,列出此式两侧的量纲有:ββξγβαξβ----+++-=KTLMMLr333两侧的幂次应该分别相等:⎪⎪⎩⎪⎪⎨⎧=-=---=-++=+003331βγβξγβαξβ解得⎪⎪⎩⎪⎪⎨⎧====1000ξγβα,8.8截面为半圆形的无限长直管中的不可压缩流体做层流运动,沿管轴方向某一长度l 上的压降为p ∆。

二维不可压缩粘性流体绕钝体流动的数值模拟

二维不可压缩粘性流体绕钝体流动的数值模拟
U /,其 中 D为 圆柱直 径 ; D ̄ , U为 自由来 流速 度 , 过 通 改变来 流速度来 得 到不 同的雷诺 数 。 3 2 计 算域及 网格 划分 .
数时 , 光滑 圆柱体的周期性尾流只是雷诺数的函数。 按 圆柱 体 直径计 算 的雷 诺数 很 小 时 , 体并 不 脱 离 流
界区, 此时柱 体表 面上 的边 界层 为层 流 , 而柱 体后 面
的涡街 已完全转变为湍流, 并按一定的频率发放漩 涡; × 0 < × 0 称为过渡区, 3 1 3 1。 此时柱体表面上的
边 界层 也 已变 为 湍 流 , 离 点 向后 移 , 力 显 著 下 分 阻
符合重新分离 , 这一类物体 的绕流在较大的雷诺数 范 围 内变化不 明显 。
3 计算 实例
3 1 控 制 方程 .
通 过 控 制 流体 由计 算域 的流 速来 控 制 雷诺 数 ,
得 到的尾 迹及漩 涡脱 落 图如 图 3所 示 。 当R e=1 , 流 中 有 一 对 稳 定 的 弗 普 尔 旋 时 尾
第8 期
曹广龙等: 二维不可压缩粘性流体绕钝体流动的数值模拟
C= p竿堕
式 中 : 和 为来 流 的静 压 和速度 ; 底 部压 强 。 P p为
对 于不 可压缩 粘性 流 体 , 直角 坐标 系下 , 在 其运 动规律 受 N—S方 程 控 制 , 续 性 方 程 和 动 量 方 程 连
分别 为 :
O t u
_

对于三维钝体绕流 , 压差阻力仍然是总阻力 的主要 部 分 。若对 物 体形 状 进 行 流线 型 处 理 ; 压 差 阻力 则
定 义是 :
. = s
决定圆柱绕 流流 态的是雷诺数 ( e 的值 , R) 当

流体力学第八章答案

流体力学第八章答案

流体力学第八章答案【篇一:流体力学第8、10、11章课后习题】>一、主要内容(一)边界层的基本概念与特征1、基本概念:绕物体流动时物体壁面附近存在一个薄层,其内部存在着很大的速度梯度和漩涡,粘性影响不能忽略,我们把这一薄层称为边界层。

2、基本特征:(1)与物体的长度相比,边界层的厚度很小;(2)边界层内沿边界层厚度方向的速度变化非常急剧,即速度梯度很大;(3)边界层沿着流体流动的方向逐渐增厚;(4)由于边界层很薄,因而可以近似地认为边界层中各截面上压强等于同一截面上边界层外边界上的压强;(5)在边界层内粘性力和惯性力是同一数量级;(6)边界层内流体的流动与管内流动一样,也可以有层流和紊流2种状态。

(二)层流边界层的微分方程(普朗特边界层方程)??v?vy?2v1?p?vy?????vx?x?y??x?y2????p??0?y???v?vy???0?x?y??其边界条件为:在y?0处,vx?vy?0 在y??处,vx?v(x)(三)边界层的厚度从平板表面沿外法线到流速为主流99%的距离,称为边界层的厚度,以?表示。

边界层的厚度?顺流逐渐加厚,因为边界的影响是随着边界的长度逐渐向流区内延伸的。

图8-1 平板边界层的厚度1、位移厚度或排挤厚度?1?1?2、动量损失厚度?2?vx1?(v?v)dy?(1?)dy x??00vv?2?1?v2???vx(v?vx)dy???vxv(1?x)dy vv(四)边界层的动量积分关系式??2???p?vdy?v?vdy?????wdx xx??00?x?x?x对于平板上的层流边界层,在整个边界层内每一点的压强都是相同的,即p?常数。

这样,边界层的动量积分关系式变为?wd?2d?vdy?vvdy?? x?x??00dxdx?二、本章难点(一)平板层流边界层的近似计算根据三个关系式:(1)平板层流边界层的动量积分关系式;(2)层流边界层内的速度分布关系式;(3)切向应力关系式。

不可压缩粘性流体的运动微分方程名师优质资料

不可压缩粘性流体的运动微分方程名师优质资料

根据边界层的特征,在边界层内惯性项和粘性项具有 同样的数量级,由方程组(8-37)可知,必须使 1 Re l 和 2
同数量级,所以 l ~ 1 Re l ,即 反比于 雷诺数越大,边界层相对厚度越小。
Re l
。这表明,
这样,将式(8-37)中的某些项略去,再变换成有量纲 量,便得到了层流边界层的微分方程(称为普朗特边界层方 程): v x v x 2vx 1 p
(8-36)
v
y
o

x
l
vx
x
图8-11
推导层流边界层的微分方程用图
可以利用边界层每一处的厚度都很小的特征,来比较方程 组(8-36)中各项的数量级,权衡主次,忽略次要项,这 样便可大大简化该方程组。
即 l或 足不等式 边界层的厚度 与平板的长度 l 相比较是很小的, ,而 l 1 y的数值限制在边界层内,并满
p xx
dy
xz
zx
fz
xy
xy
fy
xy
A
y
yx
dx x
zx zx dz z
fx
o z
第一个下标表示应力所在平面的法线方向 第二个下标表示应力本身的方向。
y
yx
yx y
dy
xy
dy
M
xy
xy x
dx
dx
yx
现在根据边界层的特征,利用不可压缩粘性流体 的运动微分方程来研究边界层内流体的运动规律。为 简单起见,只讨论流体沿平板作定常的平面流动,x轴与 壁面相重合,如图8-11所示。假定边界层内的流动全 是层流,忽略质量力,则不可压缩粘性流体平面定常 流动的微分方程和连续方程为
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
V∞ y
无滑移边界条件 y = 0时,u = 0
a=0
前缘
δ(x)
层流 x
边界层与势流衔接处, y = δ 时,u = V∞
V∞ = bδ + cδ + dδ
2 3
y = δ 时,τ = 0
流体力学
du τ =μ =0 dy
顺流平板层流边界层4
b + 2cδ + 3dδ 2 = 0
y V∞ 前缘 δ(x) 层流 x
∞ 0
δ* = ∫
u (1 − )dy U
δ*
μ=0 u=U
边界层内由粘性影响减少的流量=理想流 体流过物面时表面向外移动 δ*减少的流量
δ =∫
*
δ
0
u (1 − )dy U
流体力学
边界层内的厚度4
动量损失厚度
y dy y
u = u(y)
U
U
θ
μ≠0
μ=0 u=U
边界层内由于粘性的影响,动量流量比 理想流体流经该区域时有所减少
固体壁面上 y = 0 时
∂ u U dU =− 2 ν dx ∂y
2
dU =0 由 dx
∂ u = 2c = 0 2 ∂y
2
3V∞ V∞ a=c=0 , b= , d=− 3 2δ 2δ
流体力学
顺流平板层流边界层5
V∞ u= 2δ ⎛ y3 ⎞ ⎜3y − 2 ⎟ δ ⎠ ⎝
补充方程2-壁面切应力方程
3 ∞
= 1.771 ρμaV b
边界层内的流动状态1
V∞ 层流边界层 转捩区
紊流边界层 x
平板前缘开始
x ↑, δ ↑, Re ↑
层流边界层 转捩点 紊流边界层
流体力学
边界层内的流动状态2
边界层流动状态的判据
Re x = Ux

ν
Reδ =
ν
其中,x 为物面上一点到前缘的距离 顺流平板
流体力学
(Re x )cr = 3 × 105 ~ 3 × 106
解:沿边 a 方向拖动
1 2 D = C D ρV∞ A × 2 2
Ua b a
其中:V∞ = U a
Re l =
流体力学
A = ab
Uaa
ν
1.328 1.328 = CD = Re l Uaa ν
顺流平板层流边界层-例题1
1.328 2 3 Da = ρU a ab = 1.328 ρb U a aν Uaa ν
假设 定常不可压 二元边界层 物面曲率很小
B x y U A
C
δ
δ+ dδ
D x+ dx
x
CV所受外力之和=净流出CV的动量流率
流体力学
边界层动量积分方程2
控制体在 x 方向所 受合力
dp Fx = −δ dx − τ w dx dx
y U A
pAC
C
α
δ δ+ dδ
p
B x
τw
p+ dp
D x+ dx x
1.328 CD = Re l
3 dD = 0.664 ρμV∞
(a − ay b ) dy
流体力学
顺流平板层流边界层-例题2
由流动的对称性,以 及三角形平板两侧均 受阻力
V∞ y x o dy l x a 2b
D = 4∫
b 0
1.328 a ⎞ ⎛ 3 ρμV∞ ⎜ a − y ⎟ dy b ⎠ 2 ⎝
2
)
(∫ )
δ
0
udy = −τ w + δρ U
dU dx
边界上速度分布 U 只是 x 的函数
d ρU dx
流体力学
(∫ )
δ
0
δ d dU udy = ρ U ∫ udy − ρ 0 dx dx
(
)

δ
0
udy
边界层动量积分方程的其它形式2
d =ρ dx
(∫
δ
0
dU Uudy − ρ dx
)

第八章 粘性不可压缩流体绕物体 的流动
内流
在固壁限定的空间内流动
管流、通道流、各种动力设备内部流动
外流
流体从物体外部流过
飞机在大气中飞行,潜艇在水中航行
流体力学
概述1
粘性、不可压、定常、绕流
边界层的概念、动量积分方程、曲壁边界 层分离、绕流物体的升力、阻力
基础知识
不可压缩流体积分形式控制方程, 园管内速度的幂次分布规律,雷诺 数,理想流体圆柱绕流
控制体 x 方向动量的净流出率
d dx
流体力学
(∫
δ
0
d ρ u dy dx − U dx
2
)
(∫
δ
0
ρ udy dx
)
边界层动量积分方程3
边界层动量积分方程
d dx
(∫
δ
0
d ρ u dy − U dx
2
)
(∫
δ
0
dp ρ udy = −τ w − δ dx
)
适用条件 不可压定常 二元边界层,物面曲率很小 对层流边界层和紊流边界层均适用
⎛ du ⎞ τw = μ ⎜ ⎟ ⎝ dy ⎠ y=0
d δ u dx ∫0 V∞
流体力学
3 V∞ τw = μ 2 δ
⎛ u ⎜1 − V∞ ⎝ ⎞ τw ⎟ dy = 2 ρV∞ ⎠
顺流平板层流边界层6
1 2 140 μ x+C δ = 2 13 ρV∞
边界条件
x = 0时, δ = 0
μ x δ = 4.641 ρV∞
ρ xV∞ 由 Re x = μ
流体力学
δ = 4.641
x Re x
顺流平板层流边界层7
μ x δ = 4.641 ρV∞
δ = 4.641
x Re x
层流边界层厚度与流体性质、来流速度 及距前缘的距离有关
δ ∝ x1 2
壁面切应力
流体力学
3 V∞ τw = μ 2 δ
顺流平板层流边界层8
τ w = 0.3232 μρV
δ
0
udy
代入边界层动量积分方程
d ρ dx
(∫
δ
0
d u dy − ρ dx
2
)
(∫

δ
0
dU Uudy + ρ dx
)

δ
0
udy
dU dU = −τ w + ρ = −τ w + δρ U dx dx

δ
0
Udy
d ⎡ δ 2 ⎤ − ρ dU ⎡ δ (U − u) dy ⎤ = −τ Uu − u ) dy −ρ w ⎥ ⎢ ∫0 ( ⎥ ⎣ ⎦ ⎦ dx ⎢ ∫0 dx ⎣
∂u ∂u dU ∂ u u +v =U +ν 2 ∂x ∂y dx ∂y
2
y = 0时,u = 0,v = 0
流体力学
∂2u U dU =− 2 ν dx ∂y
速度分布在边界上应满足的条件4
∂u ∂u dU ∂ 2u u +v =U +ν 2 ∂x ∂y dx ∂y
对 y 求导
∂ u ⎛ ∂u ∂ v ⎞ ∂2u ∂2u ∂3u ⎜ ∂x + ∂y ⎟ + u ∂x ∂y + v ∂y 2 = ν ∂y 3 ∂y ⎝ ⎠
流体力学
顺流平板层流边界层10
C D = 1.293 1 Re L
层流边界层微分方程勃拉休斯精确解
边界层厚度 平板阻力系数
流体力学
μ δ = 5.0 x ρV∞
CD = 1.328 Re L
顺流平板层流边界层-例题1
例:矩形平板边长为a和b,若在静止流体中沿边 例:矩形平板边长为a和b,若在静止流体中沿边 a方向以Uaa拖动平板与沿边b方向以Ubb拖动平 a方向以U 拖动平板与沿边b方向以U 拖动平 板的阻力相等,求Uaa//Ubb,层流边界层。 板的阻力相等,求U U ,层流边界层。
从物面沿外法线到速度达到势流速度99% 处的距离 边界层厚度沿流动方向不断增大
流体力学
边界层内的厚度2
位移厚度
y dy y
u = u(y)
U
U
δ*
μ≠0
μ=0 u=U
边界层内由于粘性的影响,质量流量比 理想流体流经该区域时有所减少
流体力学
边界层内的厚度3
位移厚度 (排挤厚度)

*
U
=∫
∞ 0
(U − u)dy

流体力学
θ =∫
δ
0
u u (1 − )dy U U
δ = ∫ (1 −
∗ 0
δ
u )dy U
动量积分方程的其它形式
d dU 2 ρ (U θ ) + ρ Uδ ∗ = τ w dx dx

τw dθ 1 dU ∗ + ( 2θ + δ ) = ρU 2 dx U dx
令 H = δ ∗ θ 形状因子
y = 0时,u = 0,v = 0
流体力学
∂3u =0 3 ∂y
8.4 顺流平板层流边界层
问题
均匀来流
V∞ = C p∞ = C
V∞ 前缘 δ(x) 层流 x y
粘性、不可压、定常、二元 层流边界层,板长为L
边界层积分 方程
流体力学
τw dθ 1 dU ∗ + ( 2θ + δ ) = ρU 2 dx U dx
由 Da = Db
1.328 ρb U aν = 1.328 ρa U bν
相关文档
最新文档