Removed unnecessary return statement
[foam-extend-3.2.git] / applications / solvers / compressible / rhoPorousMRFPimpleFoam / createFields.H
blobb9a86ef9957b30923c2b3dd66606dcd8a33159e6
1     Info<< "Reading thermophysical properties\n" << endl;
3     autoPtr<basicPsiThermo> pThermo
4     (
5         basicPsiThermo::New(mesh)
6     );
7     basicPsiThermo& thermo = pThermo();
9     volScalarField& p = thermo.p();
10     volScalarField& h = thermo.h();
11     const volScalarField& psi = thermo.psi();
13     volScalarField rho
14     (
15         IOobject
16         (
17             "rho",
18             runTime.timeName(),
19             mesh,
20             IOobject::READ_IF_PRESENT,
21             IOobject::AUTO_WRITE
22         ),
23         thermo.rho()
24     );
26     Info<< "Reading field U\n" << endl;
27     volVectorField U
28     (
29         IOobject
30         (
31             "U",
32             runTime.timeName(),
33             mesh,
34             IOobject::MUST_READ,
35             IOobject::AUTO_WRITE
36         ),
37         mesh
38     );
40     #include "compressibleCreatePhi.H"
42     dimensionedScalar pMin
43     (
44         mesh.solutionDict().subDict("PIMPLE").lookup("pMin")
45     );
47     Info<< "Creating turbulence model\n" << endl;
48     autoPtr<compressible::turbulenceModel> turbulence
49     (
50         compressible::turbulenceModel::New
51         (
52             rho,
53             U,
54             phi,
55             thermo
56         )
57     );
59     //dimensionedScalar initialMass = fvc::domainIntegrate(rho);
62     Info<< "Creating field DpDt\n" << endl;
63     volScalarField DpDt =
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);