<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hello everyone,<div><br></div><div>  I've been noticing an issue in blocks that refine on the boundaries, and so I've reduced it to a simple setup that seems to recreate the issue. </div><div><br></div><div>In this simple set up I have a rectangular volume with nxb=32, nyb=16, nzb=8, so that all dimensions have equal resolution. In each block I have a density gradient in the x-direction from 2e-6 g/cc in the first cell (including guard cells), to 2e-4 g/cc in the last cell (including guard cells). I also have an x-velocity of +1e4 cm/s. I am using diode boundary conditions. I have an initial uniform pressure of 2e-2 dyne/cm^2. </div><div><br></div><div>At init, I force everything to be a level lmin=3=lmax-1. After nref=2 steps I force the entire volume to refine to level lmax=4. This results in a strip of high energy (or pressure) cells on the inner-x boundary that are approximately 3 orders of magnitude higher than the uniform value (see attached). </div><div><br></div><div>[Aside: This also leads to an increasing in the integrated internal energy of the box by more than a factor of 20. While the density appears to remain the same (but smoother) it does look like it shifts slightly 'downhill', so when density is increasing from left to right, a refinement is activated such that a band of same colour (say 3e-5 g/cc) shifts to lower x. This may be an optical illusion though.]</div><div><br></div><div> These refinements are done at the end of step=2, and while the high pressure strip is present in the images I make of the checkpoint file at step=2, they are not present in the images I make of the plotfile at step=2, as the plotfile is generated before the refinements are done. This has convinced me it is the refinements causing the incorrect pressure. </div><div><br></div><div>I've done a few tests so far:</div><div>   If there is a constant density, and no velocity, then the problem doesn't appear.</div><div>   If there is a density gradient, but no velocity, then the problem doesn't appear.</div><div>   If there is a constant density, and a positive x velocity, then it does appear. </div><div>   This is true if the gradients and velocity are in the y direction, but the pressure issue is on the inner y boundary.</div><div>   If the density gradient is negative, and the x velocity positive, the problem doesn't seem to occur.</div><div>   If the density gradient is negative, and the x velocity negative, the problem appears as before, but on the upper x boundary. </div><div><br></div><div>Attached is a tarball of the setup files found in their respective directory in SimulationMain, including a flash.par and a setup_call that was used. I've also included my sites file so you can see my optimisations/options. Currently it's -O3, but I have tested this with -debug on, with the same results. I am currently using ifort/icc/icpc (IFORT/ICC/) 15.0.2 20150121.</div><div><br></div><div>Any help in understanding why this is happening would be very useful! Is this the expected result for a prolongation on a simulation boundary with diode condition and motion away from the boundary? Or could there be a bug in gr_updateRefinement? </div><span class="HOEnZb"><font color="#888888"><div><br></div><div>-Mark<br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><br>Mark Richardson<br><div>Beecroft Postdoctoral Fellow, Department of Physics<div>Denys Wilkinson Building - 555F</div><div>University of Oxford<div><a href="mailto:Mark.Richardson@physics.ox.ac.uk" target="_blank">Mark.Richardson@physics.ox.ac.uk</a><br><div>+44 1865 280793</div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
</div></font></span></div></blockquote></div></div><font color="#888888">
</font>
<br><br>
</div>