<div dir="ltr">Hi Andrea,<div><br></div><div>That's an interesting problem you'd like to simulate! Regarding the error you're receiving, I usually see that happens when the timestep (dt) drops below 1e-10. In my case, usually a velocity gets too large. However, in your case it's interesting the error occurs immediately after mesh refinement (so possibly dx becomes quite small which also influences the CFL constraint). </div><div><br></div><div>Does that error occur even if you turn off AMR? (That is, set lrefine_min = lrefine_max = constant OR set nrefs = large).</div><div><br></div><div>Back to what you're trying to simulate, I'd think you'd want to use BDRY_VAR to simulate solids in FLASH. However, in that case you're not actually simulating the evolution of the internal structure of the solid material (instead, it applies boundary conditions).</div><div><br></div><div>Anyway, I think FLASH might be suitable for your problem since it sounds somewhat analogous to Mullen & Gammie 2020's work on the giant impact formation of the Moon, which uses Athena++ (and FLASH is a similar code).</div><div><br></div><div>Good luck!</div><div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr">--------<div>Ryan</div></div></div></div></div></div></div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Apr 12, 2021 at 2:30 PM Andrea Zorzi <<a href="mailto:zorzi@stanford.edu">zorzi@stanford.edu</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><font size="2"><span style="font-size:11pt">
<div>Dear Klaus,<br>
<br>
Thank you for your email. You are right, this is my mistake. Nonetheless, even setting:<br>
<br>
 sim_pAtm                = 1.<br>
 sim_pTarg               = 1.<br>
<br>
throws me the same error as described, even when more realistic values, e.g. sim_pAtm=101325 Pa and sim_pTarg=10^7 Pa, are input. If helpful, the new dat, log files and the command line output are attached.<br>
<br>
Andrea<br>
<br>
</div>
</span></font></div>
<div><font size="2"><span style="font-size:11pt">
<div><br>
<br>
<br>
<br>
> Il giorno 12 apr 2021, alle ore 11:11, Klaus Weide <<a href="mailto:klaus@flash.uchicago.edu" target="_blank">klaus@flash.uchicago.edu</a>> ha scritto:<br>
> <br>
> Andrea,<br>
> <br>
> Without trying to answer your questions substantially, I would just like <br>
> to point out that you have, in your flash.par,<br>
> <br>
>  sim_pAtm             = 0.<br>
>  sim_pTarg            = 0.<br>
> <br>
> so it seems you are initializing various pressures to 0.<br>
> That shouldn't really work, and I am a bit surprised your simulation got <br>
> even as far as it did.<br>
> <br>
> Klaus<br>
<br>
</div>
</span></font></div>
</div>

</blockquote></div>