alphavPsi.H
Go to the documentation of this file.
1 {
2  alphav =
3  max
4  (
5  min
6  (
7  (rho - rholSat)/(rhovSat - rholSat),
8  scalar(1)
9  ),
10  scalar(0)
11  );
12  alphal = 1.0 - alphav;
13 
14  Info<< "max-min alphav: " << max(alphav).value()
15  << " " << min(alphav).value() << endl;
16 
17  psiModel->correct();
18 
19  //Info<< "min a: " << 1.0/sqrt(max(psi)).value() << endl;
20 }
Foam::endl
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:369
rho
rho
Definition: readInitialConditions.H:88
Foam::min
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
Definition: hashSets.C:33
Foam::Info
messageStream Info
Information stream (stdout output on master, null elsewhere)
Foam::max
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
Definition: hashSets.C:47
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
alphal
alphal
Definition: alphavPsi.H:12