Saved in:
Bibliographic Details
Main Authors: Mathew, Arun, Mackey, Jonathan, Celeste, Maggie, Haworth, Thomas J., Mellema, Garrelt
Format: Preprint
Published: 2024
Subjects:
Online Access:https://arxiv.org/abs/2409.18904
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866913726374871040
author Mathew, Arun
Mackey, Jonathan
Celeste, Maggie
Haworth, Thomas J.
Mellema, Garrelt
author_facet Mathew, Arun
Mackey, Jonathan
Celeste, Maggie
Haworth, Thomas J.
Mellema, Garrelt
contents While many astrophysical plasmas can be modelled successfully assuming ionisation and thermal equilibrium, in some cases this is not appropriate and a non-equilibrium approach is required. In nebulae around evolved stars the local elemental abundances may also strongly vary in space and time. Here we present a non-equilibrium multi-ion module developed for the fluid-dynamics code \textsc{pion}, describing the physical processes included and demonstrating its capabilities with some test calculations. A non-equilbrium ionisation solver is developed that allows arbitrary elemental abundances for neutral and ionised (but not molecular) gas, for the elements H, He, C, N, O, Ne, Si, S and Fe. Collisional ionisation and recombination, photoionisation and charge-exchange reactions are included, and ion-by-ion non-equilibrium radiative cooling is calculated based on the instantaneous ion fractions of each element. Element and ion mass-fractions are advected using passive scalars, operator-split from the microphysical processes. The module is validated by comparing with equilibrium and non-equilibrium calculations in the literature. Effects of charge exchange on ion abundances in cooling plasmas are discussed. Application to modelling shocks and photoionised H~\textsc{ii} regions is demonstrated. The time-dependent expansion of a Wolf-Rayet nebula is studied, including photoionisation and collisional processes, and spectral-line luminosities calculated for non-equilibrium and equilibrium plasma states. The multi-ion module enables simulation of ionised plasmas with spatially varying elemental abundances using self-consistent ion abundances and thermal evolution. This allows prediction of spectral lines in UV, optical, IR and X-ray even in cases where the plasma is out of ionisation equilibrium.
format Preprint
id arxiv_https___arxiv_org_abs_2409_18904
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle A multi-ion non-equilibrium solver for ionised astrophysical plasmas with arbitrary elemental abundances
Mathew, Arun
Mackey, Jonathan
Celeste, Maggie
Haworth, Thomas J.
Mellema, Garrelt
Astrophysics of Galaxies
Computational Physics
While many astrophysical plasmas can be modelled successfully assuming ionisation and thermal equilibrium, in some cases this is not appropriate and a non-equilibrium approach is required. In nebulae around evolved stars the local elemental abundances may also strongly vary in space and time. Here we present a non-equilibrium multi-ion module developed for the fluid-dynamics code \textsc{pion}, describing the physical processes included and demonstrating its capabilities with some test calculations. A non-equilbrium ionisation solver is developed that allows arbitrary elemental abundances for neutral and ionised (but not molecular) gas, for the elements H, He, C, N, O, Ne, Si, S and Fe. Collisional ionisation and recombination, photoionisation and charge-exchange reactions are included, and ion-by-ion non-equilibrium radiative cooling is calculated based on the instantaneous ion fractions of each element. Element and ion mass-fractions are advected using passive scalars, operator-split from the microphysical processes. The module is validated by comparing with equilibrium and non-equilibrium calculations in the literature. Effects of charge exchange on ion abundances in cooling plasmas are discussed. Application to modelling shocks and photoionised H~\textsc{ii} regions is demonstrated. The time-dependent expansion of a Wolf-Rayet nebula is studied, including photoionisation and collisional processes, and spectral-line luminosities calculated for non-equilibrium and equilibrium plasma states. The multi-ion module enables simulation of ionised plasmas with spatially varying elemental abundances using self-consistent ion abundances and thermal evolution. This allows prediction of spectral lines in UV, optical, IR and X-ray even in cases where the plasma is out of ionisation equilibrium.
title A multi-ion non-equilibrium solver for ionised astrophysical plasmas with arbitrary elemental abundances
topic Astrophysics of Galaxies
Computational Physics
url https://arxiv.org/abs/2409.18904