From 08164161b346e5481336261b3e4dbda2f63a5ed9 Mon Sep 17 00:00:00 2001
From: David <forenkram@gmx.de>
Date: Tue, 16 Jun 2020 11:19:21 +0200
Subject: [PATCH] change parameters to make it converge
---
.../TP-TP-2-patch-same-intrinsic-perm.py | 27 +++++++------------
1 file changed, 10 insertions(+), 17 deletions(-)
diff --git a/Two-phase-Two-phase/two-patch/TP-TP-2-patch-test-case/TP-TP-2-patch-same-intrinsic-perm.py b/Two-phase-Two-phase/two-patch/TP-TP-2-patch-test-case/TP-TP-2-patch-same-intrinsic-perm.py
index 84614ce..4178c0f 100755
--- a/Two-phase-Two-phase/two-patch/TP-TP-2-patch-test-case/TP-TP-2-patch-same-intrinsic-perm.py
+++ b/Two-phase-Two-phase/two-patch/TP-TP-2-patch-test-case/TP-TP-2-patch-same-intrinsic-perm.py
@@ -36,7 +36,7 @@ thisfile = "TP-TP-2-patch-same-intrinsic-perm.py"
# GENERAL SOLVER CONFIG ######################################################
# maximal iteration per timestep
-max_iter_num = 300
+max_iter_num = 400
FEM_Lagrange_degree = 1
# GRID AND MESH STUDY SPECIFICATIONS #########################################
@@ -47,7 +47,7 @@ resolutions = {
# 4: 1e-6,
# 8: 1e-6,
# 16: 5e-6,
- 32: 3e-6,
+ 32: 1e-5,
# 64: 2e-6,
# 128: 1e-6,
# 256: 1e-6,
@@ -61,16 +61,16 @@ timestep_size = 0.001
number_of_timesteps = 1000
# LDD scheme parameters ######################################################
-Lw1 = 0.25 #/timestep_size
+Lw1 = 0.5 #/timestep_size
Lnw1= 0.25
-Lw2 = 0.25 #/timestep_size
+Lw2 = 0.5 #/timestep_size
Lnw2= 0.25
lambda_w = 4
lambda_nw = 4
-include_gravity = False
+include_gravity = True
debugflag = False
analyse_condition = False
@@ -120,16 +120,9 @@ else:
}
# OUTPUT FILE STRING #########################################################
-if mesh_study:
- output_string = "./output/{}-{}_timesteps{}_P{}".format(
- datestr, use_case, number_of_timesteps, FEM_Lagrange_degree
- )
-else:
- for tol in resolutions.values():
- solver_tol = tol
- output_string = "./output/{}-{}_timesteps{}_P{}_solver_tol{}".format(
- datestr, use_case, number_of_timesteps, FEM_Lagrange_degree, solver_tol
- )
+output_string = "./output/{}-{}_timesteps{}_P{}".format(
+ datestr, use_case, number_of_timesteps, FEM_Lagrange_degree
+ )
# DOMAIN AND INTERFACE #######################################################
@@ -243,8 +236,8 @@ lambda_param = {#
}
intrinsic_permeability = {
- 1: 0.1,
- 2: 0.1,
+ 1: 0.01,
+ 2: 0.01,
}
--
GitLab