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;
51 tmp<surfaceScalarField> talphaPhiUn
69 tmp<surfaceScalarField> talphaPhiCorr(talphaPhiUn() -
alphaPhi);
114 Info<<
"Phase-1 volume fraction = "
fvMatrix< scalar > fvScalarMatrix
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
word alpharScheme("div(phirb,alpha)")
tmp< surfaceScalarField > talphaPhiCorr0
Ostream & endl(Ostream &os)
Add newline and flush stream.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
volScalarField alpha10("alpha10", alpha1)
Info<< "Predicted p max-min : "<< max(p).value()<< " "<< min(p).value()<< endl;rho==max(psi *p+alphal *rhol0+((alphav *psiv+alphal *psil) - psi) *pSat, rhoMin);# 1 "/home/ubuntu/OpenFOAM-5.0/applications/solvers/multiphase/cavitatingFoam/alphavPsi.H" 1{ alphav=max(min((rho - rholSat)/(rhovSat - rholSat), scalar(1)), scalar(0));alphal=1.0 - alphav;Info<< "max-min alphav: "<< max(alphav).value()<< " "<< min(alphav).value()<< endl;psiModel-> correct()
Info<< "Reading field p_rgh\"<< endl;volScalarField p_rgh(IOobject("p_rgh", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Reading field U\"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Creating phaseChangeTwoPhaseMixture\"<< endl;autoPtr< phaseChangeTwoPhaseMixture > mixture
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
surfaceScalarField phir(IOobject("phir", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), mixture.cAlpha() *mag(phi/mesh.magSf()) *mixture.nHatf())
tmp< surfaceScalarField > flux(const volVectorField &vvf)
Return the face-flux field obtained from the given volVectorField.
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)