readTransportProperties.H 534 Bytes
Newer Older
shunbo's avatar
shunbo committed
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
Info<< "\nReading transportProperties" << endl;

IOdictionary transportProperties
(
    IOobject
    (
        "transportProperties",
        runTime.constant(),
        mesh,
        IOobject::MUST_READ,
        IOobject::NO_WRITE
    )
);

dimensionedScalar c0("c0", dimVelocity, transportProperties);

dimensionedScalar rho("rho", dimDensity, transportProperties);

scalar MaxCo =
    max(mesh.surfaceInterpolation::deltaCoeffs()*c0).value()
   *runTime.deltaT().value();

Info<< "Max acoustic Courant Number = " << MaxCo << endl;