<div dir="ltr">Hi Michiel,<div><br></div><div>You could be having problems for a number of reasons. Some questions that might help narrow things down:</div><div><br></div><div>1) Are you using gravity? This looks like a stellar envelope.</div><div>2) Which hydro solver are you using (split PPM, unsplit, or other)?</div><div>3) If using the split PPM solver, are you using the "hybrid_riemann" flag to deal with the odd-even decoupling instability (due to shocks aligned with the grid)?</div><div>4) Regarding the inner regions, how many cells per decade in radius are you using (i.e., how well resolved is the density gradient)?</div><div>5) What is the largest Mach number of your flow? If it is too large, you might want to use the "eint_switch" flag to evolve the internal energy instead of total energy when the Mach number exceeds a threshold value.</div><div><br></div><div>Best,</div><div><br></div><div>Rodrigo</div><div><br></div><div><br></div><div><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Mar 22, 2017 at 7:10 AM, Michiel Bustraan <span dir="ltr"><<a href="mailto:michiel.bustraan@astro.su.se" target="_blank">michiel.bustraan@astro.su.se</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" style="font-size:12pt;color:#000000;background-color:#ffffff;font-family:Calibri,Arial,Helvetica,sans-serif">
<p>Dear FLASH Users,<br>
</p>
<p><br>
</p>
<p>I've been running 2D simulations of shocks in spherical simulations.<br>
</p>
<p>The initial conditions involve a high density around the origin.<br>
</p>
<p><br>
</p>
<p>I had to introduce random temperature fluctuations in the inner region to ensure refinement.<br>
</p>
<p>I started with quite a high level of refinement, and the maximum level of refinement decreases over time (from 10 layers to 5).<br>
</p>
<p><br>
</p>
<p>I've found that ripples appear in the inner portion of the domain, as well as larger numerical artifacts.<br>
</p>
<p>I think that these might either be due to:<br>
</p>
<p><br>
</p>
<p>- The initial temperature fluctuations I added, though I feel these would be too small (spatially) to cause all these issues.<br>
</p>
<p>- Could it be that the higher pressure region in the center, combined with the lower radial velocity, cause the gas to bounce off the inner boundary?<br>
</p>
<p>This might create ripples and internal shocks.<br>
</p>
<p><br>
</p>
<p>Either way, I'm not certain how to improve my results from here.<br>
</p>
<p>If anyone has had similar problems, I would greatly appreciate some advice.<br>
</p>
<p><br>
</p>
<p>The figures attached highlight my problem.<br>
</p>
<p><br>
</p>
<p>Kind Regards,<br>
</p>
<p>Michiel Bustraan<br>
</p>
</div>

</blockquote></div><br></div></div>