From 6459a1d179cacd1af249f95dee91b2ef8ecd0f8a Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Fri, 11 Apr 2014 12:11:36 +0100
Subject: [PATCH] BUG: heat transfer Boussinesq solvers - added fvOptions MRF
 support - mantis #1267

---
 .../solvers/heatTransfer/buoyantBoussinesqPimpleFoam/pEqn.H   | 4 +++-
 .../solvers/heatTransfer/buoyantBoussinesqSimpleFoam/pEqn.H   | 4 +++-
 2 files changed, 6 insertions(+), 2 deletions(-)

diff --git a/applications/solvers/heatTransfer/buoyantBoussinesqPimpleFoam/pEqn.H b/applications/solvers/heatTransfer/buoyantBoussinesqPimpleFoam/pEqn.H
index d77eb12b..b73e5c88 100644
--- a/applications/solvers/heatTransfer/buoyantBoussinesqPimpleFoam/pEqn.H
+++ b/applications/solvers/heatTransfer/buoyantBoussinesqPimpleFoam/pEqn.H
@@ -15,13 +15,15 @@
       + phig
     );
 
+    fvOptions.makeRelative(phiHbyA);
+
     // Update the fixedFluxPressure BCs to ensure flux consistency
     setSnGrad<fixedFluxPressureFvPatchScalarField>
     (
         p_rgh.boundaryField(),
         (
             phiHbyA.boundaryField()
-          - (mesh.Sf().boundaryField() & U.boundaryField())
+          - fvOptions.relative(mesh.Sf().boundaryField() & U.boundaryField())
         )/(mesh.magSf().boundaryField()*rAUf.boundaryField())
     );
 
diff --git a/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/pEqn.H b/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/pEqn.H
index bae1bcef..ae104e8f 100644
--- a/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/pEqn.H
+++ b/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/pEqn.H
@@ -14,6 +14,8 @@
         (fvc::interpolate(HbyA) & mesh.Sf())
     );
 
+    fvOptions.makeRelative(phiHbyA);
+
     adjustPhi(phiHbyA, U, p_rgh);
 
     phiHbyA += phig;
@@ -24,7 +26,7 @@
         p_rgh.boundaryField(),
         (
             phiHbyA.boundaryField()
-          - (mesh.Sf().boundaryField() & U.boundaryField())
+          - fvOptions.relative(mesh.Sf().boundaryField() & U.boundaryField())
         )/(mesh.magSf().boundaryField()*rAUf.boundaryField())
     );