<div dir="ltr">Dear Evan,<br>I tried the new EOS file but got the same error. Then by bypassing L162 in <span style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,Helvetica,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:16px">eos_nuclear.F90 </span><span style="color:rgb(0,0,0)"><font face="arial, sans-serif">and calling EOS by the temp instead of internal energy, it could manage to get the bounce (despite a few flags pop up about internal energy pumping). </font></span><div>I also could get the bounce with LS375 with lrefine_max=10 without changing in eos_nuclear.F90.</div><div><br><font color="#000000" face="arial, sans-serif">Can I get an EOS table with higher resolution to avoid this issue by the EOS driver or EOS maker you provided in your personal website?</font><div><span style="color:rgb(0,0,0)"><font face="arial, sans-serif">Because I am now trying to evolve postbounce from checkpoints files and got the following (3ms after the bounce):<br><br></font></span> i,rl,logrho,logtemp0,lt,ye,eps,eps0,abs(eps-eps0)/eps0<br> 21 0 14.188184188353635 -0.61304754054343624 -0.61304754054343624 0.27795080232191288 19.542939200343021 19.530696364609064 6.2685095837863099E-004<br> d2= NaN<br> Tried calling bisection... didn't help... :-/<br> Bisection error: 667<br> DRIVER_ABORT: [Eos_putData] ERROR Density or Internal Energy are zero after a call to EOS!<br> EOS: Did not converge in findtemp!<br> i,rl,logrho,logtemp0,lt,ye,eps,eps0,abs(eps-eps0)/eps0<br> 21 0 14.159389155587279 -0.43346281291235161 -0.43346281291235161 0.27796184415346487 19.534894764509421 19.529483548477000 2.7707932055596459E-004<br> d2= NaN<br> Tried calling bisection... didn't help... :-/<br> Bisection error: 667<br> EOS: Did not converge in findtemp!<br> i,rl,logrho,logtemp0,lt,ye,eps,eps0,abs(eps-eps0)/eps0<br> 21 0 14.323881185734283 -0.56350764342967907 -0.56350764342967907 0.27795290984794996 19.583574339654572 19.566803819457075 8.5709042479489340E-004<br> d2= NaN<br> Tried calling bisection... didn't help... :-/<br> Bisection error: 667<br>Abort(1) on node 4 (rank 4 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 1) - process 4<br><br>real 7m37.613s<br>user 0m0.045s<br>sys 0m0.100s<br><span style="color:rgb(0,0,0)"><font face="arial, sans-serif"><br><br><br>Kind Regards, </font></span></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, May 18, 2022 at 12:43 AM Evan Patrick O'Connor <<a href="mailto:evan.oconnor@astro.su.se">evan.oconnor@astro.su.se</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div dir="ltr">
<div id="gmail-m_7261513257979920782divtagdefaultwrapper" dir="ltr" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,Helvetica,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<div id="gmail-m_7261513257979920782divtagdefaultwrapper" dir="ltr" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,Helvetica,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<p>Hi,</p>
<p><br>
</p>
<p>Sorry for the silence, I haven't had a chance to look at this till now. </p>
<p><br>
</p>
<p>The point this is failing at is not near the edge of the table (rho~1.77e14g/cm^3, temp~3.5MeV, and ye=0.278), however, in energy space it is very close to the edge. For that density and ye, I've attached a plot of the internal energy vs temperature and
there is no temperature where it equals (logged, cgs) 19.5587, so the EOS cannot be inverted (the hydro solver gives the EOS an energy and the EOS has to find the temperature that corresponds to it).</p>
<p><br>
</p>
<p>We have a different LS220 table that we use for FLASH, I've added it to the page you linked (it was there before, but I missed it when moving the data around). It has a slightly higher resolution in density and temperature. It may give better results, especially
given that we use this for some reason (I forget now, it is a decade old). <a href="https://su.drive.sunet.se/index.php/s/EGP7wFx4XCLT37G" target="_blank">https://su.drive.sunet.se/index.php/s/EGP7wFx4XCLT37G</a></p>
<p><br>
</p>
<p>If this doesn't work (we still have issues like this some times with some progenitors) you can do what Naveen suggested. With the nuclear EOS it is a bit harder since just because the minimum energy (logged, cgs) for this rho and ye is 19.5599, at other
rho, ye pairs it will be very different. We accomplish this by bypassing the abort on L162 source/physics/Eos/EosNuclear/eos_nuclear.F90 and calling the EOS with the incoming temperature instead of the incoming energy, this avoid the inversion. It essentially
pumps internal energy into the zone, so be careful with it. We have a parameter file flag to turn this on and keep track of any such calls to ensure we don't do it too often. When it does occur it is in the situation you have now, right at bounce. </p>
<p><br>
</p>
<p>Also, FYI, I had suggested changing the file io_intfTypesModule.F90, not Flash.h, you should never edit the Flash.h as it get re generated with each setup.</p>
<p><br>
</p>
<p>Evan</p>
<p><br>
</p>
</div>
<hr style="display:inline-block;width:98%">
<div id="gmail-m_7261513257979920782divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> flash-users <<a href="mailto:flash-users-bounces@flash.rochester.edu" target="_blank">flash-users-bounces@flash.rochester.edu</a>> on behalf of kiyarash Taghiniyarami <<a href="mailto:kiyarash.niyarami@gmail.com" target="_blank">kiyarash.niyarami@gmail.com</a>><br>
<b>Sent:</b> Sunday, May 15, 2022 10:03:41 AM<br>
<b>To:</b> <a href="mailto:flash-users@flash.rochester.edu" target="_blank">flash-users@flash.rochester.edu</a><br>
<b>Subject:</b> [FLASH-USERS] 2D SN Core-Collapse Leakage scheme problem</font>
<div> </div>
</div>
<div>
<div dir="ltr">Dear all FLASH users,<br>
I tried to run 2D supernova core collapse for two different nuclear EOS of Lattimer and Sewsty (<a href="https://ttt.astro.su.se/~eoco/eos.html" id="gmail-m_7261513257979920782LPlnk806626" target="_blank">https://ttt.astro.su.se/~eoco/eos.html</a>) with K=180 and K=220 but I got
the errors which I attached the images of them.<br>
for the case of k=220 it runs into problem after 231 ms (it seems before the bounce) and for k=180 just after 24 ms<br>
Could anyone tell me what I have done wrong? <br>
(changing the resolution also did not resolve the issue)<br>
<br>
P.S. Trying to compile the FLASH, I also changed the Flash.h file and replace the line as Dr. O'Connor proposed:<br>
<span style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,Helvetica,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:16px">#FLASH_HAVE_HDF5_MOD</span><br>
<div><span style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,Helvetica,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:16px">with</span></div>
<div><span style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,Helvetica,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:16px">#define IO_USE_HDF5_MOD 1<br>
</span><span style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,Helvetica,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols;font-size:16px"><br>
I also attached the sample flash.par file</span></div>
</div>
</div>
</div>
</div>
</blockquote></div>