期刊
ASTRONOMY AND COMPUTING
卷 12, 期 -, 页码 1-10出版社
ELSEVIER SCIENCE BV
DOI: 10.1016/j.ascom.2015.04.003
关键词
MHD simulation; Time-dependent ionization; Numerical calculation
资金
- NASA [NNX11AB61G, NNX12AB25G]
- National Science Foundation SHINE [AGS-1156076, AGS-1358342]
- Program 973 [2011CB811403, 2013CBA01503]
- NSFC [11273055, 11333007]
- CAS [XDB09040202]
- NASA [30900, NNX12AB25G] Funding Source: Federal RePORTER
- Directorate For Geosciences
- Div Atmospheric & Geospace Sciences [1156076] Funding Source: National Science Foundation
- Div Atmospheric & Geospace Sciences
- Directorate For Geosciences [1358342] Funding Source: National Science Foundation
Time-dependent ionization is important in astrophysical environments where the thermodynamical time scale is shorter than the ionization or recombination time scales. In this work, we report a FORTRAN program that performs fast non-equilibrium ionization calculations in post-processing based on hydrodynamics(HD) or magnetohydrodynamics(MHD) simulation results. Using HD or MHD simulation results, we track the movement of plasma in a Lagrangian framework, and obtain the evolutionary history of temperature and electron density. The time-dependent ionization equations are then solved by the Eigenvalue method. For any complex temperature and electron density histories, we introduce an adaptive time-step strategy to improve the computational efficiency. Our tests show that this program has advantages of high numerical stability and high accuracy. In addition, it is also easy to extend this solver to other HD and MHD simulations. This code is freely available for download from the Web. (C) 2015 The Authors. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据