[aspect-devel] melt in the visco plastic model

Juliane Dannberg dannberg at gfz-potsdam.de
Wed Nov 23 12:19:01 PST 2016

Hi Mohamed,

I would also like to add that combining melt migration and plasticity 
will make the problem strongly nonlinear, and probably very difficult 
for both the linear and the nonlinear solver. So if your viscosity 
contrast is high, you should expect a large number of linear iterations 
(this could be up to 1000) and a large number of nonlinear iterations.

We are working to make both parts faster, but at the moment you will 
probably need a lot of computational resources to run these kind of 
models, and we haven't tested ASPECT with this kind of models, so other 
problems might come up.


On 11/23/2016 01:15 PM, John Naliboff wrote:
> Hi Mohamed,
> With some work it is certainly possible to combine the functionality 
> of the visco plastic material model and 'melting' components of the code.
> However, do note that the visco plastic model is incompressible and 
> the melting material models are inherently compressible.  My 
> suggestion is to use the 'melt simple' material model as a starting 
> point and then add in the features in 'visco plastic'.
> However, this seems like a good time to release a general cautionary 
> note about moving to quickly on this front.  I am still in the process 
> of doing sensitivity tests with the provided extension models and 
> adding in melting is likely going to require a huge amount of 
> additional testing to obtain robust convergence behavior and results.
> I certainly encourage you to move forward with the idea, but just keep 
> in mind that this is not a trivial numerical exercise in any way, 
> shape or form.
> Cheers,
> John
> *************************************************
> John Naliboff
> Assistant Project Scientist, CIG
> Earth & Planetary Sciences Dept., UC Davis
> On 11/23/2016 08:19 AM, Mohamed Gouiza wrote:
>> Hi All,
>> Is there a way to combine the visco plastic material model with the 
>> melt simple material model?
>> I want to set up a model that simulate lithospheric extension and 
>> induce melt generation and transport at various stages of the extension.
>> Thank you
>> -------------------------------------------------
>> Mohamed Gouiza, Research Fellow
>> Basin Structure Group, Institute of Applied Geosciences
>> University of Leeds, School of Earth and Environment
>> Leeds,  LS2 9JT, UK
>> M.Gouiza at leeds.ac.uk <mailto:M.Gouiza at leeds.ac.uk>
>> +44 7985 782073
>> -------------------------------------------------
>> _______________________________________________
>> Aspect-devel mailing list
>> Aspect-devel at geodynamics.org
>> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
> _______________________________________________
> Aspect-devel mailing list
> Aspect-devel at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel

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

More information about the Aspect-devel mailing list