diff --git a/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-realistic-no-gravity-but-varying-intrinsic-perm.py b/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-realistic-no-gravity-but-varying-intrinsic-perm.py index 43645b17240e9c3f73e7b458db6ba942175d35ea..d51852eb0b811ee78b3aad39c593a54f15cd1508 100755 --- a/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-realistic-no-gravity-but-varying-intrinsic-perm.py +++ b/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-realistic-no-gravity-but-varying-intrinsic-perm.py @@ -36,7 +36,7 @@ thisfile = "TP-R-2-patch-realistic-no-gravity-but-varying-intrinsic-perm.py" # GENERAL SOLVER CONFIG ###################################################### # maximal iteration per timestep -max_iter_num = 250 +max_iter_num = 1000 FEM_Lagrange_degree = 1 # GRID AND MESH STUDY SPECIFICATIONS ######################################### @@ -47,7 +47,7 @@ resolutions = { # 4: 1e-5, # 8: 1e-5, # 16: 5e-6, - 32: 5e-6, + 32: 3e-6, # 64: 2e-6, # 128: 2e-6, # 256: 1e-6, @@ -61,18 +61,18 @@ timestep_size = 0.001 number_of_timesteps = 800 # LDD scheme parameters ###################################################### -Lw1 = 0.25 -Lnw1 = 0.25 +Lw1 = 0.025 +Lnw1 = 0.025 -Lw2 = 0.5 -Lnw2 = 0.25 +Lw2 = 0.05 +Lnw2 = 0.025 -lambda_w = 40 -lambda_nw = 40 +lambda_w = 4 +lambda_nw = 4 include_gravity = False debugflag = False -analyse_condition = True +analyse_condition = False # I/O CONFIG ################################################################# # when number_of_timesteps is high, it might take a long time to write all