set up a few layered soil examples

Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/TP-TP-layered_soil.py

@@ -31,13 +31,13 @@ datestr = date.strftime("%Y-%m-%d")
 
 
 # Name of the usecase that will be printed during simulation.
-use_case = "TP-TP-layered_soil-realistic-same-L"
+use_case = "TP-TP-layered_soil-realistic-not-VanGenuchten"
 # The name of this very file. Needed for creating log output.
 thisfile = "TP-TP-layered_soil.py"
 
 # GENERAL SOLVER CONFIG  ######################################################
 # maximal iteration per timestep
-max_iter_num = 1000
+max_iter_num = 2000
 FEM_Lagrange_degree = 1
 
 # GRID AND MESH STUDY SPECIFICATIONS  #########################################
@@ -48,7 +48,7 @@ resolutions = {
                 # 4: 1e-6,
                 # 8: 1e-5,
                 # 16: 5e-6,
-                32: 2e-6,
+                32: 1e-6,
                 # 64: 2e-6,
                 # 128: 1e-6,
                 # 256: 1e-6,
@@ -57,28 +57,28 @@ resolutions = {
 # starttimes gives a list of starttimes to run the simulation from.
 # The list is looped over and a simulation is run with t_0 as initial time
 #  for each element t_0 in starttimes.
-# starttimes = {0: 0.0, 1: 0.3, 2: 0.6, 3: 0.9}
-starttimes = {0: 0.0}
+starttimes = {0: 0.0, 1: 0.3, 2: 0.6, 3: 0.9}
+# starttimes = {0: 0.0}
 timestep_size = 0.001
-number_of_timesteps = 1500
+number_of_timesteps = 1
 
 
 # LDD scheme parameters  ######################################################
-Lw1 = 0.025
-Lnw1 = 0.025
-Lw2 = 0.025
-Lnw2 = 0.025
-Lw3 = 0.025
-Lnw3 = 0.025
-Lw4 = 0.025
-Lnw4 = 0.025
-
-lambda12_w = 1
-lambda12_nw = 0.5
-lambda23_w = 1
-lambda23_nw = 0.5
-lambda34_w = 1
-lambda34_nw = 0.5
+Lw1 = 0.1
+Lnw1 = 0.1
+Lw2 = 0.1
+Lnw2 = 0.1
+Lw3 = 0.1
+Lnw3 = 0.1
+Lw4 = 0.1
+Lnw4 = 0.1
+
+lambda12_w = 20
+lambda12_nw = 20
+lambda23_w = 20
+lambda23_nw = 20
+lambda34_w = 20
+lambda34_nw = 20
 
 include_gravity = False
 debugflag = False
@@ -88,10 +88,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 = 5
+plot_timestep_every = 1
 # 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 = 10
+number_of_timesteps_to_analyse = 1
 
 # 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.
@@ -177,9 +177,9 @@ gravity_acceleration = 9.81
 # https://www.geotechdata.info/parameter/soil-porosity.html
 # Dict of the form: { subdom_num : porosity }
 porosity = {
-    1: 0.37,  # 0.2,  # Clayey gravels, clayey sandy gravels
+    1: 0.22,  # 0.2,  # Clayey gravels, clayey sandy gravels
     2: 0.22,  # 0.22, # Silty gravels, silty sandy gravels
-    3: 0.2,  # 0.37, # Clayey sands
+    3: 0.22,  # 0.37, # Clayey sands
     4: 0.22,  # 0.2 # Silty or sandy clay
 }
 
@@ -234,11 +234,17 @@ ka_prime = rel_perm_dict["ka_prime"]
 
 # S-pc relation
 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}},
+    1: {"vanGenuchten": {"n": 3, "alpha": 0.1}},
+    2: {"vanGenuchten": {"n": 3, "alpha": 0.1}},
+    3: {"vanGenuchten": {"n": 6, "alpha": 0.1}},
+    4: {"vanGenuchten": {"n": 6, "alpha": 0.1}},
 }
+# Spc_on_subdomains = {
+#     1: {"testSpc": {"index": 1}},
+#     2: {"testSpc": {"index": 1}},
+#     3: {"testSpc": {"index": 2}},
+#     4: {"testSpc": {"index": 2}},
+# }
 Spc = fts.generate_Spc_dicts(Spc_on_subdomains)
 S_pc_sym = Spc["symbolic"]
 S_pc_sym_prime = Spc["prime_symbolic"]
@@ -251,10 +257,10 @@ x, y = sym.symbols('x[0], x[1]')  # needed by UFL
 t = sym.symbols('t', positive=True)
 
 p_e_sym_2patch = {
-    1: {'wetting': -7 - (1+t*t)*(1 + x*x + y*y),
+    1: {'wetting': -7.0 - (1.0 + t*t)*(1.0 + x*x + y*y),
         'nonwetting': -1-t*(1.1 + y + x**2)**2},
     2: {'wetting': -7.0 - (1.0 + t*t)*(1.0 + x*x),
-        'nonwetting': -1-t*(1.1 + x**2)**2 - sym.sqrt(5+t**2)*y**2},
+        'nonwetting': -1-t*(1.1 + x**2)**2},
 }