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

Differences

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

Link to this comparison view

Next revision 		Previous revision
solvers:strainrate [2012/11/20 05:25]
gag created 		solvers:strainrate [2014/08/29 15:49] (current)
mschafer [SIF contents]
Line 3: Line 3:
  
 ==== General Informations ==== ==== General Informations ====
-  * **Solver Fortran File:** //ComputeStrainRate.f90// +  * **Solver Fortran File:** ''ComputeStrainRate.f90'' 
-  * **Solver Name:** //ComputeStrainRate//  +  * **Solver Name:** ''ComputeStrainRate'' 
-  * **Required Output Variable(s):** default is //StrainRate// (else in //StrainRate Variable Name//)   +  * **Required Output Variable(s):** default is ''StrainRate'' (else in ''StrainRate Variable Name'')   
-  * **Required Input Variable(s):** A //Flow Solution// (in //Flow Solution Name//)+  * **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
Line 12: Line 12:
  
 ==== General Description ==== ==== General Description ====
-The aim of this solver is to compute the strain-rate tensor from the flow solution. For a 2D simulation there are 5 DOFs (E11, E22, E33, E12, Eii) and for a 3D simulation, 2 additional (7 DOFs) are being solved for (E11, E22, E33, E12, E23, E31, Eii). This solver uses a dummy variable and solves 5 (7 in 3D) times a 1 DOF system for each strain-rate components and the trace of the strain-rate tensor (Eii). +The aim of this solver is to compute the strain-rate tensor from the flow solution. For a 2D simulation there are 5 DOFs ''(E11, E22, E33, E12, Eii)'' and for a 3D simulation, 2 additional (7 DOFs) are being solved for ''(E11, E22, E33, E12, E23, E31, Eii)''. This solver uses a dummy variable and solves 5 (7 in 3D) times a 1 DOF system for each strain-rate components and the trace of the strain-rate tensor ''(Eii)''. 
  
 The Strain-rate tensor is computed using:\\ The Strain-rate tensor is computed using:\\
- +<m> epsilon_{ij}  =  {1/2} ({u_{i,j} + u_{j,i}}) </m> \\ 
-<m>epsilon_{ij}  =  {1/2} (u_{i,j} + u_{j,i})  </m> \\ +where <m> u </m> is the velocity vector solution of the Stokes problem.
-  +
  
 ==== SIF contents ==== ==== SIF contents ====
 The required keywords in the SIF file for this solver are: The required keywords in the SIF file for this solver are:
 +
 <code> <code>
 Solver 1 Solver 1
-  Equation = String "StressSolver" +  Equation = "Strain Rate" 
-  Procedure =  File "ComputeDevStressNS" "ComputeDevStress" +  Procedure = "ElmerIceSolvers" "ComputeStrainRate" 
-  ! this is just a dummy, hence no output is needed +! this is just a dummy, hence no output is needed 
-  !----------------------------------------------------------------------- +!-----------------------------------------------------------------------   
-  Variable = -nooutput "Sij"+  Variable = -nooutput "Eij"
   Variable DOFs = 1   Variable DOFs = 1
-  ! the name of the variable containing the flow solution (U,V,W,Pressure) + 
-  !-----------------------------------------------------------------------+  Exported Variable 1 = "StrainRate" 
 +  Exported Variable 1 DOFs = 5 !in 2D, 7 in 3D 
 + 
 +! the name of the variable containing the flow solution (U,V,W,Pressure) 
 +!-----------------------------------------------------------------------
   Flow Solver Name = String "Flow Solution"   Flow Solver Name = String "Flow Solution"
-  ! the name of the stress solution (default is 'stress') +! the name of the strain-rate solution (default is 'StrainRate') 
-  Stress Variable Name = String 'Sigma' +  StrainRate Variable Name = String "StrainRate" 
-  !----------------------------------------------------------------------- +   
-  Exported Variable 1 = "Sigma" ! [Sxx, Syy, Szz, Sxy] in 2D +  Linear System Solver = Direct 
-                                 ! [Sxx, Syy, Szz, Sxy, Syz, Szx] in 3D +  Linear System Direct Method = umfpack
-  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 End
 </code> </code>
  
 ==== Examples ==== ==== Examples ====
-Download an example using this solver. TODO+An example using this solver can be found in ''[ELMER_TRUNK]/elmerice/examples/Test_StrainRate''. 
  
-==== 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/strainrate.1353389143.txt.gz · Last modified: 2012/11/20 05:25 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