Skip to content
Snippets Groups Projects
Commit 5128ae22 authored by andy's avatar andy
Browse files

ENH: Added MRF functionality to porousSimpleFoam

parent d75fa36c
No related branches found
No related tags found
No related merge requests found
...@@ -8,6 +8,8 @@ ...@@ -8,6 +8,8 @@
sources(U) sources(U)
); );
mrfZones.addCoriolis(UEqn());
UEqn().relax(); UEqn().relax();
sources.constrain(UEqn()); sources.constrain(UEqn());
......
IOMRFZoneList mrfZones(mesh);
mrfZones.correctBoundaryVelocity(U);
IOporosityModelList pZones(mesh); IOporosityModelList pZones(mesh);
Switch pressureImplicitPorosity(false); Switch pressureImplicitPorosity(false);
......
...@@ -12,6 +12,8 @@ UEqn.clear(); ...@@ -12,6 +12,8 @@ UEqn.clear();
surfaceScalarField phiHbyA("phiHbyA", fvc::interpolate(HbyA) & mesh.Sf()); surfaceScalarField phiHbyA("phiHbyA", fvc::interpolate(HbyA) & mesh.Sf());
adjustPhi(phiHbyA, U, p); adjustPhi(phiHbyA, U, p);
mrfZones.relativeFlux(phiHbyA);
while (simple.correctNonOrthogonal()) while (simple.correctNonOrthogonal())
{ {
tmp<fvScalarMatrix> tpEqn; tmp<fvScalarMatrix> tpEqn;
...@@ -50,3 +52,4 @@ else ...@@ -50,3 +52,4 @@ else
} }
U.correctBoundaryConditions(); U.correctBoundaryConditions();
sources.correct(U);
...@@ -2,7 +2,7 @@ ...@@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | \\ / O peration |
\\ / A nd | Copyright (C) 2011 OpenFOAM Foundation \\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
...@@ -26,15 +26,17 @@ Application ...@@ -26,15 +26,17 @@ Application
Description Description
Steady-state solver for incompressible, turbulent flow with Steady-state solver for incompressible, turbulent flow with
implicit or explicit porosity treatment implicit or explicit porosity treatment and support for multiple reference
frames (MRF)
\*---------------------------------------------------------------------------*/ \*---------------------------------------------------------------------------*/
#include "fvCFD.H" #include "fvCFD.H"
#include "singlePhaseTransportModel.H" #include "singlePhaseTransportModel.H"
#include "RASModel.H" #include "RASModel.H"
#include "IOporosityModelList.H"
#include "simpleControl.H" #include "simpleControl.H"
#include "IOMRFZoneList.H"
#include "IOporosityModelList.H"
#include "IObasicSourceList.H" #include "IObasicSourceList.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
...@@ -48,7 +50,7 @@ int main(int argc, char *argv[]) ...@@ -48,7 +50,7 @@ int main(int argc, char *argv[])
simpleControl simple(mesh); simpleControl simple(mesh);
#include "createFields.H" #include "createFields.H"
#include "createPorousZones.H" #include "createZones.H"
#include "initContinuityErrs.H" #include "initContinuityErrs.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment