Journal
ASTRONOMY AND COMPUTING
Volume 12, Issue -, Pages 1-10Publisher
ELSEVIER SCIENCE BV
DOI: 10.1016/j.ascom.2015.04.003
Keywords
MHD simulation; Time-dependent ionization; Numerical calculation
Funding
- 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
Ask authors/readers for more resources
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/).
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available