! Define parameters in MPa - a - m $yearinsec = 365.25*24*60*60 $rhoi = 910.0/(1.0e6*yearinsec^2) $gravity = -9.81*yearinsec^2 $clausclap = 9.8E-08 $B = 14.46 !MPa-3a-1 T = -14°C $n = 3.0 $ function conductivity(T) { _conductivity=9.828*exp(-5.7E-03*T)} $ function capacity(T) { _capacity=146.3+(7.253*T)} Header Mesh DB "." "Data" End Constants Density Name = String "Relative Density" End Simulation Coordinate System = Cartesian 3D Simulation Type = Steady Extruded Mesh Levels = Integer 16 Steady State Min Iterations = 1 Steady State Max Iterations = 20 Output File = "porous_thermomech_slide.result" Post File = "porous_thermomech_slide.vtu" Max Output Level = 3 Initialize Dirichlet Conditions = Logical False Restart File = "porous_temp_constrained.result" Restart Position = 0 End Body 1 Equation = 1 Body Force = 1 Material = 1 Initial Condition = 1 End Initial Condition 1 Depth = Real 0.0 End Body Force 1 Porous Force 1 = Real 0.0E00 Porous Force 2 = Real 0.0E00 Porous Force 3 = Real $gravity*rhoi Heat Source = Real 0.0 End Material 1 Powerlaw Exponent = Real $n Min Second Invariant = Real 1.0E-10 Fluidity Parameter = Real $B ! MPa^{-3}a^{-1} Temperature Field Variable = String "Temp Homologous" Temp Heat Capacity = Variable Temp Real MATC "capacity(tx)*(31557600.0)^(2.0)" Temp Heat Conductivity = Variable Temp Real MATC "conductivity(tx)*31557600.0*1.0E-06" Temp Upper Limit = Variable Depth Real MATC "273.15 - clausclap* tx * 910.0 * 9.81" End Solver 1 Exec Solver = "Never" Equation = "Read DEMs" Procedure = "ElmerIceSolvers" "Grid2DInterpolator" ! Bedrock DEM Variable 1 = String "ZbDEM" Variable 1 data file = File "./Data/glacier_base.dat" Variable 1 x0 = Real 0.0d0 Variable 1 y0 = Real 0.0d0 Variable 1 lx = Real 4000.0 Variable 1 ly = Real 3000.0 Variable 1 Nx = Integer 161 Variable 1 Ny = Integer 121 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/glacier_sur2_25m.dat" Variable 2 x0 = Real 0.0d0 Variable 2 y0 = Real 0.0d0 Variable 2 lx = Real 4000.0 Variable 2 ly = Real 3000.0 Variable 2 Nx = Integer 161 Variable 2 Ny = Integer 121 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 = "HeightDepth" Procedure = "StructuredProjectToPlane" "StructuredProjectToPlane" Active Coordinate = Integer 3 Operator 1 = Depth End Solver 4 Equation = String "PorousFlow" Procedure = "ElmerIceSolvers" "PorousSolver" Variable = "Porous" Variable DOFs = 4 Optimize Bandwidth = False Exported Variable 1 = -dofs 1 -nooutput "dSdt" Exported Variable 2 = -dofs 1 -nooutput "dS" Exported Variable 3 = -dofs 1 "ZbDEM" Exported Variable 4 = -dofs 1 "ZsDEM" Exported Variable 5 = String "Relative Density" Exported variable 5 DOFs = Integer 1 Linear System Solver = "Iterative" Linear System Iterative Method = "BiCGStab" Linear System Max Iterations = 10000 !500 Linear System Convergence Tolerance = 1.0E-07 Linear System Abort Not Converged = True Linear System Preconditioning = "ILU0" Linear System Residual Output = 10 Nonlinear System Convergence Tolerance = 1.0E-05 Nonlinear System Max Iterations = 50 Nonlinear System Reset Newton = Logical True Steady State Convergence Tolerance = 1.0E-03 End Solver 5 Equation = String "Homologous Temperature Equation" Procedure = File "ElmerIceSolvers" "TemperateIceSolver" Variable = String "Temp" Variable DOFs = 1 Stabilize = True Optimize Bandwidth = Logical True Linear System Solver = "Direct" Linear System Direct Method = UMFPACK Linear System Convergence Tolerance = 1.0E-06 Linear System Abort Not Converged = True !False Linear System Preconditioning = "ILU1" Linear System Residual Output = 0 Nonlinear System Convergence Tolerance = 1.0E-07 Nonlinear System Max Iterations = 100 Nonlinear System Relaxation Factor = Real 9.999E-01 Steady State Convergence Tolerance = 1.0E-04 Apply Dirichlet = Logical True Exported Variable 1 = String -nooutput "Temp Homologous" Exported Variable 1 DOFs = 1 Exported Variable 2 = String -nooutput "Temp Residual" Exported Variable 2 DOFs = 1 End Equation 1 Active Solvers(5) = 1 2 3 4 5 Convection = "Computed" Flow Solution Name = String "Porous" End Boundary Condition 1 Name = "Sides" Target Boundaries = 1 Porous 1 = Real 0.0 Porous 2 = Real 0.0 End Boundary Condition 2 Name = "Bedrock" Bottom Surface = Equals ZbDEM Temp = Variable Depth Real MATC "273.15 - clausclap* tx * 910.0 * 9.81" Normal-Tangential Porous = Logical True Porous 1 = Real 0.0 Porous Slip Coeff 2 = Real 0.1 Porous Slip Coeff 3 = Real 0.1 End Boundary Condition 3 Name = "Surface" Top Surface = Equals ZsDEM Temp = Real 258.15 End