alphavPsi_8H_source.html

{

    alphav = clamp((rho - rholSat)/(rhovSat - rholSat), zero_one{});

    alphal = 1.0 - alphav;


    Info<< "max-min alphav: " << max(alphav).value()

        << " " << min(alphav).value() << endl;


    psiModel->correct();


    //Info<< "min a: " << 1.0/sqrt(max(psi)).value() << endl;

}

alphal
alphal
Definition alphavPsi.H:3

Foam::clamp
dimensionSet clamp(const dimensionSet &a, const dimensionSet &range)
Definition dimensionSet.C:271

rho
rho
Definition readInitialConditions.H:88

psiModel
Info<< "Reading field p\n"<< endl;volScalarField p(IOobject("p", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);volScalarField rho(IOobject("rho", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Reading field U\n"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);mesh.setFluxRequired(p.name());surfaceScalarField rhoPhi(IOobject("rhoPhi", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), fvc::interpolate(rho) *phi);Info<< "Reading transportProperties\n"<< endl;incompressibleTwoPhaseMixture mixture(U, phi);volScalarField &alphav(mixture.alpha1());alphav.oldTime();volScalarField &alphal(mixture.alpha2());Info<< "Creating compressibilityModel\n"<< endl;autoPtr< barotropicCompressibilityModel > psiModel
Definition createFields.H:72