diff --git a/applications/solvers/multiphase/interFoam/overInterDyMFoam/createDyMControls.H b/applications/solvers/multiphase/interFoam/overInterDyMFoam/createDyMControls.H index cb2dfd5ddd1e49261c47ed18216c802b08afdca0..00673cec243ac65fae2c8bf04534d24ddf8df3b2 100644 --- a/applications/solvers/multiphase/interFoam/overInterDyMFoam/createDyMControls.H +++ b/applications/solvers/multiphase/interFoam/overInterDyMFoam/createDyMControls.H @@ -19,6 +19,11 @@ bool massFluxInterpolation pimple.dict().lookupOrDefault("massFluxInterpolation", false) ); +bool adjustFringe +( + pimple.dict().lookupOrDefault("oversetAdjustPhi", false) +); + bool ddtCorr ( pimple.dict().lookupOrDefault("ddtCorr", true) diff --git a/applications/solvers/multiphase/interFoam/overInterDyMFoam/interDyMFoam.C b/applications/solvers/multiphase/interFoam/overInterDyMFoam/interDyMFoam.C index 7bb7f8af6b87112aa1b2bbc93fb68a3582f184b7..075676da7715f829c61600483e7bbb533c9eff9f 100644 --- a/applications/solvers/multiphase/interFoam/overInterDyMFoam/interDyMFoam.C +++ b/applications/solvers/multiphase/interFoam/overInterDyMFoam/interDyMFoam.C @@ -53,6 +53,7 @@ Description #include "interpolationCellPoint.H" #include "transform.H" #include "fvMeshSubset.H" +#include "oversetAdjustPhi.H" // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // diff --git a/applications/solvers/multiphase/interFoam/overInterDyMFoam/pEqn.H b/applications/solvers/multiphase/interFoam/overInterDyMFoam/pEqn.H index 0a5b0296b8b14d0b5dd4d94235636202d8ee81be..c11ce1dba5cc0c78b125a7e75d39c85124422386 100644 --- a/applications/solvers/multiphase/interFoam/overInterDyMFoam/pEqn.H +++ b/applications/solvers/multiphase/interFoam/overInterDyMFoam/pEqn.H @@ -28,6 +28,14 @@ fvc::makeAbsolute(phiHbyA, U); } + if (adjustFringe) + { + fvc::makeRelative(phiHbyA, U); + oversetAdjustPhi(phiHbyA, U); + fvc::makeAbsolute(phiHbyA, U); + } + + surfaceScalarField phig ( (