From b18108a9838a169d6c116d5be74ae2ed06567897 Mon Sep 17 00:00:00 2001
From: David Seus <david.seus@ians.uni-stuttgart.de>
Date: Tue, 7 Jul 2020 17:04:24 +0200
Subject: [PATCH] restart same intrinsic permeability with flattening
parameters
---
...R-2-patch-realistic-same-intrinsic-perm.py | 26 +++++++++----------
1 file changed, 13 insertions(+), 13 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 dcc14cc..1d41c13 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-different-L"
+use_case = "TP-R-2-patch-realistic-same-intrinsic-perm-flattening-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: 3e-6,
# 64: 2e-6,
# 128: 2e-6,
# 256: 1e-6,
@@ -57,21 +57,21 @@ 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.0}
-timestep_size = 0.01
-number_of_timesteps = 8
+timestep_size = 0.005
+number_of_timesteps = 1000
# LDD scheme parameters ######################################################
-Lw1 = 0.025
-Lnw1 = 0.025
+Lw1 = 0.01
+Lnw1 = 0.01
-Lw2 = 0.025
-Lnw2 = 0.025
+Lw2 = 0.005
+Lnw2 = 0.005
-lambda_w = 4
-lambda_nw = 4
+lambda_w = 0.5
+lambda_nw = 0.5
include_gravity = False
-debugflag = True
+debugflag = False
analyse_condition = False
# I/O CONFIG #################################################################
@@ -81,7 +81,7 @@ analyse_condition = False
plot_timestep_every = 4
# Decide how many timesteps you want analysed. Analysed means, that
# subsequent errors of the L-iteration within the timestep are written out.
-number_of_timesteps_to_analyse = 5
+number_of_timesteps_to_analyse = 10
# fine grained control over data to be written to disk in the mesh study case
# as well as for a regular simuation for a fixed grid.
@@ -334,7 +334,7 @@ if __name__ == '__main__':
# parameter=simulation_parameter
# )
- LDDsim.join()
+ # LDDsim.join()
if mesh_study:
simulation_output_dir = processQueue.get()
hlp.merge_spacetime_errornorms(isRichards=isRichards,
--
GitLab