[aspect-devel] free surface with non-zero surface pressure

John Naliboff jbnaliboff at ucdavis.edu
Fri Nov 18 13:46:18 PST 2016

Hi Carolyn,

Great question!

Currently, it is not possible to mix a free surface and traction 
boundary condition.  However, as you pointed out this feature would be 
very useful for a number of applications and I do not see any reason why 
it can't be implemented.

I think this kind of feature would likely be implemented as an added 
option to the free surface portions of the code or as a new boundary 
condition model/option.  Others will have a much better idea of exactly 
how/where this should be done.

I imagine everyone is pretty busy gearing up for AGU, but this is a 
point one or more of us could return to next year.  Alternatively, If 
you would like to try the implementation please do so! You can always 
send questions to the mailing list as you go along.

Hope this helps and if you plan to attend AGU keep in mind there will be 
a whole suite of APSECT-related presentations 


John Naliboff
Assistant Project Scientist, CIG
Earth & Planetary Sciences Dept., UC Davis

On 11/18/2016 11:29 AM, Carolyn Branecky wrote:
> Hi all,
> I'm wondering whether is any reason why the free surface option could 
> not be run with a traction boundary condition at the free surface.  
> Here, I have in mind applying a normal force at the surface, which 
> could be due to an ice sheet, sediment load, etc., and then the 
> surface would be allowed to deform under that load.
> This could be framed as two questions:
> 1. Can surface pressure be combined with a free surface in principle 
> in ASPECT?
> 2. Are there any barriers to implementing these together in either the 
> 1.4 release or the development version?
> Looking forward to hearing your thoughts.
> -Carolyn
> -- 
> Carolyn Branecky
> Ph.D. candidate
> University of California, Santa Cruz
> Earth and Planetary Sciences
> cbranecky at ucsc.edu <mailto:cbranecky at ucsc.edu>
> _______________________________________________
> 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/20161118/ba1f947e/attachment.html>

More information about the Aspect-devel mailing list