ENH: Updated sources handling in pimpleFoam

applications/solvers/incompressible/pimpleFoam/UEqn.H

@@ -5,12 +5,10 @@ tmp<fvVectorMatrix> UEqn
     fvm::ddt(U)
   + fvm::div(phi, U)
   + turbulence->divDevReff(U)
+ ==
+    sources(U)
 );
 
-mrfZones.addCoriolis(UEqn());
-
-pZones.addResistance(UEqn());
-
 UEqn().relax();
 
 sources.constrain(UEqn());
@@ -19,5 +17,5 @@ volScalarField rAU(1.0/UEqn().A());
 
 if (pimple.momentumPredictor())
 {
-    solve(UEqn() == -fvc::grad(p) + sources(U));
+    solve(UEqn() == -fvc::grad(p));
 }

applications/solvers/incompressible/pimpleFoam/createZones.H

-    IOMRFZoneList mrfZones(mesh);
-    mrfZones.correctBoundaryVelocity(U);
-
-    IOporosityModelList pZones(mesh);

applications/solvers/incompressible/pimpleFoam/pEqn.H

 volVectorField HbyA("HbyA", U);
-HbyA = rAU*(UEqn() == sources(U))().H();
+HbyA = rAU*UEqn().H();
 
 if (pimple.nCorrPISO() <= 1)
 {
@@ -15,7 +15,7 @@ surfaceScalarField phiHbyA
 
 adjustPhi(phiHbyA, U, p);
 
-mrfZones.relativeFlux(phiHbyA);
+sources.relativeFlux(phiHbyA);
 
 // Non-orthogonal pressure corrector loop
 while (pimple.correctNonOrthogonal())

applications/solvers/incompressible/pimpleFoam/pimpleFoam.C

@@ -30,8 +30,7 @@ Description
 
     Sub-models include:
     - turbulence modelling, i.e. laminar, RAS or LES
-    - porosity (explicit treatment)
-    - Multiple Reference Frame (MRF)
+    - run-time selectable sources, e.g. MRF, explicit porosity
 
 \*---------------------------------------------------------------------------*/
 
@@ -51,7 +50,6 @@ int main(int argc, char *argv[])
     #include "createTime.H"
     #include "createMesh.H"
     #include "createFields.H"
-    #include "createZones.H"
     #include "initContinuityErrs.H"
 
     pimpleControl pimple(mesh);