From 7b364c134b146ddd1af324e71ac612b486b31936 Mon Sep 17 00:00:00 2001 From: David <forenkram@gmx.de> Date: Sun, 7 Jun 2020 20:16:33 +0200 Subject: [PATCH] set up realistic example --- .../layered_soil/TP-R-layered_soil.py | 50 +++++++++++-------- 1 file changed, 28 insertions(+), 22 deletions(-) diff --git a/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil.py b/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil.py index 9cab276..2f5c4ab 100755 --- a/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil.py +++ b/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil.py @@ -23,16 +23,16 @@ datestr = date.strftime("%Y-%m-%d") # init sympy session sym.init_printing() # solver_tol = 6E-7 -use_case = "TP-R-layered-soil-realistic" -max_iter_num = 500 +use_case = "TP-R-layered-soil-realistic-new-lambda" +max_iter_num = 600 FEM_Lagrange_degree = 1 mesh_study = False resolutions = { # 1: 1e-7, # h=2 # 2: 2e-5, # h=1.1180 # 4: 1e-6, # h=0.5590 - 8: 1e-6, # h=0.2814 - # 16: 1e-6, # h=0.1412 + # 8: 1e-6, # h=0.2814 + 16: 1e-6, # h=0.1412 # 32: 1e-6, # 64: 5e-7, # 128: 5e-7 @@ -40,23 +40,29 @@ resolutions = { # GRID ####################### # mesh_resolution = 20 -timestep_size = 0.001 -number_of_timesteps = 5 -plot_timestep_every = 1 +timestep_size = 0.005 +number_of_timesteps = 500 +plot_timestep_every = 2 # 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 = 2 +number_of_timesteps_to_analyse = 5 starttimes = [0.0] -Lw = 0.025 # /timestep_size -Lnw = Lw +Lw1 = 0.25 # /timestep_size +Lnw1 = Lw1 +Lw2 = 0.025 # /timestep_size +Lnw2 = Lw2 +Lw3 = 0.025 # /timestep_size +Lnw3 = Lw3 +Lw4 = 0.025 # /timestep_size +Lnw4 = Lw4 lambda12_w = 4 lambda12_nw = 4 -lambda23_w = 4 -lambda23_nw = 4 -lambda34_w = 4 -lambda34_nw = 4 +lambda23_w = 40 +lambda23_nw = 40 +lambda34_w = 40 +lambda34_nw = 40 include_gravity = False debugflag = True @@ -298,14 +304,14 @@ porosity = { # subdom_num : subdomain L for L-scheme L = { - 1: {'wetting': Lw, - 'nonwetting': Lnw}, - 2: {'wetting': Lw, - 'nonwetting': Lnw}, - 3: {'wetting': Lw, - 'nonwetting': Lnw}, - 4: {'wetting': Lw, - 'nonwetting': Lnw} + 1: {'wetting': Lw1, + 'nonwetting': Lnw1}, + 2: {'wetting': Lw2, + 'nonwetting': Lnw2}, + 3: {'wetting': Lw3, + 'nonwetting': Lnw3}, + 4: {'wetting': Lw4, + 'nonwetting': Lnw4} } # interface_num : lambda parameter for the L-scheme on that interface. -- GitLab