set up tp example with gravity non dd

Two-phase-Two-phase/two-patch/TP-TP-2-patch-pure-dd-avoid-interface-at-origin/TP-TP-2-patch-pure-dd-horizontal-interface-avoiding-origin.py

@@ -19,7 +19,7 @@ datestr = date.strftime("%Y-%m-%d")
 # init sympy session
 sym.init_printing()
 
-use_case = "TP-TP-2-patch-really-pure-dd-horizontal-interface-avoding-origin"
+use_case = "TP-TP-2-patch-really-pure-dd-horizontal-interface-avoding-origin-finer-dt"
 # solver_tol = 5E-7
 max_iter_num = 1000
 FEM_Lagrange_degree = 1

Two-phase-Two-phase/two-patch/TP-TP-2-patch-test-case/TP-TP-2-patch-test.py

@@ -30,7 +30,7 @@ resolutions = {
                 # 4: 1e-6,  # h=0.5590
                 # 8: 1e-6,  # h=0.2814
                 # 16: 5e-7, # h=0.1412
-                32: 5e-8,
+                32: 2e-6,
                 # 64: 5e-7,
                 # 128: 5e-7
                 }
@@ -39,21 +39,21 @@ resolutions = {
 ############ GRID #######################
 # mesh_resolution = 20
 timestep_size = 0.001
-number_of_timesteps = 10
-plot_timestep_every = 1
+number_of_timesteps = 1300
+plot_timestep_every = 5
 # 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 = 1
+number_of_timesteps_to_analyse = 5
 starttime = 0.0
 
-Lw = 0.03 #/timestep_size
+Lw = 0.05 #/timestep_size
 Lnw=Lw
 
 lambda_w = 0.4
 lambda_nw = 0.4
 
-include_gravity = False
-debugflag = True
+include_gravity = True
+debugflag = False
 analyse_condition = True
 
 if mesh_study:
@@ -386,7 +386,7 @@ p_e_sym = {
     1: {'wetting': (-6 - (1+t*t)*(1 + x*x + y*y)),  #*cutoff,
         'nonwetting': (-1 -t*(1.1+ y*y))},  #*(sym.sin((1+y)/2*sym.pi)*sym.sin((1+x)/2*sym.pi))**2},
     2: {'wetting': (-6.0 - (1.0 + t*t)*(1.0 + x*x)),  #*(sym.sin((1+y)/2*sym.pi)*sym.sin((1+x)/2*sym.pi))**2,
-        'nonwetting': (-1 -t*(1.1 + y*y) - sym.sqrt(2+t**2)*(1.1+y)**2*y**2)},  #*(sym.sin((1+y)/2*sym.pi)*sym.sin((1+x)/2*sym.pi))**2},
+        'nonwetting': (-1 -t*(1.1 + y*y) - sym.sin((x*y-0.5*t)*y**2)**2)},  #*(sym.sin((1+y)/2*sym.pi)*sym.sin((1+x)/2*sym.pi))**2},
 }