[CIG-LONG] CIG-LONG Digest, Vol 42, Issue 9

[Walter Landry]

Karen Paczkowski <karen.paczkowski at yale.edu> wrote:
> Hi Gilly,
> 
> Thank you for your response.  I tried adding in the hydrostatic  
> pressure but now I am getting very high magnitude pressure values on  
> each of the 4 corners of my simulation box.
> 
> With the addition of the hydrostatic stress I now have two stress  
> boundary conditions, as shown below.  The bottom stress boundary  
> condition is equal to the weight of the material in the box.  Without  
> the hydrostatic term I do get a linear gradient, but it doesn't go to  
> zero at the top of the simulation box like I would have thought.  I  
> would think that this bottom stress boundary condition may need to be  
> altered with the addition of the hydrostatic term, but i'm not sure  
> how.  If I set the bottom boundary condition to zero then I get an  
> error and the code won't run.  Ideally I would like to set pressure  
> boundary conditions along the top, bottom and sides, so I'm hoping  
> there is a way to do this without specifying velocity boundary  
> conditions.
> 
> Do you have any idea what might be going wrong?

Can I see your entire input file?  I do not see a side StressBC.  I
have not done any simulations with StressBC's on the side, so I do not
know if they will work well.

Cheers,
Walter Landry
walter at geodynamics.org