From 125ed9136ea656ace71f40af6d94099d36cc7004 Mon Sep 17 00:00:00 2001
From: David <forenkram@gmx.de>
Date: Thu, 15 Oct 2020 23:20:55 +0200
Subject: [PATCH] set up different intrinsic TPR
---
.../layered_soil/TP-R-layered_soil.py | 34 +++++++++----------
1 file changed, 17 insertions(+), 17 deletions(-)
diff --git a/Usecases/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil.py b/Usecases/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil.py
index d6edf0f..c58def0 100755
--- a/Usecases/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil.py
+++ b/Usecases/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil.py
@@ -49,13 +49,13 @@ datestr = date.strftime("%Y-%m-%d")
# Name of the usecase that will be printed during simulation.
-use_case = "TP-R-layered-soil-realistic"
+use_case = "TP-R-layered-soil-realistic-different-intrinsic-vanG-Mualem"
# The name of this very file. Needed for creating log output.
thisfile = "TP-R-layered_soil.py"
# GENERAL SOLVER CONFIG ######################################################
# maximal iteration per timestep
-max_iter_num = 10
+max_iter_num = 2000
FEM_Lagrange_degree = 1
# GRID AND MESH STUDY SPECIFICATIONS #########################################
@@ -66,7 +66,7 @@ resolutions = {
# 4: 1e-6,
# 8: 1e-5,
# 16: 5e-6,
- 32: 5e-6,
+ 32: 3e-6,
# 64: 2e-6,
# 128: 1e-6,
# 256: 1e-6,
@@ -77,7 +77,7 @@ resolutions = {
# for each element t_0 in starttimes.
starttimes = {0: 0.0}
timestep_size = 0.001
-number_of_timesteps = 5
+number_of_timesteps = 1500
# LDD scheme parameters ######################################################
@@ -85,20 +85,20 @@ Lw1 = 0.025 # /timestep_size
Lnw1 = Lw1
Lw2 = 0.025 # /timestep_size
Lnw2 = Lw2
-Lw3 = 0.025 # /timestep_size
+Lw3 = 0.0025 # /timestep_size
Lnw3 = Lw3
-Lw4 = 0.025 # /timestep_size
+Lw4 = 0.0025 # /timestep_size
Lnw4 = Lw4
-lambda12_w = 40
-lambda12_nw = 40
-lambda23_w = 40
-lambda23_nw = 40
-lambda34_w = 40
-lambda34_nw = 40
+lambda12_w = 4
+lambda12_nw = 4
+lambda23_w = 4
+lambda23_nw = 4
+lambda34_w = 4
+lambda34_nw = 4
include_gravity = False
-debugflag = True
+debugflag = False
analyse_condition = False
# I/O CONFIG #################################################################
@@ -215,8 +215,8 @@ gravity_acceleration = 9.81
porosity = {
1: 0.37, # 0.2, # Clayey gravels, clayey sandy gravels
2: 0.22, # 0.22, # Silty gravels, silty sandy gravels
- 3: 0.2, # 0.37, # Clayey sands
- 4: 0.22, # 0.2 # Silty or sandy clay
+ 3: 0.02, # 0.37, # Clayey sands
+ 4: 0.022, # 0.2 # Silty or sandy clay
}
# subdom_num : subdomain L for L-scheme
@@ -248,8 +248,8 @@ lambda_param = {
intrinsic_permeability = {
1: 0.01, # sand
2: 0.01, # sand, there is a range
- 3: 0.01, #10e-2, # clay has a range
- 4: 0.01, #10e-3
+ 3: 0.001, #10e-2, # clay has a range
+ 4: 0.001, #10e-3
}
# relative permeabilties
--
GitLab