<div dir="ltr"><div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><div class="gmail_default" style="font-size:13px">Hi Cedric and all,</div><div class="gmail_default" style="font-size:13px"><br></div>
<div class="gmail_default" style="font-size:13px">I have also run into this and related issues.  It has to do with the fact that rigid body modes are in the nullspace of the Stokes operator, and over many timesteps the solutions can pick up power in those modes without any effect on the residuals.  You would also see rigid translations in a box with free slip BCs and periodic boundary conditions.</div>
<div class="gmail_default" style="font-size:13px"><br></div><div class="gmail_default" style="font-size:13px">Note that not only will it mess up the Vrms values, but it can really make your Courant number plummet unnecessarily.</div>
<div class="gmail_default" style="font-size:13px"><br></div><div class="gmail_default" style="font-size:13px">Citcom has the same phenomenon, and they just subtract off rigid rotations at every timestep.  I do have some code that addresses this which is not in the mainline.  I am meeting with Timo in a couple of weeks, perhaps we can add this to a todo list.</div>
<div class="gmail_default" style="font-size:13px"><br></div><div class="gmail_default" style="font-size:13px">Best,</div><div class="gmail_default" style="font-size:13px">Ian</div></div></div><div class="gmail_extra"><br>
<br><div class="gmail_quote">On Thu, Feb 27, 2014 at 6:14 AM, cedric thieulot <span dir="ltr"><<a href="mailto:geogarfield@gmail.com" target="_blank">geogarfield@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi,<br>
<br>
I am running a few simple simulations in cylindrical geometry (constant viscosity, incompressible flow, low Ra number) for a benchmark article.<br>
When I run the simulation with no-slip bc, the solution looks fine and the measurements correspond to values obtained with other codes.<br>
As soon as I use free-slip on the inner and outer bc, I observe a rigid body rotation (see attached movie). Note that the measurements in general are not very much off,<br>
but the root mean square velocity measurements will undoubtedly suffer. Is there an algorithm implemented in the ASPECT to remove it ?<br>
<br>
Cheers,<br>
<br>
Cedric.<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>_______________________________________________<br>
Aspect-devel mailing list<br>
<a href="mailto:Aspect-devel@geodynamics.org">Aspect-devel@geodynamics.org</a><br>
<a href="http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel" target="_blank">http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel</a><br></blockquote></div><br></div>