[CIG-MC] CitComCU bottom permeable boundary condition

[Eh Tan]

Hi Norm,

I don't know any implementation for permeable bottom boundary condition
available. I have tried to implement permeable top boundary condition
before. Basically, the radial component of the velocity boundary
condition becomes dVr/dr=0, instead of Vr=0. The implementation works,
but it is extremely slow for the Stokes solver to converge. If you need
it, I can dig up the implementation for you.


Cheers,

-- 
Eh Tan
Staff Scientist
Computational Infrastructure for Geodynamics
2750 E. Washington Blvd. Suite 210
Pasadena, CA 91107
(626) 395-1693
http://www.geodynamics.org




Norm Sleep wrote:
> Has anyone implemented a permeable bottom boundary condition for any
> of the CitCom codes? I know how to do this in 2-D stream function but
> I have not tried it in 3-D?  That is,fluid exits the domain
> perpendicular to the boundary, comes in at a given temperature, and
> the boundary does no work on the domain of calculations.
>
> If not, is there another public source code that does this?
>
> Moore, W. B., G. Schubert, and P. Tackley (1998) Three-dimensional
> simulations of plume-lithosphere interaction at the Hawaiian swell,/
> Science/,/ 279/, 1008-1011.
>
> thank you
>
> Norm Sleep
> ------------------------------------------------------------------------
>
> _______________________________________________
> CIG-MC mailing list
> CIG-MC at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-mc
>