FLASH 代码中 NEI 单元的特征值方法
天体物理仪器与方法
2019-03-01 v1 高能天体物理现象
摘要
我们描述了一种改进的非平衡电离（NEI）方法，该方法已作为 FLASH 磁流体力学模拟代码的一个可选模块开发。该方法采用特征值方法，而非早期的迭代常微分方程方法来求解 NEI 计算中涉及的刚性微分方程。新代码还允许从 AtomDB 数据库轻松更新原子数据。我们分别比较了更新后的原子数据和方法。新的原子数据显示在某些情况下有显著差异，尽管总体趋势保持不变。此外，新方法还允许同时计算非平衡辐射冷却，而原始方法中未包含此项。特征值方法在整体计算效率上有所提升，且不损失精度。我们通过一个示例模拟探讨了展示 NEI 状态的几种常用方式，发现使用相对于平衡态的平均离子电荷差往往是最清晰的方法。
引用
@article{arxiv.1902.10715,
title = {Eigenvalue Method for NEI Unit in FLASH Code},
author = {Gao-Yuan Zhang and Adam Foster and Randall Smith},
journal= {arXiv preprint arXiv:1902.10715},
year = {2019}
}
备注
15 pages, 8 figures