[aspect-devel] AGU meeting feedback

[Timo Heister]

> 2. Could the second peak be initiated by an instability that depends
> on the mesh? All the runs are on a mesh with an even number of
> elements. It would be easy to test this by changing the geometry model
> to call subdivided_hyper_rectangle() with an odd number of repitions
> instead of hyper_rectangle().

./ I tried that out and the time of the 2nd peak is still changing a
lot with refinement:
48^2: t=500
96^2: t=800
192^2: t=620
384^2: t=870

./ A fixed time step makes no difference.

./ Surprisingly, an adaptive computation with 3+3 levels gives the
same peak as a fixed mesh with 6 refinements, etc.

./ Is this just because of the initial condition? For example the
amount of mass in the composition at t=0? No, if I only do adaptive
refinement in the first timestep, the 2nd peak is at a different time
(similar to the coarser fixed mesh).

So, at some time t>0 at some location the size of the cells determines
the timing of the 2nd peak. I have to think about this some more. Any
ideas welcome.

-- 
Timo Heister
http://www.math.clemson.edu/~heister/