Skip to content
Snippets Groups Projects
Commit 9ffa5580 authored by David's avatar David
Browse files

setup layered soil like tpr on two domains

parent 7b364c13
No related branches found
No related tags found
No related merge requests found
...@@ -32,8 +32,8 @@ resolutions = { ...@@ -32,8 +32,8 @@ resolutions = {
# 2: 2e-5, # h=1.1180 # 2: 2e-5, # h=1.1180
# 4: 1e-6, # h=0.5590 # 4: 1e-6, # h=0.5590
# 8: 1e-6, # h=0.2814 # 8: 1e-6, # h=0.2814
16: 1e-6, # h=0.1412 # 16: 1e-6, # h=0.1412
# 32: 1e-6, 32: 1e-6,
# 64: 5e-7, # 64: 5e-7,
# 128: 5e-7 # 128: 5e-7
} }
...@@ -41,7 +41,7 @@ resolutions = { ...@@ -41,7 +41,7 @@ resolutions = {
# GRID ####################### # GRID #######################
# mesh_resolution = 20 # mesh_resolution = 20
timestep_size = 0.005 timestep_size = 0.005
number_of_timesteps = 500 number_of_timesteps = 200
plot_timestep_every = 2 plot_timestep_every = 2
# decide how many timesteps you want analysed. Analysed means, that we write # decide how many timesteps you want analysed. Analysed means, that we write
# out subsequent errors of the L-iteration within the timestep. # out subsequent errors of the L-iteration within the timestep.
...@@ -65,7 +65,7 @@ lambda34_w = 40 ...@@ -65,7 +65,7 @@ lambda34_w = 40
lambda34_nw = 40 lambda34_nw = 40
include_gravity = False include_gravity = False
debugflag = True debugflag = False
analyse_condition = False analyse_condition = False
if mesh_study: if mesh_study:
...@@ -340,7 +340,7 @@ lambda_param = { ...@@ -340,7 +340,7 @@ lambda_param = {
intrinsic_permeability = { intrinsic_permeability = {
1: 1, # sand 1: 1, # sand
2: 0.1, # sand, there is a range 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 4: 10e-4
} }
...@@ -592,10 +592,11 @@ x, y = sym.symbols('x[0], x[1]') # needed by UFL ...@@ -592,10 +592,11 @@ x, y = sym.symbols('x[0], x[1]') # needed by UFL
t = sym.symbols('t', positive=True) t = sym.symbols('t', positive=True)
p_e_sym_2patch = { 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}, 'nonwetting': 0.0*t}, # -1-t*(1.1 + y + x**2)**2},
2: {'wetting': -7.0 - (1.0 + t*t)*(1.0 + x*x), 2: {'wetting': -6.0 - (1.0 + t*t)*(1.0 + x*x),
'nonwetting': (-1-t*(1.1 + x**2)**2 - sym.sqrt(5+t**2))*y**2}, '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 = { p_e_sym = {
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment