uiae

TP-TP-layered-soil-case/mesh_study/TP-TP-layered_soil_pure_dd_mesh_study.py

@@ -34,18 +34,18 @@ mesh_study = True
 resolutions = {
                 # 1: 1e-7,  # h=2
                 2: 1e-7,  # h=1.1180
-                4: 1e-7,  # h=0.5590
-                8: 1e-7,  # h=0.2814
-                16: 1e-7, # h=0.1412
-                32: 1e-7,   # h=0.0706
-                64: 1e-7,
-                128: 1e-7,
-                256: 1e-7
+                # 4: 1e-7,  # h=0.5590
+                # 8: 1e-7,  # h=0.2814
+                # 16: 1e-7, # h=0.1412
+                # 32: 1e-7,   # h=0.0706
+                # 64: 1e-7,
+                # 128: 1e-7,
+                # 256: 1e-7
                 }
 
 ############ GRID #######################
 # mesh_resolution = 20
-timestep_size = 0.0025
+timestep_size = 0.001
 number_of_timesteps = 300
 plot_timestep_every = 2
 # decide how many timesteps you want analysed. Analysed means, that we write out
@@ -56,11 +56,11 @@ starttimes = [0.0, 0.25, 0.5]
 Lw = 0.025 #/timestep_size
 Lnw=Lw
 
-lambda_w = 40
-lambda_nw = 40
+lambda_w = 80
+lambda_nw = 80
 
 include_gravity = True
-debugflag = False
+debugflag = True
 analyse_condition = True
 
 if mesh_study:

TP-one-patch/debug_tests/R-one-patch-const-in-time.py

@@ -353,7 +353,7 @@ cutoff = gaussian/(gaussian + zero_on_shrinking)
 # }
 
 p_e_sym = {
-    0: {'wetting': -6 -(1 + x*x + y*y)),
+    0: {'wetting': -6 -(1 + x*x + y*y),
         'nonwetting': -1 -(sym.sin(3*(1+y)/2*sym.pi)*sym.sin(5*(1+x)/2*sym.pi))**2},
 }