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

[Wolfgang Bangerth]

Anne,

> On your advice, I've produced my own material models with multiple
> materials. This all works fine, thank you. I was just wondering if you
> will make it possible to output parameters in the statistics file per
> composition (other than the min/max/global composition values)?

Sure. This is currently computed in the CompositionStatistics class in 
source/postprocess/composition_statistics.cc. But you can of course 
write a similar class that computes anything else -- this data is 
computed by just one of the many postprocessor classes we have, and you 
can add your own one that computes whatever piece of data you are 
interested in.


> 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.

Best
  W.


-- 
------------------------------------------------------------------------
Wolfgang Bangerth               email:            bangerth at math.tamu.edu
                                 www: http://www.math.tamu.edu/~bangerth/