preCICE converges but the error is still too high

precice · Jan 2, 2025 · 793b0a3 · 793b0a3
1 parent 728db16
commit 793b0a3
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/tutorials/partitioned-heat-conduction/fenicsx/heat.py b/tutorials/partitioned-heat-conduction/fenicsx/heat.py
@@ -147,7 +147,7 @@ def __call__(self, x):
 if problem is ProblemType.NEUMANN:
     # modify Neumann boundary condition on coupling interface, modify weak
     # form correspondingly
-    F += coupling_expression * v  * ufl.ds
+    F += dt*coupling_expression * v * ufl.ds
 
 a = fem.form(ufl.lhs(F))
 L = fem.form(ufl.rhs(F))
@@ -216,8 +216,8 @@ def __call__(self, x):
         flux = determine_gradient(V_g, uh)
         flux_x = fem.Function(W)
         flux_x.interpolate(flux.sub(0))
-        #precice.write_data(flux_x)
-        precice.write_data(f_N)
+        precice.write_data(flux_x)
+        #precice.write_data(f_N)
     elif problem is ProblemType.NEUMANN:
         # Neumann problem reads flux and writes temperature on boundary to Dirichlet problem
         precice.write_data(uh)