From 7c63c4d759bebf92ba572f97b811ef93a2d94f40 Mon Sep 17 00:00:00 2001
From: David Seus <david.seus@ians.uni-stuttgart.de>
Date: Tue, 3 Dec 2019 14:53:39 +0100
Subject: [PATCH] setup layered soil case mesh study
---
.../TP-TP-layered_soil_pure_dd_mesh_study.py | 38 +++++++++----------
1 file changed, 19 insertions(+), 19 deletions(-)
diff --git a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/mesh_study/TP-TP-layered_soil_pure_dd_mesh_study.py b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/mesh_study/TP-TP-layered_soil_pure_dd_mesh_study.py
index c8d4c35..9b07c00 100755
--- a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/mesh_study/TP-TP-layered_soil_pure_dd_mesh_study.py
+++ b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/mesh_study/TP-TP-layered_soil_pure_dd_mesh_study.py
@@ -32,35 +32,35 @@ max_iter_num = 1000
FEM_Lagrange_degree = 1
mesh_study = True
resolutions = {
- # 1: 1e-7, # h=2
- 2: 1e-7, # h=1.1180
- # 4: 1e-7, # h=0.5590
- # 8: 1e-7, # h=0.2814
- # 16: 1e-7, # h=0.1412
- # 32: 1e-7, # h=0.0706
- # 64: 1e-7,
- # 128: 1e-7,
+ 1: 1e-6, # h=2
+ 2: 1e-6, # h=1.1180
+ 4: 1e-6, # h=0.5590
+ 8: 1e-6, # h=0.2814
+ 16: 1e-6, # h=0.1412
+ 32: 1e-6, # h=0.0706
+ 64: 1e-6,
+ 128: 1e-6,
# 256: 1e-7
}
############ GRID #######################
# mesh_resolution = 20
-timestep_size = 0.001
-number_of_timesteps = 300
+timestep_size = 0.005
+number_of_timesteps = 100
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.
number_of_timesteps_to_analyse = 0
-starttimes = [0.0, 0.25, 0.5]
+starttimes = [0.0]
Lw = 0.025 #/timestep_size
Lnw=Lw
-lambda_w = 80
-lambda_nw = 80
+lambda_w = 8
+lambda_nw = 8
-include_gravity = True
-debugflag = True
+include_gravity = False
+debugflag = False
analyse_condition = True
if mesh_study:
@@ -76,10 +76,10 @@ if mesh_study:
'space_errornorms': True,
'meshes_and_markers': True,
'L_iterations_per_timestep': False,
- 'solutions': False,
- 'absolute_differences': False,
+ 'solutions': True,
+ 'absolute_differences': True,
'condition_numbers': analyse_condition,
- 'subsequent_errors': False
+ 'subsequent_errors': True
}
else:
write_to_file = {
@@ -262,7 +262,7 @@ densities = {
'nonwetting':1}, #1.225}}
}
-gravity_acceleration = 9.81
+gravity_acceleration = 1
# porosities taken from
# https://www.geotechdata.info/parameter/soil-porosity.html
# Dict of the form: { subdom_num : porosity }
--
GitLab