diff --git a/applications/solvers/combustion/chemFoam/chemFoam.C b/applications/solvers/combustion/chemFoam/chemFoam.C index d5038d1f69..4183045a10 100644 --- a/applications/solvers/combustion/chemFoam/chemFoam.C +++ b/applications/solvers/combustion/chemFoam/chemFoam.C @@ -69,14 +69,9 @@ int main(int argc, char *argv[]) Info<< "Time = " << runTime.timeName() << nl << endl; #include "solveChemistry.H" - - { - #include "YEqn.H" - - #include "hEqn.H" - - #include "pEqn.H" - } + #include "YEqn.H" + #include "hEqn.H" + #include "pEqn.H" #include "output.H" diff --git a/applications/solvers/combustion/chemFoam/hEqn.H b/applications/solvers/combustion/chemFoam/hEqn.H index 2478e6017c..4dd593ea9e 100644 --- a/applications/solvers/combustion/chemFoam/hEqn.H +++ b/applications/solvers/combustion/chemFoam/hEqn.H @@ -9,4 +9,6 @@ { h[0] = h0 + integratedHeat; } + + thermo.correct(); } diff --git a/applications/solvers/combustion/chemFoam/pEqn.H b/applications/solvers/combustion/chemFoam/pEqn.H index 13f3d603ae..a919ff198c 100644 --- a/applications/solvers/combustion/chemFoam/pEqn.H +++ b/applications/solvers/combustion/chemFoam/pEqn.H @@ -1,5 +1,4 @@ { - thermo.correct(); rho = thermo.rho(); if (constProp == "volume") { 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