第四章 高速可压无粘流
第四章 高速可压流
作业
§4.2 一维等熵绝热流
4.2.1 一维等熵绝热流的能量方程
4.2.2 一维定常绝热流的参数间的基本关系式
27
作业
§4.2 一维等熵绝热流
一维定常流动是指流场中的物理量只是一个空间坐标的函数的 流动。 一维可压缩定常流比不可压缩流动复杂些,有四个流动参数,
v, , p, T
利用状态方程、连续方程、动量方程以及下面导出的一维等熵 绝热流的能量方程即可确定上述四个流动参数。 下面就推导一维等熵绝热流的沿流线方向的能量方程
作业
§4.1 热力学基础
定压比热容和定容比热容之间存在一定的关系,
c p cV R
绝热指数或比热比
焓用比热比来表达为,
cp cV
h
1
RT
p
1
3. 绝热过程:在该过程中d(q =0。
p
25
C
作业
§4.1 热力学基础
热力学第二定律和熵 热力学第二定律指明能量相互转化是有条件的、有方向性的, 即一个方向的变化过程可以实现,而逆方向的变化过程或者不 能实现或者只能有条件地实现。 例如热量只会从高温的物体传到低温的物体里去,而不会自发 地反向传递。因此,在热力学上有可逆过程和不可逆过程之分。 如果将变化过程一步步地倒回去,介质的一切热力学参数均回 到初始值且外界情况也都复旧,则过程可逆,否则就是不可逆 过程。 流体力学中常引用熵这个状态参数的变化来叙述热力学第二定 律。
16
作业
§4.1 热力学基础
1. 等容过程:在该过程中d(1/ =0,即该过程中所吸收的热 量 dq 都用来增加气体的内能 。
1 dq du pd
流体力学教学大纲 - 南京航空航天大学精品课程建设
6
课程教学方 式 主要适用专 业
理论讲授
流体力学等相关专业 飞行器设计等相关专业
飞行器动力工程
7
4
努利方程,基本解叠加解法,理解和掌握升力定理,了解粘流基础知识。 本章主要知识点:伯努利方程,叠加原理,升力定理,附面层,层流, 紊流。 本章重点:伯努利方程及应用,基本解叠加和库塔-儒可夫斯基升力定 理。 实验内容:不可压流水槽流动,流动现象分析 第四章 高速可压无粘流(12 学时+2 学时实验) 本章讲述高速可压无粘流。在高速时,流体运动必须考虑热力学特性。 先介绍与流体运动相关的热力学基础知识,再讲述一维等商绝热流。然后讲述 正激波、斜激波与膨胀波(波系实验演示) ;最后介绍喷管及 Δ 超音速风洞。 要求熟练掌握一维等商绝热流所满足的方程,激波前后参数关系,理解和掌握 超音速喷管原理。 本章主要知识点:状态方程,内能,焓,熵,熵增原理,总压,小 扰 动 与马赫锥,马赫波,普朗佗(Prandtl)激波公式,拉瓦尔喷管 本章重点:一维等商绝热流,正激波,拉瓦尔喷管原理 实验内容:风洞设备参观及工作原理介绍,含低速、高速、高超声速 Δ 为选讲内容,其余均为必讲内容。
1
前修课程、能力和知识结构要求:
先修课程:高等数学、工程数学、理论力学。 主要能力和知识结构:学生在学习本课程时应具备高等数学中的解析几何、微积分、 常微分方程、理论力学有关力与运动等方面的基础知识; 学生如果具有工程数学中, 矢量代数、矢量分析和数理方程等方面的知识,将有助于本课程的深入学习。作为航 空类本科生的专业基础课,还要求学生对航空航天等相关的工程背景有所了解。
课程内容共分 4 章: 第一章 流体力学基础知识(7 学时) 本章讲述流体力学的预备知识,先介绍流体力学的基本任务,研究方 向和流体力学的发展概述;然后介绍流体介质,Δ 气动力系数,Δ 矢量积分 知识;最后引入控制体,流体微团及物质导数的概念;要求掌握学习流体力学 入门的基础知识和预备知识。 本章主要知识点:连续介质假设,流体的密度、压强和温度,气体的 压缩性、粘性和传热性, 流体的模型化。 本章重点:流体的物理属性及模型化。 第二章 流体运动基本方程和基本规律(10 学时+2 学时实验) 本章讲述流体控制方程和基本规律,包括连续方程、动量方程和能量 方程的推导,方程的基本解法介绍,流体微团运动分析和 Δ 旋涡运动(实验 演示) 。要求掌握三大基本方程的推导,熟练掌握流体微团运动特点,了解流 体旋涡运动。 本章主要知识点:三大基本方程,角速度,旋度,角变形率,位函数, 流函数,环量,旋涡运动。 本章重点:三大基本方程的推导,微团运动分析结论。 实验内容:实验设备介绍及基本实验演示 第三章 不可压流(10 学时+2 学时实验) 本章先讲述不可压无粘流(理想流体) 。内容有伯努利方程及应用、 流动控制方程、方程的基本解、基本解叠加、库塔-儒可夫斯基升力定理;然 后讲述低速粘流基础;最后 Δ 介绍真实流动,给出与理想流的比较。与本章 内容相结合,实验有雷诺实验,测速实验和真实流动演示等。要求熟练掌握伯
北航水力学 第四章理想流体动力学和恒定平面势流解读
z1
p1
u12 2g
z2
p2
u22 2g
4.2.2 由动能定理推导理想流体的伯努利方程
推导过程同学们自学
z1
p1
u12 2g
z2
p2
u22 2g
本公式是由动能定理推导而得,它使伯努利方程有更加明确的 物理意义,说明伯努利方程是一能量方程。
第三节 元流伯努利方程的意义和应用
4.3.1 沿流线的伯努利方程的水力学意义
可见,在同一流线上各点的流函数为一常数,故等流函数线就是流线。
2、平面内任意两点流函数值的差等于通过这两点连线的流量。
y ABdrBnA x
d r dxi dy j
n cos i sin j dy i dx j
dr dr V ui v j
dq V
ndr
u
dy dr
v
dx dr
等 线和等Ψ线,这两族曲线互相垂直,构
成流网。
两族曲线所构成的正交网络,称为流网
流网的特征:
流网
等 线和速度矢量垂直,或者说, 等 线与等Ψ线(流线)垂直,
【例题】
已知90度角域内无粘流动,速度分布
ux kx uy ky
(k 0, x 0, y 0)
求:(1)判断该流场是否存在速度势函数, 若存在请给出并画出等势线;
流动。但粘滞性对流动 的影响很微小时,影响可以忽略。 --机械能守恒
引入势流的意义:使问题简化。
波浪运动,无分离的边界层外部的流动,多孔介质的流动(渗流) 等等可以看为势流。
4.4.1 流速势函数
以二维流动为例,根据流体运动学,它与无旋流动等价
由 ux 0 无旋流的条件→涡量 z 0
第四章 湍流流动
____ ____
uuzzuuxx
XX
xx
yy
zz
t xx
x
t yx
y
t z
x
z
——(5)
14
___
t xx ——湍流流动时x方向总法向应力。
___
r xx ——涡流粘性产生的附加法向应力。
___ ___ ___
t xx
,
t yx
,
t zx
——湍流时,总时均法向、切向应力的平均值。
表观运动粘度。
17
说明:
①涡流黏度与牛顿黏性定律中的动力黏度所表达的含义相同,但本质 不同。 ②涡流黏度不是流体的物性,而是与流道中流体所处的位置、流速及 边壁的粗糙度等因素有关的,是表示湍流中流体脉动程度的一个参数。 随时间和空间的变化很大,甚至有数量级的差别。除壁面附近外,涡 流黏度远大于分子黏度。
___ ___ ___
, , ——湍流时,法向、切向应力的时均值。
xx yx zx
(相当于层流时的应力值)
___ ___ ___
r xx
,
r yx
,
r zx
——脉动速度产生的法向、切向应力时均值。 (或附加应力时均值)
15
6.涡流粘度与混合长
宗旨:为求解上述方程,必须确立雷诺应力(脉动速度分量) 与时均速度梯度之间的关系。
7
4.湍流时的微分动量衡算方程
X方向的微分动量衡算方程
Dux X xx yx zx
D
x y z
ux
ux
ux x
uy
ux y
uz
ux z
X
xx
x
yx
第4章 高速可压缩流基础知识
的连续方程有
编
写 ρadA = (ρ + d ρ)(a − dv)dA
消去 dA 并略去二阶微量得
ad ρ = ρdv
(4-16)
再由动量方程得
dp = ρadv
(4-17)
由(4-16)和(4-17)两式消去 dv ,最后可得
西a2 = dp 北d ρ
(4-18)
工 这就是声速计算公式,它指明小扰动的传播速度与介质的压缩性大小有关, 空 业 即取决于压强变化与密度变化之比。对容易压缩的流体声速值小,对不易压缩 气 大 的流体声速值大,而完全不可压的流体声速为无限大。
气 业 基础 学 队 1 °C 一度所需的热量,单位是 J / (kg ⋅ K) 。由式(4-6),并取T = 0 时 e = 0 ,则
动 大 教 院 编 T 力 学 写 ∫ =e = 0 cV dT cVT 学 航 学 或
(4-7a)
基 天 团队 e2 − e=1 cV (T2 −T1)
(4-7b)
空 工业 基 天 团队 通过摩擦机械功可全部变成热,但热却不能百分之百地转变为功等等。据此, 气 大 础 学 编 在热力学上有可逆过程和不可逆过程之分。如果将变化过程一步步地倒回去, 动 学 教 院 写 介质的一切热力学参数均回到初始值且外界情况也都复旧,则是可逆过程,否 力 航 学 则就是不可逆过程。上面所说的高温物体向低温物体传热以及机械功通过摩擦 学 天 团 产生热均是不可逆过程。 基 学 队 热力学第二定律有许多表述方法。这里引用熵这个状态参数在不可逆过程 础 院 编 中的变化来叙述热力学第二定律。 教学 写 定义单位质量气体的熵为(非单位质量气体的熵 S 用大写)
气 业 基 学 3.焓值
动 大 础 院 在热力学,特别是气体动力学中,还常常引入另一个代表热含量的参数焓 h
空气动力学绪论PPT课件
0.3 空气动力学的发展进程
现代航空和喷气技术的迅速发展使飞行速度迅猛提高在 高速运动的情况下,必须把流体力学和热力学这两门学科 结合起来,才能正确认识和解决高速空气动力学中的问题。 1887-1896年间,奥地利科学家马赫在研究弹丸运动扰动 的传播时指出:在小于或大于声速的不同流动中,弹丸引 起的扰动传播特征是根本不同的。
高等数学计算方法大学物理理论力学绪论2学时第一章流体的基本属性和流体静力学6学时第二章流体运动学和动力学基础12学时第三章不可压缩无粘流体平面位流6学时第四章粘性流体动力学基础6学时第五章边界层理论及其近似6学时第六章可压缩高速流动基础14学时第七章高超音速流动基础4学时6学时总复习2学时陈再新刘福长鲍国华空气动力学航空工业出版社1993杨岞生俞守勤飞行器部件空气动力学航空工业出版社1987andersonjr
按速度范围分类:
低速空气动力学 (Low Aerodynamics) 亚音速空气动力学 (Subsonic Aerodynamics) 超音速空气动力学 (supersonic Aerodynamics) 高超音速空气动力学 (hypersonic Aerodynamics)
其它
36
37
38
39
21
0.3 空气动力学的发展进程
18世纪是流体力学的创建阶段。伯努利(Bernoulli) 在1738年发表“流体动力学”一书中,建立了不可压流体 的压强、高度和速度之间的关系,即伯努利公式;欧拉 (Euler)在1755年建立了理想不可压流体运动的基本方程 组,奠定了连续介质力学的基础。达朗贝尔 D'Alembert 提出著名的达朗贝尔原理:“达朗贝尔疑题”就是他在 1744年提出的。拉格朗日(Lagrange)改善了欧拉、达朗 贝尔方法,并发展了流体动力学的解析方法。关于研究气 流对物体的作用力,最早是牛顿(Newton)于1726年提出 关于流体对斜板的作用力公式,他实际上是在撞击理论的 基础上提出来的,没有考虑到流体的流动性.
流体力学II教材讲解
流体力学II(Viscous Fluid and Gas Dynamics)讲义第一章、粘性不可压缩流体运动基本方程组(学时数:6)1-1.绪论流体力学是力学的一个重要分支,主要研究流体介质(液体、气体、等离子体)的特性、状态,在各种力的作用下发生的对流、扩散、旋涡、波动现象和质量、动量、能量传输,以及同化学、生物等其他运动形式之间的相互作用。
它既是一门经典学科,又是一门现代学科,对自然科学和工程技术具有先导作用。
历史上,力学包括流体力学,曾经经历基于直观实践经验的古代力学、基于严密数学理论的经典力学、基于物理洞察能力的近代力学三个阶段。
在人类早期的生产活动过程中,力学即与数学、天文学一起发展。
17世纪,Newton基于前人的天文观测和力学实验,发明了微积分,并总结出机械运动三大定律和万有引力定律,发表了著名的《自然哲学的数学原理》一书。
由于原理是普适自然与工程领域的规律,从而使力学成为自然科学的先导。
从17世纪开始,人们逐步建立了流体力学的基本理论体系,从Pascal定律、Newton粘性定律、Pitot 管测速,到Euler方程和Bernoulli方程,标志着流体动力学正式成为力学的一个分支学科。
18世纪,人们着重发展无粘流体的位势理论。
到了19世纪,为了解决工程实际问题,开始注重粘性的影响,Navier-Stokes方程的建立为流体力学的进一步发展奠定了完整的理论基础,但该方程解的存在性与光滑性的证明至今仍是一大难题。
20世纪初,Prandtl凭借出色的物理洞察能力,提出边界层理论,从而开创了流体力学的近代发展阶段,使力学成为人类实现“飞天”梦想的重要理论先导。
60年代以来,由于超级计算机、先进测试技术的发展和应用,力学进一步凸显宏微观结合和学科交叉的特征,进入现代力学发展新阶段。
刚刚过去的2011年,人类遭遇了一系列极端事件:日本海底地震导致海啸和福岛核电站泄露事故;澳大利亚飓风;我国干旱洪水灾害等异常气候问题。
《空气动力学基础》第9章
2
1 Ma2
C py 2
2
1
1
1
4
1 Ma
2
Ma
C py
2
1
22:35
14
第九章 高超声速流动基础知识
§9-2 高超声速相仿律和马赫无关原理
•马赫高无超关声速绕流中的激波和膨胀波关系式
气流经过膨胀波后参数变化 Ma 1 若 p
tan
p
2
1
1
tan p
mz
2
M Ma
两个放射相似翼剖面
Cp c2
f
Ma
c
,
c
Cy
2
1 b
b
0
fl
fu
dx
Y
Ma
Cx
3
X
Ma
Cy c2
Y
Mac
,
c
mz c2
M
Mac
,
c
Cx c3
X
Mac
,
c
22:35
17
第九章 高超声速流动基础知识
§9-2 高超声速相仿律和马赫无关原理
•马马赫赫无无关关原理
5 激波层内高温和真实气体效应
强烈压缩导致温度剧增
P RT 不成立 cp,cV, 不为常数
T 2000K,O2 2O T 4000K,N2 2N T 9000K O O e
N N e
离解
电离
气动性能
偏离完全气体假设
真实气体效应 气动热
22:35
电磁环境
10
第九章 高超声速流动基础知识
绕翼型的空气动力系数表达式
(2)等腰三角形翼型
2c
原质论(第四章)
原质论武汉市汉阳区孟凡彬著第四章原质流体流场4. 原质流体流场定义:原质流体流场是指原质以流体状态存在和运动时所占据的空间。
释义:根据原质公设,原质流体存在的方式有两种:流体状态和演化成实物粒子的状态。
实物粒子的状态是原质流体吸收能量和信息后形成的自相封闭的流动系统,例如环形涡管等。
4.1 原质流体流场的边界条件定理:原质流体流场的外边界为整个宇宙,内边界是至密的基本粒子表面。
证明:根据原质公设和原质能够渗透任何物体的性质,原质流体流场的外边界只能是整个宇宙,内边界只能是至密的基本粒子表面。
释义:气体、液体等常见流体由分子(严格来说还有离子、质子、电子)构成,由于分子不能穿透许多种固体材料,并且总量有限,所以分子流体存在各种各样的容器和边界,液体还可以存在自由液面边界及异重流边界。
但是原质流体的边界不可能是任何固体材料的表面,因为原质是空间线度无限小的粒子,任何宏观物体或者原子都会被原质渗透,所以由微观粒子构成的任何物体均不能成为原质流体流场的边界。
推论1:原质流体流场具有唯一性,不存在多个流场之间相互流入或流出的源和汇。
证明:根据原质公设和原质的渗透性质,宇宙中的原质流体流场是一个不可分隔的整体,不可能被任何材料或力场分割。
推论2:原质流体只能做环流运动,原质流体的速度场是有旋场。
证明:因为原质流体流场是唯一的,原质流体不可压缩和膨胀,并且无源无汇,所以原质流体流体只能做环流流动。
释义:原质流体流场的性质在数学的场论中已有相关论述,这里不再详述。
值得提醒的是:场论在原质论之前只能作为一种数学抽象,在原质论中,场论所描述的对象是原质流体真实流场,排除人为设置的边界条件,欧几里德几何学和数学分析本身就是对原质及原质流体的直接描述。
推论3:原质流体没有自由表面,也不存在表面波。
证明:根据推论1,原质流体没有自由表面。
流体的表面波只能存在于液体的自由表面,所以原质流体不存在表面波。
4.2原质流体流场中的物体定理:原质流体流场中的物体不存在压力和浮力。
空气动力学第四章粘性流体动力学基础
v(x, y, z,t) (zx xz) xyx yyy zyz
w(x x, y y, z z,t) w(x, y, z,t) w x w y w z x y z
4.2、流体微团的运动形式与速度分解定理
以x方向速度分量为例,由泰勒级数展开,有
u(x x, y y, z z,t) u(x, y, z,t) u x u y u z
x y z 将上式分别加、减下列两项
1 v y , 1 w z
得到
2 x
2 x
u(x x, y y, z z,t)
1 0 0
0
2
0
0 0 3
I1 1 2 3 I2 1 2 23 13 I3 1 23
4.3、粘性流体的应力状态
1、理想流体和粘性流体作用面受力差别 流体处于静止状态,只能承受压力,几乎不能承受拉力和剪力,不具有 抵抗剪切变形的能力。理想流体在运动状态下,流体质点之间可以存在 相对运动,但不具有抵抗剪切变形的能力。因此,作用于流体内部任意 面上的力只有正向力,无切向力。 粘性流体在运动状态下,流体质点之间可以存在相对运动,流体具有 抵抗剪切变形的能力。因此,作用于流体内部任意面上力既有正向力, 也有切向力。
D ( ps cos )ds 0 2R
4.1、流体的粘性及其对流动的影响
对于粘性流体的绕流,与理想流体绕流存在很大的差别。由于流体 与固壁表面的粘附作用,在物面近区将产生边界层,受流体粘性的 阻滞作用,流体质点在由A点到B点的流程中,将消耗部分动能用之 克服摩擦阻力做功,以至使其无法满足由B点到D点压力升高的要求 ,导致流体质点在BD流程内,流经一段距离就会将全部动能消耗殆 尽(一部分转化为压能,一部分克服摩擦阻力做功),于是在壁面 某点速度变为零(S点),以后流来的流体质点将从这里离开物面进 入主流场中,这一点称为分离点。这种现象称为边界层分离。在分 离点之间的空腔内流体质点发生倒流,由下游高压区流向低压区, 从而在圆柱后面形成了旋涡区。这个旋涡涡区的出现,使得圆柱壁 面压强分布发生了变化,前后不对称(如前驻点的压强要明显大于 后驻点的压强),因此出现了阻力D。
空气动力学部分知识讲解
空气动力学及飞行原理课程空气动力学部分知识要点一、流体属性与静动力学基础1、流体与固体在力学特性上最本质的区别在于:二者承受剪应力和产生剪切变形能力上的不同。
2、静止流体在剪应力作用下(不论所加剪切应力τ多么小,只要不等于零)将产生持续不断的变形运动(流动),换句话说,静止流体不能承受剪切应力,将这种特性称为流体的易流性。
3、流体受压时其体积发生改变的性质称为流体的压缩性,而抵抗压缩变形的能力和特性称为弹性。
4、当马赫数小于0.3时,气体的压缩性影响可以忽略不计。
5、流层间阻碍流体相对错动(变形)趋势的能力称为流体的粘性,相对错动流层间的一对摩擦力即粘性剪切力。
6、流体的剪切变形是指流体质点之间出现相对运动(例如流体层间的相对运动)流体的粘性是指流体抵抗剪切变形或质点之间的相对运动的能力。
流体的粘性力是抵抗流体质点之间相对运动(例如流体层间的相对运动)的剪应力或摩擦力。
在静止状态下流体不能承受剪力;但是在运动状态下,流体可以承受剪力,剪切力大小与流体变形速度梯度有关,而且与流体种类有关7、按照作用力的性质和作用方式,可分为彻体力和表面力(面力)两类。
例如重力,惯性力和磁流体具有的电磁力等都属于彻体力,彻体力也称为体积力或质量力。
8、表面力:相邻流体或物体作用于所研究流体团块外表面,大小与流体团块表面积成正比的接触力。
由于按面积分布,故用接触应力表示,并可将其分解为法向应力和切向应力:9、理想和静止流体中的法向应力称为压强,其指向沿着表面的内法线方向,压强的量纲是[力]/[长度]210、标准大气规定在海平面上,大气温度为15℃或T0=288.15K ,压强p0 = 760 毫米汞柱= 101325牛/米2,密度ρ0 =1.225千克/米311、从基准面到11 km 的高空称为对流层,在对流层内大气密度和温度随高度有明显变化,温度随高度增加而下降,高度每增加1km,温度下降6.5 K。
从11 km 到21km 的高空大气温度基本不变,称为同温层或平流层,在同温层内温度保持为216.5 K。
《流体力学》课件
流体力学是在人类同自然界作斗争和在生产实践中逐步发展起来的。
古时中国有大禹治水疏通江河的传说;秦朝李冰父子带领劳动人民修建的都江堰,至今还在发挥着作用;大约与此同时,古罗马人建成了大规模的供水管道系统等等。
流体力学的萌芽:距今约2200年前,希腊学者阿基米德写的“论浮体”一文,他对静止时的液体力学性质作了第一次科学总结。
建立了包括物理浮力定律和浮体稳定性在内的液体平衡理论,奠定了流体静力学的基础。
此后千余年间,流体力学没有重大发展。
15世纪,意大利达·芬奇的著作才谈到水波、管流、水力机械、鸟的飞翔原理等问题;17世纪,帕斯卡阐明了静止流体中压力的概念。
但流体力学尤其是流体动力学作为一门严密的科学,却是随着经典力学建立了速度、加速度,力、流场等概念,以及质量、动量、能量三个守恒定律的奠定之后才逐步形成的。
流体力学的主要发展:17世纪,力学奠基人牛顿(英)在名著《自然哲学的数学原理》(1687年)中讨论了在流体中运动的物体所受到的阻力,得到阻力与流体密度、物体迎流截面积以及运动速度的平方成正比的关系。
他针对粘性流体运动时的内摩擦力也提出了牛顿粘性定律。
使流体力学开始成为力学中的一个独立分支。
但是,牛顿还没有建立起流体动力学的理论基础,他提出的许多力学模型和结论同实际情形还有较大的差别。
之后,皮托(法)发明了测量流速的皮托管;达朗贝尔(法)对运动中船只的阻力进行了许多实验工作,证实了阻力同物体运动速度之间的平方关系;瑞士的欧拉采用了连续介质的概念,把静力学中压力的概念推广到运动流体中,建立了欧拉方程,正确地用微分方程组描述了无粘流体的运动;伯努利(瑞士)从经典力学的能量守恒出发,研究供水管道中水的流动,精心地安排了实验并加以分析,得到了流体定常运动下的流速、压力、管道高程之间的关系——伯努利方程。
欧拉方程和伯努利方程的建立,是流体动力学作为一个分支学科建立的标志,从此开始了用微分方程和实验测量进行流体运动定量研究的阶段。
可压缩准一维管道无粘流动_欧拉方程_MacCormack
可压缩准一维无粘管道流动摘 要 本题利用一维欧拉方程求解可压缩一维无粘管道流动,并针对出口不同的条件,出口亚音速和出口超音速两种不同条件,将流道进行网格划分,利用MacCormack 进行差分求解,得到管道内的总压、马赫数、总焓、内能的分布,并给出计算过程中残差收敛的过程。
关键词 准一维;欧拉方程;MacCormack1 问题提出如Figure 1所示,流动为流经变截面流道内、理想、不可压、定常、平面的流动。
设进出口截面速度均匀分布。
来流条件为 1.5M ∞=,31.2218/kg m ρ∞=,47892.4p Pa ∞=。
计算出口超音速和亚音速(119/out u m s =)时,这个准一维流动的压强和马赫数分布,并给出残差收敛过程。
Figure 1准一维管道示意图本题的分析思路:首先,建立数学模型,之后运用边界条件定解。
接下来将模型问题转化为数值求解问题,通过网格划分、差分离散,并运用matlab 编程进行求解。
2 模型建立选定控制体,根据积分型方程导出该一维流动的欧拉方程。
控制体如Figure 1中红框所示,得到连续性、动量和能量方程如(1)、(2)和(3)()()0A dx uA dx t x ρρ∂∂+=∂∂⎰⎰(1) 2()()uA dx u A dx pAdx t x ρρ∂∂+=-∇∂∂⎰⎰⎰ (2) ()()()d EA dx EuA dx puA dx t x dx ρρ∂∂+=-∂∂⎰⎰⎰(3)控制体以上方程对于任意控制体均成立,因此可得到如下微分方程,即一维可压缩无粘流动的欧拉方程紧凑形式(4) t x∂∂+=∂∂Q FS (4)式中,20,,0dA A u A u p p dx E Hu ρρρρρρ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥==+=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦Q F S (5)其中 1.3980.347tanh(0.8 4.0)A x =+-。
对于此处的欧拉方程,有三个方程组,但是未知量分别为,,,,p u H E ρ共5个,因此为了方程组的封闭需要补充两个方程。
空气动力学前六章知识要点
空气动力学基础前六章总结第一章空气动力学一些引述1、空气动力学涉及到的物理量的定义及相应的单位①压强:是作用在单位面积上的正压力,该力是由于气体分子在单位时间内对面发生冲击(或穿过该面)而发生的动量变化,具有点属性。
p = lim i dF ,dAr 0IdA丿单位:Pa, kPa, MPa 一个标准大气压:101kPa②密度:定义为单位体积内的质量,具有点属性。
— lim dm,dv > 0dv单位:kg/tf 空气密度:1.225Kg/卅③温度:反应平均分子动能,在高速空气动力学中有重要作用。
单位:C④流速:当一个非常小的流体微元通过空间某任意一点的速度。
单位:m/s⑤剪切应力:.二卩:黏性系数dy⑥动压:4比="2空;2、空气动力及力矩的定义、来源及计算方法空气动力及力矩的来源只有两个:①物体表面的压力分布②物体表面的剪应力分布。
气动力的描述有两种坐标系:风轴系(L,D )和体轴系(A,N)。
力矩与所选的点有关系,抬头为正,低头为负。
L = N cos : - Asin :,D = Ns i n A c o s3、气动力系数的定义及其作用气动力系数是比空气动力及力矩更基本且反映本质的无量纲系数,在三维中I的力系数与二维中有差别,^口:升力系数C L=丄(3D),q =丄(2D)q(S q^cc_L Q _D Q _N Q _A CC L, C D , C N, C A, C MqSq ::Sq :Sq ::S二维:S=C (1)=C4、 压力中心的定义压力中心,作用翼剖面上的空气动力,可简化为作用于弦上某参考点的升力 L,阻力D 或法向力N ,轴向力A 及绕该点的力矩M 。
如果绕参考点的力矩为零, 则该点称为压力中心,显然压力中心就是总空气动力的作用点,气动力矩为005、 什么是量纲分析,为什么要进行量纲分析,其理论依据,具体方法在等式中,等号左边和等号右边各项的的量纲应相同, 某些物理变量可以用 一些基本量(质量,长度,时间等)来表达,据此有了量纲分析法,量纲分析可 以减少方程独立变量个数,其理论依据是白金汉n 定理。
第四章 理想流体运动基础
第四章 理想流体运动基础
流体微元沿n方向的力平衡式为
dn p dn p p dsdx p dsdx g n dsdndx an dsdndx 2 n 2 n p g n an n 式中,an是流体微元沿n方向的向心加速度,指向流线的曲 率中心,gn则是力加速度矢量在n方向的分量。对于定常 流动为 an V 2 R
或
19
式中,fr、fθ、fz分别为单位质量力在r、θ、z坐标 轴方向的分量。
7
第四章 理想流体运动基础
§4-2 自然坐标系中的欧拉方程
以流线方向为参考定义一个正交坐标系,称为流线坐 标系,或自然坐标系。 如图所示,在流线上取一点P,过P点作一个局部的正 交坐标系,其三个互相垂直的坐标方向分别为沿流线方向s、 垂直于流线的主法线方向n和副法线方向b,三个方向的单 e e e 位矢量分别表示为 、 和 s n b
12
第四章 理想流体运动基础
欧拉方程
§4-3 伯努利方程
或
fx
u u u 1 p u x ux x u y x uz x x t x y z u y u y u y 1 p u y fy ux uy uz y t x y z 1 p u z u z u z u z fz ux uy uz z t x y z
(4-1)
分别在无旋流动和有旋流动情况下求解上式 1. 欧拉积分 在无旋流动时, 0 ,式子变为
u u2 G pF 0 t 2
(4-2)
从数学分析可知,无旋的条件 0 是uxdx+uydy+uzdz成为 某一函数ψ(x, y, z, t)的全微分的必要充分条件。函数ψ( x, y, z, t)称为速度势函数,简称速度势。当以t为参变量时 ,函数ψ(x, y, z, t)的全微分可写成
第四章流体的积分关系式及其应用
第四章 流体的积分关系式及其应用众所周知,一个固体质点在保守力场中运动时,质点的动能和势能之和保持不变,这就是经典物理中的机械能守恒定律。
从数学的观点看,机械能守恒是动量方程的一次积分,称为能量积分。
有了能量积分方程,我们在处理保守场中的动力学问题时,就可通过该方程将始、末两态直接联系起来,而不必考虑中间过程的细节。
在流体力学中也有类似的积分。
前面一章建立了控制流体流动的微分方程组,原则上利用该方程组可以求解满足Stokes 假设的Newton 流体的任意流动问题。
对于理想流体流动问题,可以直接积分微分方程,得到积分方程。
利用积分方程求解流动问题显然更为简便,因而这些积分方程得到广泛应用。
什么样的流体是理想流体呢?当流体发生剪切变形时总会伴有粘性应力。
粘性应力不仅与流体的粘性性质(以粘性系数表征)有关,还依赖于速度梯度,对于低粘(μ小)流体的流动,如果速度的空间变化不太急剧,粘性应力就比较小。
如果粘性应力对所研究的流动问题影响较小,可以忽略流体的粘性,认为流体是无粘的,即理想流体。
一般常见的流体,如空气和水,粘性系数很小,在自然界和工程中遇到的这些流体的大多数流动,粘性的影响都可以忽略,都可以近似看作理想流体流动。
在流体力学发展的历史上,无粘流理论是流体力学中历史悠久,发展完善,成果辉煌,应用广泛的一个分支领域。
§4.1理想流体运动方程的进一步化简理想流体满足Euler 方程:dV pF dt ρ∇=-, (4-1)或者改写成兰姆—葛罗米柯形式2rot 2VV p V V F t ρ⎛⎫∂∇+∇+⨯=- ⎪∂⎝⎭。
(4-2)若体力有势F π=-∇,(4-3)其中π代表体力势,即单位质量流体的势能。
如果体力仅为重力,取z 轴沿g -方向,并取0z =为零势能面,则gz π=。
若流体密度是常量或仅为压强的函数,则称流体是正压流体。
若流体正压,)(p ρρ=,此时可定义压力函数()dpP p ρ=⎰(4-4-1)或()dpdP p ρ=。
4-第四讲 无粘不可压有势流动
第四讲 无粘不可压有势流动一、流函数与势函数1、势函数⎪⎪⎩⎪⎪⎨⎧∂∂=∂∂=y v x u ϕϕ 02222=∂∂+∂∂=∂∂+∂∂y x y v x u ϕϕ为调和函数,满足拉普拉斯方程。
2、流函数⎪⎪⎩⎪⎪⎨⎧∂∂-=∂∂=x v y u ψψ 022=∂∂∂-∂∂∂=∂∂+∂∂x y y x y v x u ψψ 若流动无旋,则有 02222=∂∂-∂∂-=∂∂-∂∂=yx y u x v z ψψω为调和函数,满足拉普拉斯方程。
3、科希-黎曼关系式⎪⎪⎩⎪⎪⎨⎧∂∂-=∂∂∂∂=∂∂x yy x ψϕψϕ4、等流函数线与等势线正交 对于等势线和等流函数线,有2100k yx dxdy dy y dx x d k yx dxdy dy y dx x d =∂∂∂∂-=⇒=∂∂+∂∂==∂∂∂∂-=⇒=∂∂+∂∂=ψψψψψϕϕϕϕϕ121-=⋅⋅-=∂∂⋅∂∂∂∂⋅∂∂=⋅v u v u y yx xk k ψϕ可构成正交网络。
5、等流函数线为流线vdy u dx udy vdx dy ydx x d =⇒=+-=∂∂+∂∂=0ψψψ满足流线方程。
二、基本流动及其合成1、平行来流(1)速度分布 ⎪⎪⎩⎪⎪⎨⎧∂∂-=∂∂==∂∂=∂∂==∞x y v y x V u ψϕψϕ0(2)速度势函数x V ∞=ϕ(3)流函数 y V ∞=ψ 2、点源与点汇(1)速度分布由连续方程,有 ⎪⎪⎩⎪⎪⎨⎧∂∂-=∂∂==∂∂=∂∂==⇒=r r V r r r q V q V r r r ψθϕθψϕππθ022(2)势函数 r qln 2πϕ= (3)流函数 θπψ2q = 3、点涡(1)速度环量 θπV r l d V C2=⋅=Γ⎰==〉⎪⎩⎪⎨⎧Γ==r V V r πθ20(2)势函数 θπϕ2Γ=(3)流函数 r ln 2πψΓ=4、偶极流点源和点汇配置如图所示,点源位于A 点(-a,0),点汇位于B 点(a,0),迭加后的速度势函数为B B A A r Qr Q ln 2ln 2ππϕ-=2222)()(a x y PB r a x y PA r B A -+==++==若B A Q Q =,则有2222)()(ln4ln 2)ln (ln 2a x y a x y Q r r Q r r Q B A B A -+++==-=πππϕ 流函数为P B A Q Q θπθθπψ2)(2-=-=(1) 流线方程常数=ψ 得到流线方程为 常数=P θ流线时经过点A 和点B 的圆线簇。
叶轮机械气动热力学-第 4 章
Xi’an Jiaotong University
Institute of Turbomachinery
透平中的三维无粘影响
环形面积的径向改变以及流管高度的径向改变引起径向流动; 动叶沿叶高厚度的变化引起流线的径向移动; 轴功的径向变化也引起径向二次流; 两相流引起径向密度的不连续(水轮机或水泵中气液两相流);
涡量表达式:
2011-9-28
8/52
Xi’an Jiaotong University
Institute of Turbomachinery
采用简单径向平衡方程的流型设计:
稳态、圆柱面流动(径向速度为零)时,控制方程大大简化。径向动量 方程变为简单径向平衡方程,该方程在叶片排间隙中有效。
物理解释:圆周速度引起指向轴心的向心加速度
4.1 轴对称流动分析
轴对称流动分析:
求解叶片前、后沿径向的流动参数, 再由叶栅理论预测通道中流动; 主要用于轴流叶轮机械;
2011-9-28
7/52
Xi’an Jiaotong University
Institute of Turbomachinery
轴对称理论中,假设叶片上、下游流动为轴对称; 控制方程:(无粘、稳态、圆柱坐标系)
设计实例B:h0 h0 r , s const , 给定 r 设计者拥有更大的选择范围。
径向平衡方程简化为: 将 代入,可得:
解出得: 给出中径rm处Vm,利用关系:
其中:
得:
2011-9-28
13/52
Xi’an Jiaotong University
Institute of Turbomachinery
V2 r
空气动力学部分知识要点
空气动力学及飞行原理课程空气动力学部分知识要点一、流体属性与静动力学基础1、流体与固体在力学特性上最本质的区别在于:二者承受剪应力和产生剪切变形能力上的不同。
2、静止流体在剪应力作用下(不论所加剪切应力T多么小,只要不等于零)将产生持续不断的变形运动(流动),换句话说,静止流体不能承受剪切应力,将这种特性称为流体的易流性。
3、流体受压时其体积发生改变的性质称为流体的压缩性,而抵抗压缩变形的能力和特性称为弹性。
4、当马赫数小于0.3 时,气体的压缩性影响可以忽略不计。
5、流层间阻碍流体相对错动(变形)趋势的能力称为流体的粘性,相对错动流层间的一对摩擦力即粘性剪切力。
6、流体的剪切变形是指流体质点之间出现相对运动(例如流体层间的相对运动)流体的粘性是指流体抵抗剪切变形或质点之间的相对运动的能力。
流体的粘性力是抵抗流体质点之间相对运动(例如流体层间的相对运动)的剪应力或摩擦力。
在静止状态下流体不能承受剪力;但是在运动状态下,流体可以承受剪力,剪切力大小与流体变形速度梯度有关,而且与流体种类有7、按照作用力的性质和作用方式,可分为彻体力和表面力(面力)两类。
例如重力,惯性力和磁流体具有的电磁力等都属于彻体力,彻体力也称为体积力或质量力。
8、表面力:相邻流体或物体作用于所研究流体团块外表面,大小与流体团块表面积成正比的接触力。
由于按面积分布,故用接触应力表示,并可将其分解为法向应力和切向应力:9、理想和静止流体中的法向应力称为压强,其指向沿着表面的内法线方向,压强的量纲是[力]/[长度]210、标准大气规定在海平面上,大气温度为15 C 或T o =288.15K,压强p o = 760毫米汞柱二101325牛/米2,密度p二1.225 千克/米311 、从基准面到11 km 的高空称为对流层,在对流层内大气密度和温度随高度有明显变化,温度随高度增加而下降,高度每增加1km,温度下降6.5 K。
从11 km到21km的高空大气温度基本不变,称为同温层或平流层,在同温层内温度保持为216.5 K。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
空气动力学
2 马赫数M 马赫数M
第四章
高速可压无粘流
[2] 马赫数的物理意义 M数还代表单位质量气体动能与内能之比。 数还代表单位质量气体动能与内能之比。 动能 v 2 / 2 γ (γ − 1) 2 = = M
内能 cv T 2
1. M数很小:动能相对于内能很小,速度的变 数很小:动能相对于内能很小, 化不会引起温度的显著变化, 化不会引起温度的显著变化 , 对于不可压 流体,密度、温度均可以认为是常值。 流体,密度、温度均可以认为是常值。 2. M数较高:动能相对于内能较大,速度的变 数较高:动能相对于内能较大, 化将引起温度的显著变化, 密度、 化将引起温度的显著变化 , 密度 、 温度均 是变数。 是变数。
回目录页
4.1.1 完全气体的状态方程、内能和焓 四 章 完全气体的状态方程、 第
空气动力学
高速可压无粘流
任何气体的状态方程
p = p( ρ , T )
完全气体状态方程 p = ρRT
从微观角度看, 从微观角度看,完全气体具有以下性质: ①分子体积与分子间平均距离相比可忽略不计; 分子体积与分子间平均距离相比可忽略不计; ②除弹性碰撞外,分子间没有其它相互作用。 除弹性碰撞外,分子间没有其它相互作用。
§4.2 一维等熵绝热流
4.2.1一维等熵绝热流的能量方程 4.2.1一维等熵绝热流的能量方程 4.2.2一维定常绝热流的参数间的基本关系式 4.2.2一维定常绝热流的参数间的基本关系式
回目录页
4.2.1一维等熵绝热流的能量方程 4.2.1一维等熵绝热流的能量方程 一维欧拉方程的伯努利积分: 一维欧拉方程的伯努利积分:
第四章 高速可压无粘流
静参数:流动过程中任意一点处的当地流动参数。 静参数:流动过程中任意一点处的当地流动参数。
T0 γ −1 2 = 1+ M T 2
p0 γ − 1 2 γ −1 = 1 + M p 2
γ
ρ0 γ − 1 2 = 1 + M 2 ρ
熵是一个状态参数。 熵是一个状态参数。 状态参数
∆s = s 2 − s1 = ∫ ds = ∫
1
2
2
1
T2 ρ1 dq = cV ln + R ln T T1 ρ2
下一页 回目录页
4.1.3 热力学第二定律和熵 利用状态方程: 利用状态方程:
空气动力学
第四章 高速可压无粘流
p2 ρ1 γ T2 ρ1 γ −1 ∆s = cV ln[ ( ) ] 或 ∆s = cV ln[ ( ) ] p1 ρ 2 T1 ρ 2
下一页 回目录页
4.1.1 完全气体的状态方程、内能和焓 四 章 完全气体的状态方程、 第
空气动力学
高速可压无粘流
对于完全气体,其内能只与绝对温度有关。 对于完全气体,其内能只与绝对温度有关。
e = e(T )
焓:代表热含量,定义为 代表热含量, p h = e+ ρ ① 焓表示单位质量气体的内能和压力能的和。 焓表示单位质量气体的内能和压力能的和。 ② 完全气体的焓只取决于温度,是一个状态参数。 完全气体的焓只取决于温度,
空气动力学
2 马赫数M 马赫数M
第四章
高速可压无粘流
[2] 马赫数的物理意义 高速空气动力学中, 高速空气动力学中,M数是一个反映压缩性 大小的相似准则。 大小的相似准则。
a2 = dp ∆p ρv v ~ ~ ~ dρ ∆ρ ∆ρ (∆ρ / ρ )
2 2
v 2 ∆ρ M2 = 2 ~ ρ a
显然, 越大则引起的压缩性越大。 M<0.3, 显然,M越大则引起的压缩性越大。当M<0.3, |△ρ/ ρ |<5%,认为不可压。 |<5%,认为不可压。
v2 dp +∫ = const 2 ρ
(沿流线) 沿流线)
空气动力学
第四章 高速可压无粘流
利用等熵关系: 利用等熵关系: v2 γ p + = const 2 γ −1 ρ 利用焓的定义: 利用焓的定义:
(沿流线) 流线)
v + h = const 2
微团宏观运动的动能
2
(沿流线) 沿流线) 微观的热力运动所含有的能量
空气动力学
1 音速
第四章
高速可压无粘流
a = dp / dρ
2
音速的计算公式表明小扰动的传播速度与 介质(气体、液体、固体)的压缩性有关。 介质(气体、液体、固体)的压缩性有关。 小扰动在可压气体中的传播过程非常接近 绝热可逆过程(即等熵过程) 绝热可逆过程(即等熵过程)。
dp a = dρ = γRT s =c
空气动力学
第四章 高速可压无粘流
高速流动问题中,气体的粘性内摩擦、 高速流动问题中 , 气体的粘性内摩擦 、 激波的出 现以及因温度梯度存在而引起的热传导等因素, 现以及因温度梯度存在而引起的热传导等因素 , 导致流动过程不可逆。 导致流动过程不可逆。 一般说来,在流场的大部分区域, 一般说来 , 在流场的大部分区域 , 速度梯度和温 度梯度不很大,流动过程可近似认为绝热可逆的, 度梯度不很大 , 流动过程可近似认为绝热可逆的 , 熵增△ ,称为等熵流。 熵增△s=0,称为等熵流。 在流体力学中, 在流体力学中 , 把熵值沿一条流线保持不变叫沿 流线等熵;全流场熵值相同称为均( 熵流。 流线等熵;全流场熵值相同称为均(匀)熵流。 等熵关系式
2
空气动力学
2 马赫数M 马赫数M
第四章
高速可压无粘流
[1] 马赫数的定义 流场中任一点处流速与当地音速之比定义 为马赫数M 为马赫数M: v
M = a
一般说, 一般说,流场中各点的流速和音速各不相 指的是当地值,称当地马赫数。 同,故M指的是当地值,称当地马赫数。 来流马赫数M 来流马赫数M∞= v∞/ a∞:来流速度与来流温 度所对应的音速之比。 度所对应的音速之比。
空气动力学
1 音速
第四章
高速可压无粘流
考虑一个固定扰源在静止空气中不断发生微弱扰 动的情况。采用和球面波固连在一起的动坐标: 动的情况。采用和球面波固连在一起的动坐标:
ρ 连续方程:adA = (ρ + dρ )(a − dv )dA 连续方程:
adρ = ρdv
a 2 = dp / dρ
dp 动量方程: = ρadv 动量方程:
取 T = 0 时,
e = h = 0 则有
e = cV T
h = c pT
下一页 回目录页
4.1.2 热力学第一定律和比热
空气动力学
第四章 高速可压无粘流
比热比: 比热比:等压比热和等容比热之比
γ = c p cV
由焓的定义可得
c p = cV + R
焓用比热比可写为
h=
γ γ −1
RT =
γ
p
回目录页
γ −1 ρ
4.1.3 热力学第二定律和熵
空气动力学
第四章 高速可压无粘流
热力学第二定律:指明能量相互转化是有条件、 热力学第二定律:指明能量相互转化是有条件、 有方向性的,即一个方向的变化过程可以实现, 有方向性的,即一个方向的变化过程可以实现, 而逆方向的变化过程或者不能实现或者只能有条 件地实现。 件地实现。 可逆过程、不可逆过程: 可逆过程、不可逆过程:如果将变化过程一步步 地倒回去, 地倒回去,介质的一切热力学参数均回到初始值 且外界情况也都复旧,则过程可逆, 且外界情况也都复旧,则过程可逆,否则就是不 可逆过程。 可逆过程。
结论:当沿流线速度增大时,温度、音速和焓值 结论:当沿流线速度增大时,温度、 均下降,但动能和焓的总和保持不变。 均下降,但动能和焓的总和保持不变。
回目录页
4.2.2一维定常绝热流参数间的关系 4.2.2一维定常绝热流参数间的关系 [一]驻点参数
空气动力学
第四章 高速可压无粘流
驻点是指流动速度等熵地降为零的一点。 驻点是指流动速度等熵地降为零的一点。该点 等熵地降为零的一点 的状态参数叫驻点参数,或称总参数。 的状态参数叫驻点参数,或称总参数。空气动学第四章 高速可压无粘流
§ 4.1 热力学基础 § 4.2 一维等熵绝热流 § 4.3 马赫波与膨胀波 § 4.4 正激波 § 4.5 斜激波 § 4.6 喷管及超音速风洞
§4.1热力学基础
4.1.1完全气体的状态方程、内能和焓 4.1.1完全气体的状态方程、内能和焓 4.1.2热力学第一定律和等容等压比热 4.1.2热力学第一定律和等容等压比热 4.1.3热力学第二定律和熵 4.1.3热力学第二定律和熵
v2 γ p γ p0 + = 2 γ − 1 ρ γ − 1 ρ0 v2 + h = h0 2
γ γ v2 + RT = RT0 2 γ −1 γ −1
2 v2 a2 a0 + = 2 γ −1 γ −1
下一页
回目录页
4.2.2一维定常绝热流参数间的关系 4.2.2一维定常绝热流参数间的关系
空气动力学
热力学第二定律指出, 热力学第二定律指出 , 在绝热变化过程的孤立系统 如果过程可逆熵值保持不变, 中,如果过程可逆熵值保持不变,∆s = 0 ,称为等 熵过程;如果过程不可逆,熵值必增加, 熵过程;如果过程不可逆,熵值必增加,∆s > 0 。 因此热力学第二定律亦可称为熵增原理。
回目录页
4.1.3 热力学第二定律和熵
下一页
回目录页
4.2.1一维等熵绝热流的能量方程 4.2.1一维等熵绝热流的能量方程
空气动力学
第四章 高速可压无粘流
γ v + RT = const 2 γ −1
2
(沿流线) 沿流线)
v2 a2 + = const 2 γ −1 −1