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:stress [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:stress [2012/11/20 21:39]
gag [General Informations]
solvers:stress [2016/01/20 09:34] (current)
gag [General Description]
Line 12: Line 12:
  
 ==== General Description ==== ==== 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 is computed using:\\ The Cauchy stress is computed using:\\
Line 18: Line 18:
 <m>sigma_{ij}  = 2 {eta}  {epsilon}_{ij} - p delta_{ij} </m> \\ <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. +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 Solver 1
   Equation = String "StressSolver"   Equation = String "StressSolver"
-  Procedure =  File "ComputeDevStressNS" "ComputeDevStress"+  Procedure =  File "ElmerIceSolvers" "ComputeDevStress"
   ! this is just a dummy, hence no output is needed   ! this is just a dummy, hence no output is needed
   !-----------------------------------------------------------------------   !-----------------------------------------------------------------------
Line 34: Line 36:
   !-----------------------------------------------------------------------   !-----------------------------------------------------------------------
   Flow Solver Name = String "Flow Solution"   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'   Stress Variable Name = String 'Sigma'
   !-----------------------------------------------------------------------   !-----------------------------------------------------------------------
Line 58: Line 60:
  
 ==== Examples ==== ==== Examples ====
-Download an example using this solverTODO +A 2D example can be found in ''[ELMER_TRUNK]/elmerice/Tests/ComputeDevStress''.
 ==== Reference ==== ==== Reference ====
-When used this solver can be cited using the following references:\\+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. 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.1353447588.txt.gz · Last modified: 2012/11/20 21:39 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