[FLASH-USERS] Helmhotlz EOS question

[Dean Townsley]

Hi Malia,


Usually this indicates that your abundances have gotten corrupted in 
some way.  There are various ways that this can happen, depending on 
what kind of simulation you are running.


If you are using the built-in nuclear networks, NaNs are not capture 
quite properly.  I have a patch for that if that is what you are using 
(attached).  If you aren't using the nuclear networks, then the problem 
is likely somewhere else.


Hopefully that helps some.  If you say a bit more about your simulation 
someone might be able to make a better guess at what might be going on.  
There are lots of problems that can manifest this way.


Dean


On 05/03/2017 04:02 PM, Jenks, Malia T. wrote:
> 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/20170503/d1b9c984/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: fixnanhandling.diff
Type: text/x-patch
Size: 6594 bytes
Desc: not available
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20170503/d1b9c984/attachment.bin>