Eigenvalue Method for NEI Unit in FLASH Code

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.

15 pages, 8 figures

References in corpus (11)

• A Multi-Code Analysis Toolkit for Astrophysical Simulation Data
• Supernova remnants: the X-ray perspective
• A new approach to analyzing solar coronal spectra and updated collisional ionization equilibrium calculations. II. Additional ionization rate coefficients
• Ionization Equilibrium Timescales in Collisional Plasmas
• Nonequilibrium ionization states and cooling rates of the photoionized enriched gas
• Unveiling the spatial structure of the overionized plasma in the supernova remnant W49B
• Non-equilibrium Ionization State of Warm-Hot Intergalactic Medium
• Numerical Simulations of Supernova Remnant Evolution in a Cloudy Interstellar Medium
• The Metal-enriched Thermal Composite Supernova Remnant Kesteven 41 (G337.8-0.1) in a Molecular Environment
• A Useful Approximation to the Cooling Coefficient of Trace Elements
• Non-equilibrium ionization states in galaxy clusters

Cited by in corpus (3)

• Non-equilibrium ionization in mixed-morphology supernova remnants
• Non-equilibrium ionisation and radiative transport in an evolving supernova remnant
• Non-equilibrium Ionization Modeling of Petschek-type Shocks in Reconnecting Current Sheets in Solar Eruptions