This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
solvers:sia [2012/11/19 21:31] gag [Examples] |
solvers:sia [2014/01/31 09:28] (current) ltavard [Examples] |
||
---|---|---|---|
Line 3: | Line 3: | ||
==== 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 15: | Line 15: | ||
For a flow line problem, the mesh is a plane surface, and a volume for a 3D problem. Regarding this aspect, this solver is certainly not as efficient as a classical SIA solver. But, on the other hand, it works for unstructured grid and non-constant viscosity. The SIA velocities and pressure can be use, for example, as initial conditions for the Stokes Solver. Contrary to the NS solver, the gravity must be orientated along the z-axis. | For a flow line problem, the mesh is a plane surface, and a volume for a 3D problem. Regarding this aspect, this solver is certainly not as efficient as a classical SIA solver. But, on the other hand, it works for unstructured grid and non-constant viscosity. The SIA velocities and pressure can be use, for example, as initial conditions for the Stokes Solver. Contrary to the NS solver, the gravity must be orientated along the z-axis. | ||
- | The SIA solver uses the same input parameters as the NS solver (Viscosity, Density, Viscosity Exponent, Flow BodyForce, | + | The SIA solver uses the same input parameters as the NS solver (Viscosity, Density, Viscosity Exponent, Flow BodyForce, |
The basal velocities are given as Dirichlet BC on the bedrock surface. The [[: | The basal velocities are given as Dirichlet BC on the bedrock surface. The [[: | ||
Line 28: | Line 28: | ||
Solver 2 | Solver 2 | ||
Equation = "SIA Variable" | Equation = "SIA Variable" | ||
- | Procedure = File "SIASolver" " | + | Procedure = File "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 35: | Line 35: | ||
Solver 3 | Solver 3 | ||
Equation = " | Equation = " | ||
- | Procedure = File "SIASolver" " | + | Procedure = File "ElmerIceSolvers" " |
Variable = -nooutput " | Variable = -nooutput " | ||
Variable DOFs = 1 | Variable DOFs = 1 | ||
Line 65: | Line 65: | ||
==== Examples ==== | ==== Examples ==== | ||
- | Download {{: | + | An example using the '' |
- | ==== Reference ==== | ||
- | When used this solver can be cited using the following references: | ||
- | Gagliardini O., D. Cohen, P. Råback and T. Zwinger, 2007. Finite-Element Modeling of Subglacial Cavities and Related Friction Law. J. of Geophys. Res., Earth Surface, 112, F02027. |