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
00016 Info<< "Reading field U\n" << endl;
00017 volVectorField U
00018 (
00019 IOobject
00020 (
00021 "U",
00022 runTime.timeName(),
00023 mesh,
00024 IOobject::MUST_READ,
00025 IOobject::AUTO_WRITE
00026 ),
00027 mesh
00028 );
00029
00030 # include <finiteVolume/createPhi.H>
00031
00032
00033 label pRefCell = 0;
00034 scalar pRefValue = 0.0;
00035 setRefCell(p, mesh.solutionDict().subDict("PISO"), pRefCell, pRefValue);
00036
00037
00038 singlePhaseTransportModel laminarTransport(U, phi);
00039
00040 autoPtr<incompressible::LESModel> sgsModel
00041 (
00042 incompressible::LESModel::New(U, phi, laminarTransport)
00043 );
