From 9ffa5580dd3173011c0e116ad88be9b9ff6dc99b Mon Sep 17 00:00:00 2001
From: David <forenkram@gmx.de>
Date: Sun, 7 Jun 2020 21:21:34 +0200
Subject: [PATCH] setup layered soil like tpr on two domains
---
.../layered_soil/TP-R-layered_soil.py | 17 +++++++++--------
1 file changed, 9 insertions(+), 8 deletions(-)
diff --git a/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil.py b/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil.py
index 2f5c4ab..78eb398 100755
--- a/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil.py
+++ b/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil.py
@@ -32,8 +32,8 @@ resolutions = {
# 2: 2e-5, # h=1.1180
# 4: 1e-6, # h=0.5590
# 8: 1e-6, # h=0.2814
- 16: 1e-6, # h=0.1412
- # 32: 1e-6,
+ # 16: 1e-6, # h=0.1412
+ 32: 1e-6,
# 64: 5e-7,
# 128: 5e-7
}
@@ -41,7 +41,7 @@ resolutions = {
# GRID #######################
# mesh_resolution = 20
timestep_size = 0.005
-number_of_timesteps = 500
+number_of_timesteps = 200
plot_timestep_every = 2
# decide how many timesteps you want analysed. Analysed means, that we write
# out subsequent errors of the L-iteration within the timestep.
@@ -65,7 +65,7 @@ lambda34_w = 40
lambda34_nw = 40
include_gravity = False
-debugflag = True
+debugflag = False
analyse_condition = False
if mesh_study:
@@ -340,7 +340,7 @@ lambda_param = {
intrinsic_permeability = {
1: 1, # sand
2: 0.1, # sand, there is a range
- 3: 10e-3, # clay has a range
+ 3: 10e-2, # clay has a range
4: 10e-4
}
@@ -592,10 +592,11 @@ x, y = sym.symbols('x[0], x[1]') # needed by UFL
t = sym.symbols('t', positive=True)
p_e_sym_2patch = {
- 1: {'wetting': -7 - (1+t*t)*(1 + x*x + y*y),
+ 1: {'wetting': -6 - (1+t*t)*(1 + x*x + y*y),
'nonwetting': 0.0*t}, # -1-t*(1.1 + y + x**2)**2},
- 2: {'wetting': -7.0 - (1.0 + t*t)*(1.0 + x*x),
- 'nonwetting': (-1-t*(1.1 + x**2)**2 - sym.sqrt(5+t**2))*y**2},
+ 2: {'wetting': -6.0 - (1.0 + t*t)*(1.0 + x*x),
+ 'nonwetting': (-1-t*(1.1+y + x**2))*y**2},
+ # 'nonwetting': (-1-t*(1.1 + x**2)**2 - sym.sqrt(5+t**2))*y**2},
}
p_e_sym = {
--
GitLab