Go to the documentation of this file.00001 Info<< "Reading field p\n" << endl;
00002 volScalarField p
00003 (
00004 IOobject
00005 (
00006 "p",
00007 runTime.timeName(),
00008 mesh,
00009 IOobject::MUST_READ,
00010 IOobject::AUTO_WRITE
00011 ),
00012 mesh
00013 );
00014
00015 Info<< "Reading field U\n" << endl;
00016 volVectorField U
00017 (
00018 IOobject
00019 (
00020 "U",
00021 runTime.timeName(),
00022 mesh,
00023 IOobject::MUST_READ,
00024 IOobject::AUTO_WRITE
00025 ),
00026 mesh
00027 );
00028
00029 # include <finiteVolume/createPhi.H>
00030
00031
00032 label pRefCell = 0;
00033 scalar pRefValue = 0.0;
00034 setRefCell(p, mesh.solutionDict().subDict("PISO"), pRefCell, pRefValue);
00035
00036
00037 singlePhaseTransportModel laminarTransport(U, phi);
00038
00039 autoPtr<incompressible::turbulenceModel> turbulence
00040 (
00041 incompressible::turbulenceModel::New(U, phi, laminarTransport)
00042 );
