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

update TP-R

parent c3d01fe4
No related branches found
No related tags found
No related merge requests found
......@@ -2,6 +2,8 @@ Rechungen, die noch nicht ausgewertet wurden
- TP-one-patch (purely positive pc)
- TP-TP-layered-soil-case-with-inner-patch
- TP-R-two-patch-test-case
Rechnungen, die neu gestartet werden müssen:
-TP-R-multi-patch-same-wetting-phase-as-RR-zero-nonwetting
......@@ -18,24 +18,24 @@ solver_tol = 5E-7
############ GRID #######################ü
mesh_resolution = 40
timestep_size = 0.00001
timestep_size = 0.000001
number_of_timesteps = 15
# 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 = 10
starttime = 0
Lw = 2.5 #/timestep_size
Lw = 5 #/timestep_size
Lnw=Lw
l_param_w = 40
l_param_nw = 40
l_param_w = 100
l_param_nw = 100
include_gravity = True
debugflag = True
analyse_condition = False
output_string = "./output/before_reimplementing_gravity_term_number_of_timesteps{}_".format(number_of_timesteps)
output_string = "./output/after_reimplementing_gravity_term_number_of_timesteps{}_".format(number_of_timesteps)
##### Domain and Interface ####
# global simulation domain domain
......@@ -325,7 +325,7 @@ t = sym.symbols('t', positive=True)
p_e_sym = {
1: {'wetting': (-5.0 - (1.0 + t*t)*(1.0 + x*x + y*y))}, #*(1-x)**2*(1+x)**2*(1-y)**2},
2: {'wetting': (-5.0 - (1.0 + t*t)*(1.0 + x*x)), #*(1-x)**2*(1+x)**2*(1+y)**2,
'nonwetting': (-1-t*(1.1+y + x**2) - sym.sqrt(2+t**4))*y**2}, #*(1-x)**2*(1+x)**2*(1+y)**2},
'nonwetting': (-1-t*(1.1+y + x**2))*y**2}, #*(1-x)**2*(1+x)**2*(1+y)**2},
} #-y*y*(sym.sin(-2*t+2*x)*sym.sin(1/2*y-1.2*t)) - t*t*x*(0.5-y)*y*(1-x)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment