Blatter-Pattyn Boundary Conditions

From Interactive System for Ice sheet Simulation
Revision as of 15:24, 3 August 2009 by Sprice (Talk | contribs)

Jump to: navigation, search

We will go through an approximate derivation of the boundary conditions that are implemented with Glimmer/CISM's higher-order scheme. By "approximate" we mean that some of the derivation is guided by physical intuition and what appear to be "reasonable" arguments, rather than through the application of rigorous mathematics. We take comfort in the fact that, in the end, we wind up with the same sets of equations that one ends up with from the more rigorous approach. We will look at the derivation in three parts, (1) the free surface boundary condition, (2) the specified basal traction boundary condition, and (3) lateral boundary conditions.

Stress Free Surface

At the ice surface, a stress-free boundary condition is applied. The traction vector, T, must be continuous at the ice sheet surface and, assuming that atmospheric pressure and surface tension are small, we have

  & T_{i}=-T_{i(boundary)}\approx 0 \\ 
 & T_{i}=\sigma _{ij}n_{j}=\sigma _{i1}n_{1}+\sigma _{i2}n_{2}+\sigma _{i3}n_{3}=0 \\ 

where the ni are the components of the outward facing, unit normal vector in Cartesisan coordinates.

For a function F(x,y,z) = f(x,y) - z = 0, where z = f(x,y) defines the surface, the gradient of F(x,y,z) gives the components of the surface normal vector:

\begin{align}\nabla F(x,y,z)=\left( \frac{\partial f}{\partial x},\frac{\partial f}{\partial y},-1 \right)\frac{1}{a}=n_{i}\\\ 

where a is the magnitude of the surface normal vector given by

a=\sqrt{\left( \frac{\partial f}{\partial x} \right)^{2}+\left( \frac{\partial f}{\partial y} \right)^{2}+1^{2}}