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 rho0
    (
        "rho0",
        dimDensity,
        thermodynamicProperties
    );
 
    dimensionedScalar p0
    (
        "p0",
        dimPressure,
        thermodynamicProperties
    );
 
    dimensionedScalar psi
    (
        "psi",
        dimCompressibility,
        thermodynamicProperties
    );
 
    // Density offset, i.e. the constant part of the density
    dimensionedScalar rhoO("rhoO", rho0 - psi*p0);