Skip to content
Merged
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
7 changes: 3 additions & 4 deletions CHT/flow-over-plate/buoyantPimpleFoam-fenics/Solid/heat.py
Original file line number Diff line number Diff line change
Expand Up @@ -102,9 +102,8 @@ def determine_heat_flux(V_g, u, k, flux):
# Define boundary condition
u_D = Constant('310')
u_D_function = interpolate(u_D, V)
# Define flux in y direction on coupling interface (grad(u_D) in normal direction)
f_N = Constant('0')
f_N_function = interpolate(f_N, V)
# We will only exchange flux in y direction on coupling interface. No initialization necessary.
V_flux_y = V_g.sub(1)

coupling_boundary = TopBoundary()
bottom_boundary = BottomBoundary()
Expand All @@ -116,7 +115,7 @@ def determine_heat_flux(V_g, u, k, flux):
# Adapter definition and initialization
precice = Adapter(adapter_config_filename="precice-adapter-config.json")

precice_dt = precice.initialize(coupling_boundary, read_function_space=V, write_object=f_N_function)
precice_dt = precice.initialize(coupling_boundary, read_function_space=V, write_object=V_flux_y)

# Create a FEniCS Expression to define and control the coupling boundary values
coupling_expression = precice.create_coupling_expression()
Expand Down