From b16e058221bb024f0395af6f1838b75de0699f92 Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Wed, 21 Dec 2011 10:16:33 +0000
Subject: [PATCH] MRFZoneTemplates: Call internalField on phi to ensure the
old-time value is stored for ddtPhiCorr
---
.../cfdTools/general/MRF/MRFZoneTemplates.C | 12 ++++++++++--
1 file changed, 10 insertions(+), 2 deletions(-)
diff --git a/src/finiteVolume/cfdTools/general/MRF/MRFZoneTemplates.C b/src/finiteVolume/cfdTools/general/MRF/MRFZoneTemplates.C
index 55821dbfde6..23cae40603c 100644
--- a/src/finiteVolume/cfdTools/general/MRF/MRFZoneTemplates.C
+++ b/src/finiteVolume/cfdTools/general/MRF/MRFZoneTemplates.C
@@ -44,11 +44,15 @@ void Foam::MRFZone::relativeRhoFlux
const vector& origin = origin_.value();
const vector& Omega = Omega_.value();
+ const vectorField& Cfi = Cf.internalField();
+ const vectorField& Sfi = Sf.internalField();
+ scalarField& phii = phi.internalField();
+
// Internal faces
forAll(internalFaces_, i)
{
label facei = internalFaces_[i];
- phi[facei] -= rho[facei]*(Omega ^ (Cf[facei] - origin)) & Sf[facei];
+ phii[facei] -= rho[facei]*(Omega ^ (Cfi[facei] - origin)) & Sfi[facei];
}
// Included patches
@@ -91,11 +95,15 @@ void Foam::MRFZone::absoluteRhoFlux
const vector& origin = origin_.value();
const vector& Omega = Omega_.value();
+ const vectorField& Cfi = Cf.internalField();
+ const vectorField& Sfi = Sf.internalField();
+ scalarField& phii = phi.internalField();
+
// Internal faces
forAll(internalFaces_, i)
{
label facei = internalFaces_[i];
- phi[facei] += rho[facei]*(Omega ^ (Cf[facei] - origin)) & Sf[facei];
+ phii[facei] += rho[facei]*(Omega ^ (Cfi[facei] - origin)) & Sfi[facei];
}
// Included patches
--
GitLab