From 76008bc22d58150f32c670d4795d017e0a477d05 Mon Sep 17 00:00:00 2001 From: IskanderI <123633876+IskanderI@users.noreply.github.com> Date: Sun, 5 Nov 2023 22:42:48 +0100 Subject: [PATCH 1/7] Update Volcano.dat Fixed running parameters and added magma chamber as a thermal anomaly --- input_models/ScalingTests/Volcano.dat | 15 +++++++++++---- 1 file changed, 11 insertions(+), 4 deletions(-) diff --git a/input_models/ScalingTests/Volcano.dat b/input_models/ScalingTests/Volcano.dat index d11a18ce..8c0d4d35 100644 --- a/input_models/ScalingTests/Volcano.dat +++ b/input_models/ScalingTests/Volcano.dat @@ -97,10 +97,10 @@ unit_stress = 1e9 #=============================================================================== SolverType = multigrid # solver [direct or multigrid] - MGLevels = 6 # number of MG levels [default=3] - MGSweeps = 5 # number of MG smoothening steps per level [default=10] - MGSmoother = jacobi # type of smoothener used [chebyshev or jacobi] - MGJacobiDamp = 0.5 # Dampening parameter [only employed for Jacobi smoothener; default=0.6] + MGLevels = 5 # number of MG levels [default=3] + MGSweeps = 25 # number of MG smoothening steps per level + MGSmoother = jacobi # type of smoothener used [chebyshev or jacobi] + MGJacobiDamp = 0.5 # Dampening parameter [only employed for Jacobi smoothener; default=0.6] #=============================================================================== @@ -125,6 +125,13 @@ unit_stress = 1e9 + + phase = 1 + axes = 14.0 14.0 1.0 # semi-axes of ellipsoid in x, y and z + center = 0.0 0.0 -4.0 + Temperature = constant # optional: Temperature of the sphere. possibilities: [constant] + cstTemp = 1000 # required in case of [constant]: temperature value [in Celcius in case of GEO units] + #=============================================================================== # Output #=============================================================================== From f12215178441a8a096aa22f8a4c05b1ed9651c22 Mon Sep 17 00:00:00 2001 From: IskanderI <123633876+IskanderI@users.noreply.github.com> Date: Mon, 6 Nov 2023 10:33:10 +0100 Subject: [PATCH 2/7] Update Volcano.dat activated steady-state diffusion --- input_models/ScalingTests/Volcano.dat | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/input_models/ScalingTests/Volcano.dat b/input_models/ScalingTests/Volcano.dat index 8c0d4d35..a34f5e16 100644 --- a/input_models/ScalingTests/Volcano.dat +++ b/input_models/ScalingTests/Volcano.dat @@ -80,9 +80,9 @@ unit_stress = 1e9 gravity = 0.0 0.0 -9.81 # gravity vector act_temp_diff = 1 # temperature diffusion activation flag - #act_steady_temp = 0 # steady-state temperature initial guess activation flag - #steady_temp_t = 0.1 # time for (quasi-)steady-state temperature initial guess - #nstep_steady = 50 # number of steps for (quasi-)steady-state temperature initial guess (default = 1) + act_steady_temp = 1 # steady-state temperature initial guess activation flag + steady_temp_t = 0.1 # time for (quasi-)steady-state temperature initial guess + nstep_steady = 50 # number of steps for (quasi-)steady-state temperature initial guess (default = 1) act_heat_rech = 0 # heat recharge activation flag init_lith_pres = 1 # initial pressure with lithostatic pressure init_guess = 1 # initial guess flag From 33dc3fbc14c1471d79d4493a834e5d4464be2f32 Mon Sep 17 00:00:00 2001 From: IskanderI <123633876+IskanderI@users.noreply.github.com> Date: Thu, 9 Nov 2023 10:57:47 +0100 Subject: [PATCH 3/7] Update Volcano.dat Changed more parameters --- input_models/ScalingTests/Volcano.dat | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/input_models/ScalingTests/Volcano.dat b/input_models/ScalingTests/Volcano.dat index a34f5e16..46c50f9a 100644 --- a/input_models/ScalingTests/Volcano.dat +++ b/input_models/ScalingTests/Volcano.dat @@ -82,7 +82,7 @@ unit_stress = 1e9 act_temp_diff = 1 # temperature diffusion activation flag act_steady_temp = 1 # steady-state temperature initial guess activation flag steady_temp_t = 0.1 # time for (quasi-)steady-state temperature initial guess - nstep_steady = 50 # number of steps for (quasi-)steady-state temperature initial guess (default = 1) + nstep_steady = 150 # number of steps for (quasi-)steady-state temperature initial guess (default = 1) act_heat_rech = 0 # heat recharge activation flag init_lith_pres = 1 # initial pressure with lithostatic pressure init_guess = 1 # initial guess flag @@ -97,7 +97,7 @@ unit_stress = 1e9 #=============================================================================== SolverType = multigrid # solver [direct or multigrid] - MGLevels = 5 # number of MG levels [default=3] + MGLevels = 5 # number of MG levels MGSweeps = 25 # number of MG smoothening steps per level MGSmoother = jacobi # type of smoothener used [chebyshev or jacobi] MGJacobiDamp = 0.5 # Dampening parameter [only employed for Jacobi smoothener; default=0.6] From 77ec62c600a49804d628573649364f454cbb25ee Mon Sep 17 00:00:00 2001 From: IskanderI <123633876+IskanderI@users.noreply.github.com> Date: Thu, 23 Nov 2023 13:54:48 +0100 Subject: [PATCH 4/7] Update Volcano.dat Changed plastic properties --- input_models/ScalingTests/Volcano.dat | 27 +++++++++++++++------------ 1 file changed, 15 insertions(+), 12 deletions(-) diff --git a/input_models/ScalingTests/Volcano.dat b/input_models/ScalingTests/Volcano.dat index 46c50f9a..1b925116 100644 --- a/input_models/ScalingTests/Volcano.dat +++ b/input_models/ScalingTests/Volcano.dat @@ -17,14 +17,14 @@ unit_stress = 1e9 # Time stepping parameters #=============================================================================== - dt = 5e-4 # time step - dt_min = 4.5e-4 # minimum time step (declare divergence if lower value is attempted) + dt = 1e-9 # time step + dt_min = 1e-12 # minimum time step (declare divergence if lower value is attempted) dt_max = 1e-3 # maximum time step inc_dt = 1 # time step increment per time step (fraction of unit) CFL = 0.5 # CFL (Courant-Friedrichs-Lewy) criterion CFLMAX = 0.8 # CFL criterion for elasticity nstep_ini = 0 # save output for n initial steps - nstep_max = 5 # maximum allowed number of steps (lower bound: time_end/dt_max) + nstep_max = 10 # maximum allowed number of steps (lower bound: time_end/dt_max) nstep_out = 1 # save output every n steps #=============================================================================== @@ -64,8 +64,8 @@ unit_stress = 1e9 # noslip = 0 0 0 0 0 0 - temp_top = 0 # Temperature @ top - temp_bot = 1350 # Temperature @ bottom; side BC's are flux-free + temp_top = 0 # Temperature @ top + temp_bot = 800 # Temperature @ bottom; side BC's are flux-free # Background strain rate parameters # exx_num_periods = 1 # number intervals of constant strain rate (x-axis) @@ -86,21 +86,24 @@ unit_stress = 1e9 act_heat_rech = 0 # heat recharge activation flag init_lith_pres = 1 # initial pressure with lithostatic pressure init_guess = 1 # initial guess flag - eta_min = 1e19 # viscosity upper bound + eta_min = 1e17 # viscosity upper bound eta_max = 1e23 # viscosity lower limit eta_ref = 1e22 # reference viscosity (initial guess) T_ref = 20 # reference temperature p_lim_plast = 1 - + min_cohes = 1e6 # cohesion lower bound [Pa] + min_fric = 1.0 # friction lower bound [degree] + tau_ult = 1e9 # ultimate yield stress [Pa] + #=============================================================================== # Solver options #=============================================================================== SolverType = multigrid # solver [direct or multigrid] - MGLevels = 5 # number of MG levels + MGLevels = 4 # number of MG levels MGSweeps = 25 # number of MG smoothening steps per level MGSmoother = jacobi # type of smoothener used [chebyshev or jacobi] - MGJacobiDamp = 0.5 # Dampening parameter [only employed for Jacobi smoothener; default=0.6] + MGJacobiDamp = 0.5 # Dampening parameter [only employed for Jacobi smoothener; default=0.6] #=============================================================================== @@ -224,8 +227,8 @@ unit_stress = 1e9 # Plastic parameters ch = 10e6 - fr = 30 - eta_st = 5e20 # stabilization viscosity (default is eta_min) + fr = 20 + eta_st = 5e20 # stabilization viscosity (default is eta_min) @@ -251,7 +254,7 @@ unit_stress = 1e9 # Plastic parameters ch = 1e6 - fr = 10 + fr = 5 eta_st = 5e20 # stabilization viscosity (default is eta_min) From 079ababc2f69b31f332fdfa0c858605f39002308 Mon Sep 17 00:00:00 2001 From: IskanderI <123633876+IskanderI@users.noreply.github.com> Date: Fri, 24 Nov 2023 14:26:36 +0100 Subject: [PATCH 5/7] Update Volcano.dat Fixed running parameters: plastic strain in 1 timestep --- input_models/ScalingTests/Volcano.dat | 24 ++++++++++++------------ 1 file changed, 12 insertions(+), 12 deletions(-) diff --git a/input_models/ScalingTests/Volcano.dat b/input_models/ScalingTests/Volcano.dat index 1b925116..34641ed5 100644 --- a/input_models/ScalingTests/Volcano.dat +++ b/input_models/ScalingTests/Volcano.dat @@ -8,23 +8,23 @@ units = geo # Always in SI units!! -unit_temperature = 500 +unit_temperature = 1000 unit_length = 1000 -unit_viscosity = 1e19 +unit_viscosity = 1e20 unit_stress = 1e9 #=============================================================================== # Time stepping parameters #=============================================================================== - dt = 1e-9 # time step - dt_min = 1e-12 # minimum time step (declare divergence if lower value is attempted) + dt = 1e-7 # time step + dt_min = 1e-9 # minimum time step (declare divergence if lower value is attempted) dt_max = 1e-3 # maximum time step - inc_dt = 1 # time step increment per time step (fraction of unit) + inc_dt = 0,5 # time step increment per time step (fraction of unit) CFL = 0.5 # CFL (Courant-Friedrichs-Lewy) criterion CFLMAX = 0.8 # CFL criterion for elasticity nstep_ini = 0 # save output for n initial steps - nstep_max = 10 # maximum allowed number of steps (lower bound: time_end/dt_max) + nstep_max = 1 # maximum allowed number of steps (lower bound: time_end/dt_max) nstep_out = 1 # save output every n steps #=============================================================================== @@ -100,7 +100,7 @@ unit_stress = 1e9 #=============================================================================== SolverType = multigrid # solver [direct or multigrid] - MGLevels = 4 # number of MG levels + MGLevels = 3 # number of MG levels MGSweeps = 25 # number of MG smoothening steps per level MGSmoother = jacobi # type of smoothener used [chebyshev or jacobi] MGJacobiDamp = 0.5 # Dampening parameter [only employed for Jacobi smoothener; default=0.6] @@ -187,7 +187,7 @@ unit_stress = 1e9 Name = air ID = 0 - rho = 100 + rho = 1 alpha = 3e-5 # Linear Viscosity @@ -202,7 +202,7 @@ unit_stress = 1e9 Cp = 1000 # Plastic parameters - ch = 10e6 + ch = 20e6 fr = 30 @@ -226,8 +226,8 @@ unit_stress = 1e9 Cp = 1000 # Plastic parameters - ch = 10e6 - fr = 20 + ch = 5e6 + fr = 10 eta_st = 5e20 # stabilization viscosity (default is eta_min) @@ -286,7 +286,7 @@ unit_stress = 1e9 -snes_max_it 200 -snes_rtol 1e-5 -snes_atol 1e-6 - -snes_PicardSwitchToNewton_rtol 5e-2 #-7 + -snes_PicardSwitchToNewton_rtol 1e-3 #-7 #-snes_NewtonSwitchToPicard_it 1 # number of Newton iterations after which we switch back to Picard # -snes_monitor -snes_linesearch_monitor From efe4f7f545b6b44b2205f33c8ddb38c8387a1eb6 Mon Sep 17 00:00:00 2001 From: IskanderI <123633876+IskanderI@users.noreply.github.com> Date: Fri, 24 Nov 2023 15:06:05 +0100 Subject: [PATCH 6/7] Update Volcano.dat Fixed for dot instead of comma --- input_models/ScalingTests/Volcano.dat | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/input_models/ScalingTests/Volcano.dat b/input_models/ScalingTests/Volcano.dat index 34641ed5..b81cadd8 100644 --- a/input_models/ScalingTests/Volcano.dat +++ b/input_models/ScalingTests/Volcano.dat @@ -20,7 +20,7 @@ unit_stress = 1e9 dt = 1e-7 # time step dt_min = 1e-9 # minimum time step (declare divergence if lower value is attempted) dt_max = 1e-3 # maximum time step - inc_dt = 0,5 # time step increment per time step (fraction of unit) + inc_dt = 0.5 # time step increment per time step (fraction of unit) CFL = 0.5 # CFL (Courant-Friedrichs-Lewy) criterion CFLMAX = 0.8 # CFL criterion for elasticity nstep_ini = 0 # save output for n initial steps From 49277df257eaf496eed8bb7b23552ab5e8f18c2f Mon Sep 17 00:00:00 2001 From: Boris Kaus <61824822+boriskaus@users.noreply.github.com> Date: Tue, 28 Nov 2023 08:30:44 +0100 Subject: [PATCH 7/7] Update README.md --- README.md | 1 + 1 file changed, 1 insertion(+) diff --git a/README.md b/README.md index 2788c330..495d4753 100644 --- a/README.md +++ b/README.md @@ -1,6 +1,7 @@ # LaMEM [![Build Status](https://github.com/UniMainzGeo/LaMEM/workflows/CI/badge.svg)](https://github.com/UniMainzGeo/LaMEM/actions) [![Dev](https://img.shields.io/badge/docs-dev-blue.svg)](https://unimainzgeo.github.io/LaMEM/dev) +[![DOI](https://zenodo.org/badge/641438500.svg)](https://zenodo.org/doi/10.5281/zenodo.10211634) LaMEM (Lithosphere and Mantle Evolution Model) is a parallel 3D numerical code that can be used to simulate various thermo-mechanical geodynamical processes such as mantle-lithosphere interaction for rocks