Microphysics

Microphysics copied to clipboard

The maximum order in VODE is [hardcoded to 5](https://github.com/starkiller-astro/Microphysics/blob/21.10/integration/VODE/vode_type.H#L50). We should experiment to see what happens if we use a lower limit on maximum order. On GPUs thread-local memory is...

maxpkatz

add the Potekhin and Chabrier (2010,2013) EOS

3

The Potekhin and Chabrier (2010,2013) EOS is freely available (in Fortran) here: http://www.ioffe.ru/astro/EIP/ This EOS seems to be preferred by MESA: https://docs.mesastar.org/en/r15140/eos/overview.html There is also support here for a magnetized...

zingale

other 3-alpha rates?

1

We are still using an old rate for 3-alpha in most nets. Some updates are found in: Angulo, C., Arnould, M., Rayet, M., Descouvemont, P., Baye, D., et al. (1999):...

zingale

this addresses #762

zingale

we should be able to query an EOS to see if an output is supported

1

Not all EOSes fill all the thermodynamic data. We should be able to query whether something is supported? Like `is_filled("ne")`.

zingale

analytic Jacobian species derivatives are missing terms

1

We want to express our Jacobian in terms of (rho, e, X_k), but the reaction system is more naturally in terms of (ρ, T, X_k). Consider dX/dt = f(ρ, T,...

zingale

what if we pass VODE failure message out of burn through burn_t

4

Currently we don't dump out the burn failure when we are using CUDA, which makes debugging hard on GPUs. What if we add a `std::string` to the `burn_t` that can...

zingale

It would be good to have a script that ran a single zone burn for a given network using all of the integrators and plotted the mass fraction evolution vs....

zingale

use C++ autodiff to do Jacobians

4

Instead of hand-coding the Jacobians, we should explore using a C++ autodiff package, e.g., https://autodiff.github.io/ to construct the analytic Jacobians.

zingale

for simplified SDC, if T is outside of range do we still need to include advective terms?

1

In `vode_rhs_simplified_sdc.H`, we have a check: ``` if (state.T = MAX_TEMP) { for (int n = 1; n

zingale