[aspect-devel] Morency-Dorin with free surface

Payman Janbakhsh payman.janbakhsh at mail.utoronto.ca
Fri Oct 21 13:48:58 PDT 2016

Hi John
2D , delamination, geometry aspect ration of 2x1 , free top surface, other boundaries zero velocity.
Linear temperature profile within the upper 100km and lower 100km and a constant T in between.
25km upper crust, 30km lower crust, a weak zone within the lower crust( with a lower activation energy compare to the enclosing lower crust, so can set off the delamination).
Have tried different reasonable densities and activation energies for every layer except mantle.
Also have put Boundary composition model as : Initial composition  . not sure if this is correct.
Running with IMPES or iterated Stokes made no difference
Also velocity polynomial set to 2 and other polynomials set to 1
also have turned off Local conservative discretization



From: Aspect-devel [mailto:aspect-devel-bounces at geodynamics.org] On Behalf Of John Naliboff
Sent: October 21, 2016 3:36 PM
To: aspect-devel at geodynamics.org
Subject: Re: [aspect-devel] Morency-Dorin with free surface

Hi Payman,

The Morency-Dorin material model should not have intrinsic problem with handling a free surface.  The issue is likely related to some other aspect of your model design.  What type of processes are you modeling?



John Naliboff

Assistant Project Scientist, CIG

Earth & Planetary Sciences Dept., UC Davis
On 10/21/2016 12:17 PM, Payman Janbakhsh wrote:
Does Morency-Dorin material model have difficulty with free surface?
My models get stuck at solving Stokes system even at step zero when I switch top boundary to free surface.
I have set CFL to 0.1 and theta to 0.7 even chanh3d max time step to 1000 but still same issue. I have all other boundaries set to no slip boundaries


Sent from my Samsung Galaxy smartphone.


Aspect-devel mailing list

Aspect-devel at geodynamics.org<mailto:Aspect-devel at geodynamics.org>


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

More information about the Aspect-devel mailing list