diff --git a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/mesh_study/TP-TP-layered_soil_pure_dd_mesh_study.py b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/mesh_study/TP-TP-layered_soil_pure_dd_mesh_study.py index c8d4c35b13ff1739e6e29a4874aa7705e2d0f3e5..9b07c00cf0bebf0d55e6606d15009a998ab3f9ad 100755 --- a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/mesh_study/TP-TP-layered_soil_pure_dd_mesh_study.py +++ b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/mesh_study/TP-TP-layered_soil_pure_dd_mesh_study.py @@ -32,35 +32,35 @@ max_iter_num = 1000 FEM_Lagrange_degree = 1 mesh_study = True resolutions = { - # 1: 1e-7, # h=2 - 2: 1e-7, # h=1.1180 - # 4: 1e-7, # h=0.5590 - # 8: 1e-7, # h=0.2814 - # 16: 1e-7, # h=0.1412 - # 32: 1e-7, # h=0.0706 - # 64: 1e-7, - # 128: 1e-7, + 1: 1e-6, # h=2 + 2: 1e-6, # h=1.1180 + 4: 1e-6, # h=0.5590 + 8: 1e-6, # h=0.2814 + 16: 1e-6, # h=0.1412 + 32: 1e-6, # h=0.0706 + 64: 1e-6, + 128: 1e-6, # 256: 1e-7 } ############ GRID ####################### # mesh_resolution = 20 -timestep_size = 0.001 -number_of_timesteps = 300 +timestep_size = 0.005 +number_of_timesteps = 100 plot_timestep_every = 2 # decide how many timesteps you want analysed. Analysed means, that we write out # subsequent errors of the L-iteration within the timestep. number_of_timesteps_to_analyse = 0 -starttimes = [0.0, 0.25, 0.5] +starttimes = [0.0] Lw = 0.025 #/timestep_size Lnw=Lw -lambda_w = 80 -lambda_nw = 80 +lambda_w = 8 +lambda_nw = 8 -include_gravity = True -debugflag = True +include_gravity = False +debugflag = False analyse_condition = True if mesh_study: @@ -76,10 +76,10 @@ if mesh_study: 'space_errornorms': True, 'meshes_and_markers': True, 'L_iterations_per_timestep': False, - 'solutions': False, - 'absolute_differences': False, + 'solutions': True, + 'absolute_differences': True, 'condition_numbers': analyse_condition, - 'subsequent_errors': False + 'subsequent_errors': True } else: write_to_file = { @@ -262,7 +262,7 @@ densities = { 'nonwetting':1}, #1.225}} } -gravity_acceleration = 9.81 +gravity_acceleration = 1 # porosities taken from # https://www.geotechdata.info/parameter/soil-porosity.html # Dict of the form: { subdom_num : porosity }