From 352d62aa48cae279e11de72c12d765d6bf5553b6 Mon Sep 17 00:00:00 2001
From: David Seus <david.seus@ians.uni-stuttgart.de>
Date: Wed, 19 Aug 2020 14:18:41 +0200
Subject: [PATCH] set up TPTP example
---
...layered_soil_with_inner_patch-realistic.py | 8 +-
.../TP-TP-layered_soil_with_inner_patch.py | 82 +++++++++----------
2 files changed, 45 insertions(+), 45 deletions(-)
diff --git a/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/TP-R-layered_soil_with_inner_patch-realistic.py b/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/TP-R-layered_soil_with_inner_patch-realistic.py
index 94efbe4..dcca4a3 100755
--- a/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/TP-R-layered_soil_with_inner_patch-realistic.py
+++ b/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/TP-R-layered_soil_with_inner_patch-realistic.py
@@ -47,8 +47,8 @@ resolutions = {
# 4: 2e-6, # h=0.5590
# 8: 2e-6, # h=0.2814
# 16: 1e-6, # h=0.1412
- 32: 1e-6,
- # 64: 2e-6,
+ # 32: 1e-6,
+ 64: 1e-7,
# 128: 2e-6
}
@@ -58,7 +58,7 @@ resolutions = {
starttimes = {0: 0.0}
# starttimes = {0: 0.0, 1:0.3, 2:0.6, 3:0.9}
timestep_size = 0.001
-number_of_timesteps = 5000
+number_of_timesteps = 2000
# LDD scheme parameters ######################################################
Lw1 = 0.01 # /timestep_size
@@ -117,7 +117,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 = 6
+number_of_timesteps_to_analyse = 5
# 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.
diff --git a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch.py b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch.py
index 6e81215..27663b2 100755
--- a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch.py
+++ b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch.py
@@ -45,9 +45,9 @@ resolutions = {
# 2: 5e-5,
# 4: 5e-5,
# 8: 5e-5,
- 16: 3e-5,
- # 32: 8e-6,
- # 64: 2e-6,
+ # 16: 3e-5,
+ 32: 1e-6,
+ 64: 1e-7,
# 128: 1e-6,
# 256: 1e-6,
}
@@ -55,55 +55,55 @@ 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.3}
+starttimes = {0: 0.0}
timestep_size = 0.001
-number_of_timesteps = 5
+number_of_timesteps = 2000
# LDD scheme parameters ######################################################
-Lw1 = 0.25 # /timestep_size
-Lnw1 = 0.25
+Lw1 = 0.01 # /timestep_size
+Lnw1 = 0.001
-Lw2 = 0.25 # /timestep_size
-Lnw2 = 0.25
+Lw2 = 0.01 # /timestep_size
+Lnw2 = 0.001
-Lw3 = 0.25 # /timestep_size
-Lnw3 = 0.25
+Lw3 = 0.01 # /timestep_size
+Lnw3 = 0.001
-Lw4 = 0.25 # /timestep_size
-Lnw4 = 0.25
+Lw4 = 0.01 # /timestep_size
+Lnw4 = 0.001
-Lw5 = 0.25 # /timestep_size
-Lnw5 = 0.25
+Lw5 = 0.01 # /timestep_size
+Lnw5 = 0.001
-Lw6 = 0.25 # /timestep_size
-Lnw6 = 0.25
+Lw6 = 0.01 # /timestep_size
+Lnw6 = 0.001
-lambda12_w = 40
-lambda12_nw = 40
+lambda12_w = 0.75
+lambda12_nw = 0.75
-lambda23_w = 40
-lambda23_nw = 40
+lambda23_w = 0.75
+lambda23_nw = 0.75
-lambda24_w = 40
-lambda24_nw= 40
+lambda24_w = 0.75
+lambda24_nw= 0.75
-lambda25_w= 40
-lambda25_nw= 40
+lambda25_w= 0.75
+lambda25_nw= 0.75
-lambda34_w = 40
-lambda34_nw = 40
+lambda34_w = 0.75
+lambda34_nw = 0.75
-lambda36_w = 40
-lambda36_nw = 40
+lambda36_w = 0.75
+lambda36_nw = 0.75
-lambda45_w = 40
-lambda45_nw = 40
+lambda45_w = 0.75
+lambda45_nw = 0.75
-lambda46_w = 40
-lambda46_nw = 40
+lambda46_w = 0.75
+lambda46_nw = 0.75
-lambda56_w = 40
-lambda56_nw = 40
+lambda56_w = 0.75
+lambda56_nw = 0.75
include_gravity = False
debugflag = True
@@ -113,7 +113,7 @@ analyse_condition = False
# when number_of_timesteps is high, it might take a long time to write all
# timesteps to disk. Therefore, you can choose to only write data of every
# plot_timestep_every timestep to disk.
-plot_timestep_every = 1
+plot_timestep_every = 3
# 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
@@ -337,13 +337,13 @@ p_e_sym = {
2: {'wetting': -6.0 - (1.0 + t*t)*(1.0 + x*x + y*y),
'nonwetting': (-1 -t*(1.1 + y + x**2)) },
3: {'wetting': (-6.0 - (1.0 + t*t)*(1.0 + x*x)),
- 'nonwetting': (-1 -t*(1.0 + x**2) - sym.sin(2+t**2)*y**2) },
+ 'nonwetting': (-1 -t*(1.0 + x**2) - sym.sin(2+t**2)**2*y**2) },
4: {'wetting': (-6.0 - (1.0 + t*t)*(1.0 + x*x)),
- 'nonwetting': (-1 -t*(1.0 + x**2) - sym.sin(2+t**2)*y**2) },
+ 'nonwetting': (-1 -t*(1.0 + x**2) - sym.sin(2+t**2)**2*y**2) },
5: {'wetting': (-6.0 - (1.0 + t*t)*(1.0 + x*x)),
- 'nonwetting': (-1 -t*(1.0 + x**2) - sym.sin(2+t**2)*y**2) },
+ 'nonwetting': (-1 -t*(1.0 + x**2) - sym.sin(2+t**2)**2*y**2) },
6: {'wetting': (-6.0 - (1.0 + t*t)*(1.0 + x*x)),
- 'nonwetting': (-1 -t*(1.0 + x**2) - sym.sin(2+t**2)*y**2) },
+ 'nonwetting': (-1 -t*(1.0 + x**2) - sym.sin(2+t**2)**2*y**2) },
}
pc_e_sym = hlp.generate_exact_symbolic_pc(
@@ -463,7 +463,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