Model help:GIPL: Difference between revisions
No edit summary |
No edit summary |
||
Line 61: | Line 61: | ||
|width="20%"|<span> Max distance between next temperature | |width="20%"|<span> Max distance between next temperature | ||
iterations (E1) (0.001-0.1)</span> | iterations (E1) (0.001-0.1)</span> | ||
|width="60%"|<span>If the absolute difference between current and next temperature iteration is greater than E1 | |width="60%"|<span>If the absolute difference between current and next temperature iteration is greater than E1 then apply sweep method (see Ismayil-zadeh and Tackley 2010) over again. Used in Stefan subroutine.</span> | ||
|width="20%"|<span>-</span> | |width="20%"|<span>-</span> | ||
Line 67: | Line 67: | ||
|-valign="top" | |-valign="top" | ||
|width="20%"|<span> Error of saturation for next iteration (UWK)</span> | |width="20%"|<span> Error of saturation for next iteration (UWK)</span> | ||
|width="60%"|<span>If the absolute difference between current and next saturated unfrozen water coefficient is greater than UWK | |width="60%"|<span>If the absolute difference between current and next saturated unfrozen water coefficient is greater than UWK then apply sweep method over again. Used in Stefan subroutine.</span> | ||
|width="20%"|<span>-</span> | |width="20%"|<span>-</span> | ||
Line 229: | Line 229: | ||
==References== | ==References== | ||
<span | Ismail-Zadeh A, Tackley P. 2010. Computational Methods for Geodynamics. Cambridge. | ||
<span>Key papers</span> | |||
==Links== | ==Links== |
Revision as of 16:44, 17 August 2011
GIPL
GIPL is a heat flow with phase change model use to model permafrost thermal state and active layer depth, can be used for site specific and regional ground temperature distribution modeling.
Model introduction
GIPL(Geophysical Institute Permafrost Laboratory) is an implicit finite difference transient one-dimensional heat flow model. The model simulates ground temperature dynamics and the depth of the active layer by solving non-linear heat equation with phase change numerically. The model employs the Enthalpy method which does not require explicit treatment of the freeze/thaw moving boundary. In this model the process of freezing or thawing is occurring in accordance with unfrozen water content and soil thermal properties, and depends on the degree of soil saturation.
Model parameters
Uses ports
This will be something that the CSDMS facility will add
Provides ports
This will be something that the CSDMS facility will add
Main equations
1) Heat Flow Equation
[math]\displaystyle{ \frac{\partial H(x,t)}{\partial \tau}=\texttt{div}(k(x,t)\nabla t(x,\tau)) }[/math] (1)
[math]\displaystyle{ H(x,t)=\int\limits_0^tC(x,s)ds+L\Theta(x,t) }[/math] (2)
[math]\displaystyle{ \frac{\partial t(l_2,\tau)}{\partial x}=g }[/math] (3)
[math]\displaystyle{ t(x,0)=t_0(x) }[/math] (4)
[math]\displaystyle{ \Theta(x,t)=\eta(x)\cdot\begin{cases} 1 , & t\ge t_* \\ a|t|^{-b}, & t\lt t_* \end{cases} }[/math] (5)
Notes
Any notes, comments, you want to share with the user
Numerical scheme
Examples
An example run with input parameters, BLD files, as well as a figure / movie of the output
Follow the next steps to include images / movies of simulations:
- Upload file: http://csdms.colorado.edu/wiki/Special:Upload
- Create link to the file on your page: [[Image:<file name>]].
See also: Help:Images or Help:Movies
Developer(s)
Name of the module developer(s)
References
Ismail-Zadeh A, Tackley P. 2010. Computational Methods for Geodynamics. Cambridge. Key papers
Links
Any link, eg. to the model questionnaire, etc.