Differences

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

--- solvers:stress [2012/11/13 21:19]
gag [SIF contents]
+++ solvers:stress [2016/01/20 09:34] (current)
gag [General Description]
@@ Line 1: / Line 1: @@
+~~NOTOC~~
 ===== Solver ComputeDevStress =====
 ==== General Informations ====
-  * **Solver Fortran File:** //ComputeDevStressNS.f90//
+  * **Solver Fortran File:** ''ComputeDevStressNS.f90''
-  * **Solver Name:** //ComputeDevStress//
+  * **Solver Name:** ''ComputeDevStress''
-  * **Required Output Variable(s):** default is //Stress//
+  * **Required Output Variable(s):** default is ''Stress'' (else in ''Stress Variable Name'')
-  * **Required Input Variable(s):** //A Flow Solution//
+  * **Required Input Variable(s):** A ''Flow Solution'' (in ''Flow Solution Name'')
   * **Optional Output Variable(s):** None
   * **Optional Input Variable(s):** None
-==== 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.
 ==== General Description ====
-The aim of this solver is to compute deviatoric or Cauchy stress from flow solution. For a 2D simulation there are 4 DOFs (S11, S22, S33, S12), for a 3D simulation, 2 additional are being solved for (S11, S22, S33, S12, S23, S31). This solver uses a dummy variable and solves 4 (6 in 3D) times a 1 DOF system for each stress components.
+The aim of this solver is to compute deviatoric or Cauchy stress from flow solution. For a 2D simulation there are 4 DOFs ''(S11, S22, S33, S12)'', for a 3D simulation, 2 additional are being solved for ''(S11, S22, S33, S12, S23, S31)''. This solver uses a dummy variable and solves 4 (6 in 3D) times a 1 DOF system for each stress components.
-The Cauchy stress are computed using:\\
+The Cauchy stress is computed using:\\
-<m>sigma_{ij}  = 2 {eta}  {epsilon}_{ij} - p delta_{ij} </m>
-where <m>epsilon</m> is directly evaluated from the velocity field and <m>p</m> is the isotropic pressure.
+<m>sigma_{ij}  = 2 {eta}  {epsilon}_{ij} - p delta_{ij} </m> \\
+where <m>epsilon</m> is directly evaluated from the velocity field and <m>p</m> is the isotropic pressure.The convention is that a positive stress corresponds to a tensile stress (opposite to the isotropic pressure convention).
+This solver doesn't work for the GOLF anisotropic ([[solvers:aiflow|AIFlow Solver]]) and the snow/firn ([[solvers:porous|Porous Solver]]) rheologies. Nevertheless, these two solvers have intrinsic functions which allow to compute the stress directly.
@@ Line 26: / Line 28: @@
 Solver 1
   Equation = String "StressSolver"
-  Procedure =  File "ComputeDevStressNS" "ComputeDevStress"
+  Procedure =  File "ElmerIceSolvers" "ComputeDevStress"
   ! this is just a dummy, hence no output is needed
   !-----------------------------------------------------------------------
@@ Line 34: / Line 36: @@
   !-----------------------------------------------------------------------
   Flow Solver Name = String "Flow Solution"
-  ! the name of the stress solution (default is 'stress')
+  ! no default value anymore for "Stress Variable Name"
+  Stress Variable Name = String 'Sigma'
   !-----------------------------------------------------------------------
   Exported Variable 1 = "Sigma" ! [Sxx, Syy, Szz, Sxy] in 2D
@@ Line 57: / Line 60: @@
 ==== Examples ====
-Download an example using this solver. TODO
+A 2D example can be found in ''[ELMER_TRUNK]/elmerice/Tests/ComputeDevStress''.
+==== Reference ====
+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/stress.1352841598.txt.gz · Last modified: 2012/11/13 21:19 by gag

Show page Old revisions

Media Manager Back to top