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

reduce timestep size to see if solver behaviour is better

parent a8fc3a9a
No related branches found
No related tags found
No related merge requests found
...@@ -102,19 +102,19 @@ lambda56_w = 4 ...@@ -102,19 +102,19 @@ lambda56_w = 4
lambda56_nw = 4 lambda56_nw = 4
include_gravity = False include_gravity = False
debugflag = True debugflag = False
analyse_condition = False analyse_condition = False
if mesh_study: # if mesh_study:
output_string = "./output/{}-{}_timesteps{}_P{}".format( output_string = "./output/{}-{}_timesteps{}_P{}".format(
datestr, use_case, number_of_timesteps, FEM_Lagrange_degree datestr, use_case, number_of_timesteps, FEM_Lagrange_degree
) )
else: # else:
for tol in resolutions.values(): # # for tol in resolutions.values():
solver_tol = tol # # solver_tol = tol
output_string = "./output/{}-{}_timesteps{}_P{}_solver_tol{}".format( # output_string = "./output/{}-{}_timesteps{}_P{}_solver_tol{}".format(
datestr, use_case, number_of_timesteps, FEM_Lagrange_degree, solver_tol # datestr, use_case, number_of_timesteps, FEM_Lagrange_degree, solver_tol
) # )
# toggle what should be written to files # toggle what should be written to files
......
...@@ -30,26 +30,27 @@ sym.init_printing() ...@@ -30,26 +30,27 @@ sym.init_printing()
# solver_tol = 6E-7 # solver_tol = 6E-7
use_case = "TP-R-layered-soil-with-inner-patch-all-params-one" use_case = "TP-R-layered-soil-with-inner-patch-all-params-one"
# name of this very file. Needed for log output. # name of this very file. Needed for log output.
thisfile = "TP-R-layered_soil_with_inner_patch-all-params-one-mesh-study.py" thisfile = "TP-R-layered_soil_with_inner_patch-all-params-one.py"
max_iter_num = 300 max_iter_num = 150
FEM_Lagrange_degree = 1 FEM_Lagrange_degree = 1
mesh_study = True mesh_study = False
resolutions = { resolutions = {
1: 5e-6, # h=2 # 1: 2e-6, # h=2
2: 5e-6, # h=1.1180 # 2: 2e-6, # h=1.1180
4: 5e-6, # h=0.5590 # 4: 2e-6, # h=0.5590
8: 5e-6, # h=0.2814 # 8: 2e-6, # h=0.2814
16: 5e-6, # h=0.1412 # 16: 7e-6, # h=0.1412
32: 5e-6, 32: 5e-6,
64: 5e-6, 64: 2e-6,
128: 5e-6 128: 2e-6,
256: 2e-6
} }
# GRID ####################### # GRID #######################
# mesh_resolution = 20 # mesh_resolution = 20
timestep_size = 0.001 timestep_size = 0.00025
number_of_timesteps = 1000 number_of_timesteps = 4000
plot_timestep_every = 4 plot_timestep_every = 4
# decide how many timesteps you want analysed. Analysed means, that we write # decide how many timesteps you want analysed. Analysed means, that we write
# out subsequent errors of the L-iteration within the timestep. # out subsequent errors of the L-iteration within the timestep.
...@@ -102,19 +103,19 @@ lambda56_w = 4 ...@@ -102,19 +103,19 @@ lambda56_w = 4
lambda56_nw = 4 lambda56_nw = 4
include_gravity = False include_gravity = False
debugflag = True debugflag = False
analyse_condition = True analyse_condition = False
if mesh_study: # if mesh_study:
output_string = "./output/{}-{}_timesteps{}_P{}".format( output_string = "./output/{}-{}_timesteps{}_P{}".format(
datestr, use_case, number_of_timesteps, FEM_Lagrange_degree datestr, use_case, number_of_timesteps, FEM_Lagrange_degree
) )
else: # else:
for tol in resolutions.values(): # # for tol in resolutions.values():
solver_tol = tol # # solver_tol = tol
output_string = "./output/{}-{}_timesteps{}_P{}_solver_tol{}".format( # output_string = "./output/{}-{}_timesteps{}_P{}_solver_tol{}".format(
datestr, use_case, number_of_timesteps, FEM_Lagrange_degree, solver_tol # datestr, use_case, number_of_timesteps, FEM_Lagrange_degree, solver_tol
) # )
# toggle what should be written to files # toggle what should be written to files
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment