readThermodynamicProperties.H

Go to the documentation of this file.

Info<< "Reading thermodynamicProperties\n" << endl;
 
IOdictionary thermodynamicProperties
(
    IOobject
    (
        "thermodynamicProperties",
        runTime.constant(),
        mesh,
        IOobject::MUST_READ_IF_MODIFIED,
        IOobject::NO_WRITE
    )
);
 
dimensionedScalar psil("psil", dimCompressibility, thermodynamicProperties);
 
dimensionedScalar rholSat("rholSat", dimDensity, thermodynamicProperties);
 
dimensionedScalar psiv("psiv", dimCompressibility, thermodynamicProperties);
 
dimensionedScalar pSat("pSat", dimPressure, thermodynamicProperties);
 
dimensionedScalar rhovSat("rhovSat", psiv*pSat);
 
dimensionedScalar rhol0("rhol0", rholSat - pSat*psil);
 
const dimensionedScalar rhoMin("rhoMin", dimDensity, thermodynamicProperties);