From 12a0ba4dec87abdfd2bf524ca3e23e56abf1084a Mon Sep 17 00:00:00 2001
From: David Seus <david.seus@ians.uni-stuttgart.de>
Date: Fri, 16 Aug 2019 17:09:23 +0200
Subject: [PATCH] update TP-R
---
Rechnungen_auszuwerten.txt | 2 ++
TP-R-two-patch-test-case/TP-R-2-patch-test.py | 12 ++++++------
2 files changed, 8 insertions(+), 6 deletions(-)
diff --git a/Rechnungen_auszuwerten.txt b/Rechnungen_auszuwerten.txt
index b4fed71..486461d 100644
--- a/Rechnungen_auszuwerten.txt
+++ b/Rechnungen_auszuwerten.txt
@@ -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
diff --git a/TP-R-two-patch-test-case/TP-R-2-patch-test.py b/TP-R-two-patch-test-case/TP-R-2-patch-test.py
index f0710f0..168870d 100755
--- a/TP-R-two-patch-test-case/TP-R-2-patch-test.py
+++ b/TP-R-two-patch-test-case/TP-R-2-patch-test.py
@@ -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)
--
GitLab