diff --git a/Two-phase-Two-phase/two-patch/TP-TP-2-patch-test-case/TP-TP-2-patch-test.py b/Two-phase-Two-phase/two-patch/TP-TP-2-patch-test-case/TP-TP-2-patch-test.py index 2957da4b6c95167c1fbbbb377a0faba3d75d2ab6..d892719b99c484791ba52f70344f6e38ebe2de52 100755 --- a/Two-phase-Two-phase/two-patch/TP-TP-2-patch-test-case/TP-TP-2-patch-test.py +++ b/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: 2e-6, + 32: 1e-6, # 64: 5e-7, # 128: 5e-7 } @@ -39,7 +39,7 @@ resolutions = { ############ GRID ####################### # mesh_resolution = 20 timestep_size = 0.001 -number_of_timesteps = 1300 +number_of_timesteps = 1500 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. @@ -49,8 +49,8 @@ starttime = 0.0 Lw = 0.05 #/timestep_size Lnw=Lw -lambda_w = 0.4 -lambda_nw = 0.4 +lambda_w = 4 +lambda_nw = 4 include_gravity = True debugflag = False @@ -184,7 +184,7 @@ densities = { 'nonwetting': 1}, #1225}, } -gravity_acceleration = 9.81 +gravity_acceleration = 1#9.81 L = {# # subdom_num : subdomain L for L-scheme