[FLASH-USERS] Negative Internal Energy in 3T Simulation

[Michiel Bustraan]

??Hello everyone,


I've only recently started using FLASH, and have been trying to run spherical shock waves.

I used one of the existing test simulations, in radflaHD/RadBlastWave, as a basis for this.


It seemed to work fine in 1D, but when I attempted to adapt the same initial conditions to 2D, either Cartesian or spherical,

the code inevitably crashes due to a negative 3T energy. Looking at the log files, the code crashes due to a negative advected internal energy.

However, even before the crash, the log shows warnings regarding a negative internal energy. Plots made right before the crash don't indicate a problem.


I wanted to ask if anyone else had similar problems, and what the cause might be. Would switching to a 1T setup be the easiest solution?


Kind Regards,

Michiel Bustraan
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20160816/c5c661c6/attachment.htm>