Skip to content
Snippets Groups Projects
Commit a339015f authored by David Seus's avatar David Seus
Browse files

setup layered soil examples

parent a41a6494
No related branches found
No related tags found
No related merge requests found
...@@ -47,8 +47,8 @@ resolutions = { ...@@ -47,8 +47,8 @@ resolutions = {
# 2: 1e-6, # 2: 1e-6,
# 4: 1e-6, # 4: 1e-6,
# 8: 1e-5, # 8: 1e-5,
# 16: 5e-6, 16: 5e-6,
32: 2e-6, # 32: 2e-6,
# 64: 2e-6, # 64: 2e-6,
# 128: 1e-6, # 128: 1e-6,
# 256: 1e-6, # 256: 1e-6,
...@@ -65,12 +65,12 @@ number_of_timesteps = 1 ...@@ -65,12 +65,12 @@ number_of_timesteps = 1
# LDD scheme parameters ###################################################### # LDD scheme parameters ######################################################
Lw1 = 0.007 # /timestep_size Lw1 = 0.007 # /timestep_size
Lnw1 = 0.005 Lnw1 = 0.001
Lw2 = 0.007 # /timestep_size Lw2 = 0.007 # /timestep_size
Lnw2 = 0.005 Lnw2 = 0.001
Lw3 = 0.0007 # /timestep_size Lw3 = 0.007 # /timestep_size
Lnw3 = 0.0005 Lnw3 = 0.0005
Lw4 = 0.0007 # /timestep_size Lw4 = 0.007 # /timestep_size
Lnw4 = 0.0005 Lnw4 = 0.0005
lambda12_w = 0.5 lambda12_w = 0.5
......
...@@ -37,7 +37,7 @@ thisfile = "TP-TP-layered_soil.py" ...@@ -37,7 +37,7 @@ thisfile = "TP-TP-layered_soil.py"
# GENERAL SOLVER CONFIG ###################################################### # GENERAL SOLVER CONFIG ######################################################
# maximal iteration per timestep # maximal iteration per timestep
max_iter_num = 10 max_iter_num = 1000
FEM_Lagrange_degree = 1 FEM_Lagrange_degree = 1
# GRID AND MESH STUDY SPECIFICATIONS ######################################### # GRID AND MESH STUDY SPECIFICATIONS #########################################
...@@ -47,8 +47,8 @@ resolutions = { ...@@ -47,8 +47,8 @@ resolutions = {
# 2: 1e-6, # 2: 1e-6,
# 4: 1e-6, # 4: 1e-6,
# 8: 1e-5, # 8: 1e-5,
# 16: 5e-6, 16: 5e-6,
32: 5e-6, # 32: 2e-6,
# 64: 2e-6, # 64: 2e-6,
# 128: 1e-6, # 128: 1e-6,
# 256: 1e-6, # 256: 1e-6,
...@@ -57,30 +57,31 @@ resolutions = { ...@@ -57,30 +57,31 @@ resolutions = {
# starttimes gives a list of starttimes to run the simulation from. # starttimes gives a list of starttimes to run the simulation from.
# The list is looped over and a simulation is run with t_0 as initial time # The list is looped over and a simulation is run with t_0 as initial time
# for each element t_0 in starttimes. # for each element t_0 in starttimes.
starttimes = {0: 0.0} starttimes = {0: 0.0, 1: 0.3, 2: 0.6, 3: 0.9}
# starttimes = {0: 0.0}
timestep_size = 0.001 timestep_size = 0.001
number_of_timesteps = 5 number_of_timesteps = 1
# LDD scheme parameters ###################################################### # LDD scheme parameters ######################################################
Lw1 = 0.025 # /timestep_size Lw1 = 0.007
Lnw1 = Lw1 Lnw1 = 0.005
Lw2 = 0.025 # /timestep_size Lw2 = 0.007
Lnw2 = Lw2 Lnw2 = 0.005
Lw3 = 0.025 # /timestep_size Lw3 = 0.007
Lnw3 = Lw3 Lnw3 = 0.005
Lw4 = 0.025 # /timestep_size Lw4 = 0.007
Lnw4 = Lw4 Lnw4 = 0.005
lambda12_w = 40 lambda12_w = 0.05
lambda12_nw = 40 lambda12_nw = 0.05
lambda23_w = 40 lambda23_w = 0.05
lambda23_nw = 40 lambda23_nw = 0.05
lambda34_w = 40 lambda34_w = 0.05
lambda34_nw = 40 lambda34_nw = 0.05
include_gravity = False include_gravity = False
debugflag = True debugflag = False
analyse_condition = False analyse_condition = False
# I/O CONFIG ################################################################# # I/O CONFIG #################################################################
...@@ -385,7 +386,7 @@ if __name__ == '__main__': ...@@ -385,7 +386,7 @@ if __name__ == '__main__':
# parameter=simulation_parameter # parameter=simulation_parameter
# ) # )
LDDsim.join() # LDDsim.join()
if mesh_study: if mesh_study:
simulation_output_dir = processQueue.get() simulation_output_dir = processQueue.get()
hlp.merge_spacetime_errornorms(isRichards=isRichards, hlp.merge_spacetime_errornorms(isRichards=isRichards,
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment