Impurities in the ETS-6
Impurity actor needed. Irena’s FORTRAN exists but needs to be wrapped so it can be integrated into the ETS-6.
Lackner scheme is NOT implemented in solver at present.
DC - See paper about where to implement Lackner scheme, http://zfn.mpdl.mpg.de/data/Reihe_A/37/ZNA-1982-37a-0931.pdf
DECISION: The temperature of the impurities should be calculated from the ratio between the total main ion pressure and the total main ion density. This option is already implemented.
DECISION: The present code parameters for the impurities in the ets_init are sufficient for the time being. But make minor changes
- Rename the option 'renormalisation' to something that clarifies that this refers to Corona equilibrium
- There are more than one interpretative mode; we need implicit for each state and implicit total density with coronal-split between charge state.
- Add documentation of these options.
DC - If we want to add bundling later on it'll be straight forward, but it's not a priority.
DC - Note that an edge models might change the boundary conditions
Lackner Scheme
Alexei - Why are neutrals not part of the Lackner scheme.
- Answer DC - neutrals treated by separate models - pushed to the future
Strategy as proposed by David Coster:
From the source terms for ionisation and recombination the rates can be calculated in the solver be dividing by the relevant densities. Since the rates do not depend on the ion densities the Lackner scheme can be evaluated using these rates. However, the top priority should be to get something working. The implementation of the Lackner scheme is a later concern.