Skip to content
Snippets Groups Projects
Commit a731fb1e authored by David Seus's avatar David Seus
Browse files

setup new layered soil calculation

parent 4236119b
Branches
Tags
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch-realistic.py
+ 8
− 8
View file @ a731fb1e
...@@ -38,18 +38,18 @@ resolutions = { ...@@ -38,18 +38,18 @@ resolutions = {
# 8: 1e-6, # h=0.2814 # 8: 1e-6, # h=0.2814
# 16: 5e-7, # h=0.1412 # 16: 5e-7, # h=0.1412
# 32: 4e-7, # h=0.0706 # 32: 4e-7, # h=0.0706
64: 7e-7, 64: 1e-6,
# 128: 5e-7 # 128: 5e-7
} }
############ GRID ####################### ############ GRID #######################
# mesh_resolution = 20 # mesh_resolution = 20
timestep_size = 0.005 timestep_size = 0.001
number_of_timesteps = 240 number_of_timesteps = 1100
plot_timestep_every = 1 plot_timestep_every = 1
# decide how many timesteps you want analysed. Analysed means, that we write out # decide how many timesteps you want analysed. Analysed means, that we write out
# subsequent errors of the L-iteration within the timestep. # subsequent errors of the L-iteration within the timestep.
number_of_timesteps_to_analyse = 8 number_of_timesteps_to_analyse = 10
starttime = 0.0 starttime = 0.0
Lw = 0.025 #/timestep_size Lw = 0.025 #/timestep_size
...@@ -655,13 +655,13 @@ p_e_sym = { ...@@ -655,13 +655,13 @@ p_e_sym = {
2: {'wetting': -5.0 - (1.0 + t*t)*(1.0 + x*x + y*y), 2: {'wetting': -5.0 - (1.0 + t*t)*(1.0 + x*x + y*y),
'nonwetting': (-1 -t*(1.1 + y + x**2)) }, 'nonwetting': (-1 -t*(1.1 + y + x**2)) },
3: {'wetting': (-5.0 - (1.0 + t*t)*(1.0 + x*x)), 3: {'wetting': (-5.0 - (1.0 + t*t)*(1.0 + x*x)),
'nonwetting': (-1 -t*(1 + x**2) - sym.sqrt(2+t**2)*(1+y)*y**2) }, 'nonwetting': (-1 -t*(1 + x**2) - sym.sin(2+t**2)**2*(1+y)*y**2) },
4: {'wetting': (-5.0 - (1.0 + t*t)*(1.0 + x*x)), 4: {'wetting': (-5.0 - (1.0 + t*t)*(1.0 + x*x)),
'nonwetting': (-1 -t*(1 + x**2) - sym.sqrt(2+t**2)*(1+y)*y**2) }, 'nonwetting': (-1 -t*(1 + x**2) - sym.sin(2+t**2)**2*(1+y)*y**2) },
5: {'wetting': (-5.0 - (1.0 + t*t)*(1.0 + x*x)), 5: {'wetting': (-5.0 - (1.0 + t*t)*(1.0 + x*x)),
'nonwetting': (-1 -t*(1 + x**2) - sym.sqrt(2+t**2)*(1+y)*y**2) }, 'nonwetting': (-1 -t*(1 + x**2) - sym.sin(2+t**2)**2*(1+y)*y**2) },
6: {'wetting': (-5.0 - (1.0 + t*t)*(1.0 + x*x)), 6: {'wetting': (-5.0 - (1.0 + t*t)*(1.0 + x*x)),
'nonwetting': (-1 -t*(1 + x**2) - sym.sqrt(2+t**2)*(1+y)*y**2) }, 'nonwetting': (-1 -t*(1 + x**2) - sym.sin(2+t**2)**2*(1+y)*y**2) },
# 2: {'wetting': 1.0 - (1.0 + t*t)*(10.0 + x*x + (y-5.0)*(y-5.0)), # 2: {'wetting': 1.0 - (1.0 + t*t)*(10.0 + x*x + (y-5.0)*(y-5.0)),
# 'nonwetting': - 2 - t*(1 + (y-5.0) + x**2)**2 -sym.sqrt(2+t**2)*(1 + (y-5.0))}, # 'nonwetting': - 2 - t*(1 + (y-5.0) + x**2)**2 -sym.sqrt(2+t**2)*(1 + (y-5.0))},
# 3: {'wetting': 1.0 - (1.0 + t*t)*(10.0 + x*x + (y-5.0)*(y-5.0)*3*sym.sin(-2*t+2*x)*sym.sin(1/2*y-1.2*t)), # 3: {'wetting': 1.0 - (1.0 + t*t)*(10.0 + x*x + (y-5.0)*(y-5.0)*3*sym.sin(-2*t+2*x)*sym.sin(1/2*y-1.2*t)),
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment