From a339015f55e0db1680bff747e01125caf516211b Mon Sep 17 00:00:00 2001
From: David Seus <david.seus@ians.uni-stuttgart.de>
Date: Fri, 10 Jul 2020 13:21:41 +0200
Subject: [PATCH] setup layered soil examples
---
.../TP-TP-layered_soil-different-intrinsic.py | 12 ++---
.../TP-TP-layered_soil.py | 45 ++++++++++---------
2 files changed, 29 insertions(+), 28 deletions(-)
diff --git a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/TP-TP-layered_soil-different-intrinsic.py b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/TP-TP-layered_soil-different-intrinsic.py
index fd21276..bfb754e 100755
--- a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/TP-TP-layered_soil-different-intrinsic.py
+++ b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/TP-TP-layered_soil-different-intrinsic.py
@@ -47,8 +47,8 @@ resolutions = {
# 2: 1e-6,
# 4: 1e-6,
# 8: 1e-5,
- # 16: 5e-6,
- 32: 2e-6,
+ 16: 5e-6,
+ # 32: 2e-6,
# 64: 2e-6,
# 128: 1e-6,
# 256: 1e-6,
@@ -65,12 +65,12 @@ number_of_timesteps = 1
# LDD scheme parameters ######################################################
Lw1 = 0.007 # /timestep_size
-Lnw1 = 0.005
+Lnw1 = 0.001
Lw2 = 0.007 # /timestep_size
-Lnw2 = 0.005
-Lw3 = 0.0007 # /timestep_size
+Lnw2 = 0.001
+Lw3 = 0.007 # /timestep_size
Lnw3 = 0.0005
-Lw4 = 0.0007 # /timestep_size
+Lw4 = 0.007 # /timestep_size
Lnw4 = 0.0005
lambda12_w = 0.5
diff --git a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/TP-TP-layered_soil.py b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/TP-TP-layered_soil.py
index a176766..64ccc3c 100755
--- a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/TP-TP-layered_soil.py
+++ b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/TP-TP-layered_soil.py
@@ -37,7 +37,7 @@ thisfile = "TP-TP-layered_soil.py"
# GENERAL SOLVER CONFIG ######################################################
# maximal iteration per timestep
-max_iter_num = 10
+max_iter_num = 1000
FEM_Lagrange_degree = 1
# GRID AND MESH STUDY SPECIFICATIONS #########################################
@@ -47,8 +47,8 @@ resolutions = {
# 2: 1e-6,
# 4: 1e-6,
# 8: 1e-5,
- # 16: 5e-6,
- 32: 5e-6,
+ 16: 5e-6,
+ # 32: 2e-6,
# 64: 2e-6,
# 128: 1e-6,
# 256: 1e-6,
@@ -57,30 +57,31 @@ resolutions = {
# starttimes gives a list of starttimes to run the simulation from.
# 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}
+starttimes = {0: 0.0, 1: 0.3, 2: 0.6, 3: 0.9}
+# starttimes = {0: 0.0}
timestep_size = 0.001
-number_of_timesteps = 5
+number_of_timesteps = 1
# LDD scheme parameters ######################################################
-Lw1 = 0.025 # /timestep_size
-Lnw1 = Lw1
-Lw2 = 0.025 # /timestep_size
-Lnw2 = Lw2
-Lw3 = 0.025 # /timestep_size
-Lnw3 = Lw3
-Lw4 = 0.025 # /timestep_size
-Lnw4 = Lw4
-
-lambda12_w = 40
-lambda12_nw = 40
-lambda23_w = 40
-lambda23_nw = 40
-lambda34_w = 40
-lambda34_nw = 40
+Lw1 = 0.007
+Lnw1 = 0.005
+Lw2 = 0.007
+Lnw2 = 0.005
+Lw3 = 0.007
+Lnw3 = 0.005
+Lw4 = 0.007
+Lnw4 = 0.005
+
+lambda12_w = 0.05
+lambda12_nw = 0.05
+lambda23_w = 0.05
+lambda23_nw = 0.05
+lambda34_w = 0.05
+lambda34_nw = 0.05
include_gravity = False
-debugflag = True
+debugflag = False
analyse_condition = False
# I/O CONFIG #################################################################
@@ -385,7 +386,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