Microphysics
Microphysics copied to clipboard
AMReX-Astro
allow for multiple networks
It would be nice to have two different reaction networks built in and be able to chose which net we use in a zone based on thermodynamic conditions.
E.g., for a sub-Ch simulation, have a He burning net (like subch) and also a carbon-burning net, and use the He burning net in the He layer and the carbon burning once we ignite the core.
An alternative is to employ a mask that determines which species of a network we are currently integrating.
If we wanted to use the mask and also make it automated, perhaps we could specify some tolerance levels below which a given isotope is automatically not integrated (e.g. X < 1.0e-4).
To do this we need to relax the requirement that unevolved species are ordered after evolved species and generalize the unevolved species update. We should also avoid evaluating rates that are only needed to evolve those species.
We can more easily automate tracking the map between species and rates with pynucastro, for other networks we will need to construct such maps.
I think that for now, the only way we are going to do this is with the tabular NSE solver.