carre2D.sif

! Square 10x10 submitted to horizontal shear                     

Check Keywords Warn

Header
  Mesh DB "." "mesh"
End


Constants
End


Simulation
  Coordinate System = Cartesian 2D
  Simulation Type = Steady State
  Steady State Max Iterations  = 1
  Output Intervals = 1
  Post File   = "carre2D.vtu"
  Output File = "carre2D.result"
  max output level = 4
End  


Body 1
   Equation = 1
   Material = 1
   Body Force = 1
   Initial Condition = 1
End

Body Force 1
Flow BodyForce 1 = Real 0.0
Flow BodyForce 2 = Real 0.0
End


Material 1
  Density = Real 1.0
  Viscosity = Real 1.0
  Viscosity Model = String "power law"
  Viscosity Exponent = Real $1.0/3.0
  Critical Shear Rate = Real 1.0e-10
  Cauchy = Logical True  
End

Initial Condition 1
  velocity 1 = real 0.0
  velocity 2 = real 0.0
  Pressure = real 0.0
End

Solver 1
   Equation = Navier-Stokes

   Linear System Solver = Direct
   Linear System Direct Method = umfpack

   Nonlinear System Max Iterations = 100
   Nonlinear System Convergence Tolerance  = 1.0e-6
   Nonlinear System Newton After Iterations = 100
   Nonlinear System Newton After Tolerance = 1.0e-06
   Nonlinear System Relaxation Factor = 1.0
    
   Stabilization Method = String Bubbles
   Flow Model = Stokes
   Steady state Convergence tolerance = Real 1.0e-5

   Exported Variable 1 = Flow Solution Loads[Stress Vector:2 CEQ Residual:1] 
   Calculate Loads = Logical True
End

Solver 2
  Equation = "Sij"          
  Variable = -nooutput "Sij"     
  Variable DOFs = 1

  Exported Variable 1 = "Stresses"
  Exported Variable 1 DOFs = 4

  Procedure = "ElmerIceSolvers" "ComputeDevStress"
  Flow Solver Name = String "Flow Solution"
  Stress Variable Name = String "Stresses"

  Linear System Solver = Direct         
  Linear System Direct Method = umfpack
End 

Solver 3
   Equation = "EigenStresses"
   Variable = -nooutput dumy
   Variable DOFs = 1

   Procedure = "ElmerIceSolvers" "ComputeEigenValues"

   ! 3 Eigenvalues    
   Exported Variable 1 = "EigenS"
   EigenValue Variable Name = String "EigenS"
   Tensor Variable Name = String "Stresses"
   Exported Variable 1 DOFS = 3

  ! Principal vectors (optional) 
   Exported Variable 2 = EigenVector1
   Exported Variable 2 DOFS = 3
   Exported Variable 3 = EigenVector2
   Exported Variable 3 DOFS =  3
   Exported Variable 4 = EigenVector3
   Exported Variable 4 DOFS = 3
End

Solver 4
  Exec Solver = After All
  Equation = SaveScalars
  Procedure = File "SaveData" "SaveScalars"
  Filename = "results.dat"

  Variable 1 = String "Stresses"
  Operator 1 = String "int"
  Variable 2 = String "Stresses"
  Operator 2 = String "max"
  Variable 3 = String "Stresses"
  Operator 3 = String "min"

  Operator 4 = String "cpu time"
End

Equation 1
    Active Solvers(4) = Integer 1 2 3 4
End

Boundary Condition 1
   Target Boundaries  = 1
   Velocity 2 = Real 0.0
End

Boundary Condition 2
   Target Boundaries = 2
End

Boundary Condition 3
    Target Boundaries  = 3
    Pressure 2 = Real 1.0  
End
 

! Periodic Left
Boundary Condition 4
  Target Boundaries = 4
  Velocity 1 = Real 0.0
End

Solver 3 :: Reference Norm = Real 0.0000000
Solver 3 :: Reference Norm Tolerance = Real 1E-06