<div dir="ltr">Dear Kiya,<div><br></div><div>You might have better luck lowering the <b>cfl</b> number as this more directly limits the timestep. <b>dtinit</b> and <b>tstep_change_factor </b>only limits the initial timestep and subsequent steps by <b>tstep_change_factor </b>until you reach the hydro time step determined by the cfl condition. <a href="https://flash.rochester.edu/site/flashcode/user_support/flash4_ug_4p62/node51.html#SECTION04115000000000000000">Here</a> is the relevant section from the users guide on how FLASH determines the timestep.</div><div><br></div><div><br></div><br clear="all"><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div dir="ltr" style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px">*************************************************************************</div><div dir="ltr" style="color:rgb(0,0,0);font-family:Helvetica;font-size:12px">Adam Reyes<div><img src="https://flash.rochester.edu/site/common/images/email_signature/FLASH.jpg"><br></div><div>Code Group Leader, Flash Center for Computational Science <span> </span><br></div><div>Research Scientist, Dept. of Physics and Astronomy</div><div>University of Rochester</div><div>River Campus: Bausch and Lomb Hall, 369 <span> </span><br></div><div>500 Wilson Blvd. PO Box 270171, Rochester, NY 14627</div><div dir="ltr">Email<span> </span><a href="mailto:adam.reyes@rochester.edu" target="_blank">adam.reyes@rochester.edu</a></div><div>Web<span> </span><a href="https://flash.rochester.edu/" target="_blank">https://flash.rochester.edu</a><br></div><div> (he / him / his)</div><div><img src="https://flash.rochester.edu/site/common/images/email_signature/FLASH-pride-sml.png"><br></div>*************************************************************************</div></div></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, May 16, 2022 at 4:18 AM Naveen Yadav <<a href="mailto:naveen.phys@gmail.com">naveen.phys@gmail.com</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 class="gmail_default" style="font-family:monospace,monospace;font-size:small">Dear Haakon and Kiyarash,</div><div class="gmail_default" style="font-family:monospace,monospace;font-size:small"><br></div><div class="gmail_default" style="font-family:monospace,monospace;font-size:small">I think the second error is caused by your thermodynamic variables moving out of the range of the EoS you are using. You can probably fix it by setting limits on the thermodynamic variables (floor values). It seems you are running into a low temperature, low Ye regime. You can plot the values of temp and Ye on a map of EoS in the T-Ye plane and see if you are outside the boundaries. In this case, you might need to patch up with a low-temperature EoS. I hope this helps along with a lower value of CFL condition.</div><div class="gmail_default" style="font-family:monospace,monospace;font-size:small"><br></div><div class="gmail_default" style="font-family:monospace,monospace;font-size:small">Regards<br clear="all"></div><div><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div><div><font size="2"><span style="font-family:monospace,monospace">Naveen Yadav<br></span></font></div><font size="2"><span style="font-family:monospace,monospace">Post-doctoral Researcher,<br></span></font></div><font size="2"><span style="font-family:monospace,monospace">Max Planck Institute For Astrophysics,<br></span></font></div><font size="2"><span style="font-family:monospace,monospace">Garching, Germany.</span></font><br></div></div><div dir="ltr"><font size="2"><span style="font-family:monospace,monospace"><br></span></font></div><div dir="ltr"><font face="monospace">"The purpose of (scientific) computing is insight, not numbers" -- Richard Hamming.</font><font size="2"><span style="font-family:monospace,monospace"><br></span></font></div></div></div></div></div></div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, May 16, 2022 at 8:44 AM Haakon Andresen <<a href="mailto:haakon.andresen@astro.su.se" target="_blank">haakon.andresen@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>
<div id="gmail-m_5694376117704111609gmail-m_5930247685054064909divtagdefaultwrapper" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif" dir="ltr">
<p>Dear Kiya,</p>
<p><br>
</p>
<p>I do not have a particular parfile nor a concrete FLASH setup to </p>
<p>recommend, I am still somewhat of a FLASH novice. I am not sure</p>
<p>what exact spatial resolution your setup results in, what are the smallest</p>
<p>and largest grid cells you have in the core during bounce? </p>
<p><br>
</p>
<p>Changing the timestep factor might not have solved your problem, but it</p>
<p>improved it. Your simulation now progresses further and I am willing to bet</p>
<p>that a greater constraint on the timestep will allow even further progress.</p>
<p>I suspect that the problem is indeed related to the rapid changes occurring during</p>
<p>bounce. Plotting the central density, temperature and pressure might help you</p>
<p>check if is indeed the case. You can also continue to increase your temporal and</p>
<p>spatial resolution and check if your simulation progresses further and further.</p>
<p><br>
</p>
<p>Best,<br>
Haakon</p>
<p><br>
</p>
<p><br>
</p>
<p><br>
</p>
<p><br>
</p>
<p></p>
</div>
<hr style="display:inline-block;width:98%">
<div id="gmail-m_5694376117704111609gmail-m_5930247685054064909divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>From:</b> kiyarash Taghiniyarami <<a href="mailto:kiyarash.niyarami@gmail.com" target="_blank">kiyarash.niyarami@gmail.com</a>><br>
<b>Sent:</b> 15 May 2022 19:09:59<br>
<b>To:</b> Haakon Andresen; <a href="mailto:flash-users@flash.rochester.edu" target="_blank">flash-users@flash.rochester.edu</a><br>
<b>Subject:</b> Re: [FLASH-USERS] 2D SN Core-Collapse Leakage scheme problem</font>
<div> </div>
</div>
<div>
<div dir="ltr">Dear Haakon Andresen,<br>
<br>
Thanks for your reply. for the second error do you have any parfile or time and resolution available to suggest?<br>
I changed the resolution to:<br>
# Grid / Refinement<br>
nblockx = 15<br>
nblocky = 30<br>
nblockz = 1<br>
<br>
gr_lrefineMaxRedDoByLogR = .true.<br>
gr_lrefineMaxRedDoByTime = .FALSE.<br>
gr_lrefineMaxRedLogBase = 10.0<br>
gr_lrefineMaxRedRadiusFact = 0.15<br>
gr_lrefineMaxRedTRef = 2.0 <br>
gr_lrefineMaxRedTimeScale = 0.5<br>
lrefine_max = 11 <br>
nrefs = 2<br>
lrefine_min = 1<br>
refine_var_1 = "dens"<br>
refine_var_2 = "pres"<br>
refine_var_3 = "entr"<br>
refine_var_4 = "none"<br>
refine_cutoff_1 = 0.8 <br>
refine_cutoff_2 = 0.8<br>
refine_cutoff_3 = 0.8<br>
refine_cutoff_4 = 0.8<br>
<br>
but this time instead of 231 ms it gave me the same error at later time (260 ms). changing the
<b>tstep_change_factor </b> , and <b>dtinit</b> did not resolve the issue.<br>
<br>
Kind regards,<br>
Kiya</div>
<br>
<div class="gmail_quote">
<div dir="ltr" class="gmail_attr">On Sun, May 15, 2022 at 1:42 PM Haakon Andresen <<a href="mailto:haakon.andresen@astro.su.se" target="_blank">haakon.andresen@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_5694376117704111609gmail-m_5930247685054064909gmail-m_1816844387868435849divtagdefaultwrapper" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif" dir="ltr">
<div id="gmail-m_5694376117704111609gmail-m_5930247685054064909gmail-m_1816844387868435849divtagdefaultwrapper" dir="ltr" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<p>Dear Kiyarash Taghiniyarami,</p>
<p><br>
</p>
<p>The first error (for LS180) seem to indicate a lack of storage space on your drive. </p>
<p>You see this written by the error handler in your first screenshot. Make sure that there is space</p>
<p>and that you are not trying to create files larger than the maximum filesize allowed by your system. </p>
<p><br>
</p>
<p>The second errors seems to be related to physics. During bounce things change very fast and it </p>
<p>might be that you need to reduce your timestep during this time. However, I am not sure about this</p>
<p>because I do not know exactly how flash regulates the timestep. The error is in the EOS module where the</p>
<p>routine fails to converge, I have seen this happen around bounce and for me it always helped to restrict dt or</p>
<p>the resolution. However, I was working with an other code at the time.</p>
<p><br>
</p>
<p>Best,<br>
Haakon Andresen</p>
</div>
<hr style="display:inline-block;width:98%">
<div id="gmail-m_5694376117704111609gmail-m_5930247685054064909gmail-m_1816844387868435849divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><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> 15 May 2022 10:03:41<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_5694376117704111609gmail-m_5930247685054064909gmail-m_1816844387868435849LPlnk532912" 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>
</div>
</div>
_______________________________________________<br>
flash-users mailing list<br>
<a href="mailto:flash-users@flash.rochester.edu" target="_blank">flash-users@flash.rochester.edu</a><br>
<br>
For list info, including unsubscribe:<br>
<a href="https://flash.rochester.edu/mailman/listinfo/flash-users" rel="noreferrer" target="_blank">https://flash.rochester.edu/mailman/listinfo/flash-users</a><br>
</blockquote></div>
_______________________________________________<br>
flash-users mailing list<br>
<a href="mailto:flash-users@flash.rochester.edu" target="_blank">flash-users@flash.rochester.edu</a><br>
<br>
For list info, including unsubscribe:<br>
<a href="https://flash.rochester.edu/mailman/listinfo/flash-users" rel="noreferrer" target="_blank">https://flash.rochester.edu/mailman/listinfo/flash-users</a><br>
</blockquote></div>