[论文解读] The Structure of Radiative Shock Waves. I. The method of global iterations

Yu. A. Fadeyev, Denis Gillet|arXiv (Cornell University)|Feb 6, 1998

Gas Dynamics and Kinetic Theory参考文献 4被引用 4

一句话总结

本文提出了一种全局迭代方法，用于自洽求解部分电离氢气中稳态平面平行辐射冲击波的流体动力学、辐射转移和原子动力学耦合方程。该方法收敛至一个稳定解，能够准确捕捉冲击流与辐射场之间的强反馈作用，对上游速度为15至60 km/s、马赫数在2.3至9.3之间的冲击波表现出鲁棒性。

ABSTRACT

The structure of steady plane-parallel radiative shock waves propagating through the hydrogen gas undergoing partial ionization and excitation of bound atomic states is investigated in terms of the self-consistent solution of the equations of fluid dynamics, radiation transfer and atomic kinetics. The shock wave model is represented by a flat finite slab with no incoming radiation from external sources at both its boundaries. The self-consistent solution is obtained using the global iteration procedure each step of which involves (1) integration of the fluid dynamics and rate equations for the preshock and postshock regions, consecutively, both solutions being fitted by the Rankine-Hugoniot relations at the discontinuous jump; (2) solution of the radiation transfer equation for the whole slab. The global iteration procedure is shown to converge to the stable solution which allows for the strong coupling of the gas flow and the radiation field produced by this flow. Application of the method is demonstrated for the shock waves with upstream velocities from 15 km/s to 60 km/s (that is with upstream Mach numbers from 2.3 to 9.3) and the hydrogen gas of unperturbed temperature T=3000K and density rho = 1e-10 gm/cm^3.

研究动机与目标

• 开发一种自洽的数值方法，用于建模流体动力学、辐射转移和原子动力学完全耦合的辐射冲击波。
• 解决部分电离氢气中冲击波产生的辐射场与气体结构之间强反馈的挑战。
• 实现对上游速度范围（15–60 km/s）和马赫数范围（2.3–9.3）内冲击波结构的精确计算。
• 验证全局迭代过程在模拟复杂辐射冲击系统时的收敛性和稳定性。
• 为未来在H II区和超新星遗迹等天体物理环境中的辐射冲击波研究提供基础。

提出的方法

• 该方法采用全局迭代过程，交替求解预激波区和激波后区的流体动力学方程与速率方程，通过Rankine-Hugoniot跃变条件实现解的匹配。
• 在每次迭代步骤中，对整个有限厚度平板区域求解辐射转移方程，考虑自洽的发射与吸收过程。
• 在激波前区和后区依次求解流体与原子动力学方程，确保在激波间断面上的连续性。
• 将前一次迭代中计算得到的辐射场作为下一次迭代的输入，逐步优化辐射与气体结构之间的耦合。
• 采用有限厚度平板几何结构，边界无外部辐射，模拟孤立的冲击系统。
• 当迭代之间辐射场与气体结构的变化低于指定容差时，认为达到收敛，确保获得稳定自洽的解。

实验结果

研究问题

• RQ1 如何在辐射冲击波中实现流体动力学、辐射转移和原子动力学的自洽耦合？
• RQ2 何种数值方法可确保平面平行辐射冲击波中完整方程组的稳定与收敛解？
• RQ3 辐射场如何影响马赫数在2.3至9.3之间的冲击波的结构与热力学特性？
• RQ4 氢的局部电离与激发在塑造冲击波结构与辐射输出中起到何种作用？
• RQ5 该全局迭代方法能否在无外部辐射源的条件下，可靠地模拟上游速度范围为15–60 km/s的冲击波？

主要发现

• 全局迭代方法成功收敛至部分电离氢气中辐射冲击波的稳定自洽解。
• 该方法准确捕捉了冲击波辐射场与气体结构之间的强反馈作用，尤其在高马赫数冲击波中表现显著。
• 对于上游速度15–60 km/s、马赫数2.3–9.3的条件，该方法生成了具有正确能量重分布的物理解释性冲击波剖面。
• 辐射场主要由氢的谱线发射和连续过程主导，尤其在激波后区，其中束缚-束缚与束缚-自由跃迁贡献显著。
• 即使在强耦合条件下，解仍保持稳定与收敛，证明了全局迭代方法的鲁棒性。
• 计算得到的冲击波结构与非辐射模型相比存在显著差异，特别是在激波后冷却层，这是由于高效的辐射能量损失所致。