alphavPsi.H
Go to the documentation of this file.
1 {
2  alphav = clamp((rho - rholSat)/(rhovSat - rholSat), zero_one{});
3  alphal = 1.0 - alphav;
4 
5  Info<< "max-min alphav: " << max(alphav).value()
6  << " " << min(alphav).value() << endl;
7 
8  psiModel->correct();
9 
10  //Info<< "min a: " << 1.0/sqrt(max(psi)).value() << endl;
11 }
Foam::max
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
Definition: hashSets.C:40
Foam::endl
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:531
alphal
alphal
Definition: alphavPsi.H:3
Foam::min
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
Definition: hashSets.C:26
psiModel
Info<< "Reading field p\"<< 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\"<< 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\"<< endl;incompressibleTwoPhaseMixture mixture(U, phi);volScalarField &alphav(mixture.alpha1());alphav.oldTime();volScalarField &alphal(mixture.alpha2());Info<< "Creating compressibilityModel\"<< endl;autoPtr< barotropicCompressibilityModel > psiModel
Definition: createFields.H:72
Foam::Info
messageStream Info
Information stream (stdout output on master, null elsewhere)
Foam::clamp
dimensionSet clamp(const dimensionSet &a, const dimensionSet &range)
Definition: dimensionSet.C:271