!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !! !! !! Elmer/Ice Course 2011 !! !! !! !! Olivier Gagliardini May 2011 !! !! !! !! Tete Rousse Step 1 - No cavity !! !! !! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ! Updated April 2013 to use internal Mesh Extrusion ! Updated Oct 2014 to use the Grid2dInterpolator solver check keywords warn echo on ! Name of output files $Step = "step1_" ! Define the parameter in MPa - a - m $yearinsec = 365.25*24*60*60 $rhoi = 900.0/(1.0e6*yearinsec^2) $rhow = 1000.0/(1.0e6*yearinsec^2) $gravity = -9.81*yearinsec^2 ! Prefactor from Cuffey and Paterson (2010) in MPa^{-3} a^{-1} $A1 = 2.89165e-13*yearinsec*1.0e18 $A2 = 2.42736e-02*yearinsec*1.0e18 $Q1 = 60.0e3 $Q2 = 115.0e3 Header Mesh DB "." "teterousse" End Constants Water Density = Real $rhow End !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Simulation Coordinate System = Cartesian 3D Simulation Type = Steady Extruded Mesh Levels = Integer 16 Steady State Min Iterations = 1 Steady State Max Iterations = 1 Output File = "teterousse_$Step".result" Post File = "teterousse_$Step".vtu" max output level = 3 End !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ! The ice 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 Velocity 3 = Real 0.0 End !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Body Force 1 Flow BodyForce 1 = Real 0.0 Flow BodyForce 2 = Real 0.0 Flow BodyForce 3 = Real $gravity End !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Material 1 Density = Real $rhoi Viscosity Model = String "glen" Viscosity = 1.0 ! Dummy but avoid warning output Glen Exponent = Real 3.0 Limit Temperature = Real -10.0 Rate Factor 1 = Real $A1 Rate Factor 2 = Real $A2 Activation Energy 1 = Real $Q1 Activation Energy 2 = Real $Q2 Glen Enhancement Factor = Real 1.0 Critical Shear Rate = Real 1.0e-10 Constant Temperature = Real -1.0 Cauchy = Logical True End !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Solver 1 Exec Solver = "Before Simulation" Equation = "Read DEMs" Procedure = "ElmerIceSolvers" "Grid2DInterpolator" ! Bedrock DEM Variable 1 = String "bedDEM" Variable 1 data file = File "../Data/DEM_TR_bed.dat" Variable 1 x0 = Real -1666000.0d0 Variable 1 y0 = Real 841000.0d0 Variable 1 lx = Real 28000.0 Variable 1 ly = Real 33000.0 Variable 1 Nx = Integer 281 Variable 1 Ny = Integer 331 Variable 1 Invert = Logical False Variable 1 Fill = Logical False Variable 1 Position Tol = Real 1.0e-1 Variable 1 No Data = Real -9999.0 Variable 1 No Data Tol = Real 1.0 ! Surface DEM Variable 2 = String "ZsDEM" Variable 2 data file = File "../Data/DEM_TR_surf.dat" Variable 2 x0 = Real -1666000.0d0 Variable 2 y0 = Real 841000.0d0 Variable 2 lx = Real 28000.0 Variable 2 ly = Real 33000.0 Variable 2 Nx = Integer 281 Variable 2 Ny = Integer 331 Variable 2 Invert = Logical False Variable 2 Fill = Logical False Variable 2 Position Tol = Real 1.0e-1 Variable 2 No Data = Real -9999.0 Variable 2 No Data Tol = Real 1.0 End Solver 2 Equation = "MapCoordinate" Procedure = "StructuredMeshMapper" "StructuredMeshMapper" Active Coordinate = Integer 3 Mesh Velocity Variable = String "dSdt" Mesh Update Variable = String "dS" Mesh Velocity First Zero = Logical True Displacement Mode = Logical False Correct Surface = Logical True Minimum Height = Real 1.0 End Solver 3 Equation = "Navier-Stokes" Stabilization Method = String Stabilized Flow Model = Stokes Exported Variable 1 = -dofs 1 "dSdt" Exported Variable 2 = -dofs 1 "dS" Exported Variable 3 = -dofs 1 "BedDEM" Exported Variable 4 = -dofs 1 "ZsDEM" Linear System Solver = Iterative Linear System Iterative Method = BICGStab Linear System Preconditioning = ILU0 Linear System Convergence Tolerance = Real 1.0e-8 Linear System Max Iterations = 500 Linear System Residual Output = 500 Linear System Abort Not Converged = True Nonlinear System Max Iterations = 50 Nonlinear System Convergence Tolerance = 1.0e-5 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 4 Equation = Sij Procedure = "ElmerIceSolvers" "ComputeDevStress" Variable = -nooutput "Sij" Variable DOFs = 1 Exported Variable 1 = -dofs 6 Stress Stress Variable Name = String "Stress" Flow Solver Name = String "Flow Solution" Linear System Solver = Direct Linear System Direct Method = umfpack End Solver 5 Equation = "EigenStresses" Procedure = "ElmerIceSolvers" "ComputeEigenValues" Variable = -nooutput dumy Variable DOFs = 1 ! The 3 eigenvalues Exported Variable 1 = EigenStress Exported Variable 1 DOFS = 3 ! The 3 eigenvectors (Option) 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 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Equation 1 Active Solvers(5) = 1 2 3 4 5 End !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ! lateral side of the glacier Boundary Condition 1 Target Boundaries = 1 Velocity 1 = real 0.0 Velocity 2 = real 0.0 End ! cavity roof and Bedrock Boundary Condition 2 Bottom Surface = Equals "BedDEM" !-------------------------------------------- ! --------- No sliding Velocity 1 = Real 0.0 Velocity 2 = Real 0.0 Velocity 3 = Real 0.0 !----------- OR ----------------------------- ! --------- Sliding ! Flow Force BC = Logical True ! Normal-Tangential Velocity = Logical True ! Velocity 1 = Real 0.0e0 ! Slip Coefficient 2 = Real 0.1 ! Slip Coefficient 3 = Real 0.1 !-------------------------------------------- End ! Upper Surface Boundary Condition 3 Top Surface = Equals "ZsDEM" End