This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
solvers:strainrate [2012/11/20 05:26] gag [General Description] |
solvers:strainrate [2014/08/29 15:49] (current) mschafer [SIF contents] |
||
---|---|---|---|
Line 3: | Line 3: | ||
==== General Informations ==== | ==== General Informations ==== | ||
- | * **Solver Fortran File: | + | * **Solver Fortran File: |
- | * **Solver Name: | + | * **Solver Name: |
- | * **Required Output Variable(s): | + | * **Required Output Variable(s): |
- | * **Required Input Variable(s): | + | * **Required Input Variable(s): |
* **Optional Output Variable(s): | * **Optional Output Variable(s): | ||
* **Optional Input Variable(s): | * **Optional Input Variable(s): | ||
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 '' |
The Strain-rate tensor is computed using:\\ | The Strain-rate tensor is computed using:\\ | ||
<m> epsilon_{ij} | <m> epsilon_{ij} | ||
- | + | 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: | ||
+ | |||
< | < | ||
Solver 1 | Solver 1 | ||
- | Equation = String | + | Equation = "Strain Rate" |
- | Procedure = | + | 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 | ||
- | | + | |
- | !----------------------------------------------------------------------- | + | Exported Variable 1 = " |
+ | Exported Variable 1 DOFs = 5 !in 2D, 7 in 3D | ||
+ | |||
+ | ! the name of the variable containing the flow solution (U, | ||
+ | !----------------------------------------------------------------------- | ||
Flow Solver Name = String "Flow Solution" | Flow Solver Name = String "Flow Solution" | ||
- | | + | ! the name of the strain-rate |
- | | + | |
- | !----------------------------------------------------------------------- | + | |
- | Exported Variable 1 = "Sigma" | + | Linear System Solver = Direct |
- | ! [Sxx, Syy, Szz, Sxy, Syz, Szx] in 3D | + | Linear System |
- | | + | |
- | Linear System Solver = " | + | |
- | Linear System | + | |
- | Linear System Max Iterations = 300 | + | |
- | Linear System Convergence Tolerance = 1.0E-09 | + | |
- | Linear System Abort Not Converged = True | + | |
- | Linear System Preconditioning = " | + | |
- | Linear System Residual Output = 1 | + | |
- | End | + | |
- | + | ||
- | Material 1 | + | |
- | ... | + | |
- | ! we want to have the Cauchy stress | + | |
- | !---------------------------------- | + | |
- | Cauchy = Logical True | + | |
End | End | ||
</ | </ | ||
==== Examples ==== | ==== Examples ==== | ||
- | Download an example using this solver. | + | An example using this solver |
- | ==== 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. |