!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !! !! !! Elmer/Ice Course 2011 !! !! !! !! Olivier Gagliardini May 2011 !! !! !! !! Tete Rousse Step 1 - No cavity !! !! !! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !! Mod by Norbert Bader to extrude measured surface !! Test stage! ! Updated April 2013 to use internal Mesh Extrusion ! Updated Oct 2014 to use the Grid2dInterpolator solver ! Updated Oct 2018 to use LUA language instead of MATC 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 "." "dimple2d" 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 = 40 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 "DEM_TR_bed.dat" Variable 1 x0 = Real 0.0 Variable 1 y0 = Real 0.0 Variable 1 lx = Real 6512.347200 Variable 1 ly = Real 2022.023200 Variable 1 Nx = Integer 237 Variable 1 Ny = Integer 74 Variable 1 Invert = Logical False Variable 1 Fill = Logical False Variable 1 Position Tol = Real 1.0 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 "DEM_TR_surf.dat" Variable 2 x0 = Real 0.0 Variable 2 y0 = Real 0.0 Variable 2 lx = Real 6512.347200 Variable 2 ly = Real 2022.023200 Variable 2 Nx = Integer 237 Variable 2 Ny = Integer 74 Variable 2 Invert = Logical False Variable 2 Fill = Logical False Variable 2 Position Tol = Real 1.0 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 Mapped Mesh Name = String my_mapped_mesh End Solver 3 Equation = "Navier-Stokes" Exec Solver = "never" 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" Exec Solver = "never" 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" Exec Solver = "never" Procedure = "ElmerIceSolvers" "ComputeEigenValues" Variable = -nooutput dumy Variable DOFs = 1 Tensor Variable Name = String "Stress" ! 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 = 0 Bottom Surface = Equals BedDEM !-------------------------------------------- ! --------- No sliding ! Velocity 1 = Real 0.0 ! Velocity 2 = Real 0.0 ! Velocity 3 = Real 0.0 !----------- OR ----------------------------- ! --------- Sliding ! 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 !Boundary Condition 4 ! Target Boundaries = 2 ! Velocity 1 = 0.0 !End !Boundary Condition 5 ! Target Boundaries = 3 ! Velocity 1 = 0.0 !End !Boundary Condition 6 ! Target Boundaries = 4 ! Velocity 1 = 0.0 !End