[aspect-devel] Aspect melt_fraction

Magali Billen mibillen at ucdavis.edu
Fri Mar 14 13:00:13 PDT 2014

I think it would be great to have a module for calculating melt fraction in the main branch.
This doesn't come up for every problem, but there's often a problem that you're working 
on and then part way through you realize that melt fraction might be helpful. If we are looking
ahead to more realistic theologic models, for example, being able to have rheology depend
on melt fraction would require calculating melt fraction first.

On Mar 14, 2014, at 12:04 PM, Jonathan Perry-Houts wrote:

> Hash: SHA1
> Hi all,
> I'm starting work on melt migration in ASPECT, and the first step was
> to make melt fraction a more easily tracked quantity.
> In December Juliane made a nice postprocessor to visualize the amount
> of melt generated in the model based on Katz, 2003 paper. The
> algorithm is great but I don't like the fact that it's hard-coded in
> to the postprocessor. I think it would be more useful to make melt
> fraction a quantity calculated by the material model and have the
> postprocessor display it in the same way viscosity and density are
> displayed.
> This also allows material properties to be easily derived from the
> melt fraction without calculating the same thing twice for each
> timestep. If anyone has a better suggestion for dealing with melt I'd
> be happy to hear it.
> I'm not sure that other people would be interested in these changes
> enough to put them in the main branch, but I figured I'd share them
> with the mailing list in case anyone does find them useful. The
> changes are relatively unobtrusive, except for breaking Juliane's
> melt_fraction visualization postprocessor, but that functionality can
> be easily replicated by including the same algorithm in the material
> model.
> I can also send out an example material model which calculates melt
> based on Juliane's algorithm if people are interested.
> The modified files are:
> include/aspect/material_model/interface.h
> include/aspect/postprocess/visualization/melt_fraction.h
> source/material_model/interface.cc
> source/postprocess/visualization/melt_fraction.cc
> Version: GnuPG v1.4.14 (GNU/Linux)
> Comment: Using GnuPG with Thunderbird - http://www.enigmail.net/
> iQEcBAEBAgAGBQJTI1KfAAoJEGe6xJ1FYRpR6ykH/3aGtp1GotuXjxbYlGdu+LhF
> 17Pi3z4DBfGW3s2P+wx7W05nsLHmad3cSFtLbwmdpPq1gZJ41nrweOordcFy8W/v
> kVQryawlBW+29LV5aqy8k5T4eB0FizzssU5yTjc/0jOJQeANQkmdeB1yFTLVlkiX
> 8f0i5owJb0ZeuDp9FElQPOQaNTJf5U/GKdCXk5oENz6w2yEkaGmloMLNnf/1uNK9
> lT+yLp8DL74n2p9sBiMlZUD2hyHrdhu0fQ7p32OyoseZ5aLlvRGcGMD81lUfWT8=
> =c9+n
> <interface.h><melt_fraction.h><interface.cc><melt_fraction.cc>_______________________________________________
> Aspect-devel mailing list
> Aspect-devel at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel

Professor of Geophysics & UCD Chancellor Fellow
Chair, Geology Graduate Program
Earth & Planetary Sciences Dept., UC Davis
Davis, CA 95616
2129 Earth & Physical Sciences Bldg.
Office Phone: (530) 752-4169

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20140314/b3755a57/attachment.html>

More information about the Aspect-devel mailing list