<div dir="ltr">Dear Jeremy,<div><br></div><div>Do you see this instability issue only with FLASH4, but NOT with FLASH3? More specifically, I wonder what happens if you run the exact same problem with FLASH3. Have you tried this type of FLASH3 and FLASH4 comparison tests?</div><div><br></div><div>Cheers,</div><div>Dongwook</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Feb 18, 2018 at 10:41 AM, Jeremy S  Ritter <span dir="ltr"><<a href="mailto:jritter@mail.utexas.edu" target="_blank">jritter@mail.utexas.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hello Flash Users,<div><br></div><div>I have been using FLASH since version 3 to run simulations of supernova blastwaves in cosmological simulations. Since updating to version 4 I have had a lot of stability problems with the unsplit solver and strong shocks. The setup is similar to a Sedov example, with a spherical blastwave inserted into a realistic cosmological environment. We zoom in to high resolution to insert the blastwaves with radii of ~1 pc and 16-32 grid blocks across the diameter. The background density/temperature/pressure in the vicinity of the blastwaves is effectively constant when zoomed in to this scale. We refine a large volume around the blastwave so that the entirety of it is kept at the same refinement level while expanding.</div><div><br></div><div>The blastwaves will evolve as they should for many years, doubling or more in radius, and then suddenly within 1 step the solution will break down. This happens in the free expansion phase, before the reverse shock. The timestep dt will drop from ~1 year to 10^-7 seconds, tracer particles in the blastwave will be flung outside of the simulation box, gas density in some cells will drop by 20 orders of magnitude, etc. </div><div><br></div><div>I have tried both the HLLC and the hybrid HLLC solver with shock detection and shock lower CFL. By lowering the CFL even further, to 0.1, I was able to push the simulation further but eventually it breaks down too. I never had any similar issues in FLASH3. I was wondering if anybody has dealt with anything similar in FLASH4?</div><div><br></div><div>Cheers,</div><div>-Jeremy</div><div><br></div></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div style="font-size:12.8px"><font color="#000000" size="2" face="arial, helvetica, sans-serif"><br>=========================================</font></div><div style="font-size:12.8px"><font color="#000000" size="2" face="arial, helvetica, sans-serif">Dongwook Lee, Ph.D., Assistant Professor</font></div><div style="font-size:12.8px"><font color="#000000" size="2" face="arial, helvetica, sans-serif">Applied Mathematics and Statistics<br></font></div><div style="font-size:12.8px"><font color="#000000" size="2" face="arial, helvetica, sans-serif">University of California, Santa Cruz</font></div><div style="font-size:12.8px"><font color="#000000" size="2" face="arial, helvetica, sans-serif">Baskin Engineering, Room 353C</font></div><div style="font-size:12.8px"><font color="#000000" size="2" face="arial, helvetica, sans-serif">1156 High Street, Santa Cruz, CA 95064</font></div><div style="font-size:12.8px"><font color="#000000" size="2" face="arial, helvetica, sans-serif"><a href="https://users.soe.ucsc.edu/~dongwook/" target="_blank">https://users.soe.ucsc.edu/~dongwook/</a></font></div><div><br></div></div></div>
</div>