diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/alphaEqn.H b/applications/solvers/multiphase/twoPhaseEulerFoam/alphaEqn.H index ddb923cf87..ae21059446 100644 --- a/applications/solvers/multiphase/twoPhaseEulerFoam/alphaEqn.H +++ b/applications/solvers/multiphase/twoPhaseEulerFoam/alphaEqn.H @@ -1,7 +1,4 @@ { - label nAlphaCorr(readLabel(pimple.dict().lookup("nAlphaCorr"))); - label nAlphaSubCycles(readLabel(pimple.dict().lookup("nAlphaSubCycles"))); - word alphaScheme("div(phi,alpha1)"); word alpharScheme("div(phir,alpha1)"); @@ -24,7 +21,7 @@ !(++alphaSubCycle).end(); ) { - surfaceScalarField phiAlpha + surfaceScalarField alphaPhic1 ( fvc::flux ( @@ -44,7 +41,7 @@ ( alpha1, phi, - phiAlpha, + alphaPhic1, (g0.value() > 0 ? alphaMax : 1), 0 ); @@ -54,11 +51,6 @@ { surfaceScalarField alpha1f(fvc::interpolate(alpha1)); - // ppMagf = rAU1f*fvc::interpolate - // ( - // (1.0/(rho1*(alpha1 + scalar(0.0001)))) - // *g0*min(exp(preAlphaExp*(alpha1 - alphaMax)), expMax) - // ); ppMagf = rAU1f/(alpha1f + scalar(0.0001)) *(g0/rho1)*min(exp(preAlphaExp*(alpha1f - alphaMax)), expMax); diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/readTwoPhaseEulerFoamControls.H b/applications/solvers/multiphase/twoPhaseEulerFoam/readTwoPhaseEulerFoamControls.H index 4e6847debe..1abe3ef10d 100644 --- a/applications/solvers/multiphase/twoPhaseEulerFoam/readTwoPhaseEulerFoamControls.H +++ b/applications/solvers/multiphase/twoPhaseEulerFoam/readTwoPhaseEulerFoamControls.H @@ -1,3 +1,5 @@ #include "readTimeControls.H" + int nAlphaCorr(readInt(pimple.dict().lookup("nAlphaCorr"))); + int nAlphaSubCycles(readInt(pimple.dict().lookup("nAlphaSubCycles"))); Switch correctAlpha(pimple.dict().lookup("correctAlpha"));