激波管模拟器

激波前后满足质量、动量、能量守恒。这些条件汇总成的兰金-雨贡尼奥关系式是激波管问题的核心。

$$ \frac{p_2}{p_1}= 1 + \frac{2\gamma}{\gamma+1}(M_s^2 - 1) $$ $$ \frac{\rho_2}{\rho_1}= \frac{(\gamma+1)M_s^2}{(\gamma-1)M_s^2 + 2}$$

其中，$p$ 是压力，$\rho$ 是密度，$\gamma$ 是比热比，$M_s$ 是激波马赫数。下标 1、2 分别表示激波前（低压侧）和后（高压侧）。压力比与马赫数的平方成正比。

激波管整体解由激波、接触不连续、稀疏波满足的条件决定，通过求解区域间的压力 $p^*$ 和速度 $u^*$ 的非线性方程得到。

$$ \frac{p_4}{p_1}= \frac{p_2}{p_1}\left[ 1 - \frac{(\gamma_4-1)(\frac{a_1}{a_4})(\frac{p_2}{p_1}-1)}{\sqrt{2\gamma_1}\sqrt{2\gamma_1+(\gamma_1+1)(\frac{p_2}{p_1}-1)}}\right]^{-\frac{2\gamma_4}{\gamma_4-1}} $$

这是激波关系式与稀疏波的等熵关系通过接触不连续面的压力、速度连续条件联系起来的方程。$a$ 是音速，下标 1、4 分别表示初始状态的右侧（低压）和左侧（高压）区域。数值求解该方程可确定所有波的强度和位置。

CFD 求解器的验证与基准测试：激波管问题（索德问题）是验证压缩性流体数值计算代码精度的最基本、最重要的测试用例。开发新的 CFD 求解器时，首先要与该精确解比较，检查激波的锐度和稀疏波的分辨率。

超音速与极超音速空气动力学研究：激波管虽然作用时间短，但能产生高焓气流。该特性用于研究再入舱或极超音速导弹头部产生的强激波及伴随的高温气体流。

爆轰与爆炸现象建模：炸药爆轰和气体爆炸产生极强激波（爆轰波）。激波管理论是预测爆轰波传播速度、诱发流动及压力分布的物理模型基础。

航空航天发动机进气口设计：超音速飞行时，发动机进气口内形成复杂激波系统。激波与边界层干涉不稳定化流动。激波管研究成果被用于进气口形状优化和冲击控制技术开发。

首先，很容易误认为"激波总是超音速传播"。这指的是相对激波的相对速度。模拟器中观测的波传播速度由初始态音速相对的马赫数决定。例如，压力差很小的"弱激波"马赫数仅 1.05 左右，几乎以音速传播；"强激波"预设的马赫数可超过 3，接近爆炸现象。

其次，容易忽视接触不连续面的重要性。该面压力、速度连续，乍看图表容易漏掉。但密度、温度、熵在此处不连续，刻录了流体的"历史"。例如左侧为氦气、右侧为空气时，该面前后密度差异巨大，对后续流动和混合现象影响深远。实际 CFD 中，如何抑制该面的数值扩散是计算精度的关键。

最后，比热比 γ 的设置不可轻视。从默认空气值（γ=1.4）改为其他气体（二原子气体、单原子气体、高温气体）时，结果变化远超预期。例如单原子气体氩（γ=1.67）在相同压力比下，激波后温升比空气明显。这是因为比热比反映分子自由度，决定内能分配。养成确认目标流体物性参数的习惯，而非盲目套用默认值。