English

Eigenvalue Method for NEI Unit in FLASH Code

Instrumentation and Methods for Astrophysics 2019-03-01 v1 High Energy Astrophysical Phenomena

Abstract

We describe an improved nonequilibrium ionization (NEI) method that we have developed as an optional module for the FLASH magnetohydrodynamic simulation code. The method employs an eigenvalue approach rather than the earlier iterative ordinary differential equation approach to solve the stiff differential equations involved in NEI calculations. The new code also allows the atomic data to be easily updated from the AtomDB database. We compare both the updated atomic data and the methods separately. The new atomic data are shown to make a significant difference in some circumstances, although the general trends remain the same. Additionally, the new method also allows simultaneous calculation of the nonequilibrium radiative cooling, which is not included in the original method. The eigenvalue method improves the calculation efficiency overall with no loss of accuracy. We explore some common ways to present the NEI state with a sample simulation and find that using the average ionic charge difference from the equilibrium tends to be the clearest method.

Cite

@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}
}

Comments

15 pages, 8 figures

R2 v1 2026-06-23T07:53:24.405Z