diff --git a/applications/solvers/compressible/rhoPimpleFoam/rhoPorousMRFLTSPimpleFoam/rhoPorousMRFLTSPimpleFoam.C b/applications/solvers/compressible/rhoPimpleFoam/rhoPorousMRFLTSPimpleFoam/rhoPorousMRFLTSPimpleFoam.C index 82ba5cbac1..18a32a1791 100644 --- a/applications/solvers/compressible/rhoPimpleFoam/rhoPorousMRFLTSPimpleFoam/rhoPorousMRFLTSPimpleFoam.C +++ b/applications/solvers/compressible/rhoPimpleFoam/rhoPorousMRFLTSPimpleFoam/rhoPorousMRFLTSPimpleFoam.C @@ -37,7 +37,7 @@ Description #include "fvCFD.H" #include "psiThermo.H" #include "turbulenceModel.H" -#include "MRFZones.H" +#include "IOMRFZoneList.H" #include "IOporosityModelList.H" #include "IObasicSourceList.H" #include "fvcSmooth.H" diff --git a/applications/solvers/compressible/rhoPimpleFoam/rhoPorousMRFPimpleFoam/createZones.H b/applications/solvers/compressible/rhoPimpleFoam/rhoPorousMRFPimpleFoam/createZones.H index d2522b8c75..f6c7e4dddf 100644 --- a/applications/solvers/compressible/rhoPimpleFoam/rhoPorousMRFPimpleFoam/createZones.H +++ b/applications/solvers/compressible/rhoPimpleFoam/rhoPorousMRFPimpleFoam/createZones.H @@ -1,4 +1,4 @@ - MRFZones mrfZones(mesh); + IOMRFZoneList mrfZones(mesh); mrfZones.correctBoundaryVelocity(U); IOporosityModelList pZones(mesh); diff --git a/applications/solvers/compressible/rhoPimpleFoam/rhoPorousMRFPimpleFoam/rhoPorousMRFPimpleFoam.C b/applications/solvers/compressible/rhoPimpleFoam/rhoPorousMRFPimpleFoam/rhoPorousMRFPimpleFoam.C index ac5dc1eadc..0975b0de2b 100644 --- a/applications/solvers/compressible/rhoPimpleFoam/rhoPorousMRFPimpleFoam/rhoPorousMRFPimpleFoam.C +++ b/applications/solvers/compressible/rhoPimpleFoam/rhoPorousMRFPimpleFoam/rhoPorousMRFPimpleFoam.C @@ -37,7 +37,7 @@ Description #include "psiThermo.H" #include "turbulenceModel.H" #include "bound.H" -#include "MRFZones.H" +#include "IOMRFZoneList.H" #include "IOporosityModelList.H" #include "IObasicSourceList.H" #include "pimpleControl.H" diff --git a/applications/solvers/compressible/rhoSimpleFoam/rhoPorousMRFSimpleFoam/createZones.H b/applications/solvers/compressible/rhoSimpleFoam/rhoPorousMRFSimpleFoam/createZones.H index 4eb2cb193c..614460f82b 100644 --- a/applications/solvers/compressible/rhoSimpleFoam/rhoPorousMRFSimpleFoam/createZones.H +++ b/applications/solvers/compressible/rhoSimpleFoam/rhoPorousMRFSimpleFoam/createZones.H @@ -1,4 +1,4 @@ - MRFZones mrfZones(mesh); + IOMRFZoneList mrfZones(mesh); mrfZones.correctBoundaryVelocity(U); IOporosityModelList pZones(mesh); diff --git a/applications/solvers/compressible/rhoSimpleFoam/rhoPorousMRFSimpleFoam/rhoPorousMRFSimpleFoam.C b/applications/solvers/compressible/rhoSimpleFoam/rhoPorousMRFSimpleFoam/rhoPorousMRFSimpleFoam.C index b552cadfce..e906c56f80 100644 --- a/applications/solvers/compressible/rhoSimpleFoam/rhoPorousMRFSimpleFoam/rhoPorousMRFSimpleFoam.C +++ b/applications/solvers/compressible/rhoSimpleFoam/rhoPorousMRFSimpleFoam/rhoPorousMRFSimpleFoam.C @@ -34,7 +34,7 @@ Description #include "fvCFD.H" #include "rhoThermo.H" #include "RASModel.H" -#include "MRFZones.H" +#include "IOMRFZoneList.H" #include "IObasicSourceList.H" #include "IOporosityModelList.H" #include "simpleControl.H" diff --git a/applications/solvers/heatTransfer/buoyantSimpleFoam/UEqn.H b/applications/solvers/heatTransfer/buoyantSimpleFoam/UEqn.H index ca28910aaf..e1593c3691 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 4a9387c188..fc47fc90a1 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 0000000000..08480c68d2 --- /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 f63e12f363..66ae75b650 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; diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionFoam.C b/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionFoam.C index 557017ccad..a0bd72c3b4 100644 --- a/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionFoam.C +++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionFoam.C @@ -60,14 +60,9 @@ int main(int argc, char *argv[]) #include "createFluidMeshes.H" #include "createSolidMeshes.H" - #include "createPorousFluidRegions.H" - #include "createPorousSolidMeshes.H" #include "createFluidFields.H" #include "createSolidFields.H" - #include "createPorousFluidFields.H" - #include "createPorousSolidFields.H" - #include "initContinuityErrs.H" #include "readTimeControls.H" @@ -116,24 +111,6 @@ int main(int argc, char *argv[]) #include "solveFluid.H" } - forAll(porousFluidRegions, i) - { - Info<< "\nSolving for fluid porous region " - << porousFluidRegions[i].name() << endl; - #include "setPorousFluidFields.H" - #include "readPorousFluidRegionPIMPLEControls.H" - #include "solvePorousFluid.H" - } - - forAll(porousSolidRegions, i) - { - Info<< "\nSolving for porous solid region " - << porousSolidRegions[i].name() << endl; - #include "setPorousRegionSolidFields.H" - #include "readPorousSolidMultiRegionPIMPLEControls.H" - #include "solvePorousSolid.H" - } - forAll(solidRegions, i) { Info<< "\nSolving for solid region " diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/Make/options b/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/Make/options index d5c1e47394..7418bfd694 100644 --- a/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/Make/options +++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/Make/options @@ -1,8 +1,7 @@ EXE_INC = \ -Ifluid \ -Isolid \ - -I./porousFluid \ - -I./porousSolid \ + -I../solid \ -I$(LIB_SRC)/finiteVolume/lnInclude \ -I$(LIB_SRC)/meshTools/lnInclude \ -I$(LIB_SRC)/finiteVolume/cfdTools \ diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/chtMultiRegionSimpleFoam.C b/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/chtMultiRegionSimpleFoam.C index 2f4b801cb3..4586122006 100644 --- a/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/chtMultiRegionSimpleFoam.C +++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/chtMultiRegionSimpleFoam.C @@ -50,13 +50,9 @@ int main(int argc, char *argv[]) #include "createFluidMeshes.H" #include "createSolidMeshes.H" - #include "createPorousFluidRegions.H" - #include "createPorousSolidMeshes.H" #include "createFluidFields.H" #include "createSolidFields.H" - #include "createPorousFluidFields.H" - #include "createPorousSolidFields.H" #include "initContinuityErrs.H" @@ -74,24 +70,6 @@ int main(int argc, char *argv[]) #include "solveFluid.H" } - forAll(porousFluidRegions, i) - { - Info<< "\nSolving for fluid porous region " - << porousFluidRegions[i].name() << endl; - #include "setPorousFluidFields.H" - #include "readPorousFluidRegionSIMPLEControls.H" - #include "solvePorousFluid.H" - } - - forAll(porousSolidRegions, i) - { - Info<< "\nSolving for porous solid region " - << porousSolidRegions[i].name() << endl; - #include "setPorousRegionSolidFields.H" - #include "readPorousSolidMultiRegionSIMPLEControls.H" - #include "solvePorousSolid.H" - } - forAll(solidRegions, i) { Info<< "\nSolving for solid region " diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/porousFluid/UPorousFluidEqn.H b/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/porousFluid/UPorousFluidEqn.H deleted file mode 100644 index 4c1081f5d4..0000000000 --- a/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/porousFluid/UPorousFluidEqn.H +++ /dev/null @@ -1,11 +0,0 @@ - // Solve the Momentum equation - tmp porousUEqn - ( - fvm::div(porousPhi, porousU) - + turbPorous.divDevRhoReff(porousU) - + porousSources(porousRho, porousU) - ); - - porousUEqn().relax(); - - solve(porousUEqn() == -fvc::grad(porousP)); diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/porousFluid/createPorousFluidFields.H b/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/porousFluid/createPorousFluidFields.H deleted file mode 100644 index 00762909f3..0000000000 --- a/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/porousFluid/createPorousFluidFields.H +++ /dev/null @@ -1,181 +0,0 @@ - // Initialise porous field pointer lists - PtrList thermoPorous(porousFluidRegions.size()); - PtrList rhoPorous(porousFluidRegions.size()); - PtrList kappaPorous(porousFluidRegions.size()); - PtrList UPorous(porousFluidRegions.size()); - PtrList phiPorous(porousFluidRegions.size()); - PtrList turbulencePorous - ( - porousFluidRegions.size() - ); - PtrList pPorous(porousFluidRegions.size()); - - List initialMassFluidPorous(porousFluidRegions.size()); - List