[aspect-devel] Compositional fields with different densities and viscosities

[Glerum, A.C. (Anne)]

Hi all,

Thank you for your suggestions. I'll have a go at producing a new postprocessor by adapting the velocity statistics one.

I did run into a problem with a model considering 8 compositional fields; it got stopped because it was exceeding the allowed virtual and physical memory usage.
I restarted the same run, but with a CLF number of 1.0e-6. Hence, the number of elements, degrees of freedom, global mass etc is constant. However, the memory usage keeps increasing (see snapshot).
A similar increase is seen for a high resolution, uniform and constant grid Rayleigh-Taylor instability using 1 compositional field and material model simple.cc, only then the usage doesn't exceed the limit.
Do you have an idea about the cause of this memory runaway?

Thank you,

Anne

[cid:753E56BE-4AD5-4D60-A4FE-583EA32227C1 at geo.uu.nl]





Anne Glerum, MSc | PhD candidate | Department of Earth Sciences | Utrecht University|  Budapestlaan 4, 3584 CD Utrecht | Room Z.204 | A.C.Glerum at uu.nl<mailto:D.J.J.vanHinsbergen at uu.nl> |



On Jan 27, 2013, at 11:38 PM, Timo Heister wrote:

For
example, I would be interested in the Vrms and dissipation of energy of
only the downgoing slab (where the slab is represented by one
compositional field). This would greatly help with the benchmark model I
am working on.

Yes, in a case like this I would look at how the velocity statistics
postprocessor is written. What you want to do is clone it to compute
what you describe. Essentially, what you'll have to do is only add to
the integral and area whenever you find that a particular quadrature
point satisfies the criteria by which you define whether it is part of
the slab.

To extend on that, if your compositional field is 1 for "interesting"
and 0 for "not interesting", you can just multiply the integral by the
value of the compositional field at that quadrature point.

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

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://geodynamics.org/pipermail/aspect-devel/attachments/20130129/db60e639/attachment-0001.htm 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screen Shot 2013-01-29 at 2.55.37 PM.png
Type: image/png
Size: 45384 bytes
Desc: Screen Shot 2013-01-29 at 2.55.37 PM.png
Url : http://geodynamics.org/pipermail/aspect-devel/attachments/20130129/db60e639/attachment-0001.png