Use beta consistent with literature

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

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)

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

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)

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;

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'

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;

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'

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