左状态（高压侧）
右状态（低压侧）
预设
激波参数
—
激波马赫数
—
后方压力 P*
—
激波速度
—
接触面速度
精确计算一维黎曼解。实时可视化激波、接触间断和稀疏扇的压力、密度和速度分布，支持x-t特征线图。
整个模拟的核心是求解一个被称为“黎曼问题”的初值问题。它由描述气体运动的一维欧拉方程组控制，包含了质量、动量和能量守恒：
$$ \frac{\partial \mathbf{U}}{\partial t}+ \frac{\partial \mathbf{F}(\mathbf{U})}{\partial x}= 0 $$其中，$\mathbf{U}= [\rho, \rho u, E]^T$ 是守恒变量向量，分别代表密度、动量和总能量。$\mathbf{F}$ 是对应的通量。给定初始时刻左右两侧不同的均匀状态，求解这个方程组就得到了激波管的精确解。
求解的关键步骤是确定接触间断两侧的均匀压力 $p^*$ 和速度 $u^*$。这需要通过迭代求解一个由两侧黎曼不变量连接起来的压力函数 $f(p)=0$：
$$ u_L + \frac{2a_L}{\gamma-1}\left[ \left( \frac{p}{p_L}\right)^{(\gamma-1)/(2\gamma)}- 1 \right] = u_R - \frac{2a_R}{\gamma-1}\left[ \left( \frac{p}{p_R}\right)^{(\gamma-1)/(2\gamma)} - 1 \right] $$这里，$p_L, p_R$ 和 $u_L, u_R$ 是初始左右状态的压力和速度，$a$ 是声速，$\gamma$ 是比热比。解出 $p^*$ 后，所有波（激波、接触间断、稀疏扇）的速度和中间状态都能解析求出。
计算流体力学（CFD）软件验证：激波管问题是CFD求解器的“标尺”或“基准测试”。任何新开发的可压缩流求解器，都必须先能准确模拟这个简单的一维问题，其结果需要与这里的精确解高度吻合，工程师才会信任它去模拟更复杂的飞机超音速飞行或发动机燃烧。
航空航天工程：在超音速进气道、火箭发动机喷管的设计中，内部流动会产生复杂的激波系。激波管原理帮助工程师理解激波与边界层相互作用、激波反射等基础现象，为设计提供理论依据。
爆炸与冲击波研究：激波管本身就是一个产生可控冲击波的实验装置。它被用于研究材料在冲击载荷下的性能、爆炸冲击波的传播规律，以及建筑物的抗爆设计。
汽车安全测试：在安全气囊点爆的瞬间，气体高速充入气囊，这个过程类似于一个微型激波管。理解其中的气体动力学，有助于优化气囊的展开速度和形状，从而更有效地保护乘员。
首先，人们常误以为“冲击波总是以超音速传播”，但这其实是相对于冲击波的速度。在模拟器中观测到的波传播速度，是由初始状态声速对应的马赫数决定的。例如，左右压差仅略微增大的“弱冲击波”，其马赫数会接近1（如1.05），几乎以声速传播。反之，在“强冲击波”预设中，马赫数可能超过3，这已接近爆炸现象。
其次，接触间断面的重要性容易被忽视。由于该界面上的压力和速度连续，乍看图线时很容易忽略。然而，这里正是密度、温度乃至熵发生间断的界面，承载着流体的“历史印记”。例如，若初始时左右气体不同（左：氦气，右：空气），该界面两侧的密度将显著差异，从而对后续流动和混合现象产生巨大影响。在实际CFD应用中，如何抑制该界面的数值扩散是计算精度的关键。
最后，需注意比热比γ的设定不可轻视。从默认的空气（γ=1.4）改为双原子气体、单原子气体或高温气体时，结果的变化可能超乎预期。例如，对于单原子气体氩气（γ=1.67），即使压力比相同，冲击波后的温升也会比空气更为显著。这是因为比热比反映了气体分子自由度，决定了内能分配比例。建议养成习惯：勿盲目沿用默认值，务必确认目标流体的物性参数。
本工具涉及的一维黎曼问题解法，正是超音速/高超音速气动设计的基础。例如，预测超音速飞行器前端产生的冲击波角度及其后方压力/温度分布时，便应用了此概念。尤其在超燃冲压发动机进气道内部，为通过控制复杂冲击波系实现高效空气压缩，以黎曼求解器为核心的CFD分析不可或缺。
此外，该方法也直接关联爆炸安全工程与炸药爆轰现象分析。TNT等炸药爆轰时，会产生极高压力的“爆轰波”——一种强冲击波。激波管理论作为理解其传播的基础模型，已被广泛应用和拓展。在模拟器中观察“强冲击波”，正是直观理解爆轰波导致超压骤增的第一步。
更出人意料的是，在天体物理学领域（特别是超新星爆发残骸的膨胀模拟中），该理论经多维化与相对论化改造后亦被使用。冲击波在星际介质中的传播机制、以及其加速宇宙线的原理阐释，都离不开描述冲击波前后突跃条件的根本物理规律。
建议下一步学习“二维黎曼问题”。即从一维管道拓展至平面，研究四种不同状态的气体区块相互作用的问题。由此可观察到冲击波反射、不同族冲击波相交的“激波-激波相互作用”，以及涡旋生成等更为丰富的现象。这正是实用CFD代码中复杂流动问题的原型。
若想深化数学背景，必须理解“特征线法”。本工具的后台正是运用该方法构建稀疏波区域的连续解。通过学习双曲型偏微分方程（欧拉方程）的解在特征线上保持定值（黎曼不变量）的特性，便能从数学结构上理解膨胀波出现的必然性。
最后，可探索“数值黎曼求解器”领域。这是对当前“精确解”进行快速近值的各类数值算法（如Roe求解器、HLL求解器）。实际CFD应用出于计算成本与鲁棒性平衡的考虑，并不使用精确解，而是以这些近似求解器为核心。了解它们近似了什么、简化了何处，是工程师批判性验证CFD计算结果的必备素养。