第七章_气体的一维流动
第7章 一元气体动力学基础

第7章一元气体动力学基础本章目录§7.1 理想气体一元恒定流动的运动方程§7.2 音速、滞止参数、马赫数§7.3 气体一元恒定流动的连续性方程§7.4 等温管路中的流动§7.5 绝热管路中的流动本章概述气体动力学研究可压缩流体运动规律及其工程应用。
气体的密度随着压强和温度变化而变化,此时必须考虑气体的可压缩性。
气体动力学不仅研究流速、压强问题,而且包含密度和温度问题,不仅需要流体力学知识,还需要热力学知识。
进行气体动力学计算时,压强和温度只能用绝对压强和热力学温度。
理想气体状态方程:定容过程:热力学中,定容过程系指气体在容积不变或比容不变条件下进行的热力过程。
定温过程:热力学中,定温过程系指气体在温度不变条件下进行的热力过程。
绝热过程:热力学中,在无能量损失且与外界无能量交换的条件下进行的热力过程称为可逆的绝热过程,又称为等熵过程。
§7.1 理想气体一元恒定流动的运动方程§7.1.1 一元理想流体欧拉运动微分方程此即欧拉运动微分方程,也称为微分形式的伯努利方程。
§ 7.1.2 气体一元定容流动该方程的物理意义:沿流各断面上单位重量理想气体的压能和动能之和守恒,二者可以互相转换。
§7.1.3 气体一元定温流动定温流动也就是气体在温度保持不变情况下的流动。
§7.1.4 气体一元绝热流动绝热条件下的流动就是绝热流动,又称为等熵流动。
在绝热条件下,气体参数变化服从等熵过程方程理想气体绝热流动(等熵流动),沿流任意断面上,单位重量的气体所具有的内能、压能和动能之和为一常量。
§7.1.5 例题§7.1.6 关于气体一元绝热流动方程使用理想气体一元绝热流动方程,不仅适用于无摩阻的绝热流动中,也适用实际气流。
由于流动系统与外界无热量交换,摩擦产生的热量保存在管路中,所消耗的机械能转化为内能,其总和将保持不变。
工程流体力学课件-气体一维高速流动

由于气体一维流动中,气体参数 不随位置变化,因此流动是线性 的,可以应用一维流动方程进行 描述。
气体一维流动的分类
等熵流动
气体在流动过程中,熵值保持不变的 流动。等熵流动中,气体压力和密度 随速度增加而减小。
等温流动
气体在流动过程中,温度保持不变的 流动。等温流动中,气体压力和密度 随速度增加而增加。
火箭发动机喷管中的气体一维流动特性研究
总结词
火箭发动机喷管中的气体一维流动特性研究对于喷管 设计和火箭性能优化至关重要。
详细描述
火箭发动机喷管中的气体流动具有极高的速度和压力变 化,直接模拟三维流场非常困难且计算量大。因此,采 用一维流动模型进行研究和分析是常用的方法。一维流 动模型可以模拟喷管中气体的流动、加速和膨胀过程, 分析喷管的性能和特性。通过研究喷管中气体的流动特 性,可以优化喷管设计,提高火箭发动机的推力和效率 ,为火箭设计和发射提供重要的理论支持和技术保障。
动量守恒方程
表示动量在流动过程中的 变化,即动量在流场中不 增加也不减少。
能量守恒方程
表示能量在流动过程中的 变化,即能量在流场中不 增加也不减少。
初始条件和边界条件
初始条件
表示流动开始时流场中各物理量的值 。
边界条件
表示流场边界上各物理量的值或其变 化规律。
控制方程的离散化
有限差分法
将控制方程中的偏导数用差分近似代替 ,将连续的物理量离散为离散的数值。
有限差分法的优点是简单直观,易于编程实现,适用于各种类型的偏微分方程,特别是对波动问题和 稳定性问题有较好的处理能力。
有限元法
有限元法是一种将连续的物理量离散化为有限个单元,并在 每个单元上设置节点,通过节点上的等效源代替单元内的源 ,从而将偏微分方程离散化为线性方程组的方法。这种方法 在气体一维流动数值模拟中也有应用。
气体一维流动

7.4 正激波与斜激波
激波的定义
当超声速气流流过大的障碍物(或超声速飞机、炮弹和火箭等在空 中飞行)时,气流在障碍物前将受到急剧的压缩,它的压强、温度 和密度都将突跃地升高,而速度则突跃地降低。这种使流动参数发 生突跃变化的强压缩波叫做激波。
各种超声速飞行器飞行时 超声速气流绕过叶片、叶栅或其它物体流动时 超声速风洞启动时 缩放喷管在非设计工况运行时,在喷管的超声速流中也可能出 现激波。 原子弹、氢弹爆炸时产生的破坏力很大的高压强锋面是激波, 又称冲击波.
2
7.4 正激波与斜激波
完全气体蓝金-许贡纽公式
经过激波的密度突跃和温度突跃只决定于压强突跃
7.4 正激波与斜激波
7.4 正激波与斜激波
正激波前后气流参数的变化
连续方程: 动量方程: 能量方程:
7.4 正激波与斜激波
v 1v 2 cc2r M * 1M * 2 1
普朗特激波公式
M *1 1
7.4 正激波与斜激波
激波的分类
正激波:波面与气流方向相垂直的平面激波。
斜激波:波面与气流方向不垂直的平面激波。
曲激波:波形是弯曲的。
7.4 正激波与斜激波
正激波的形成 (0 t1)
7.4 正激波与斜激波
后面的微弱压缩波总比它前面的微弱压缩波传播得快
7.4 正激波与斜激波
( t1
t 2)
7.4 正激波与斜激波
斜激波的形成
超声速气流沿内折转一微小角度 dθ的内凹壁面流动时产生的马赫 波是微弱压缩波。气流越过微弱压缩波的流动为等熵的压缩过程。
7.4 正激波与斜激波
超声速气流经过每一个扰动点,都要产生一道微弱压缩波。气流越过 这一系列的微弱压缩渡时,其速度逐渐降低,而压强、密度和温度逐 渐升高,气流的马赫数逐渐减小,而马赫角逐渐增大。
第七章气体一维高速流动

Ma V c
(7-5)
Ma称为马赫数,是一个无量纲数,也是气体动力学中一
个重要参数。
我们常根据马赫数的大小,把气流分为亚声速流Ma<1,
跨声速流Ma ≈1,超声速流1< M<a3和高超声速流 M>a3等
几类。亚声速流动和超声速流动有许多显著的差别,我
们将在以后各节中逐一介绍。
1.静止流场(V=0) 在静止流场中,扰动源产生的微弱扰动波以声速c向四周
传播,形成以扰动源所在位置为中心的同心球面波,微弱 扰动波在3s末的传播情况如图7-2(a)所示。如果不考虑微弱 扰动波在传播过程中的损失,随着时间的延续,扰动必将 传遍整个流场。也就是说,微弱扰动波在静止气体中的传 播是无界的。
上述关系也适用于气流流过一静止微小障碍 物时的情况。假如气体以与上述扰动源的运动速 度数值相等而方向相反的速度作等速直线运动, 则扰动源就成为静止微小障碍物,即图7-2中的3 点就是静止扰动源,而扰动源所发出的扰动波 (图中的各圆)不断地被气流以速度-V带走。很 明显,在 V c(即Ma 1)的亚声速流动时,带 走的各扰动波在一定时间后可达到空间中的任何 一点。也就是说,扰动波不仅能顺流传播,而且 也能逆流传播。但在 V c(即 Ma 1)的超声速 流动时,带走的各扰动波只能在马赫锥内顺流传 播,不能逆流传播,也就是说在超声速流动中的 微弱扰动不能传播到整个空间。这就是超声速流 动和亚声速流动的一个重要差别,从而使这两种 流动的图形有着根本的不同。
第一节 微弱扰动波的传播
一. 微弱扰动波的一维传播
如图7-1所示,在一个截面积为A、足够长的直圆管中充满 了静止的气体,将圆管左端的活塞以微小速度 dV 向右轻微地 推动一下,使活塞右侧的气体压强升高一个微小增量 dp ,dp 所产生的微弱压强扰动向右传播。活塞将首先压缩紧贴活塞的 那一层气体,这层气体受压后,又传及下一层气体,这样依次 一层一层地传下去,就在圆管中形成一个不连续的微弱的压强 突跃,就是压缩波mn,它以速度 向右推进。压缩波面mn是 受活塞微小推移的影响而被扰动过的气体与未被扰动过的静止 气体的分界面。设在压缩波前未被扰动过的静止气体的压强
流体流动的一维流分析

流体流动的一维流分析引言流体流动的研究是流体力学中的重要内容之一。
流体流动可以分为一维流和多维流两种情况。
在一维流动中,流体在一个方向上的速度变化可以忽略不计,因此可以简化为一维问题进行分析。
本文将对流体流动的一维流进行分析,包括基本概念、数学模型、基本方程、解析方法及其应用等方面进行阐述。
一维流的基本概念一维流是指流体在一个方向上的速度变化可以忽略不计的流动。
在实际情况中,一维流动可以近似地描述一些特定的流体流动现象,如河流、管道流动等。
一维流动的速度场和压力场只与流动方向有关,与流动方向垂直的任意截面上的速度和压力分布均相等。
一维流动可以看作是通过管道或河道等局部的流动现象,对于整个系统来说,仍然是三维空间中的流动。
一维流的数学模型一维流动可以通过一维流动方程进行描述。
一维流动方程包括质量守恒方程和动量守恒方程两个基本方程。
质量守恒方程描述了单位时间内单位截面积内的流体质量守恒,动量守恒方程描述了单位时间内单位截面积内的动量守恒。
在一维流动中,流体的密度通常是恒定的,因此可以简化为不可压缩流动的数学模型。
质量守恒方程质量守恒方程可以写成以下形式:$$\\frac{\\partial \\rho}{\\partial t} + \\frac{\\partial (\\rho u)}{\\partial x} = 0$$其中,$\\rho$是流体密度,t是时间,x是坐标轴方向,u是流体速度。
动量守恒方程动量守恒方程可以写成以下形式:$$\\frac{\\partial (\\rho u)}{\\partial t} + \\frac{\\partial (\\rhou^2)}{\\partial x} = - \\frac{\\partial P}{\\partial x} + \\frac{\\partial}{\\partial x}\\left(\\mu\\frac{\\partial u}{\\partial x}\\right)$$其中,P是压力,$\\mu$是动力粘度。
第七章气体的一维流动

第七章 气体的一维流动
三、马赫数
马赫数:气体在某点的流速与当地声速之比定义为该点气流的马 赫数,用Ma表示。
Ma v / c
马赫数是零量纲速度 完全气体
Ma2 v2
RT
过程装备与控制工程教研室
16
第七章 气体的一维流动
Ma v / c
马赫数代表的是气体的宏观运动动能与气体内分子运动动能之比。 在气体流动的分析和计算中,将以马赫数作为判断气体压缩性的影响
波后气体以和活塞同样的微小速度dv运动。
过程装备与控制工程教研室
8
第七章 气体的一维流动
微弱压强波在圆管中的传播速度c
假定微弱压强波的波面已传到A-A,右侧尚未传 到,速度为零,压强为p,密度为ρ;
A-A的左侧是已受扰区,气体速度为dv,压强为 p+dp,密度为ρ+dρ;
对静止观测者,流动是非定常的; 如果取以波速c同步运动的坐标观测该流场,则流
动是定常的。
过程装备与控制工程教研室
9
第七章 气体的一维流动
取图中虚线部分为控制面 流体始终以速度c流向波面,压强和密度分别为p、ρ 流体又始终以c-dv的速度离开波面,其压强和密度分别为p+dp,
ρ+dρ 由连续方程 (ρ+dρ)(c-dv)A- ρcA=0
ρcA-ρdvA+ cdρA- dρdv A -ρcA=0 略去二阶微量 cdρ=ρdv
过程装备与控制工程教研室
7
第七章 气体的一维流动
倘若管内的活塞突然以微小的速度dv向左运动
首先使紧靠活塞右侧的一层气体膨胀,这层气体膨胀后,接着又 使下一层气体膨胀,一层一层地依次传下去,便在管内形成一道 以速度c向左传播的微弱膨胀波。
气体的一维流动共43页PPT

16、自己选择的路、跪着也要把它走 完。 17、一般情况下)不想三年以后的事, 只想现 在的事 。现在 有成就 ,以后 才能更 辉煌。
18、敢于向黑暗宣战的人,心里必须 充满光 明。 19、学习的关键--重复。
20、懦弱的人只会裹足不前,莽撞的 人只能 引为烧 身,只 有真正 勇敢的 人才能 所向披 靡。
Hale Waihona Puke 16、业余生活要有意义,不要越轨。——华盛顿 17、一个人即使已登上顶峰,也仍要自强不息。——罗素·贝克 18、最大的挑战和突破在于用人,而用人最大的突破在于信任人。——马云 19、自己活着,就是为了使别人过得更美好。——雷锋 20、要掌握书,莫被书掌握;要为生而读,莫为读而生。——布尔沃
END
气体的一维定常流动

速度系数
马赫数与速度系数的关系
2 2 M* 1 Ma 2 1 2 1 M* 1
M* < 1 M* = 1 M* > 1
1 1
亚声速流动 声速流动 超声速流动
M*
1
M*
2
1
1
2
2
Ma 2 Ma
2
1
Ma 1 §6-4 气体流动的三种状态和速度系数
0
§6-1 气体一维流动的基本概念
声速
取虚线所围控制体。 连续方程:
d c dvA cA
动量方程:
cd dv
cAc dv c p p dp A
p+dp ρ+dρ T+dT
cdv dp
p ρ T
§6-1 气体一维流动的基本概念
§6-2 微弱扰动在气体中的传播
气体静止不动
v0 Ma 0
扰动波是球形波,向所有方向传遍全部空间。
§6-2 微弱扰动在气体中的传播
气流亚声速流动
vc Ma 1
扰动波可以逆流传播,向所有方向传遍全部空 间。
§6-2 微弱扰动在气体中的传播
气流以声速流动
vc Ma 1
扰动波不能逆流传播,传播限制在下游半个空 间。
§6-2 微弱扰动在气体中的传播
气流超声速流动
扰动面为一系列 与圆锥面相切的 马赫锥 球面,该圆锥称 马赫角 为马赫锥(弱扰 动锥),圆锥面 称为马赫波。
ct c 1 sin vt v Ma
马赫波
vc Ma 1
扰动波不能逆流传播,传播限制在马赫锥空间 内。
§6-2 微弱扰动在气体中的传播
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
液体: 液体:
dp K=ρ →c= dρ
K
ρ
完全气体: 完全气体:忽略粘性和传热视作等熵过程
微分: 微分:
dp = γ
p
ρ
dp → c = γ
p
ρ
= γRT
第一节 微弱压强波的一维传播
讨论: 讨论: (1)声速是状态参数
(2)
c=
1 dρ dp
dρ / dp
越大,越易压缩,c越小
音速是反映流体压缩性大小的物理参数
γ
γ
当T=0时,v = vmax 则
vmax = 2 γRTT γ 1
第二节 气体特定状态和参考速度
(3)临界状态 当气流速度等于当地音速的状态,便是临界状态。以
vcr ,Tcr , pcr , ρ cr
分别表示临界速度、临界温度、临界压力和临界密度,则
2 γ 1 2γ 12 vcr = ccr = cT = vmax = (γRTcr ) = RTT γ +1 γ +1 γ +1
c1 ρ 2 vs = 1 2 γ ρ1
p2 1 p1 ρ2 1 ρ1
12
气流的温度突跃、 气流的温度突跃、密度突跃与压强突跃 一一对应; 一一对应;
12
p2 ρ 2 1 1 c1 p1 ρ1 vg = 1 2 ρ2 γ ρ1
第三节 正激波
二、激波的形成和厚度
由于活塞先后发出的压缩波并不以相同的绝对速度向前传播, 由于活塞先后发出的压缩波并不以相同的绝对速度向前传播,因为后面 的波是在前面的波已扰动的基础上发出的, 的波是在前面的波已扰动的基础上发出的,而后面的扰动波的速度比前 面波的速度要快, 面波的速度要快,故后面的波最终将追赶上前面的波而形成一道强的压 缩波即激波。 缩波即激波。
1 P1 ρ1 T1 v1 2 v2 P2 ρ2 T2 连续性方程 动量方程 能量方程
ρ1v1 = ρ 2 v2
波阻
ρ 2 v22 ρ1v12 = p1 p2
1 1 h1 + v12 = h2 + v22 2 2
2 γ + 1 ccr γ 1 v1 v2 = (v1 v2 ) 2γ v v + 2γ 1 2
M * max γ +1 v = max = ccr γ 1
12
第二节 气体特定状态和参考速度
速度系数与马赫数的关系图 3
M *max
2
T c2 γ 1 2 = 2 = 1 γ + 1 M* TT cT p γ 1 2 = 1 γ +1 M* pT
γ γ 1
1 2 v2 1 1 1 γ +1 2 γ +1 1 c + = ccr + = 2 2 2( γ 1 ) γ 1 γ 1 Ma 2 2( γ 1 ) M *2
可以得到 Ma 2 =
2M *2 (γ + 1) (γ 1)M *2
作用:绝能流中,ccr=常数, M c = v * cr 当v=vmax时,
激波与微弱压缩波有着本质的区别,其主要表现为: 激波与微弱压缩波有着本质的区别,其主要表现为: (1)激波是强压缩波 经过激波气流参数变化是突跃的; 激波是强压缩波, (1)激波是强压缩波,经过激波气流参数变化是突跃的; (2)气体经过激波受到突然地 强烈地压缩, 气体经过激波受到突然地, (2)气体经过激波受到突然地,强烈地压缩,必然在气体内 部造成强烈的摩擦和热传导,因此气流经过激波是绝能不等熵流动; 部造成强烈的摩擦和热传导,因此气流经过激波是绝能不等熵流动; (3)激波的强弱与气流受压缩的程度 或扰动的强弱)有直接关系。 激波的强弱与气流受压缩的程度( (3)激波的强弱与气流受压缩的程度(或扰动的强弱)有直接关系。
第三节 正激波
一、激波的分类
正激波:气流方向与波面垂直, ⒈ 正激波:气流方向与波面垂直,如图 (a); 缩放喷管的超声速段 ; 斜激波:气流方向与波面不垂直,如图(b); ⒉ 斜激波:气流方向与波面不垂直,如图 ; 超声速气流绕叶片 曲线激波:波形为曲线形,如当超声速气流流过钝头物体时, ⒊ 曲线激波:波形为曲线形,如当超声速气流流过钝头物体时,在物体 前面往往产生脱体激波,这种激波就是曲线激波,如图(c)所示 所示。 前面往往产生脱体激波,这种激波就是曲线激波,如图 所示。
(γ + 1)Ma1 ρ2 = ρ1 2 + (γ 1)Ma12
2
正激波前后气流参数只与Ma和 正激波前后气流参数只与 和 完全气体比热比有关系。 完全气体比热比有关系。
2γ γ 1 2 γ + 1 Ma1 γ + 1
γ γ 1
2 pT 2 (γ + 1)Ma1 = 2 pT 1 2 + (γ 1)Ma1
ρ dp 2 dv = vdv = γMa p p v
p dp dρ = C, = γ ρ ρ p dp dρ dT + p / ρ = RT , = p ρ T
γ
第三节 正激波
当超声速气流绕过大的障碍物时,气流 当超声速气流绕过大的障碍物时, 在障碍物前将受到急剧的压缩, 在障碍物前将受到急剧的压缩,它的压 温度、密度都将突跃地升高, 强、温度、密度都将突跃地升高,而 速度则突跃地降低, 速度则突跃地降低,这种流动参数发 生突跃变化的强压缩波叫激波 激波。 生突跃变化的强压缩波叫激波。激波 是超声速气流中经常出现的现象。 是超声速气流中经常出现的现象。
2γ γ 1 2 γ + 1 Ma1 γ + 1
1 γ 1
标志激波强度的压强几乎与波 前马赫数的平方成正比,说明来 前马赫数的平方成正比, 流马赫数的高低是激波强弱的主 要标志。 要标志。 激波强度越强,经过激波,机 激波强度越强,经过激波, 械能损失越大。 械能损失越大。
2 + (γ 1)Ma 2 γ 2γ γ 1 s 2 1 = ln 2 γ + 1 Ma1 γ + 1 Cv (γ + 1)Ma1
第三节 正激波
二、激波的形成和厚度
由于速度、温度等参数是连续变化的, 由于速度、温度等参数是连续变化的,实际的激波 是有厚度的。 是有厚度的。 Ma=2时,激波厚度为2.54×10-4mm,只有几个分 时 激波厚度为 × , 子平均自由行程。 子平均自由行程。
第三节 正激波
三、正激波的传播速度
12 12
12
2 Tcr = TT γ + 1 pcr 2 γ 1 = γ +1 pT
γ
ρ cr 2 = γ +1 ρT
1 γ 1
第二节 气体特定状态和参考速度
速度系数
M* =
气流速度与临界音速之比称为速度系数,用M* 表示,即
cT =
v ccr
γ +1
2
ccr
1 2 v2 1 2 c + = cT γ 1 2 γ 1
第三节 正激波
四、正激波前后气流参数的关系 波阻
气流经过激波,速度降低,动量减小,熵值增加, 气流经过激波,速度降低,动量减小,熵值增加,因而 必有作用在气流上与来流方向相反的力; 必有作用在气流上与来流方向相反的力; 引起激波的物体,必然受到与上述作用力大小相等、 引起激波的物体,必然受到与上述作用力大小相等、与 来流方向相同的反作用力,即气流作用在物体上的阻力, 来流方向相同的反作用力,即气流作用在物体上的阻力, 这种阻力因激波的存在而产生。 这种阻力因激波的存在而产生。 激波越强,波阻越大。 激波越强,波阻越大。
(3)
c = f (T ) = f ( p ,V ,T )
当地声速
(4)空气
c = 1.4 × 287T
T = 288K
c = 340.3m / s
第一节 微弱压强波的一维传播
马赫数
流场中任一点的流速与当地音速之比,称为该点处气流的马赫数,以符号Ma表示
v v Ma = ; Ma = c γRT
根据马赫数的大小把流动分为: Ma < 1:亚音速流 Ma = 1:音速流 Ma > 1
2
3
Ma
4
5
6
7
ρ γ 1 2 = 1 γ + 1 M* ρT
1 γ 1
当Ma = 0时,M * = 0 当Ma < 0时,M * < 0 当Ma = 1时,M * = 1 当Ma > 1时,M * > 1
例 题
视空气为 视空气为γ=1.4的完全气体,在一无摩擦的 的完全气体, 的完全气体 渐缩管道中流动,在位置1处的平均流速为 渐缩管道中流动,在位置 处的平均流速为 152.4m/s,气温为 ,气温为333.3K,气压为 , 2.086×105Pa,在管道出口 处达到临界状 × ,在管道出口2处达到临界状 试计算2处气流的温度 压强、 处气流的温度、 态。试计算 处气流的温度、压强、密度和 平均流速。 平均流速。
第四节 变截面管流
前提: 定比热容;完全气体;一维;定常;绝 前提: 定比热容;完全气体;一维;定常; 能等熵流。 能等熵流。 激波出现时,另当考虑。 激波出现时,另当考虑。
第四节 变截面管流
一、气流速度与通道截面的关系
dA dv dρ + + =0 A v ρ
动量方程
ρvdv = dp
c = γp / ρ
第七章 气体的一维流动
应用范围 气体流速较高,气体压缩性将明显地影响它 气体流速较高 气体压缩性将明显地影响它 许多热力学和动力学特性。 许多热力学和动力学特性。 各类气体输送管道、喷管、扩压管、 各类气体输送管道、喷管、扩压管、斜切 喷管、叶栅流道等。 喷管、叶栅流道等。
第一节 微弱压强波的一维传播
第二节 气体特定状态和参考速度