From 220639c8d7a0f00dfbf343d84d6dbe20b14faa0b Mon Sep 17 00:00:00 2001 From: Benjamin Rodenberg Date: Wed, 4 Oct 2023 19:04:29 +0200 Subject: [PATCH] Update beta from 1.3 to 1.2. --- partitioned-heat-conduction-complex/fenics/heat.py | 2 +- partitioned-heat-conduction-direct/nutils/heat.py | 2 +- partitioned-heat-conduction/fenicsx/heat.py | 2 +- partitioned-heat-conduction/nutils/heat.py | 2 +- partitioned-heat-conduction/openfoam-dirichlet/0.orig/T | 2 +- .../openfoam-dirichlet/setInitialField.sh | 2 +- partitioned-heat-conduction/openfoam-neumann/0.orig/T | 2 +- partitioned-heat-conduction/openfoam-neumann/setInitialField.sh | 2 +- partitioned-heat-conduction/openfoam-solver/heatTransfer.C | 2 +- 9 files changed, 9 insertions(+), 9 deletions(-) diff --git a/partitioned-heat-conduction-complex/fenics/heat.py b/partitioned-heat-conduction-complex/fenics/heat.py index cd44241cc..766755ec5 100644 --- a/partitioned-heat-conduction-complex/fenics/heat.py +++ b/partitioned-heat-conduction-complex/fenics/heat.py @@ -72,7 +72,7 @@ def determine_gradient(V_g, u, flux): # Error is bounded by coupling accuracy. In theory we can obtain the analytical solution. error_tol = 10 ** -6 alpha = 3 # parameter alpha -beta = 1.3 # parameter beta +beta = 1.2 # parameter beta gamma = args.gamma # parameter gamma, dependence of heat flux on time # Create mesh and separate mesh components for grid, boundary and coupling interface diff --git a/partitioned-heat-conduction-direct/nutils/heat.py b/partitioned-heat-conduction-direct/nutils/heat.py index 2d74fae33..92bf3ae33 100644 --- a/partitioned-heat-conduction-direct/nutils/heat.py +++ b/partitioned-heat-conduction-direct/nutils/heat.py @@ -7,7 +7,7 @@ import precice -def main(side='Dirichlet', n=10, degree=1, timestep=.1, alpha=3., beta=1.3): +def main(side='Dirichlet', n=10, degree=1, timestep=.1, alpha=3., beta=1.2): if side == 'Dirichlet': x_grid = np.linspace(0, 1, n) diff --git a/partitioned-heat-conduction/fenicsx/heat.py b/partitioned-heat-conduction/fenicsx/heat.py index 4f9d7ba7f..7ecef9fc4 100644 --- a/partitioned-heat-conduction/fenicsx/heat.py +++ b/partitioned-heat-conduction/fenicsx/heat.py @@ -68,7 +68,7 @@ def determine_gradient(V_g, u): error_tol = args.error_tol alpha = 3 # parameter alpha -beta = 1.3 # parameter beta +beta = 1.2 # parameter beta if args.dirichlet and not args.neumann: problem = ProblemType.DIRICHLET diff --git a/partitioned-heat-conduction/nutils/heat.py b/partitioned-heat-conduction/nutils/heat.py index c0ce77f3d..edea5b5be 100644 --- a/partitioned-heat-conduction/nutils/heat.py +++ b/partitioned-heat-conduction/nutils/heat.py @@ -7,7 +7,7 @@ import precice -def main(side='Dirichlet', n=10, degree=1, timestep=.1, alpha=3., beta=1.3): +def main(side='Dirichlet', n=10, degree=1, timestep=.1, alpha=3., beta=1.2): if side == 'Dirichlet': x_grid = np.linspace(0, 1, n) diff --git a/partitioned-heat-conduction/openfoam-dirichlet/0.orig/T b/partitioned-heat-conduction/openfoam-dirichlet/0.orig/T index fed3ddbc4..f90cb60d5 100644 --- a/partitioned-heat-conduction/openfoam-dirichlet/0.orig/T +++ b/partitioned-heat-conduction/openfoam-dirichlet/0.orig/T @@ -22,7 +22,7 @@ boundaryField DirichletBoundary { type groovyBC; - variables "val=1+pow(pos().x,2)+(3*pow(pos().y,2))+1.3*time();"; + variables "val=1+pow(pos().x,2)+(3*pow(pos().y,2))+1.2*time();"; valueExpression "val"; value uniform 0; evaluateDuringConstruction 1; diff --git a/partitioned-heat-conduction/openfoam-dirichlet/setInitialField.sh b/partitioned-heat-conduction/openfoam-dirichlet/setInitialField.sh index bf3dfd29c..02d282160 100755 --- a/partitioned-heat-conduction/openfoam-dirichlet/setInitialField.sh +++ b/partitioned-heat-conduction/openfoam-dirichlet/setInitialField.sh @@ -5,4 +5,4 @@ set -e -u rm -rf ./0 cp -r ./0.orig 0 # Initialize the new field -funkySetFields -keepPatches -field T -expression '1+pow(pos().x,2)+(3*pow(pos().y,2))+1.3*time()' -time '0' +funkySetFields -keepPatches -field T -expression '1+pow(pos().x,2)+(3*pow(pos().y,2))+1.2*time()' -time '0' diff --git a/partitioned-heat-conduction/openfoam-neumann/0.orig/T b/partitioned-heat-conduction/openfoam-neumann/0.orig/T index 01f8fc1af..753212b67 100644 --- a/partitioned-heat-conduction/openfoam-neumann/0.orig/T +++ b/partitioned-heat-conduction/openfoam-neumann/0.orig/T @@ -22,7 +22,7 @@ boundaryField DirichletBoundary { type groovyBC; - variables "val=1+pow(pos().x,2)+(3*pow(pos().y,2))+1.3*time();"; + variables "val=1+pow(pos().x,2)+(3*pow(pos().y,2))+1.2*time();"; valueExpression "val"; value uniform 0; evaluateDuringConstruction 1; diff --git a/partitioned-heat-conduction/openfoam-neumann/setInitialField.sh b/partitioned-heat-conduction/openfoam-neumann/setInitialField.sh index bf3dfd29c..02d282160 100755 --- a/partitioned-heat-conduction/openfoam-neumann/setInitialField.sh +++ b/partitioned-heat-conduction/openfoam-neumann/setInitialField.sh @@ -5,4 +5,4 @@ set -e -u rm -rf ./0 cp -r ./0.orig 0 # Initialize the new field -funkySetFields -keepPatches -field T -expression '1+pow(pos().x,2)+(3*pow(pos().y,2))+1.3*time()' -time '0' +funkySetFields -keepPatches -field T -expression '1+pow(pos().x,2)+(3*pow(pos().y,2))+1.2*time()' -time '0' diff --git a/partitioned-heat-conduction/openfoam-solver/heatTransfer.C b/partitioned-heat-conduction/openfoam-solver/heatTransfer.C index 6194f6fc1..65269836f 100644 --- a/partitioned-heat-conduction/openfoam-solver/heatTransfer.C +++ b/partitioned-heat-conduction/openfoam-solver/heatTransfer.C @@ -62,7 +62,7 @@ int main(int argc, char *argv[]) Info<< "\nCalculating temperature distribution\n" << endl; const double alpha = 3; - const double beta = 1.3; + const double beta = 1.2; const double rhs = beta - 2 - 2 * alpha; volScalarField f