zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

BUG: interIsoFoam nAlphaSubCycles, nOuterCorrectors issue (fixes #1300)

Browse Source 
- failed when nAlphaSubCycles > 1 and nOuterCorrectors > 1.
  Fixed thanks to Henning Scheufler (DLR Bremen).
This commit is contained in:
Johan Roenby
2019-04-30 10:59:32 +02:00
committed by Andrew Heather
parent ed6b26ba45
commit d8e8306bbc
2 changed files with 36 additions and 37 deletions
Download Patch File Download Diff File Expand all files Collapse all files

37
applications/solvers/multiphase/interIsoFoam/alphaEqn.H
View File

@ -1,48 +1,11 @@
{ {
if (pimple.nCorrPIMPLE() > 1)
{
// If nOuterCorrectors > 1 then for all but the first loop the advection
// of alpha is done using an average, 0.5*phi+0.5*phiNew where phi is
// the flux at the beginning of the time step and phiNew is the flux
// estimate at the end of the time step from the previous outer
// iteration. Similarly we use 0.5*U + 0.5*UNew in later iterations.
if (pimple.firstIter())
{
// To recalculate the alpha1 update in subsequent iterations, we
// must store its current value before overwriting with the new
// value
alpha1.storePrevIter();
// Storing initial phi and U for use in later outer iterations.
phi.storePrevIter();
U.storePrevIter();
}
else
{
// Resetting alpha1 to value before advection in first PIMPLE
// iteration.
alpha1 = alpha1.prevIter();
// Setting U and phi with which to advect interface.
U = 0.5*U.prevIter() + 0.5*U;
phi = 0.5*phi.prevIter() + 0.5*phi;
}
}
// Updating alpha1 // Updating alpha1
advector.advect(); advector.advect();
#include "rhofs.H" #include "rhofs.H"
rhoPhi = advector.getRhoPhi(rho1f, rho2f); rhoPhi = advector.getRhoPhi(rho1f, rho2f);
if (!pimple.firstIter())
{
// Resetting U and phi to value at latest iteration.
U = 2.0*U - U.prevIter();
phi = 2.0*phi - phi.prevIter();
}
alpha2 = 1.0 - alpha1; alpha2 = 1.0 - alpha1;
mixture.correct(); mixture.correct();
} }
Info<< "Phase-1 volume fraction = " Info<< "Phase-1 volume fraction = "

36
applications/solvers/multiphase/interIsoFoam/alphaEqnSubCycle.H
View File

@ -1,3 +1,32 @@
if (pimple.nCorrPIMPLE() > 1)
{
// If nOuterCorrectors > 1 then for all but the first loop the advection
// of alpha is done using an average, 0.5*phi+0.5*phiNew where phi is
// the flux at the beginning of the time step and phiNew is the flux
// estimate at the end of the time step from the previous outer
// iteration. Similarly we use 0.5*U + 0.5*UNew in later iterations.
if (pimple.firstIter())
{
// To recalculate the alpha1 update in subsequent iterations, we
// must store its current value before overwriting with the new
// value
alpha1.storePrevIter();
// Storing initial phi and U for use in later outer iterations.
phi.storePrevIter();
U.storePrevIter();
}
else
{
// Resetting alpha1 to value before advection in first PIMPLE
// iteration.
alpha1 = alpha1.prevIter();
// Temporarily setting U and phi with which to advect interface.
U = 0.5*U.prevIter() + 0.5*U;
phi = 0.5*phi.prevIter() + 0.5*phi;
}
}
if (nAlphaSubCycles > 1) if (nAlphaSubCycles > 1)
{ {
dimensionedScalar totalDeltaT = runTime.deltaT(); dimensionedScalar totalDeltaT = runTime.deltaT();
@ -30,4 +59,11 @@ else
#include "alphaEqn.H" #include "alphaEqn.H"
} }
if (!pimple.firstIter())
{
// Resetting U and phi to value at latest iteration.
U = 2.0*U - U.prevIter();
phi = 2.0*phi - phi.prevIter();
}
rho == alpha1*rho1 + alpha2*rho2; rho == alpha1*rho1 + alpha2*rho2;