2 word alphaScheme(
"div(phi,alpha)");
10 + fv::gaussConvectionScheme<scalar>
20 Info<<
"Phase-1 volume fraction = " 26 tmp<surfaceScalarField> talphaPhiUD(alpha1Eqn.flux());
31 Info<<
"Applying the previous iteration correction flux" <<
endl;
39 UniformField<scalar>(
mixture.alphaMax()),
52 tmp<surfaceScalarField> talphaPhiUn
70 tmp<surfaceScalarField> talphaPhiCorr(talphaPhiUn() -
alphaPhi);
79 UniformField<scalar>(
mixture.alphaMax()),
104 UniformField<scalar>(
mixture.alphaMax()),
117 Info<<
"Phase-1 volume fraction = " fvMatrix< scalar > fvScalarMatrix
tmp< surfaceScalarField > flux(const volVectorField &vvf)
Return the face-flux field obtained from the given volVectorField.
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
word alpharScheme("div(phirb,alpha)")
tmp< surfaceScalarField > talphaPhiCorr0
Info<< "Creating temperaturePhaseChangeTwoPhaseMixture\"<< endl;autoPtr< temperaturePhaseChangeTwoPhaseMixture > mixture
Ostream & endl(Ostream &os)
Add newline and flush stream.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
volScalarField alpha10("alpha10", alpha1)
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
surfaceScalarField phir(IOobject("phir", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), mixture.cAlpha() *mag(phi/mesh.magSf()) *mixture.nHatf())
messageStream Info
Information stream (stdout output on master, null elsewhere)
surfaceScalarField alphaPhi(phi.name()+alpha1.name(), fvc::flux(phi, alpha1, alphaScheme))
void explicitSolve(const RdeltaTType &rDeltaT, const RhoType &rho, volScalarField &psi, const surfaceScalarField &phiPsi, const SpType &Sp, const SuType &Su)
const volScalarField & alpha1