<div dir="ltr"><div><div><div><div>Happy New Year!<br><br></div>  I was part of the AGU Aspect work meeting and wanted to clarify what I'm working on so we don't have redundant work :<br><br></div>1) I've corrected the non-linear solver in Aspect and am testing it against CitcomS using the non-linear rheology from Billen and Hirth, 2007 ("Rheologic controls on slab dynamics"). The nonlinear solver now solves for the temperature ONCE at the beginning of a timestep then iterates the viscosity and velocity fields until the change from the prior iteration is small (like CitcomS).<br>
<br></div>2) I'm building subduction zones from the same paper and will compare them to identical models in CitcomS.<br><br></div> - Katrina <br><div><div><div><br></div></div></div></div><div class="gmail_extra"><br>
<br><div class="gmail_quote">On Tue, Dec 31, 2013 at 6:50 AM, Timo Heister <span dir="ltr"><<a href="mailto:heister@clemson.edu" target="_blank">heister@clemson.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Okay, sometimes it helps to look at the output:<br>
<a href="http://i.imgur.com/0o9ojFa.png" target="_blank">http://i.imgur.com/0o9ojFa.png</a><br>
these are two runs, on the left with 5 global refinements and 6 on the right.<br>
<br>
So, not only does the second peak occur at a different time, it is a<br>
plume that rises from a different location (see t=600).<br>
<br>
Why does this happen? It looks like the plume on the left comes from<br>
the first staircase (on the right) in the initial condition. Now it<br>
makes sense that the solution is dependent on the mesh size. We<br>
obviously do a poor job with a good initial condition for the<br>
compositional field.<br>
<br>
I guess we could solve this by:<br>
- asking the user to only input "smooth" initial conditions<br>
- implementing smoothing/reinitialization for fields<br>
<div class="HOEnZb"><div class="h5"><br>
<br>
<br>
On Tue, Dec 31, 2013 at 5:52 AM, Timo Heister <<a href="mailto:heister@clemson.edu">heister@clemson.edu</a>> wrote:<br>
>> 2. Could the second peak be initiated by an instability that depends<br>
>> on the mesh? All the runs are on a mesh with an even number of<br>
>> elements. It would be easy to test this by changing the geometry model<br>
>> to call subdivided_hyper_rectangle() with an odd number of repitions<br>
>> instead of hyper_rectangle().<br>
><br>
> ./ I tried that out and the time of the 2nd peak is still changing a<br>
> lot with refinement:<br>
> 48^2: t=500<br>
> 96^2: t=800<br>
> 192^2: t=620<br>
> 384^2: t=870<br>
><br>
> ./ A fixed time step makes no difference.<br>
><br>
> ./ Surprisingly, an adaptive computation with 3+3 levels gives the<br>
> same peak as a fixed mesh with 6 refinements, etc.<br>
><br>
> ./ Is this just because of the initial condition? For example the<br>
> amount of mass in the composition at t=0? No, if I only do adaptive<br>
> refinement in the first timestep, the 2nd peak is at a different time<br>
> (similar to the coarser fixed mesh).<br>
><br>
> So, at some time t>0 at some location the size of the cells determines<br>
> the timing of the 2nd peak. I have to think about this some more. Any<br>
> ideas welcome.<br>
><br>
> --<br>
> Timo Heister<br>
> <a href="http://www.math.clemson.edu/~heister/" target="_blank">http://www.math.clemson.edu/~heister/</a><br>
<br>
<br>
<br>
--<br>
Timo Heister<br>
<a href="http://www.math.clemson.edu/~heister/" target="_blank">http://www.math.clemson.edu/~heister/</a><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>
</div></div></blockquote></div><br><br clear="all"><br>-- <br><br><div>Katrina Arredondo<br><span><span>Ph</span>.D. Candidate</span><br>Geodynamics/Subduction Zones<br>University of California, Davis<br>
Department of Geology<br>
One Shield Avenue, Davis CA 95616<br>
<br></div>
</div>