Difference between revisions of "Solving the equation for thickness evolution"

From Interactive System for Ice sheet Simulation
Jump to: navigation, search
(A generic conservation equation)
(A generic conservation equation)
Line 7: Line 7:
  
  
On the left-hand side we have the change in ''&phi;'' per unit time, which is given by the negative of the divergence of the flux, ''Q'' (a vector), plus any contribution from sources, <math>S_{\varphi }</math>. The negative sign in front of the divergence operator (the <math>\nabla</math>) insures that, for example in the case of heat flow, a negative temperature gradient causes an ''increase'' in temperature over time in the absence of sources (intuitively, we expect heat to flow from warmer to colder regions rather than vice versa). This is more obvious if we expand the equation above in terms of partial derivatives and the vector-valued flux components,  
+
On the left-hand side we have the change in ''&phi;'' per unit time, which is given by the negative of the divergence of the flux, ''Q'' (a vector), plus any contribution from sources, <math>S_{\varphi }</math>. The negative sign in front of the divergence operator (the <math>\nabla</math>) insures that, for example in the case of heat flow, a negative temperature gradient causes an ''increase'' in temperature over time in the absence of source. Intuitively, we expect heat to flow from warmer to colder regions rather than vice versa. This is, perhaps, more obvious if we expand the equation above in terms of partial derivatives and the vector-valued flux components,  
  
  
 
<math>\frac{\partial \varphi }{\partial t}=-\left( \frac{\partial Q_{x}}{\partial x}+\frac{\partial Q_{y}}{\partial y} \right)+S_{\varphi }</math>.
 
<math>\frac{\partial \varphi }{\partial t}=-\left( \frac{\partial Q_{x}}{\partial x}+\frac{\partial Q_{y}}{\partial y} \right)+S_{\varphi }</math>.
 
  
 
==Conservation of mass for an ice sheet==
 
==Conservation of mass for an ice sheet==

Revision as of 11:36, 5 August 2009

A generic conservation equation

Evolution of the ice sheet geometry is simply the result of mass conservation. Consider a generic conservation equation for some scalar value φ


\frac{\partial \varphi }{\partial t}=-\nabla \cdot \vec{Q}+S_{\varphi }


On the left-hand side we have the change in φ per unit time, which is given by the negative of the divergence of the flux, Q (a vector), plus any contribution from sources, S_{\varphi }. The negative sign in front of the divergence operator (the \nabla) insures that, for example in the case of heat flow, a negative temperature gradient causes an increase in temperature over time in the absence of source. Intuitively, we expect heat to flow from warmer to colder regions rather than vice versa. This is, perhaps, more obvious if we expand the equation above in terms of partial derivatives and the vector-valued flux components,


\frac{\partial \varphi }{\partial t}=-\left( \frac{\partial Q_{x}}{\partial x}+\frac{\partial Q_{y}}{\partial y} \right)+S_{\varphi }.

Conservation of mass for an ice sheet

For an ice sheet (or any depth-integrated mass flow for that matter), the equation becomes


\begin{align}
  & \varphi =H,\quad \vec{Q}=\vec{U}H,\quad \vec{U}=\left( U,V \right),\quad S_{H}=\dot{b} \\ 
 &  \\ 
 & \frac{\partial H}{\partial t}=-\left( \frac{\partial \left( UH \right)}{\partial x}+\frac{\partial \left( VH \right)}{\partial y} \right)+\dot{b} \\ 
\end{align}