test_glenDistrib.sif

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!                                       !!
!!  Compare MATC Glen / Glen law         !!
!!                                       !!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

check keywords warn
echo on
$L = 20.0e3
$Slope = -tan(0.5 * pi / 180.0)

Header
  Mesh DB "." "mesh_B"
End

Constants
! No constant needed
End

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Simulation
  Coordinate System  = Cartesian 2D 
  Simulation Type = Steady State

  Steady State Min Iterations = 1
  Steady State Max Iterations = 1

!  Output File = "test_glen_dist.result"
  Post File = "test_glen_dist.vtu"
  max output level = 3
End

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

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Initial Condition 1
  Pressure = Real 0.0 
  Velocity 1 = Real 0.0 
  Velocity 2 = Real 0.0 
End

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Body Force 1
  Flow BodyForce 1 = Real 0.0                          
  Flow BodyForce 2 = Real -9.7696e15  !MPa - a - m
End

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Material 1
  Density = Real 9.1376e-19  ! MPa - a - m 
  
  Viscosity = 1.0 ! Dummy but avoid warning output

  Viscosity Model = String "glen"
  Glen Exponent = Real 3.0
  Critical Shear Rate = Real 3.0e-11
  Limit Temperature = Real -10.0
  Rate Factor 1 = Real 1.258e13  
  Rate Factor 2 = Real 6.046e28
  Activation Energy 1 = Real 60e3
  Activation Energy 2 = Real 139e3  
  Glen Enhancement Factor = Real 2.5

  Constant Temperature = Variable Depth
   Real  
   0.0 -25.0
   1500.0 -5.0 
   End 
End
 
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Solver 1
  Exec Solver = Before all    

  Equation = Depth
  Variable = Depth
  Variable DOFs = 1

  Procedure =  "ElmerIceSolvers" "FlowdepthSolver"

  Linear System Solver = direct
  Linear System Direct Method = umfpack

  Gradient = Real -1.0E00
  Calc Free Surface = Logical False
End

Solver 2
  Equation = "Navier-Stokes"
  
  Stabilization Method = String Bubbles
  Flow model = Stokes 

  Linear System Solver = Direct         
  Linear System Direct Method = umfpack

  Nonlinear System Max Iterations = 100
  Nonlinear System Convergence Tolerance  = 1.0e-8
  Nonlinear System Newton After Iterations = 5   
  Nonlinear System Newton After Tolerance = 1.0e-02
  Nonlinear System Relaxation Factor = 1.00

  Steady State Convergence Tolerance = Real 1.0e-3
End

 
Solver 3
  Exec Solver =  After TimeStep  ! For transient simulation
  Procedure = File "SaveData" "SaveScalars"
  Filename =  "results.dat"
  Show Norm = True
  Show Norm Index = 1

  Variable 1 = String "Depth"
  Operator 1 = String "boundary max"

  Variable 2 = String "Flow Solution"
  Operator 2 = String "Volume"
  
  Variable 3 = String "Velocity 1"
  Operator 3 = String "Max Abs"

  Operator 4 = String "cpu time"
End

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

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

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! Bedrock
Boundary Condition 1
  Target Boundaries = 1
  Velocity 1 = Real 0.0e0
  Velocity 2 = Real 0.0e0

  Save Scalars = logical True
End

! Periodic Right
Boundary Condition 2
  Target Boundaries = 2
  Periodic BC = 4 
! Periodic BC Translate(2) = Real $ (L) (L*Slope)  
  Periodic BC Translate(2) = Real 2.00000000e+04 -1.74537356e+02
  Periodic BC Velocity 1  = Logical True
  Periodic BC Velocity 2  = Logical True
  Periodic BC Pressure  = Logical True
End

! Upper Surface
Boundary Condition 3
  Target Boundaries = 3
  
  Depth = Real 0.0
End

! Periodic Left
Boundary Condition 4
  Target Boundaries = 4
End

Solver 2 :: Reference Norm = Real 6.8808207
Solver 2 :: Reference Norm Tolerance = Real 1E-05

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