FreeFOAM The Cross-Platform CFD Toolkit
createFields.H
Go to the documentation of this file.
1  Info<< "Reading thermophysical properties\n" << endl;
2 
3  autoPtr<basicPsiThermo> pThermo
4  (
5  basicPsiThermo::New(mesh)
6  );
7  basicPsiThermo& thermo = pThermo();
8 
9  volScalarField& p = thermo.p();
10  volScalarField& h = thermo.h();
11  const volScalarField& psi = thermo.psi();
12 
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  );
25 
26  Info<< "Reading field U\n" << endl;
28  (
29  IOobject
30  (
31  "U",
32  runTime.timeName(),
33  mesh,
34  IOobject::MUST_READ,
35  IOobject::AUTO_WRITE
36  ),
37  mesh
38  );
39 
41 
43  (
44  mesh.solutionDict().subDict("PIMPLE").lookup("pMin")
45  );
46 
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  );
58 
59  //dimensionedScalar initialMass = fvc::domainIntegrate(rho);
60 
61 
62  Info<< "Creating field DpDt\n" << endl;
65 
66 // ************************ vim: set sw=4 sts=4 et: ************************ //