[FLASH-USERS] Helmhotlz EOS question

[Jenks, Malia T.]

I am not using nuclear burning for the simulation. I do call physics/sourceTerms/Burn/BurnMain/nuclearBurn and Simulation/SimulationComposition in my Config file to get the species information from the burn unit but I have useBurn set to false. I use the same mass fractions for all points on my initial grid, so all the abar and zbar values should be the same. I am not sure that the error isn't appearing in the t=0 data. I know it is after successfully reading in my data file and then I get the following outputs before the error message:

  iteration, no. not moved =            0           9
  iteration, no. not moved =            1           2
  iteration, no. not moved =            2           0
 refined: total leaf blocks =           64
 refined: total blocks =           80
 [amr_morton_process]: Initializing surr_blks using standard orrery implementation

Thanks

Malia


________________________________
From: Jenks, Malia T.
Sent: Wednesday, May 3, 2017 3:02:52 PM
To: flash-users at flash.uchicago.edu
Subject: Helmhotlz EOS question


I am getting the following error from the helmholtz eos about negative pressure:

 [eos_helm] Negative total pressure.
   values: dens,temp:   2.321442471007820E-004   1549.01220000000

   values: abar,zbar:                 Infinity                     NaN

   coulomb coupling parameter Gamma:                      NaN

  Prad   1.45198E-02 Pion  0.00000E+00 Pele           NaN Pcoul  0.00000E+00

  Ptot            NaN df_d    8.31451E+21


I know that the initial mass fraction supplied to the eos is correct. The initial call of the eos sets abar and zbar correctly. I am not sure what is causing them to error. I also don't know why I am getting a NaN for the electron Pressure. Could anyone help point me to what is going wrong and how I can solve it?

Thanks,

Malia Jenks

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