可压缩流动计算器
返回
Compressible Flow

可压缩流动与马赫数关系计算器

实时计算等熵流动关系式 · 正激波 · 斜激波。调节M和γ,在P/P₀、T/T₀、ρ/ρ₀、A/A*曲线上标记当前工作点。

参数设置
马赫数 M₁
M=1:声速   M>1:超声速
比热比 γ
空气:1.40 · 氦气:1.67 · CO₂:1.30
预设

等熵流动关系:

$$\frac{T_0}{T}=1+\frac{\gamma-1}{2}M^2$$ $$\frac{P_0}{P}=\left(1+\frac{\gamma-1}{2}M^2\right)^{\!\gamma/(\gamma-1)}$$ $$\frac{A}{A^*}=\frac{1}{M}\left[\frac{2}{\gamma+1}\left(1+\frac{\gamma-1}{2}M^2\right)\right]^{\!(\gamma+1)/[2(\gamma-1)]}$$

正激波(Rankine-Hugoniot):

$$M_2=\sqrt{\frac{M_1^2(\gamma-1)+2}{2\gamma M_1^2-(\gamma-1)}}$$ $$\frac{P_2}{P_1}=\frac{2\gamma M_1^2-(\gamma-1)}{\gamma+1}$$
计算结果
P/P₀
T/T₀
ρ/ρ₀
A/A*
等熵流动关系式 vs 马赫数

什么是可压缩流动与马赫数

🙋
老师,总听人说“马赫数”,它到底是个啥?和速度有啥区别?
🎓
简单来说,马赫数(M)就是物体速度与当地声速的比值。它衡量的是空气的“可压缩性”。比如飞机速度是340米/秒,当地声速也是340米/秒,那马赫数就是1。你可以在模拟器里把马赫数M₁从0.3拖到3试试,你会发现当M<1时是亚音速,M>1时就是超音速,流动性质会发生剧变!
🙋
诶,真的吗?那为什么超音速飞机前面会有一道“白雾墙”(激波)?这和马赫数有什么关系?
🎓
问得好!那堵“墙”就是激波。当物体以超音速飞行时,前面的空气来不及“躲开”,被剧烈压缩和加热,形成一道很薄的波面。在实际工程中,比如导弹或航天飞机再入大气层时,头部就会产生强烈的弓形激波。你试着在模拟器里把马赫数调到2以上,然后看看“正激波”相关的参数,比如激波后的压力P₂会突然跳升好几倍,这就是激波压缩的效果。
🙋
原来如此!那旁边还有个叫“比热比γ”的参数是干嘛的?空气不都是一样的吗?
🎓
γ(伽马)是气体的一个关键属性,等于定压比热容除以定容比热容。它决定了气体被压缩或膨胀时,温度和压力变化的“剧烈程度”。比如,对于常温空气,γ≈1.4;但对于火箭发动机喷出的高温燃气,或者氦气这样的单原子气体,γ就不同。你可以在模拟器里把γ从1.4改成1.67(模拟氦气),再看看同样的马赫数下,温度比T/T₀和压力比P/P₀的变化曲线,会发现曲线形状都变了!这对设计不同工质的发动机至关重要。

物理模型与关键公式

等熵流动关系(核心):描述气流在绝热、无摩擦(熵不变)条件下,马赫数如何影响温度、压力和密度。这是分析喷管、风洞的基础。

$$\frac{T_0}{T}=1+\frac{\gamma-1}{2}M^2$$ $$\frac{P_0}{P}=\left(1+\frac{\gamma-1}{2}M^2\right)^{\!\gamma/(\gamma-1)}$$ $$\frac{\rho_0}{\rho}=\left(1+\frac{\gamma-1}{2}M^2\right)^{\!1/(\gamma-1)}$$

其中,$T_0, P_0, \rho_0$ 是气流速度降为零时的“总温”、“总压”、“总密度”(也叫驻点参数)。$T, P, \rho$ 是当地的实际参数。$\gamma$ 是比热比,$M$ 是马赫数。

正激波关系(Rankine-Hugoniot方程):描述气流垂直通过一道无限薄的激波时,参数发生的突跃变化。这是一个不可逆过程,总压会下降(熵增)。

$$M_2 = \sqrt{\frac{M_1^2(\gamma-1)+2}{2\gamma M_1^2 - (\gamma-1)}}$$ $$\frac{P_2}{P_1}= \frac{2\gamma M_1^2 - (\gamma-1)}{\gamma+1}$$

其中,下标 $1$ 和 $2$ 分别代表激波前和激波后的状态。$M_1$ 必须大于1(超声速),而 $M_2$ 必然小于1(亚声速)。

现实世界中的应用

航空航天发动机喷管设计:火箭和喷气式发动机的尾喷管是收缩-扩张形的(拉瓦尔喷管)。工程师利用等熵流动中的 $A/A^*$(面积比)公式,精确设计喉部($M=1$)和扩张段,确保气流能加速到设计马赫数(如 $M>5$),产生最大推力。

超音速飞行器与再入舱热防护:飞船或导弹以高超音速再入大气层时,前方会产生极强的弓形激波。激波后的气体温度极高(可达数千度)。通过正激波关系计算压力与温度跃升,是设计防热大底材料厚度和冷却系统的关键输入。

风洞试验与CFD仿真设置:在风洞中,通过测量静压和总压,利用等熵公式反推出实验段的马赫数。在CFD软件(如Fluent、OpenFOAM)中,这些公式用于设置超声速入口、出口和远场的边界条件,确保计算收敛且符合物理。

压气机与涡轮叶片流动分析:在航空发动机的压气机和涡轮中,气流速度很高,属于可压缩流动。叶片表面的局部区域可能达到跨音速甚至超音速,形成复杂的激波系。分析这些激波的位置和强度,对于提高效率、防止失速和喘振至关重要。

常见误解与注意事项

开始使用此工具时,有几个需要特别注意的要点。首先是“比热比γ并非固定值”。默认值1.4对应常温空气,但在例如喷气发动机内部的高温环境下γ会变小,而氩气则为1.67。若此值设置错误,特别是激波后的温度和压力计算结果会产生显著偏差。在实际工程中,务必先确认所处理流体及其温度范围再设定γ值,这是基本原则。

其次,请牢记“等熵流动终究是一种理想状态”。实际流道必然存在摩擦和热传导,因此计算得出的完美温降比(T/T₀)不会完全实现。例如在拉瓦尔喷管设计中,通常会在基于理想计算得到的形状上,额外增加考虑摩擦损失的“安全余量”。

最后,要正确理解斜激波“最大偏转角”的含义。在工具中将马赫数设为3、γ=1.4并逐渐增加偏转角θ时,大约在34度附近计算会失效。这表示“在此条件下物体无法弯曲得更尖锐”的极限。但实际工程中可能需要更大的弯曲角度吧?此时激波会从物体表面“分离”,产生更复杂的干涉波系。工具报错本身正是重要物理现象的信号。