[论文解读] Incorporation of cosmic ray transport into the ZEUS MHD code. Application for studies of Parker instability in the ISM

M. Hanasz, H. Lesch|ArXiv.org|Sep 24, 2003

Solar and Space Plasma Dynamics参考文献 31被引用 40

一句话总结

本文提出了一种数值算法，将宇宙射线输运通过扩散-平流方程整合进 ZEUS-3D MHD 代码中，利用平行和垂直扩散系数实现沿磁场和跨磁场的各向异性扩散。该方法在真实的银河系扩散系数范围内（例如，$3 \times 10^{28}\ \text{cm}^2\text{s}^{-1}$）保持数值稳定，并表明有限扩散会降低帕克不稳定性增长率，挑战线性近似，凸显宇宙射线驱动的星际介质动力学中的非线性效应。

ABSTRACT

We present a numerical algorithm for the incorporation of the active cosmic ray transport, into the ZEUS-3D magnetohydrodynamical code. The cosmic ray transport is described by the diffusion-advection equation. The applied form of the diffusion tensor allows for anisotropic diffusion of cosmic rays along and across the magnetic field direction, which is controlled by two parameters: the parallel and perpendicular diffusion coefficients. The implemented numerical algorithm is tested by comparison of the diffusive transport of cosmic rays to analytical solutions of the diffusion equation. Our method is numerically stable for a wide range of diffusion coefficients, including the realistic values inferred from the observational data for the Milky Way of about $ 6 imes 10^{28} \cm^2 \s^{-1}$. The presented algorithm is applied for exemplary simulations of the the Parker instability triggered by cosmic rays injected by a single SN remnant.

研究动机与目标

• 开发一种在 ZEUS-3D MHD 代码中实现活跃宇宙射线输运的数值稳定方法。
• 通过分别使用平行和垂直扩散系数，对沿磁场和跨磁场的宇宙射线扩散进行建模。
• 通过与扩散方程的解析解对比，对算法进行验证。
• 将扩展后的代码应用于研究由单个超新星遗迹引发的帕克不稳定性。
• 研究有限扩散系数如何影响帕克不稳定性非线性演化的特性，挑战线性分析中的假设。

提出的方法

• 将 ZEUS-3D MHD 代码修改为在交错网格上使用算子分裂和有限差分方法求解宇宙射线扩散-平流方程。
• 通过具有独立平行（$K_\parallel$）和垂直（$K_\perp$）扩散系数的张量实现各向异性扩散。
• 采用约束传输方法，在宇宙射线输运模拟中保持磁场散度为零。
• 将该方法应用于模拟来自单个超新星遗迹的宇宙射线注入磁化且具有引力分层的星际介质。
• 通过将模拟的宇宙射线扩散与扩散方程在广泛 $K_\parallel$ 值范围内的解析解进行比较，对算法进行验证。
• 对帕克不稳定性进行三维模拟，变化 $K_\parallel$ 值，包括真实值（$3 \times 10^{28}\ \text{cm}^2\text{s}^{-1}$）和较低值（$3 \times 10^{27}\ \text{cm}^2\text{s}^{-1}$）以作对比。

实验结果

研究问题

• RQ1 有限宇宙射线扩散，特别是具有真实 $K_\parallel \approx 3 \times 10^{28}\ \text{cm}^2\text{s}^{-1}$ 的情况，如何影响帕克不稳定性增长率？
• RQ2 当宇宙射线扩散受限制时，非线性效应在帕克不稳定性演化中占据多大程度的主导地位？
• RQ3 随着扩散系数的变化，宇宙射线压力的空间分布如何随时间演变？
• RQ4 该数值算法能否准确模拟广泛范围内的扩散系数下的宇宙射线输运，包括从观测数据推断出的系数？
• RQ5 有限扩散的引入是否会改变帕克不稳定的模态结构，例如改变交换模态与波动模态的相对贡献？

主要发现

• 该数值算法在广泛扩散系数范围内（包括银河系真实值 $3 \times 10^{28}\ \text{cm}^2\text{s}^{-1}$）成功模拟了宇宙射线输运，且保持高度数值稳定性。
• 在真实扩散系数（$K_\parallel = 3 \times 10^{28}\ \text{cm}^2\text{s}^{-1}$）下，宇宙射线扩散更均匀，峰值宇宙射线压力降低，初始浮力力弱于低扩散情况。
• 与低扩散情况（$K_\parallel = 3 \times 10^{27}\ \text{cm}^2\text{s}^{-1}$）相比，真实扩散情况下帕克不稳定性增长更慢，与线性分析的预期相矛盾。
• 在真实扩散情况下，波动模态在演化后期（约 $t = 150$）成为主导，而交换模态初始占主导，表明非线性发展被延迟。
• 在 $t = 100$ 时，真实扩散情况下的最大垂直速度低于低扩散情况，但随后因波动模态的出现而加速。
• 结果表明，对于低扩散系数，线性近似的适用性存疑，因为强非线性效应（如垂直速度接近声速）出现，需要进行完整的非线性模拟。