<div dir="ltr"><div><div><div><div><div><div><div><div>Hi Michiel,<br><br></div>If you suspect that your temperature fluctuations are part of the problem, try removing them. You mentioned that you need high resolution; in order to get full refinement, you can set the minimum refinement level equal to the maximum refinement level:<br><br>lrefine_min = 10<br></div>lrefine_max = 10<br><br></div>Of course, this effectively turns off the AMR feature, so you may have to add processors in order to get the resolution you need across the entire domain. You might also try reducing your domain dimensions for this test run, in order to reduce the cluster resource requirements. <br><br></div>Also, you mention that the gas might be bouncing off the inner boundary: are your boundary conditions reflective, then? If not, there will be no bouncing; if so, turn the boundary conditions to 'outflow' and see what happens.<br><br></div>Cheers,<br></div>Jason<br><br>------<br></div>Jason Galyardt <br></div>University of Georgia<br><div><div><div><div><div><div><div><div><br></div></div></div></div></div></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Mar 22, 2017 at 9: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>