[FLASH-USERS] Crashing with a sudden jump in kinetic energy FLASH4.4

[Jenks, Malia T.]

I am using FLASH4.4. I am running a 2D simulation in cylindrical coordinates using the MHD staggered mesh solver and a uniform grid. My simulation will run for about 30 seconds (simulation time) then it crashes. The error I get can vary based on the cfl and minimum time step I give it, but the issue is that the y momentum begins to sharply decrease and then the Kinetic energy jumps significantly and the run crashes.


Here is part of an output file for one of the runs:


This run had a cfl condition of .2, and a minimum time step of 1e-6 seconds.

  time                                           y-momentum                              E_kinetic                                     E_internal                          MagEnergy

  3.165470637854052782E+01  5.587253776952739325E+30  4.753600339710994969E+49  8.743916647261655525E+49  1.754228010830377560E+42
  3.165493812157636455E+01  4.754893533206718301E+30  4.754122794831045632E+49  8.743926772696261717E+49  1.754446232246603526E+42
  3.165493813157636538E+01  2.596391772160589512E+26  4.979926949539627564E+49  8.743926773126580588E+49  1.754485014199999339E+42
  3.165493814157636621E+01  1.616370861886859215E+26  7.629586798455982783E+54  8.743926773725260724E+49  1.778462439566404185E+42

The error given for this is:

[eos_helm] Negative total pressure.
   values: dens,temp:    471500237.463805       5.857646036122962E+015

   values: abar,zbar:    14.1176470588276        7.05882352941380

   coulomb coupling parameter Gamma:   6.231188937874524E-007

  Prad   2.96915E+48 Pion  1.62660E+31 Pele  -4.17392E+49 Pcoul  0.00000E+00

  Ptot   -3.87700E+49 df_d   -7.51000E+32

 DRIVER_ABORT: [eos_helm] ERROR: negative total pressure.


On similar runs I have gotten errors suggesting lowering the cfl condition and something about needing a finer grid.

I am not sure what is going on to cause this issue. I have also tried increasing cvisc but it still crashes. Any help figuring out what is causing this issue would be greatly appreciated.

Malia Jenks

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20180710/b74651cd/attachment.htm>