Commit 0daa3af2 authored by sergio's avatar sergio
Browse files

ENH: adding oversetAdjustPhi to overInterFoam plus reconstructing U from phi in pEq.

parent ea173eea
...@@ -23,3 +23,8 @@ bool ddtCorr ...@@ -23,3 +23,8 @@ bool ddtCorr
( (
pimple.dict().lookupOrDefault("ddtCorr", true) pimple.dict().lookupOrDefault("ddtCorr", true)
); );
bool adjustFringe
(
pimple.dict().lookupOrDefault("oversetAdjustPhi", false)
);
...@@ -53,6 +53,7 @@ Description ...@@ -53,6 +53,7 @@ Description
#include "interpolationCellPoint.H" #include "interpolationCellPoint.H"
#include "transform.H" #include "transform.H"
#include "fvMeshSubset.H" #include "fvMeshSubset.H"
#include "oversetAdjustPhi.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
...@@ -64,7 +65,7 @@ int main(int argc, char *argv[]) ...@@ -64,7 +65,7 @@ int main(int argc, char *argv[])
#include "createTime.H" #include "createTime.H"
#include "createDynamicFvMesh.H" #include "createDynamicFvMesh.H"
#include "initContinuityErrs.H" #include "initContinuityErrs.H"
#include "createControl.H" pimpleControl pimple(mesh);
#include "createTimeControls.H" #include "createTimeControls.H"
#include "createDyMControls.H" #include "createDyMControls.H"
#include "createFields.H" #include "createFields.H"
......
...@@ -38,6 +38,13 @@ ...@@ -38,6 +38,13 @@
phiHbyA += phig; phiHbyA += phig;
if (adjustFringe)
{
fvc::makeRelative(phiHbyA, U);
oversetAdjustPhi(phiHbyA, U);
fvc::makeAbsolute(phiHbyA, U);
}
// Update the pressure BCs to ensure flux consistency // Update the pressure BCs to ensure flux consistency
constrainPressure(p_rgh, U, phiHbyA, rAUf, MRF); constrainPressure(p_rgh, U, phiHbyA, rAUf, MRF);
...@@ -64,7 +71,8 @@ ...@@ -64,7 +71,8 @@
"minGradP", "minGradP",
fvc::reconstruct((phig - p_rghEqn.flux())/rAUf) fvc::reconstruct((phig - p_rghEqn.flux())/rAUf)
); );
U = HbyA + rAU*cellMask*minGradP; //U = HbyA + rAU*cellMask*minGradP;
U = fvc::reconstruct(phi);
U.correctBoundaryConditions(); U.correctBoundaryConditions();
fvOptions.correct(U); fvOptions.correct(U);
} }
......
...@@ -11,5 +11,6 @@ moveMeshOuterCorrectors = ...@@ -11,5 +11,6 @@ moveMeshOuterCorrectors =
massFluxInterpolation = massFluxInterpolation =
pimple.dict().lookupOrDefault("massFluxInterpolation", false); pimple.dict().lookupOrDefault("massFluxInterpolation", false);
ddtCorr = ddtCorr = pimple.dict().lookupOrDefault("ddtCorr", true);
pimple.dict().lookupOrDefault("ddtCorr", true);
\ No newline at end of file adjustFringe = pimple.dict().lookupOrDefault("oversetAdjustPhi", false);
...@@ -85,6 +85,8 @@ PIMPLE ...@@ -85,6 +85,8 @@ PIMPLE
ddtCorr yes; ddtCorr yes;
correctPhi no; correctPhi no;
oversetAdjustPhi no;
moveMeshOuterCorrectors no; moveMeshOuterCorrectors no;
turbOnFinalIterOnly no; turbOnFinalIterOnly no;
} }
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment