[CIG-LONG] Alleviating pressure singularity in corner-flow mantle-wedge subduction zone model

[Bradford H. Hager]

Instead of imposing a velocity gradient, try allowing a few nodes to be free - that way they can determine the appropriate velocity gradient.

Brad Hager

On Feb 11, 2011, at 3:06 PM, Karen Paczkowski wrote:

> Hi,
> 
> I am trying to write a simple corner-flow mantle-wedge subduction zone  
> model using Gale.  Since, in the corner of the wedge the velocity  
> abruptly changes from v=vslab along the slab to v=0 along the  
> overriding plate a pressure singularity develops in the corner of the  
> wedge.
> 
> I was wondering if anyone else has encountered this issue and what  
> they did to alleviate the pressure singularity in the corner?
> 
> I have tried adding a small shape along the slab-mantle interface  
> right near the corner that has a velocity with a gradient going from  
> v=0 to v=slab over a couple of resolution points.  This mainly seems  
> to move my pressure singularity around, but not elimiate it.
> 
> Any advice people have on alleviating the pressure singularity in a  
> corner-flow subduction model would be greatly appreciated.
> 
> Thank you,
> Karen
> 
> 
> 
> _______________________________________________
> CIG-LONG mailing list
> CIG-LONG at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-long