Skip to content
Snippets Groups Projects
Commit 03c8cf06 authored by Henry's avatar Henry
Browse files

twoPhaseEulerFoam: Update flux naming convention

parent 0156c654
Branches
Tags
No related merge requests found
{ {
label nAlphaCorr(readLabel(pimple.dict().lookup("nAlphaCorr")));
label nAlphaSubCycles(readLabel(pimple.dict().lookup("nAlphaSubCycles")));
word alphaScheme("div(phi,alpha1)"); word alphaScheme("div(phi,alpha1)");
word alpharScheme("div(phir,alpha1)"); word alpharScheme("div(phir,alpha1)");
...@@ -24,7 +21,7 @@ ...@@ -24,7 +21,7 @@
!(++alphaSubCycle).end(); !(++alphaSubCycle).end();
) )
{ {
surfaceScalarField phiAlpha surfaceScalarField alphaPhic1
( (
fvc::flux fvc::flux
( (
...@@ -44,7 +41,7 @@ ...@@ -44,7 +41,7 @@
( (
alpha1, alpha1,
phi, phi,
phiAlpha, alphaPhic1,
(g0.value() > 0 ? alphaMax : 1), (g0.value() > 0 ? alphaMax : 1),
0 0
); );
...@@ -54,11 +51,6 @@ ...@@ -54,11 +51,6 @@
{ {
surfaceScalarField alpha1f(fvc::interpolate(alpha1)); surfaceScalarField alpha1f(fvc::interpolate(alpha1));
// ppMagf = rAU1f*fvc::interpolate
// (
// (1.0/(rho1*(alpha1 + scalar(0.0001))))
// *g0*min(exp(preAlphaExp*(alpha1 - alphaMax)), expMax)
// );
ppMagf = ppMagf =
rAU1f/(alpha1f + scalar(0.0001)) rAU1f/(alpha1f + scalar(0.0001))
*(g0/rho1)*min(exp(preAlphaExp*(alpha1f - alphaMax)), expMax); *(g0/rho1)*min(exp(preAlphaExp*(alpha1f - alphaMax)), expMax);
......
#include "readTimeControls.H" #include "readTimeControls.H"
int nAlphaCorr(readInt(pimple.dict().lookup("nAlphaCorr")));
int nAlphaSubCycles(readInt(pimple.dict().lookup("nAlphaSubCycles")));
Switch correctAlpha(pimple.dict().lookup("correctAlpha")); Switch correctAlpha(pimple.dict().lookup("correctAlpha"));
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment