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 2f5c4ab2db754d047f1817df56da6b4087ea73b4..78eb3980abd5789359e37cad4dc5f9608fe4e294 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 @@ -32,8 +32,8 @@ resolutions = { # 2: 2e-5, # h=1.1180 # 4: 1e-6, # h=0.5590 # 8: 1e-6, # h=0.2814 - 16: 1e-6, # h=0.1412 - # 32: 1e-6, + # 16: 1e-6, # h=0.1412 + 32: 1e-6, # 64: 5e-7, # 128: 5e-7 } @@ -41,7 +41,7 @@ resolutions = { # GRID ####################### # mesh_resolution = 20 timestep_size = 0.005 -number_of_timesteps = 500 +number_of_timesteps = 200 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. @@ -65,7 +65,7 @@ lambda34_w = 40 lambda34_nw = 40 include_gravity = False -debugflag = True +debugflag = False analyse_condition = False if mesh_study: @@ -340,7 +340,7 @@ lambda_param = { intrinsic_permeability = { 1: 1, # sand 2: 0.1, # sand, there is a range - 3: 10e-3, # clay has a range + 3: 10e-2, # clay has a range 4: 10e-4 } @@ -592,10 +592,11 @@ x, y = sym.symbols('x[0], x[1]') # needed by UFL t = sym.symbols('t', positive=True) p_e_sym_2patch = { - 1: {'wetting': -7 - (1+t*t)*(1 + x*x + y*y), + 1: {'wetting': -6 - (1+t*t)*(1 + x*x + y*y), 'nonwetting': 0.0*t}, # -1-t*(1.1 + y + x**2)**2}, - 2: {'wetting': -7.0 - (1.0 + t*t)*(1.0 + x*x), - 'nonwetting': (-1-t*(1.1 + x**2)**2 - sym.sqrt(5+t**2))*y**2}, + 2: {'wetting': -6.0 - (1.0 + t*t)*(1.0 + x*x), + 'nonwetting': (-1-t*(1.1+y + x**2))*y**2}, + # 'nonwetting': (-1-t*(1.1 + x**2)**2 - sym.sqrt(5+t**2))*y**2}, } p_e_sym = {