add layered soil-case-with-inner-patch

Rechnungen_auszuwerten.txt

+Rechungen, die noch nicht ausgewertet wurden 
+
+- TP-one-patch (purely positive pc)
+- TP-TP-layered-soil-case-with-inner-patch

TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch.py

@@ -22,26 +22,26 @@ import helpers as hlp
 sym.init_printing()
 
 use_case = "TP-TP-layered-soil-with-inner-patch"
-solver_tol = 5E-7
+solver_tol = 1E-6
 
 ############ GRID #######################ü
-mesh_resolution = 40
+mesh_resolution = 50
 timestep_size = 0.001
 number_of_timesteps = 1500
 # 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 = 10
 starttime = 0
 
 Lw = 0.25  #/timestep_size
 Lnw=Lw
 
-l_param_w = 40
-l_param_nw = 40
+lambda_w = 41
+lambda_nw = 41
 
-include_gravity = False
+include_gravity = True
 debugflag = False
-analyse_condition = True
+analyse_condition = False
 
 output_string = "./output/nondirichlet_number_of_timesteps{}_".format(number_of_timesteps)
 
@@ -370,18 +370,18 @@ L = {
 
 # subdom_num : lambda parameter for the L-scheme
 lambda_param = {
-    1: {'wetting': l_param_w,
-         'nonwetting': l_param_nw},#
-    2: {'wetting': l_param_w,
-         'nonwetting': l_param_nw},#
-    3: {'wetting': l_param_w,
-         'nonwetting': l_param_nw},#
-    4: {'wetting': l_param_w,
-         'nonwetting': l_param_nw},#
-    5: {'wetting': l_param_w,
-         'nonwetting': l_param_nw},#
-    6: {'wetting': l_param_w,
-         'nonwetting': l_param_nw},#
+    1: {'wetting': lambda_w,
+         'nonwetting': lambda_nw},#
+    2: {'wetting': lambda_w,
+         'nonwetting': lambda_nw},#
+    3: {'wetting': lambda_w,
+         'nonwetting': lambda_nw},#
+    4: {'wetting': lambda_w,
+         'nonwetting': lambda_nw},#
+    5: {'wetting': lambda_w,
+         'nonwetting': lambda_nw},#
+    6: {'wetting': lambda_w,
+         'nonwetting': lambda_nw},#
 }
 
 
@@ -595,18 +595,18 @@ t = sym.symbols('t', positive=True)
 
 
 p_e_sym = {
-    1: {'wetting': -3.0 - (1.0 + t*t)*(1.0 + x*x + y*y),
-        'nonwetting': (-1 -t*(1-y - x**2)**2) },
-    2: {'wetting': -3.0 - (1.0 + t*t)*(1.0 + x*x + y*y),
-        'nonwetting': (-1 -t*(1-y - x**2)**2) },
-    3: {'wetting': (-3.0 - (1.0 + t*t)*(1.0 + x*x)),
-        'nonwetting': (-1 -t*(1- x**2)**2 - sym.sqrt(2+t**2)*(1+y)*y**2) },
-    4: {'wetting': (-3.0 - (1.0 + t*t)*(1.0 + x*x)),
-        'nonwetting': (-1 -t*(1- x**2)**2 - sym.sqrt(2+t**2)*(1+y)*y**2) },
-    5: {'wetting': (-3.0 - (1.0 + t*t)*(1.0 + x*x)),
-        'nonwetting': (-1 -t*(1- x**2)**2 - sym.sqrt(2+t**2)*(1+y)*y**2) },
-    6: {'wetting': (-3.0 - (1.0 + t*t)*(1.0 + x*x)),
-        'nonwetting': (-1 -t*(1- x**2)**2 - sym.sqrt(2+t**2)*(1+y)*y**2) },
+    1: {'wetting': -5.0 - (1.0 + t*t)*(1.0 + x*x + y*y),
+        'nonwetting': (-1 -t*(1.1 + y + x**2)) },
+    2: {'wetting': -5.0 - (1.0 + t*t)*(1.0 + x*x + y*y),
+        'nonwetting': (-1 -t*(1.1 + y + x**2)) },
+    3: {'wetting': (-5.0 - (1.0 + t*t)*(1.0 + x*x)),
+        'nonwetting': (-1 -t*(1 + x**2) - sym.sqrt(2+t**2)*(1+y)*y**2) },
+    4: {'wetting': (-5.0 - (1.0 + t*t)*(1.0 + x*x)),
+        'nonwetting': (-1 -t*(1 + x**2) - sym.sqrt(2+t**2)*(1+y)*y**2) },
+    5: {'wetting': (-5.0 - (1.0 + t*t)*(1.0 + x*x)),
+        'nonwetting': (-1 -t*(1 + x**2) - sym.sqrt(2+t**2)*(1+y)*y**2) },
+    6: {'wetting': (-5.0 - (1.0 + t*t)*(1.0 + x*x)),
+        'nonwetting': (-1 -t*(1 + x**2) - sym.sqrt(2+t**2)*(1+y)*y**2) },
     # 2: {'wetting': 1.0 - (1.0 + t*t)*(10.0 + x*x + (y-5.0)*(y-5.0)),
     #     'nonwetting': - 2 - t*(1 + (y-5.0) + x**2)**2 -sym.sqrt(2+t**2)*(1 + (y-5.0))},
     # 3: {'wetting': 1.0 - (1.0 + t*t)*(10.0 + x*x + (y-5.0)*(y-5.0)*3*sym.sin(-2*t+2*x)*sym.sin(1/2*y-1.2*t)),