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]

This is an old revision of the document!


Solver Shallow Ice Approximation (SIA)

General Informations

  • Solver Fortran File: SIASolver.f90
  • Solver Name: SIAVariable and SIASolver
  • Required Output Variable(s): SIAFlow
  • Required Input Variable(s): Depth, FreeSurfGrad1 and FreeSurfGrad2
  • Optional Output Variable(s): None
  • Optional Input Variable(s): None

General Description

This SIA solver is not classical in that sense that the equations are not solved on a grid of dimension lower than the problem dimension itself. The geometry (H, B and S) is here given by the mesh. 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 basal velocities are given as Dirichlet BC on the bedrock surface. The SSA Solver can be used to this purpose.

SIF contents

The required keywords in the SIF file for this solver are:

Solver 1
  Equation = String "StressSolver"
  Procedure =  File "ComputeDevStressNS" "ComputeDevStress"
  ! this is just a dummy, hence no output is needed
  !-----------------------------------------------------------------------
  Variable = -nooutput "Sij"
  Variable DOFs = 1
  ! the name of the variable containing the flow solution (U,V,W,Pressure)
  !-----------------------------------------------------------------------
  Flow Solver Name = String "Flow Solution"
  ! the name of the stress solution (default is 'stress')
  Stress Variable Name = String 'Sigma'
  !-----------------------------------------------------------------------
  Exported Variable 1 = "Sigma" ! [Sxx, Syy, Szz, Sxy] in 2D
                                 ! [Sxx, Syy, Szz, Sxy, Syz, Szx] in 3D
  Exported Variable 1 DOFs = 6   ! 4 in 2D, 6 in 3D
  Linear System Solver = "Iterative"
  Linear System Iterative Method = "BiCGStab"
  Linear System Max Iterations = 300
  Linear System Convergence Tolerance = 1.0E-09
  Linear System Abort Not Converged = True
  Linear System Preconditioning = "ILU0"
  Linear System Residual Output = 1
End

Material 1
  ...
  ! we want to have the Cauchy stress
  !----------------------------------
  Cauchy = Logical True
End

Examples

Download an example using this solver. TODO

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.1353348806.txt.gz · Last modified: 2012/11/19 18:13 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