set up different intrinsic TPR

Usecases/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil.py

@@ -49,13 +49,13 @@ datestr = date.strftime("%Y-%m-%d")
 
 
 # Name of the usecase that will be printed during simulation.
-use_case = "TP-R-layered-soil-realistic"
+use_case = "TP-R-layered-soil-realistic-different-intrinsic-vanG-Mualem"
 # The name of this very file. Needed for creating log output.
 thisfile = "TP-R-layered_soil.py"
 
 # GENERAL SOLVER CONFIG  ######################################################
 # maximal iteration per timestep
-max_iter_num = 10
+max_iter_num = 2000
 FEM_Lagrange_degree = 1
 
 # GRID AND MESH STUDY SPECIFICATIONS  #########################################
@@ -66,7 +66,7 @@ resolutions = {
                 # 4: 1e-6,
                 # 8: 1e-5,
                 # 16: 5e-6,
-                32: 5e-6,
+                32: 3e-6,
                 # 64: 2e-6,
                 # 128: 1e-6,
                 # 256: 1e-6,
@@ -77,7 +77,7 @@ resolutions = {
 #  for each element t_0 in starttimes.
 starttimes = {0: 0.0}
 timestep_size = 0.001
-number_of_timesteps = 5
+number_of_timesteps = 1500
 
 
 # LDD scheme parameters  ######################################################
@@ -85,20 +85,20 @@ Lw1 = 0.025  # /timestep_size
 Lnw1 = Lw1
 Lw2 = 0.025  # /timestep_size
 Lnw2 = Lw2
-Lw3 = 0.025  # /timestep_size
+Lw3 = 0.0025  # /timestep_size
 Lnw3 = Lw3
-Lw4 = 0.025  # /timestep_size
+Lw4 = 0.0025  # /timestep_size
 Lnw4 = Lw4
 
-lambda12_w = 40
-lambda12_nw = 40
-lambda23_w = 40
-lambda23_nw = 40
-lambda34_w = 40
-lambda34_nw = 40
+lambda12_w = 4
+lambda12_nw = 4
+lambda23_w = 4
+lambda23_nw = 4
+lambda34_w = 4
+lambda34_nw = 4
 
 include_gravity = False
-debugflag = True
+debugflag = False
 analyse_condition = False
 
 # I/O CONFIG  #################################################################
@@ -215,8 +215,8 @@ gravity_acceleration = 9.81
 porosity = {
     1: 0.37,  # 0.2,  # Clayey gravels, clayey sandy gravels
     2: 0.22,  # 0.22, # Silty gravels, silty sandy gravels
-    3: 0.2,  # 0.37, # Clayey sands
-    4: 0.22,  # 0.2 # Silty or sandy clay
+    3: 0.02,  # 0.37, # Clayey sands
+    4: 0.022,  # 0.2 # Silty or sandy clay
 }
 
 # subdom_num : subdomain L for L-scheme
@@ -248,8 +248,8 @@ lambda_param = {
 intrinsic_permeability = {
     1: 0.01,  # sand
     2: 0.01,  # sand, there is a range
-    3: 0.01,  #10e-2,  # clay has a range
-    4: 0.01,  #10e-3
+    3: 0.001,  #10e-2,  # clay has a range
+    4: 0.001,  #10e-3
 }
 
 # relative permeabilties