Warning: Undefined array key 1 in /home/np29546/public_html/elmerice/wiki/inc/auth.php on line 78

Warning: Cannot modify header information - headers already sent by (output started at /home/np29546/public_html/elmerice/wiki/inc/auth.php:78) in /home/np29546/public_html/elmerice/wiki/inc/auth.php on line 431

Warning: Cannot modify header information - headers already sent by (output started at /home/np29546/public_html/elmerice/wiki/inc/auth.php:78) in /home/np29546/public_html/elmerice/wiki/inc/actions.php on line 38
solvers:sia [Elmer/Ice Wiki]

Warning: Undefined array key -1 in /home/np29546/public_html/elmerice/wiki/inc/html.php on line 1458

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

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:** //SIASolver.f90// +  * **Solver Fortran File:** ''SIASolver.f90'' 
-  * **Solver Name:** //SIAVariable// and //SIASolver// +  * **Solver Name:** ''SIAVariable'' and ''SIASolver'' 
-  * **Required Output Variable(s):** //SIAFlow//    +  * **Required Output Variable(s):** ''SIAFlow''    
-  * **Required Input Variable(s):** //Depth////FreeSurfGrad1// and //FreeSurfGrad2//+  * **Required Input Variable(s):** ''Depth''''FreeSurfGrad1'' and ''FreeSurfGrad2''
   * **Optional Output Variable(s):** None     * **Optional Output Variable(s):** None  
   * **Optional Input Variable(s):** None   * **Optional Input Variable(s):** None
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,...). Doesn't work with the build-in Glen's flow law (TODO). Needs to run first the [[solvers:depth|FlowDepth solver]] +The SIA solver uses the same input parameters as the NS solver (Viscosity, Density, Viscosity Exponent, Flow BodyForce,...). Doesn't work with the build-in Glen's flow law (TODO). Needs to run first the [[solvers:depth|FlowDepth solver]].
  
 The basal velocities are given as Dirichlet BC on the bedrock surface. The [[:solvers:ssa|SSA Solver]] can be used to this purpose.  The basal velocities are given as Dirichlet BC on the bedrock surface. The [[:solvers:ssa|SSA Solver]] can be used to this purpose. 
Line 28: Line 28:
 Solver 2 Solver 2
   Equation = "SIA Variable"   Equation = "SIA Variable"
-  Procedure = File "SIASolver" "SIAVariable"+  Procedure = File "ElmerIceSolvers" "SIAVariable"
   Variable = "SIAFlow"   Variable = "SIAFlow"
   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 = "SIA"   Equation = "SIA"
-  Procedure = File "SIASolver" "SIASolver"+  Procedure = File "ElmerIceSolvers" "SIASolver"
   Variable = -nooutput "SIAvar"   Variable = -nooutput "SIAvar"
   Variable DOFs = 1   Variable DOFs = 1
Line 65: Line 65:
  
 ==== Examples ==== ==== Examples ====
-Download {{:solvers:test_sia.tar.gz|here}} an example using the //SIASolver// applied to experiment A160 of ISMIP-HOM benchamrks. +An example using the ''SIASolver'' applied to experiment A160 of ISMIP-HOM benchamrks can be found in ''[ELMER_TRUNK]/elmerice/Tests/SIA'' 
  
-==== 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. 
solvers/sia.1353360709.txt.gz · Last modified: 2012/11/19 21:31 by gag
CC Attribution-Share Alike 4.0 International
www.chimeric.de Valid CSS Driven by DokuWiki do yourself a favour and use a real browser - get firefox!! Recent changes RSS feed Valid XHTML 1.0