[FLASH-USERS] 2D SN Core-Collapse Leakage scheme problem

[Evan Patrick O'Connor]

Hi,


Sorry for the silence, I haven't had a chance to look at this till now.


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).


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). https://su.drive.sunet.se/index.php/s/EGP7wFx4XCLT37G


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.


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.


Evan


________________________________
From: flash-users <flash-users-bounces at flash.rochester.edu> on behalf of kiyarash Taghiniyarami <kiyarash.niyarami at gmail.com>
Sent: Sunday, May 15, 2022 10:03:41 AM
To: flash-users at flash.rochester.edu
Subject: [FLASH-USERS] 2D SN Core-Collapse Leakage scheme problem

Dear all FLASH users,
I tried to run 2D supernova core collapse for two different nuclear EOS of Lattimer and Sewsty (https://ttt.astro.su.se/~eoco/eos.html)  with K=180 and K=220 but I got the errors which I attached the images of them.
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
Could anyone tell me what I have done wrong?
(changing the resolution also did not resolve the issue)

P.S. Trying to compile the FLASH, I also changed the Flash.h file and replace the line  as Dr. O'Connor  proposed:
#FLASH_HAVE_HDF5_MOD
with
#define IO_USE_HDF5_MOD 1

I also attached the sample flash.par file
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20220517/6b7061a2/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: eint_vs_temp_ls220.png
Type: image/png
Size: 68226 bytes
Desc: eint_vs_temp_ls220.png
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20220517/6b7061a2/attachment-0001.png>