setup another gravity example with better parameters

Usecases/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil-g-but-same-perm.py

@@ -49,7 +49,7 @@ datestr = date.strftime("%Y-%m-%d")
 
 
 # Name of the usecase that will be printed during simulation.
-use_case = "TP-R-layered-soil-realistic-same-perm-vanG-Mualem-g-test"
+use_case = "TP-R-layered-soil-realistic-same-perm-vanG-Mualem"
 # The name of this very file. Needed for creating log output.
 thisfile = "TP-R-layered_soil-g-but-same-perm.py"
 
@@ -65,8 +65,8 @@ resolutions = {
                 # 2: 1e-6,
                 # 4: 1e-6,
                 # 8: 1e-5,
-                16: 5e-6,
-                # 32: 4e-6,
+                # 16: 5e-6,
+                32: 4e-5,
                 # 64: 2e-6,
                 # 128: 1e-6,
                 # 256: 1e-6,
@@ -78,32 +78,32 @@ resolutions = {
 # that this initial timestep should be counted as.
 # Format:
 # timestep_num: initial_time
-# starttimes = {0: 0.0}
-starttimes = {
-    0: 0.0,
-    1: 0.3,
-    2: 0.6,
-    3: 0.9,
-    }
-timestep_size = 0.0001
-number_of_timesteps = 1
+starttimes = {0: 0.0}
+# starttimes = {
+#     0: 0.0,
+#     1: 0.3,
+#     2: 0.6,
+#     3: 0.9,
+#     }
+timestep_size = 0.001
+number_of_timesteps = 1500
 
 # LDD scheme parameters  ######################################################
-Lw1 = 0.025  # /timestep_size
+Lw1 = 0.1  # /timestep_size
 Lnw1 = Lw1
-Lw2 = 0.005  # /timestep_size
+Lw2 = 0.9  # /timestep_size
 Lnw2 = Lw2
-Lw3 = 0.025  # /timestep_size
+Lw3 = 0.9  # /timestep_size
 Lnw3 = 0.025
-Lw4 = 0.025  # /timestep_size
+Lw4 = 0.08  # /timestep_size
 Lnw4 = 0.025
 
-lambda12_w = 1
-lambda12_nw = 1
-lambda23_w = 1
+lambda12_w = 0.75
+lambda12_nw = 2
+lambda23_w = 0.3
 lambda23_nw = 1
-lambda34_w = 1
-lambda34_nw = 1
+lambda34_w = 0.4
+lambda34_nw = 2
 
 include_gravity = True
 debugflag = False
@@ -113,10 +113,10 @@ analyse_condition = False
 # when number_of_timesteps is high, it might take a long time to write all
 # timesteps to disk. Therefore, you can choose to only write data of every
 # plot_timestep_every timestep to disk.
-plot_timestep_every = 1
+plot_timestep_every = 5
 # Decide how many timesteps you want analysed. Analysed means, that
 # subsequent errors of the L-iteration within the timestep are written out.
-number_of_timesteps_to_analyse = 2
+number_of_timesteps_to_analyse = 10
 
 # fine grained control over data to be written to disk in the mesh study case
 # as well as for a regular simuation for a fixed grid.
@@ -269,10 +269,10 @@ rel_perm_definition = {
         "nonwetting": {"vanGenuchtenMualem": {"n": 3}} },
     2: {"wetting": {"vanGenuchtenMualem": {"n": 3}},
         "nonwetting": {"vanGenuchtenMualem": {"n": 3}} },
-    3: {"wetting": {"vanGenuchtenMualem": {"n": 6}},
-        "nonwetting": {"vanGenuchtenMualem": {"n": 6}} },
-    4: {"wetting": {"vanGenuchtenMualem": {"n": 6}},
-        "nonwetting": {"vanGenuchtenMualem": {"n": 6}} },
+    3: {"wetting": {"vanGenuchtenMualem": {"n": 5}},
+        "nonwetting": {"vanGenuchtenMualem": {"n": 5}} },
+    4: {"wetting": {"vanGenuchtenMualem": {"n": 5}},
+        "nonwetting": {"vanGenuchtenMualem": {"n": 5}} },
 }
 
 rel_perm_dict = fts.generate_relative_permeability_dicts(rel_perm_definition)
@@ -283,8 +283,8 @@ ka_prime = rel_perm_dict["ka_prime"]
 Spc_on_subdomains = {
     1: {"vanGenuchten": {"n": 3, "alpha": 0.001}},
     2: {"vanGenuchten": {"n": 3, "alpha": 0.001}},
-    3: {"vanGenuchten": {"n": 6, "alpha": 0.001}},
-    4: {"vanGenuchten": {"n": 6, "alpha": 0.001}},
+    3: {"vanGenuchten": {"n": 5, "alpha": 0.001}},
+    4: {"vanGenuchten": {"n": 5, "alpha": 0.001}},
 }
 Spc = fts.generate_Spc_dicts(Spc_on_subdomains)
 S_pc_sym = Spc["symbolic"]