From 286971f65e89d1b1057026f1aea8c83fa95ba114 Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Mon, 29 Oct 2012 17:11:08 +0000
Subject: [PATCH] ENH: Added MRF functionality to buoyantSimpleFoam
---
applications/solvers/heatTransfer/buoyantSimpleFoam/UEqn.H | 2 ++
.../solvers/heatTransfer/buoyantSimpleFoam/buoyantSimpleFoam.C | 2 ++
.../solvers/heatTransfer/buoyantSimpleFoam/createZones.H | 3 +++
applications/solvers/heatTransfer/buoyantSimpleFoam/pEqn.H | 2 ++
4 files changed, 9 insertions(+)
create mode 100644 applications/solvers/heatTransfer/buoyantSimpleFoam/createZones.H
diff --git a/applications/solvers/heatTransfer/buoyantSimpleFoam/UEqn.H b/applications/solvers/heatTransfer/buoyantSimpleFoam/UEqn.H
index ca28910aaf3..e1593c3691c 100644
--- a/applications/solvers/heatTransfer/buoyantSimpleFoam/UEqn.H
+++ b/applications/solvers/heatTransfer/buoyantSimpleFoam/UEqn.H
@@ -8,6 +8,8 @@
UEqn().relax();
+ mrfZones.addCoriolis(rho, UEqn());
+
if (simple.momentumPredictor())
{
solve
diff --git a/applications/solvers/heatTransfer/buoyantSimpleFoam/buoyantSimpleFoam.C b/applications/solvers/heatTransfer/buoyantSimpleFoam/buoyantSimpleFoam.C
index 4a9387c1882..fc47fc90a1c 100644
--- a/applications/solvers/heatTransfer/buoyantSimpleFoam/buoyantSimpleFoam.C
+++ b/applications/solvers/heatTransfer/buoyantSimpleFoam/buoyantSimpleFoam.C
@@ -34,6 +34,7 @@ Description
#include "RASModel.H"
#include "fixedGradientFvPatchFields.H"
#include "simpleControl.H"
+#include "IOMRFZoneList.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@@ -44,6 +45,7 @@ int main(int argc, char *argv[])
#include "createMesh.H"
#include "readGravitationalAcceleration.H"
#include "createFields.H"
+ #include "createZones.H"
#include "initContinuityErrs.H"
simpleControl simple(mesh);
diff --git a/applications/solvers/heatTransfer/buoyantSimpleFoam/createZones.H b/applications/solvers/heatTransfer/buoyantSimpleFoam/createZones.H
new file mode 100644
index 00000000000..08480c68d2d
--- /dev/null
+++ b/applications/solvers/heatTransfer/buoyantSimpleFoam/createZones.H
@@ -0,0 +1,3 @@
+ IOMRFZoneList mrfZones(mesh);
+ mrfZones.correctBoundaryVelocity(U);
+
diff --git a/applications/solvers/heatTransfer/buoyantSimpleFoam/pEqn.H b/applications/solvers/heatTransfer/buoyantSimpleFoam/pEqn.H
index f63e12f3637..66ae75b650e 100644
--- a/applications/solvers/heatTransfer/buoyantSimpleFoam/pEqn.H
+++ b/applications/solvers/heatTransfer/buoyantSimpleFoam/pEqn.H
@@ -17,6 +17,8 @@
fvc::interpolate(rho)*(fvc::interpolate(HbyA) & mesh.Sf())
);
+ mrfZones.relativeFlux(fvc::interpolate(rho), phiHbyA);
+
bool closedVolume = adjustPhi(phiHbyA, U, p_rgh);
phiHbyA += phig;
--
GitLab