瑞利流(加热管道流)模拟器 返回
可压缩流体

瑞利流(加热管道流)模拟器

在等截面、无摩擦管道中加热时气体如何变化?调节入口马赫数、温度、压力与加热量,可实时看到出口马赫数与离热壅塞的余量,帮助建立对喷气发动机燃烧室与冲压发动机的直觉。

参数设置
入口马赫数 M₁
进入加热段的气体马赫数(亚音速)
入口温度 T₁
K
入口压力 P₁
kPa
比热比 γ
空气 γ=1.40,高温燃气 γ≈1.33
加热量 q
kJ/kg
单位质量的加热量;约等于燃料发热量 × 燃料/空气比
计算结果
入口滞止温度 T₀₁ (K)
加热后滞止温度 T₀₂ (K)
加热极限 q_max (kJ/kg)
出口马赫数 M₂
加热率 q/q_max (%)
热壅塞判定
加热管道流动 — 加速动画

亚音速气体进入等截面管道,加热使其密度下降并加速。颜色表示温度(蓝→橙→红),箭头长度表示当地速度(马赫数)。M 到达 1 时显示热壅塞标记。

瑞利函数 T₀/T₀* vs 马赫数
出口马赫数 vs 加热量 q
理论与主要公式

$$\frac{T_{0}}{T_{0}^{*}}=\frac{(\gamma+1)M^{2}\left(2+(\gamma-1)M^{2}\right)}{\left(1+\gamma M^{2}\right)^{2}}$$

瑞利函数:滞止温度 T₀ 与热壅塞滞止温度 T₀* 之比,是 M 的函数,在 M=1 处等于 1。亚音速下加热使 T₀ 上升并把 M 推向 1。

$$T_{0}=T\left(1+\frac{\gamma-1}{2}M^{2}\right), \qquad q_{\max}=c_{p}\,(T_{0}^{*}-T_{0,1})$$

滞止温度定义与最大加热量。c_p = γR/(γ−1),R = 287 J/(kg·K)。当 q 达到 q_max 时出口 M=1(热壅塞)。

什么是瑞利流(加热管道流)

🙋
"瑞利流"这个名字第一次听到。它和普通的管道流有什么区别?
🎓
好问题。瑞利流是一个理想化模型:理想气体在等截面管道里流动,沿途加入(或带走)热量,而忽略摩擦。这正是喷气发动机燃烧室、冲压发动机加热段的"教科书模型"。真实燃烧室里燃烧释热远比摩擦剧烈,所以用瑞利流就能抓住直觉。
🙋
原来如此。但是加热气体应该只是"温度上升"吧?可在模拟器里把加热量 q 调大,出口马赫数 M₂ 也跟着上升,让我很意外。
🎓
这正是瑞利流的核心。亚音速流加热后气体膨胀,管截面又不变,所以由质量守恒 ρuA=常数,密度 ρ 下降必然导致速度 u 上升。马赫数 M=u/√(γRT) 里温度上升也会让音速上升,但在亚音速区"速度上升"压倒"音速上升",结果 M 朝 1 攀升。"加热使流动加速"是可压缩流力学最反直觉、最标志性的结果。
🙋
那只要不停加热,流动就会一直加速吗?
🎓
不行——存在一个硬性天花板叫做"热壅塞(thermal choking)"。当加入特定热量时,出口马赫数刚好到 1,再想加热,物理上同一根管子已经塞不下了。本工具用 q_max 表示这个极限。默认条件(M₁=0.3、空气)下 q_max≈578 kJ/kg。超过 q_max 会怎样?上游压力会自发上升、流量下降,或出现激波,使工作点回到 q ≤ q_max 的范围内。实机燃烧室通常把 q 限制在 q_max 的 60~70% 以下。
🙋
出口马赫数 M₂ 能用解析公式算吗?瑞利方程看起来是按 T₀ 整理的…
🎓
观察得好。瑞利函数 T₀/T₀* = (γ+1)M²(2+(γ−1)M²)/(1+γM²)² 的分子、分母都是 M 的四次多项式,无法把 M 显式解出来。所以本工具先算加热后的 T₀₂/T₀* 作为目标值,再在 [M₁, 1] 区间用二分法迭代 60 次求 M₂。亚音速区瑞利函数单调上升,二分法可靠收敛。经典气体动力学教材里要查 Rayleigh 表,现在几行代码就能瞬时得到答案。
🙋
这种流动在实际机型里用在哪里?
🎓
最典型的就是涡喷、涡扇发动机的主燃烧室,冲压发动机、超燃冲压发动机、加力燃烧室,以及固体火箭尾喷管前的加热段。它们都符合"等截面、瞬间大量加热"的设定。地面工业炉里长管道加电热丝控制流量的场合,也用瑞利流来估算"再加热就会被流量限制——即达到热壅塞"。在可压缩流体力学教材里,瑞利流是最贴近真实硬件的一维理想模型之一。

常见问题

瑞利流是一种理想化模型:理想气体在等截面管道内流动,过程中加入或带走热量,但忽略摩擦。它最适合用于摩擦次要、热释放占主导的喷气发动机燃烧室等场合。范诺流则是它的孪生模型——绝热(无热量交换)但只有摩擦,常用于分析长管中的流动损失。两者都是一维等截面的可压缩流,区别在于驱动能量项是热还是摩擦。
加热使气体膨胀,密度ρ下降;由于管截面积 A 不变,根据质量守恒 ρ·u·A = 常数,速度 u 必须上升。马赫数 M = u/√(γRT) 中,温度上升固然会使音速上升,但在亚音速区,速度上升的贡献超过音速上升,结果马赫数会朝 1 攀升。这正是瑞利流最反直觉的标志性结果。
持续对亚音速瑞利流加热,出口马赫数最终会精确等于 1,再多的热量都无法通过同一根管道。这就是热壅塞。本工具用 q_max(最大加热量)表示。避免方法有三:(1) 降低上游质量流量;(2) 扩大管截面积;(3) 接受上游压力上升使运行点改变。实机燃烧室通常限制喷油量,并加宽燃烧室以降低入口马赫数。
瑞利函数 T0/T0* = (γ+1)M²(2+(γ−1)M²)/(1+γM²)² 在亚音速区是 M 的单调函数,且在 M=1 处等于 1。本工具先计算加热后的滞止温度 T02,然后在区间 [M1, 1] 上对方程 T02/T0* = rayleighFunc(M2) 用 60 次二分法求 M2。当加热量超过 q_max 时 target>1,将 M2 钳位为 1 并显示热壅塞警告。

实际应用

喷气发动机燃烧室:涡喷、涡扇发动机的主燃烧室就是瑞利流的直接对象。压气机将空气增压,喷油燃烧,向气流释放数百到上千 kJ/kg 的热量。设计时把燃烧室做得足够宽,使入口马赫数维持在 0.2~0.3,让 q 远低于 q_max,从而避免压力损失与出口温度失控。在本工具里把 M₁ 调到 0.5 以上、再增大 q,就能直观感受 q_max 变得多么严格。

冲压发动机与超燃冲压发动机:这两类发动机在超音速进气后,于燃烧室加热以产生推力,燃烧室内的流速管理关系生死。冲压发动机把燃烧室设为亚音速(进气道减速),严格遵循瑞利流;超燃冲压发动机则保持燃烧室超音速以避免热壅塞,利用超音速瑞利流"加热反而减速"且 M 仍趋向 1 的镜像行为。

加力燃烧室(Afterburner):战机瞬间增加 50% 推力的"蓝色火焰",就是在涡轮排气后再次喷油加热的瑞利流装置。它常工作在接近热壅塞的极限,燃油控制把 q 维持在 q_max 的 80~90% 之间,需要严密的闭环反馈以避免热壅塞与熄火。

CFD 解析的前期评估:本工具这类一维瑞利模型,可在投入详细 CFD(化学反应+湍流+三维形状)解析之前做"量级估算"。若 CFD 给出的出口马赫数或压比与瑞利估算差出一个数量级,应怀疑边界条件或燃烧模型的设定问题。本科和研究生的可压缩流体力学课程也总会把瑞利流与 h-s 图中的瑞利线一起讲解。

常见误解与注意事项

最大的误解是"加热气体只会让温度升高"。在封闭容器内的确如此,但在开口的管道流中,加热会使气体膨胀→密度下降→速度上升,亚音速时马赫数趋向 1。超音速瑞利流则恰好相反:加热使流动减速,但 M 仍朝 1 移动。"加热总把 M 推向 1"才是统一两侧的正确说法。"加热让流动变快"只适用于亚音速分支。

其次是"加热量超过 q_max 会怎样"这一疑问。物理上不存在稳态解:上游压力会自动上升(背压传到上游)、流量下降,或者出现激波,使工作点回到满足 q ≤ q_max 的状态。本工具在 q>q_max 时把 M₂ 钳位到 1 并显示热壅塞警告,但在实机中这种工况会表现为流量振荡(pumping)与燃烧振荡,是极其危险的。设计准则就是 q ≤ 0.6~0.7 × q_max。

最后请勿忘记瑞利流忽略了摩擦。实际较长的燃烧室与加热管道中,摩擦与加热同时起作用,实际出口马赫数会略高于纯瑞利预测(亚音速下摩擦同样把 M 推向 1)。短而高热负荷的燃烧室用瑞利流即可,长而低热负荷的管道则需要 Fanno+Rayleigh 综合分析或更一般的一维方法。本工具应作为"加热主导的短段"理想模型来使用。