Differences
This shows you the differences between two versions of the page.
| Both sides previous revision Previous revision Next revision | Previous revision | ||
|
solvers:porous [2012/11/14 19:41] gag [General Description] |
solvers:porous [2015/12/02 12:49] (current) gag [SIF contents] |
||
|---|---|---|---|
| Line 1: | Line 1: | ||
| ~~NOTOC~~ | ~~NOTOC~~ | ||
| - | ===== Solver Porous Solver ===== | + | ===== Snow/firn rheology - Solver Porous Solver ===== |
| ==== General Informations ==== | ==== General Informations ==== | ||
| - | * **Solver Fortran File: | + | * **Solver Fortran File: |
| - | * **Solver Name: | + | * **Solver Name: |
| - | * **Required Output Variable(s): | + | * **Required Output Variable(s): |
| - | * **Required Input Variable(s): | + | * **Required Input Variable(s): |
| - | * **Optional Output Variable(s): | + | * **Optional Output Variable(s): |
| * **Optional Input Variable(s): | * **Optional Input Variable(s): | ||
| Line 13: | Line 13: | ||
| ==== General Description ==== | ==== General Description ==== | ||
| This solver computes the flow of snow/firn material (i.e., porous incompressible ice) using the snow/firn law proposed by Gagliardini and Meyssonnier (1997). The snow/ice rheological law is function of the relative density, which is a required input variable for this solver. The law depends on two function, a(D) and b(D), which are parametrized functions of the relative density D. | This solver computes the flow of snow/firn material (i.e., porous incompressible ice) using the snow/firn law proposed by Gagliardini and Meyssonnier (1997). The snow/ice rheological law is function of the relative density, which is a required input variable for this solver. The law depends on two function, a(D) and b(D), which are parametrized functions of the relative density D. | ||
| - | More details about the snow/firn law can be found {{: | + | |
| + | As optional output variables, strain-rate, | ||
| + | |||
| + | More details about the snow/firn law can be found here{{:solvers: | ||
| Line 22: | Line 25: | ||
| $yearinsec = 365.25*24*60*60 | $yearinsec = 365.25*24*60*60 | ||
| $rhoi = 900.0/ | $rhoi = 900.0/ | ||
| - | $B = 20.0 ! MPa-3a-1 T = -2?C | + | ! B = 2 A, where A is the classical Glen's fluidity |
| + | $B = 20.0 ! MPa{^-3}a{^-1} T = -10°C | ||
| $n = 3.0 | $n = 3.0 | ||
| $gravity = -9.81*yearinsec^2 | $gravity = -9.81*yearinsec^2 | ||
| Line 35: | Line 39: | ||
| Solver 1 | Solver 1 | ||
| Equation = String " | Equation = String " | ||
| - | Procedure = "./ | + | Procedure = "ElmerIceSolvers" " |
| Variable = " | Variable = " | ||
| Variable DOFs = 4 ! 4 in 3D (u,v,w,p) ; 3 in 2D (u,v,p) | Variable DOFs = 4 ! 4 in 3D (u,v,w,p) ; 3 in 2D (u,v,p) | ||
| Line 49: | Line 53: | ||
| ! switch that in for post-processing issues only | ! switch that in for post-processing issues only | ||
| | | ||
| - | | + | |
| | | ||
| - | | + | |
| + | | ||
| + | | ||
| Linear System Solver = ' | Linear System Solver = ' | ||
| Line 67: | Line 73: | ||
| End | End | ||
| - | ! Gravity force | + | ! Gravity force is directly the ice density time the gravity |
| + | ! It is further multiplied by the relative density in the Porous solver | ||
| Body Force 1 | Body Force 1 | ||
| Porous Force 1 = Real 0.0E00 | Porous Force 1 = Real 0.0E00 | ||
| Line 78: | Line 85: | ||
| Min Second Invariant = Real 1.0E-10 | Min Second Invariant = Real 1.0E-10 | ||
| Fluidity Parameter = Real $B ! MPa^{-3}a^{-1} | Fluidity Parameter = Real $B ! MPa^{-3}a^{-1} | ||
| + | |||
| + | ! Just for output purpose, not needed by the Porous solver | ||
| ! Density as a function of relative density | ! Density as a function of relative density | ||
| Density = Variable Relative Density | Density = Variable Relative Density | ||
| Line 105: | Line 114: | ||
| ==== Examples ==== | ==== Examples ==== | ||
| - | Download an example using the Porous Solver. | + | An example using the //Porous// Solver |
| ==== Reference ==== | ==== Reference ==== | ||
| Line 111: | Line 120: | ||
| Gagliardini O. and J. Meyssonnier, | Gagliardini O. and J. Meyssonnier, | ||
| - | Its implementation within Elmer/Ice and an application are presented in this reference: | + | Its implementation within Elmer/Ice and an application are presented in this reference:\\ |
| Zwinger T. , R. Greve, O. Gagliardini , T. Shiraiwa and M. Lyly, 2007. A full Stokes-flow thermo-mechanical model for firn and ice applied to the Gorshkov crater glacier, Kamchatka. Annals of Glaciol., 45, p. 29-37. | Zwinger T. , R. Greve, O. Gagliardini , T. Shiraiwa and M. Lyly, 2007. A full Stokes-flow thermo-mechanical model for firn and ice applied to the Gorshkov crater glacier, Kamchatka. Annals of Glaciol., 45, p. 29-37. | ||
