Field ionization rates in AtomicPhysics(FLYonPIC) do not reduce with consecutive sub-steps since we do not update the local field yet.
These high rates require a large number o atomicPhyiscs sub steps to resolve and therefore lead to a slow down of the simulation.
To avoid this this PR adds a new sub-stage checking for states with time step lengths above a user defined limit and instantly applies field ionization transitions to ions in a state found to have a field ionization loss rate above the limit
[x] requires PR #5164 to be merged first
[ ] needs to be rebased to dev
This PR will break existing setups using a custom atomicPhysics.param-file. To fix a setup add the following lines inside the RateSolverParam -struct
// maximum number of atomicPhysics sub-steps per PIC time step for fieldIonization
static constexpr uint32_t maximumNumberSubStepsFieldIonization = 2000;
BrianMarre
commented
1 week ago
Field ionization rates in AtomicPhysics(FLYonPIC) do not reduce with consecutive sub-steps since we do not update the local field yet. These high rates require a large number o atomicPhyiscs sub steps to resolve and therefore lead to a slow down of the simulation. To avoid this this PR adds a new sub-stage checking for states with time step lengths above a user defined limit and instantly applies field ionization transitions to ions in a state found to have a field ionization loss rate above the limit
This PR will break existing setups using a custom
atomicPhysics.param-file. To fix a setup add the following lines inside theRateSolverParam-struct