2 word alphaScheme(
"div(phi,alpha)");
7 fvScalarMatrix alpha1Eqn
10 + fv::gaussConvectionScheme<scalar>
20 Info<<
"Phase-1 volume fraction = "
22 <<
" Min(" <<
alpha1.name() <<
") = " << min(
alpha1).value()
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 = "
119 <<
" Min(" <<
alpha1.name() <<
") = " << min(
alpha1).value()
const volScalarField & alpha1
surfaceScalarField phir(fvc::flux(UdmModel.Udm()))
word alpharScheme("div(phirb,alpha)")
volScalarField alpha10("alpha10", alpha1)
tmp< surfaceScalarField > talphaPhiCorr0
Info<< "Creating temperaturePhaseChangeTwoPhaseMixture\n"<< endl;autoPtr< temperaturePhaseChangeTwoPhaseMixture > mixture
surfaceScalarField alphaPhi(phi.name()+alpha1.name(), fvc::flux(phi, alpha1, alphaScheme))