[aspect-devel] Aspect melt_fraction

[Jonathan Perry-Houts]

-----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-----
-------------- next part --------------
A non-text attachment was scrubbed...
Name: interface.h
Type: text/x-chdr
Size: 37961 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20140314/6f637dbc/attachment-0002.h>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: melt_fraction.h
Type: text/x-chdr
Size: 2380 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20140314/6f637dbc/attachment-0003.h>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: interface.cc
Type: text/x-c++src
Size: 15741 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20140314/6f637dbc/attachment-0002.cc>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: melt_fraction.cc
Type: text/x-c++src
Size: 4069 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20140314/6f637dbc/attachment-0003.cc>