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:45] gag [SIF contents] |
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 69: | 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 80: | 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 107: | Line 114: | ||
==== Examples ==== | ==== Examples ==== | ||
- | Download an example using the Porous Solver. | + | An example using the //Porous// Solver |
==== Reference ==== | ==== Reference ==== |