[[solvers:stress]]
Trace:
Show page
AdminRecent ChangesSitemap

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
Last revision Both sides next revision
solvers:stress [2012/11/13 21:20]
gag [General Description] 		solvers:stress [2016/01/20 09:34]
gag [General Description]
Line 1: Line 1:
 +~~NOTOC~~
 ===== Solver ComputeDevStress ===== ===== Solver ComputeDevStress =====
  
 ==== General Informations ==== ==== 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 Output Variable(s):** None  
   * **Optional Input 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 ==== ==== 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 are computed using:\\ 
 <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 tensiile 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 27: 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 35: 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'
   !-----------------------------------------------------------------------   !-----------------------------------------------------------------------
   Exported Variable 1 = "Sigma" ! [Sxx, Syy, Szz, Sxy] in 2D   Exported Variable 1 = "Sigma" ! [Sxx, Syy, Szz, Sxy] in 2D
Line 58: Line 60:
  
 ==== Examples ==== ==== 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.txt · Last modified: 2016/01/20 09:34 by gag
Show pageOld revisions
Media ManagerBack to top
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