Go to the documentation of this file.00001 Info<< "Reading thermophysical properties\n" << endl;
00002
00003 autoPtr<basicPsiThermo> pThermo
00004 (
00005 basicPsiThermo::New(mesh)
00006 );
00007 basicPsiThermo& thermo = pThermo();
00008
00009 volScalarField& p = thermo.p();
00010 volScalarField& h = thermo.h();
00011 const volScalarField& psi = thermo.psi();
00012
00013 volScalarField rho
00014 (
00015 IOobject
00016 (
00017 "rho",
00018 runTime.timeName(),
00019 mesh,
00020 IOobject::NO_READ,
00021 IOobject::AUTO_WRITE
00022 ),
00023 thermo.rho()
00024 );
00025
00026 Info<< "\nReading field U\n" << endl;
00027 volVectorField U
00028 (
00029 IOobject
00030 (
00031 "U",
00032 runTime.timeName(),
00033 mesh,
00034 IOobject::MUST_READ,
00035 IOobject::AUTO_WRITE
00036 ),
00037 mesh
00038 );
00039
00040 # include <finiteVolume/compressibleCreatePhi.H>
00041
00042
00043 Info<< "Creating turbulence model\n" << endl;
00044 autoPtr<compressible::turbulenceModel> turbulence
00045 (
00046 compressible::turbulenceModel::New
00047 (
00048 rho,
00049 U,
00050 phi,
00051 thermo
00052 )
00053 );
00054
00055
00056 Info<< "Creating field DpDt\n" << endl;
00057 volScalarField DpDt =
00058 fvc::DDt(surfaceScalarField("phiU", phi/fvc::interpolate(rho)), p);