<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Arial,Helvetica,sans-serif;" dir="ltr">
<p></p>
<div>
<div>
<div>
<div>I am not using nuclear burning for the simulation. I do call physics/sourceTerms/Burn/<wbr>BurnMain/nuclearBurn and Simulation/<wbr>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: <br>
<br>
iteration, no. not moved = 0 9<br>
iteration, no. not moved = 1 2<br>
iteration, no. not moved = 2 0<br>
refined: total leaf blocks = 64<br>
refined: total blocks = 80<br>
[amr_morton_process]: Initializing surr_blks using standard orrery implementation<br>
</div>
<br>
</div>
Thanks<br>
<br>
</div>
Malia </div>
<br>
<p></p>
</div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> Jenks, Malia T.<br>
<b>Sent:</b> Wednesday, May 3, 2017 3:02:52 PM<br>
<b>To:</b> flash-users@flash.uchicago.edu<br>
<b>Subject:</b> Helmhotlz EOS question</font>
<div> </div>
</div>
<div>
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Arial,Helvetica,sans-serif;" dir="ltr">
<p></p>
<div>
<div>
<div>
<div>I am getting the following error from the helmholtz eos about negative pressure:<br>
<br>
[eos_helm] Negative total pressure.<br>
values: dens,temp: 2.321442471007820E-004 1549.01220000000 <br>
<br>
values: abar,zbar: Infinity NaN <br>
<br>
coulomb coupling parameter Gamma: NaN <br>
<br>
Prad 1.45198E-02 Pion 0.00000E+00 Pele NaN Pcoul 0.00000E+00 <br>
<br>
Ptot NaN df_d 8.31451E+21 <wbr>
<br>
<br>
<br>
</div>
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?<br>
<br>
</div>
Thanks,<br>
</div>
<br>
Malia Jenks</div>
<br>
<p></p>
</div>
</div>
</body>
</html>