createAlphaFluxes.H

Go to the documentation of this file.

IOobject alphaPhi10Header
(
    IOobject::groupName("alphaPhi0", alpha1.group()),
    runTime.timeName(),
    mesh,
    IOobject::READ_IF_PRESENT,
    IOobject::AUTO_WRITE
);
 
const bool alphaRestart =
    alphaPhi10Header.typeHeaderOk<surfaceScalarField>(true);
 
if (alphaRestart)
{
    Info << "Restarting alpha" << endl;
}
 
// MULES flux from previous time-step
surfaceScalarField alphaPhi10
(
    alphaPhi10Header,
    phi*fvc::interpolate(alpha1)
);
 
// MULES Correction
tmp<surfaceScalarField> talphaPhi1Corr0;