[aspect-devel] Aspect melt_fraction

[Louise Kellogg]

I agree, this seems like a useful addition to the main branch.

On Mar 14, 2014, at 1:00 PM, Magali Billen <mibillen at ucdavis.edu> wrote:

> 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.
> Magali
> 
> On Mar 14, 2014, at 12:04 PM, Jonathan Perry-Houts wrote:
> 
>> -----BEGIN PGP SIGNED MESSAGE-----
>> 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
>> -----BEGIN PGP SIGNATURE-----
>> 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
>> vCEXIQ/aKKXY+VP4Dq2RVwoo5Kw0/60COxCELWJUGUuijX7GoT6rVzjNouFlrXBT
>> 8f0i5owJb0ZeuDp9FElQPOQaNTJf5U/GKdCXk5oENz6w2yEkaGmloMLNnf/1uNK9
>> lT+yLp8DL74n2p9sBiMlZUD2hyHrdhu0fQ7p32OyoseZ5aLlvRGcGMD81lUfWT8=
>> =c9+n
>> -----END PGP SIGNATURE-----
>> <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
> http://mygeologypage.ucdavis.edu/billen/
> --------------------------------------------------
> 
> 
> 
> 
> 
> 
> 
> 
> 
> _______________________________________________
> Aspect-devel mailing list
> Aspect-devel at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel

********************************************
Louise Kellogg
Professor, Dept. of Earth & Planetary Sciences
One Shields Avenue, University of California, Davis, CA 95616
Phone: (530)752-3690
********************************************




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