FreeFOAM The Cross-Platform CFD Toolkit
createFields.H
Go to the documentation of this file.
1  Info<< "Reading field p\n" << endl;
3  (
4  IOobject
5  (
6  "p",
7  runTime.timeName(),
8  mesh,
9  IOobject::MUST_READ,
10  IOobject::AUTO_WRITE
11  ),
12  mesh
13  );
14 
15 
16  Info<< "Reading field U\n" << endl;
18  (
19  IOobject
20  (
21  "U",
22  runTime.timeName(),
23  mesh,
24  IOobject::MUST_READ,
25  IOobject::AUTO_WRITE
26  ),
27  mesh
28  );
29 
30 # include <finiteVolume/createPhi.H>
31 
32 
33  label pRefCell = 0;
34  scalar pRefValue = 0.0;
35  setRefCell(p, mesh.solutionDict().subDict("PIMPLE"), pRefCell, pRefValue);
36 
37  singlePhaseTransportModel laminarTransport(U, phi);
38 
39  autoPtr<incompressible::turbulenceModel> turbulence
40  (
41  incompressible::turbulenceModel::New(U, phi, laminarTransport)
42  );
43 
44  Info<< "Reading field rAU if present\n" << endl;
45  volScalarField rAU
46  (
47  IOobject
48  (
49  "rAU",
50  runTime.timeName(),
51  mesh,
52  IOobject::READ_IF_PRESENT,
53  IOobject::AUTO_WRITE
54  ),
55  mesh,
56  runTime.deltaT(),
57  zeroGradientFvPatchScalarField::typeName
58  );
59 
60 // ************************ vim: set sw=4 sts=4 et: ************************ //