From e1580cabb250b96f29db8b245f26ab073f551af7 Mon Sep 17 00:00:00 2001
From: David <forenkram@gmx.de>
Date: Wed, 17 Jun 2020 13:57:59 +0200
Subject: [PATCH] adjust L for non g example
---
...R-2-patch-realistic-same-intrinsic-perm.py | 24 +++++++++----------
1 file changed, 12 insertions(+), 12 deletions(-)
diff --git a/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-realistic-same-intrinsic-perm.py b/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-realistic-same-intrinsic-perm.py
index 4a0087b..e2a7086 100755
--- a/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-realistic-same-intrinsic-perm.py
+++ b/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-realistic-same-intrinsic-perm.py
@@ -30,7 +30,7 @@ date = datetime.datetime.now()
datestr = date.strftime("%Y-%m-%d")
# Name of the usecase that will be printed during simulation.
-use_case = "TP-R-2-patch-realistic-same-intrinsic-perm"
+use_case = "TP-R-2-patch-realistic-same-intrinsic-perm-different-L"
# The name of this very file. Needed for creating log output.
thisfile = "TP-R-2-patch-realistic-same-intrinsic-perm.py"
@@ -47,7 +47,7 @@ resolutions = {
# 4: 1e-5,
# 8: 1e-5,
# 16: 5e-6,
- 32: 5e-6,
+ 32: 4e-6,
# 64: 2e-6,
# 128: 2e-6,
# 256: 1e-6,
@@ -57,22 +57,22 @@ resolutions = {
# The list is looped over and a simulation is run with t_0 as initial time
# for each element t_0 in starttimes.
starttimes = [0.0]
-timestep_size = 0.01
+timestep_size = 0.001
number_of_timesteps = 800
# LDD scheme parameters ######################################################
-Lw1 = 0.25
-Lnw1 = 0.25
+Lw1 = 0.025
+Lnw1 = 0.025
-Lw2 = 0.5
-Lnw2 = 0.25
+Lw2 = 0.025
+Lnw2 = 0.025
-lambda_w = 40
-lambda_nw = 40
+lambda_w = 4
+lambda_nw = 4
include_gravity = False
debugflag = False
-analyse_condition = False
+analyse_condition = True
# I/O CONFIG #################################################################
# when number_of_timesteps is high, it might take a long time to write all
@@ -198,9 +198,9 @@ isRichards = {
viscosity = {
# subdom_num : viscosity
- 1: {'wetting': 1,
+ 1: {'wetting': 1.0,
'nonwetting': 1/50},
- 2: {'wetting': 1,
+ 2: {'wetting': 1.0,
'nonwetting': 1/50}
}
--
GitLab