setup new TPTP example

Two-phase-Two-phase/two-patch/mesh_studies/TP-TP-2-patch-mesh-study.py

@@ -21,32 +21,32 @@ sym.init_printing()
 
 use_case = "TP-TP-2-patch-realistic"
 # solver_tol = 6E-7
-max_iter_num = 1000
+max_iter_num = 300
 FEM_Lagrange_degree = 1
 mesh_study = True
 resolutions = {
-                1: 5e-6,
-                2: 5e-6,
-                4: 5e-6,
-                8: 5e-6,
-                16: 2e-6,
-                32: 2e-6,
-                64: 2e-6,
-                128: 2e-6,
-                #256: 7e-7,
+                1: 1e-6,
+                2: 1e-6,
+                4: 1e-6,
+                8: 1e-6,
+                16: 1e-6,
+                32: 1e-6,
+                64: 1e-6,
+                128: 1e-6,
+                256: 1e-6,
                 }
 
 ############ GRID #######################
 # mesh_resolution = 20
 timestep_size = 0.001
-number_of_timesteps = 500
-plot_timestep_every = 1
+number_of_timesteps = 800
+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 = 5
+number_of_timesteps_to_analyse = 6
 starttimes = [0.0]
 
-Lw = 0.5 #/timestep_size
+Lw = 0.025 #/timestep_size
 Lnw=Lw
 
 lambda_w = 4
@@ -179,10 +179,10 @@ densities = {
 }
 
 intrinsic_permeability = {
-    1: {"wetting": 10e-6,
-        "nonwetting": 10e-6},
-    2: {"wetting": 10e-8,
-        "nonwetting": 10e-8},
+    1: {"wetting": 1,
+        "nonwetting": 1},
+    2: {"wetting": 1,
+        "nonwetting": 1},
 }
 
 gravity_acceleration = 9.81