Go to the source code of this file.
|
Info<< "Reading field U\n"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Creating face flux\n"<< endl;surfaceScalarField phi(IOobject("phi", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), mesh, dimensionedScalar("zero", mesh.Sf().dimensions()*U.dimensions(), 0.0));singlePhaseTransportModel laminarTransport(U, phi);autoPtr< incompressible::turbulenceModel > | turbulence (incompressible::turbulenceModel::New(U, phi, laminarTransport)) |
|
dimensionedVector | Ubar ("Ubar", dimVelocity, laminarTransport) |
|
dimensionedVector | gradP ("gradP", dimensionSet(0, 1,-2, 0, 0), Zero) |
|
Info<< "Reading field U\n" << endl;volVectorField U( IOobject ( "U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh);Info<< "Creating face flux\n" << endl;surfaceScalarField phi( IOobject ( "phi", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE ), mesh, dimensionedScalar("zero", mesh.Sf().dimensions()*U.dimensions(), 0.0));singlePhaseTransportModel laminarTransport(U, phi);autoPtr<incompressible::turbulenceModel> turbulence |
( |
incompressible::turbulenceModel::New(U, phi, laminarTransport) |
| ) |
|
dimensionedVector gradP |
( |
"gradP" |
, |
|
|
dimensionSet(0, 1,-2, 0, 0) |
, |
|
|
Zero |
|
|
) |
| |
vector flowDirection = (Ubar/mag(Ubar)).value() |