diff --git a/applications/solvers/heatTransfer/buoyantSimpleFoam/UEqn.H b/applications/solvers/heatTransfer/buoyantSimpleFoam/UEqn.H index ca28910aaf3c51186663b2fb5919104d4fd07398..e1593c3691cda41ce18ded23b5925b5813bcb9ab 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 4a9387c1882b17967d218a95e68553f29d52b64c..fc47fc90a1cc25556e6e52077db6976588ccf1a8 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 0000000000000000000000000000000000000000..08480c68d2d6f4c9efc89409a46ed1abe1042ad9 --- /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 f63e12f363774a88972459f748bd5162501fc0a4..66ae75b650e22a69f980c41e5b1365658c05d43a 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;