1 Info<< "Reading thermophysical properties\n" << endl;
3 autoPtr<basicPsiThermo> pThermo
5 basicPsiThermo::New(mesh)
7 basicPsiThermo& thermo = pThermo();
9 volScalarField& p = thermo.p();
10 volScalarField& h = thermo.h();
11 const volScalarField& psi = thermo.psi();
20 IOobject::READ_IF_PRESENT,
26 Info<< "Reading field U\n" << endl;
40 #include "compressibleCreatePhi.H"
42 dimensionedScalar pMin
44 mesh.solutionDict().subDict("PIMPLE").lookup("pMin")
47 Info<< "Creating turbulence model\n" << endl;
48 autoPtr<compressible::turbulenceModel> turbulence
50 compressible::turbulenceModel::New
59 //dimensionedScalar initialMass = fvc::domainIntegrate(rho);
62 Info<< "Creating field DpDt\n" << endl;
64 fvc::DDt(surfaceScalarField("phiU", phi/fvc::interpolate(rho)), p);
66 MRFZones mrfZones(mesh);
67 mrfZones.correctBoundaryVelocity(U);
69 porousZones pZones(mesh);
70 Switch pressureImplicitPorosity(false);