diff --git a/applications/solvers/coupled/conjugateHeatFoam/conjugateHeatFoam.C b/applications/solvers/coupled/conjugateHeatFoam/conjugateHeatFoam.C
index 7366b5618..46dc0d1f9 100644
--- a/applications/solvers/coupled/conjugateHeatFoam/conjugateHeatFoam.C
+++ b/applications/solvers/coupled/conjugateHeatFoam/conjugateHeatFoam.C
@@ -109,6 +109,7 @@ int main(int argc, char *argv[])
 
         // Update density according to Boussinesq approximation
         rhok = 1.0 - beta*(T - TRef);
+        rhok.correctBoundaryConditions();
 
         runTime.write();
 
diff --git a/applications/solvers/coupled/conjugateHeatSimpleFoam/conjugateHeatSimpleFoam.C b/applications/solvers/coupled/conjugateHeatSimpleFoam/conjugateHeatSimpleFoam.C
index 5e7478deb..9a1e44be9 100644
--- a/applications/solvers/coupled/conjugateHeatSimpleFoam/conjugateHeatSimpleFoam.C
+++ b/applications/solvers/coupled/conjugateHeatSimpleFoam/conjugateHeatSimpleFoam.C
@@ -97,6 +97,7 @@ int main(int argc, char *argv[])
 
         // Update density according to Boussinesq approximation
         rhok = 1.0 - beta*(T - TRef);
+        rhok.correctBoundaryConditions();
 
         runTime.write();