[论文解读] A radiation-hydrodynamics scheme valid from the transport to the diffusion limit

E. Audit, Pierre Charrier+2 more|arXiv (Cornell University)|Jun 17, 2002

Computational Fluid Dynamics and Aerodynamics参考文献 12被引用 21

一句话总结

本文提出了一种基于M¹闭合模型的新型辐射流体动力学格式，能够准确捕捉辐射传输（自由流）和扩散极限。通过渐近保持方式修改HLLE黎曼求解器，并与移动网格、时间隐式框架耦合，该方法在极端光学厚度梯度下仍保持鲁棒性和准确性，已在辐射冲击波和原恒星形成模拟中得到验证，能够高分辨率解析强烈的激波和扩散前峰。

ABSTRACT

We present in this paper the numerical treatment of the coupling between hydrodynamics and radiative transfer. The fluid is modeled by classical conservation laws (mass, momentum and energy) and the radiation by the grey moment $M_1$ system. The scheme introduced is able to compute accurate numerical solution over a broad class of regimes from the transport to the diffusive limits. We propose an asymptotic preserving modification of the HLLE scheme in order to treat correctly the diffusion limit. Several numerical results are presented, which show that this approach is robust and have the correct behavior in both the diffusive and free-streaming limits. In the last numerical example we test this approach on a complex physical case by considering the collapse of a gas cloud leading to a proto-stellar structure which, among other features, exhibits very steep opacity gradients.

研究动机与目标

• 开发一种单一数值格式，能够在无需区域分解的情况下，准确模拟从自由流到扩散的所有辐射流体动力学状态。
• 克服针对不同状态使用独立模型的局限性，避免在过渡区域引入误差。
• 确保在强光学厚度梯度问题中，正确表现流体和扩散极限的行为。
• 实现对复杂天体物理现象（如原恒星形成）的高分辨率模拟，准确处理辐射激波。
• 在移动网格、时间隐式框架中实现该格式，以应对天体物理模拟中时间与空间尺度的巨大差异。

提出的方法

• 流体通过质量、动量和能量的经典守恒定律建模。
• 辐射通过灰度M¹矩系统描述，该系统利用最小熵原理闭合矩方程，以捕捉各向异性。
• 黎曼求解器通过解耦流体与辐射子系统构建，假设耦合项与源项可忽略，分别求解两个子系统。
• 通过在辐射方程中使用流体速度，重新引入非守恒耦合项，从而在自由流状态下提高精度。
• 通过修改HLLE格式以处理辐射子系统中的刚性弛豫项，实现渐近保持，确保正确表现扩散极限行为。
• 将完整求解器集成到移动网格、时间隐式框架中，以应对天体物理模拟中时间与空间尺度的巨大差异。

实验结果

研究问题

• RQ1 能否设计一种单一数值格式，在无需区域分解的情况下，准确模拟从自由流到扩散的整个辐射流体动力学范围？
• RQ2 Godunov型格式应如何修改，以在保持自由流状态精度的同时，确保正确表现扩散极限？
• RQ3 高空间分辨率与自适应网格对解析原恒星形成中尖锐的辐射激波与光学厚度梯度有何影响？
• RQ4 M¹闭合模型在真实天体物理流中，对扩散与自由流状态之间的过渡捕捉效果如何？
• RQ5 所提出的格式能否在具有强光度与温度梯度的坍缩气体云长期模拟中保持稳定与准确？

主要发现

• 该格式在自由流与扩散极限下均成功捕捉了正确的物理行为，未出现人为的过渡误差。
• 在辐射激波测试中，方法在激波处高分辨率解析出尖锐的温度峰，证实其能够准确模拟超临界激波。
• 光度分布显示在吸积激波处出现尖锐跃迁，与能量耗散和辐射冷却一致，而辐射温度则平滑变化。
• 在中心核心区域，Eddington因子降至三分之一，确认处于扩散区；在激波附近上升，随后在大半径区域再次下降，反映出向自由流状态的过渡。
• 网格分辨率动态自适应，激波与核心区域实现高分辨率，空间分辨率跨越超过10个数量级，以捕捉陡峭梯度。
• 原恒星形成模拟显示密度对比超过18个数量级，速度分布表明外区膨胀、内区坍缩，且吸积激波被良好解析。