It looks like numerical viscosity is not turned in the default configuration. This seems counter-intuitive to me, as it should help stabilize simulations. What is the reason?
Additionally, we should update the beam heating configuration example to turn it on. It certainly should help with stability issues in flare simulations.