编辑推荐：
研究人员提出了一种基于逆向时间积分的激波管非平衡流动建模方法，以解决传统方法在处理非均匀激波速度及驱动气体波系传播时的局限性。该方法结合准一维欧拉方程与Mirels边界层源项模型，通过引入实验测量的激波轨迹与最终压力分布作为边界条件，实现了对测试段内流场时空演
研究人员提出了一种基于逆向时间积分的激波管非平衡流动建模方法，以解决传统方法在处理非均匀激波速度及驱动气体波系传播时的局限性。该方法结合准一维欧拉方程与Mirels边界层源项模型，通过引入实验测量的激波轨迹与最终压力分布作为边界条件，实现了对测试段内流场时空演化的高精度重构。研究采用两温度（平动-转动温度T与振动-电子温度TV）非平衡热力学模型，并耦合详细化学反应动力学机制，以模拟高焓流动中的能量传递与物种演化。数值求解采用有限差分格式与牛顿迭代法，通过特征线分解处理前向与后向传播波系，显著提升了计算效率。验证结果表明，该方法与二维轴对称Navier-Stokes模拟及现有准一维非平衡求解器NESS在氩气与空气测试案例中均表现出良好的一致性。应用于月球返回与土卫六进入条件下的激波管辐射实验时，该方法成功量化了激波速度不确定性对辐射峰值强度的贡献，并揭示了反应动力学参数在远离激波区域的支配作用。研究成果为激波管实验数据的精确解读与热化学模型优化提供了新范式。
研究背景与意义
激波管实验是高焓非平衡流动研究的核心手段，但其流场常受非均匀激波速度、边界层增长及驱动气体波系干扰的影响，导致传统基于恒定激波速度与平衡假设的模型难以准确重构测试段内的热力学状态。尤其在辐射预测中，流场参数的微小偏差会引发显著的辐射强度误差，制约了对热化学模型的验证精度。针对上述问题，研究人员开发了LAgrangian Shock Tube Analysis 2.0（LASTA 2.0）软件，通过融合逆向时间积分策略与两温度非平衡模型，首次实现了对激波管全流场时空演化的自洽模拟，为消除实验数据解读中的系统误差提供了新方法。该成果发表于流体力学领域权威期刊《Journal of Fluid Mechanics》。
关键技术方法
研究采用准一维欧拉方程框架，嵌入Mirels边界层质量引射模型以表征粘性效应。核心创新在于逆向时间积分技术：将实验获得的激波轨迹通过Rankine-Hugoniot关系转化为激波边界条件，同时将测试段末端压力分布作为后向特征线（w+）的约束条件，通过求解嵌套三对角Jacobian矩阵的牛顿方程组获得全流场解。非平衡热力学采用Park两温度模型，耦合OCEAN热力学库计算输运性质与反应速率。网格收敛性研究表明，200个空间点与60个时间步可满足计算精度需求。验证阶段分别对比了FROSST（二维轴对称Navier-Stokes求解器）与NESS（准一维非平衡求解器）的结果，并应用于NASA EAST与牛津T6激波管的发射光谱实验数据。
研究结果
2.2 边界层形状因子
基于Mirels分析，引入形状因子β的修正公式，其中CM项用于修正自由流物性变化对流场的影响。研究指出，尽管Mirels模型针对稳态激波发展，但在典型激波管条件下瞬时边界层重塑假设仍具合理性。
2.3 双温非平衡系统的特征线
定义w-（源自激波）与w+（源自驱动气体）两组特征线。传统前向积分因缺乏接触面压力边界条件而无法处理w+波系，本研究通过逆向时间积分，利用最终压力剖面反推历史状态，成功捕捉了w+波对测试段后部流场的显著影响。
2.4 数值格式
建立拉格朗日坐标系下的有限差分格式，采用特征变量离散控制方程。通过构造包含测试段长度η的Schur补矩阵，解决了流场全域耦合问题。稳定性分析表明，后向特征线离散格式的系数大于1，保证了数值稳定性。
2.5 边界条件
激波边界采用高阶多项式拟合实验轨迹，通过Rankine-Hugoniot关系计算初态；压力边界通过时间插值关联空间位置，并讨论了壁面测量与中心线差异的修正方法。
2.6 求解方法
采用牛顿迭代法求解非线性代数方程组，利用Jacobian矩阵的嵌套三对角结构实现高效求逆。通过分区求解策略（???B/???B与???T/???T）加速收敛，松弛因子ω<1确保迭代稳定性。
3.1 流体模型验证
在氩气激波管案例中，LASTA 2.0预测的8 m、8.5 m处压力分布与FROSST模拟偏差小于1%，温度分布在初始200 μs内偏差低于4%。粒子轨迹分析显示，边界层压缩导致流线在测试段后部聚集，证实了粘性效应的空间非均匀性。
3.2 双温模型验证
在4.72 km/s空气激波案例中，与NESS求解器对比显示振动-电子温度TV演化趋势一致，近激波区差异源于NESS包含的粘性耗散效应。
3.3 网格独立性
空间网格加密至500点时TV收敛于307 K以内；时间步增至60步时结果波动小于0.5%，验证了瞬时边界层重塑假设的适用性。
4.3 月球返回条件应用
9.4 km/s N2/O2激波实验中，结合激波轨迹与压力边界的模拟显著改善了辐射预测：在760–785 nm波段，Boltzmann态布居模型在平衡区（>18 mm）与实验吻合，而非平衡区偏差归因于Rankine-Hugoniot初态的高温近似（>40,000 K）及反应速率不确定性。激波速度±100 m/s扰动可导致平衡区辐射强度±18%的变化，凸显了精密测速的重要性。
4.4 土卫六进入条件应用
4.8 km/s N2/CH4激波中，新方法成功捕捉了CN辐射的非平衡演化。与Satchell方法相比，其避免了平衡假设导致的CN数密度高估，揭示了激波减速历史对产物分布的调控作用。
讨论与结论
研究通过Bowman准则与Hooker边界层厚度评估，证实所有案例均处于边界层主导区，准一维假设有效。不确定性量化表明：在近激波区（<15 mm），辐射预测误差主要来自激波速度测量（可达±100 m/s）；而在远场区（>15 mm），反应动力学参数（如C+N2?CN+N与N2+M?2N+M）的不确定性成为主导因素。因此，高精度激波速度测量与两温度非平衡模型的协同优化是提升激波管数据解读可靠性的关键。该方法消除了传统模型中非物理的压力约束，为热化学模型验证提供了迄今最完备的流场重构工具。
