Radial recycling causes crash on cold start in 2D with core sources

[mikekryjak]

I have just started a simulation on a new grid and wanted to do it from scratch. I was only able to start the simulation without core sources and with everything turned down/off, as per the last time I did this with the MAST-U cases.

I just did some more testing and found that the case can be started from cold with all the bells and whistles (core sources present, final density and pressure settings, etc) if radial recycling is disabled.

This suggests there is likely a bug in recycling.cxx which causes an unphysical situation only at the very start of the solution.

[mikekryjak]

The error message is as below.

====== Back trace ======
-> virtual void Recycling::transform(Options&) on line 148 of '/users/mjk557/scratch/hermes-3/src/recycling.cxx'

====== Exception thrown ======
Setting invalid value for 'species:d:density_source'

There is a workaround: set sol_recycle = false, pfr_recycle = false, neutral_pump = false for the first iteration. You can then re-enable them upon restart.