Difference between revisions of "Stokes equations"

From Interactive System for Ice sheet Simulation
Jump to: navigation, search
(Stokes Equations)
Line 19: Line 19:
 
where ''<big>&sigma;</big>'' is the full stress tensor and ''<big><math>\delta _{ij}</math></big>'' is the Kronecker delta (or identity matrix).
 
where ''<big>&sigma;</big>'' is the full stress tensor and ''<big><math>\delta _{ij}</math></big>'' is the Kronecker delta (or identity matrix).
  
The "nonlinearity" is not obvious here, but comes in when we invoke a constitutive law that links stresses and strain rates. A standard Newtonian (constant viscosity) constitutive law for ice looks like
+
The "nonlinearity" is not obvious here, but becomes clear when we invoke a constitutive law that links stresses and strain rates. A standard Newtonian (constant viscosity) constitutive law for ice looks like
  
  
Line 26: Line 26:
  
  
where ''&eta;'' is the viscosity. We can write the constitutive law the same way, but define the viscosity instead as an "effective" viscosity such that
+
where ''<big>&eta;</big>'' is the viscosity. We can write the constitutive law the same way, but define the viscosity instead as an "effective" viscosity such that
  
  
 
<math>\begin{align}\quad \eta \equiv \frac{1}{2}B\dot{\varepsilon }_{e}^{\frac{1-n}{n}} \\\end{align}</math>
 
<math>\begin{align}\quad \eta \equiv \frac{1}{2}B\dot{\varepsilon }_{e}^{\frac{1-n}{n}} \\\end{align}</math>
 +
 +
 +
where ''<big>B</big>'' is the temperature dependent pre-factor found in Glen's (inverse) flow law and <math>\begin{align}\quad\dot{\varepsilon }_{e}^{\frac{1-n}{n}} \\\end{align}</math> is the "effective strain rate", a norm of the strain-rate tensor. With this definition, as the strain rate at any point increases, the effective viscosity of the ice decreases (i.e. the ice becomes "softer" and deforms more easily). The stress depends on the strain rate, but in this case the coefficient that links the two also depends on the strain rate (hence the nonlinearity in the equations).

Revision as of 12:36, 3 August 2009

Stokes Equations

The conservation equations that describe the balance of momentum for an ice sheet (assuming that accelerations are negligible) are the non-linear Stokes equations. In Cartesian coordinates, these are given by,


\begin{align}
  & x:\quad \frac{\partial \tau _{xx}}{\partial x}-\frac{\partial P}{\partial x}+\frac{\partial \tau _{xy}}{\partial y}+\frac{\partial \tau _{xz}}{\partial z}=0 \\ 
 & y:\quad \frac{\partial \tau _{yy}}{\partial y}-\frac{\partial P}{\partial y}+\frac{\partial \tau _{xy}}{\partial x}+\frac{\partial \tau _{yz}}{\partial z}=0 \\ 
 & z:\quad \frac{\partial \tau _{zz}}{\partial z}-\frac{\partial P}{\partial z}+\frac{\partial \tau _{zy}}{\partial y}+\frac{\partial \tau _{xz}}{\partial x}=\rho g \\ 
\end{align},


where P is the pressure and τ is the deviatoric stress tensor. The latter is given by


\quad\tau _{ij}=\sigma _{ij}+P\delta _{ij},


where σ is the full stress tensor and \delta _{ij} is the Kronecker delta (or identity matrix).

The "nonlinearity" is not obvious here, but becomes clear when we invoke a constitutive law that links stresses and strain rates. A standard Newtonian (constant viscosity) constitutive law for ice looks like


\begin{align}\quad \tau _{ij}=2\eta \dot{\varepsilon }_{ij} \\ 
\end{align}


where η is the viscosity. We can write the constitutive law the same way, but define the viscosity instead as an "effective" viscosity such that


\begin{align}\quad \eta \equiv \frac{1}{2}B\dot{\varepsilon }_{e}^{\frac{1-n}{n}} \\\end{align}


where B is the temperature dependent pre-factor found in Glen's (inverse) flow law and \begin{align}\quad\dot{\varepsilon }_{e}^{\frac{1-n}{n}} \\\end{align} is the "effective strain rate", a norm of the strain-rate tensor. With this definition, as the strain rate at any point increases, the effective viscosity of the ice decreases (i.e. the ice becomes "softer" and deforms more easily). The stress depends on the strain rate, but in this case the coefficient that links the two also depends on the strain rate (hence the nonlinearity in the equations).