[论文解读] Pencil-Code special edition: modelling supernova driven turbulence

Frederick A. Gent, Mordecai‐Mark Mac Low+3 more|arXiv (Cornell University)|Jun 5, 2018

Solar and Space Plasma Dynamics参考文献 24被引用 1

一句话总结

本文通过引入人工粘性、热导率、质量扩散以及广义CFL条件，扩展了Pencil Code以模拟超新星驱动的湍流，从而准确捕捉高马赫数激波。该方法成功模拟了一维和三维激波管以及超新星遗迹演化，其结果在绝热与冷却状态下与解析解和半解析解高度一致。

ABSTRACT

High Mach number shocks are ubiquitous in interstellar turbulence. The Pencil Code is particularly well suited to the study of magnetohydrodynamics in weakly compressible turbulence and the numerical investigation of dynamos because of its high-order advection and time evolution algorithms. However, the high-order algorithms and lack of Riemann solver to follow shocks make it less well suited to handling high Mach number shocks, such as those produced by supernovae (SNe). Here, we outline methods required to enable the code to efficiently and accurately model SNe, using parameters that allow stable simulation of SN-driven turbulence, in order to construct a physically realistic galactic dynamo model. These include the resolution of shocks with artificial viscosity, thermal conductivity, and mass diffusion; the correction of the mass diffusion terms; and generalization of the Courant condition to include all source terms in the momentum and energy equations. We test our methods with the numerical solution of the 1D Riemann shock tube (Sod 1978), also extended to a 1D adiabatic shock with parameters and Mach number relevant to SN shock evolution, including shocks with radiative losses. We extend our test with the 3D numerical simulation of individual SN remnant evolution for a range of ambient gas densities typical of the interstellar medium and relate these to the analytical solutions of Sedov-Taylor (adiabatic) and the snowplough and Cioffi et al. (1988) results incorporating cooling and heating processes.

研究动机与目标

• 使Pencil Code能够模拟超新星（SN）爆发中典型的高马赫数激波，此前由于缺乏Riemann求解器和高阶格式，该代码在处理此类激波时存在困难。
• 开发数值方法，以在具有强激波的弱可压缩磁流体动力学流动中实现激波的稳定与精确解析。
• 引入人工粘性、热导率和质量扩散等物理过程，以真实地模拟超新星驱动的湍流。
• 将CFL条件推广至包含动量和能量方程中的所有源项，以提升稳定性和精度。
• 通过一维激波管和三维超新星遗迹演化在不同星际介质密度下的解析解，对增强后的代码进行验证。

提出的方法

• 引入人工粘性以解析传统高阶格式失效时的强激波。
• 实现热导率和质量扩散项，以模拟激波前缘的能量传递与混合过程。
• 修正质量扩散项，以确保与守恒定律及物理行为的一致性。
• 将CFL条件推广至包含动量和能量方程中的所有源项，提升时间步长的稳定性。
• 利用一维Riemann激波管（Sod 1978）和绝热激波测试验证数值方法。
• 将验证扩展至不同环境气体密度下的一系列三维单个超新星遗迹模拟，与Sedov-Taylor模型及Cioffi et al. (1988)模型进行比较。

实验结果

研究问题

• RQ1 能否通过引入人工粘性和改进的数值格式，将Pencil Code扩展以准确模拟超新星产生的高马赫数激波？
• RQ2 热导率和质量扩散对Pencil Code中激波模拟的稳定性和精度有何影响？
• RQ3 广义CFL条件在强源项存在的情况下，能在多大程度上提升数值稳定性？
• RQ4 三维超新星遗迹模拟在绝热与冷却状态下，对解析解的再现程度如何？
• RQ5 增强后的Pencil Code能否在星际介质密度变化范围内真实模拟超新星驱动的湍流？

主要发现

• 增强后的Pencil Code通过引入人工粘性、热导率和校正的质量扩散，成功解析了高马赫数激波，实现了超新星驱动湍流的稳定模拟。
• 一维激波管测试（Sod 1978）与解析解高度一致，验证了数值实现的正确性。
• 采用超新星相关参数的一维绝热激波测试证实了该方法在高精度下模拟强激波的能力。
• 在一维测试中引入辐射损失后，物理真实性得到提升，且激波结构与预期相符。
• 三维超新星遗迹模拟在绝热状态下重现了Sedov-Taylor解，并在包含冷却与加热时与雪plough模型及Cioffi et al. (1988)模型一致。
• 广义CFL条件显著提升了在不同源项强度下时间步长的稳定性和模拟鲁棒性。