UrelEqn.H
Go to the documentation of this file.
1  // Relative momentum predictor
2 
3  tmp<fvVectorMatrix> UrelEqn
4  (
5  fvm::div(phi, Urel)
6  + turbulence->divDevReff(Urel)
7  + SRF->Su()
8  ==
9  fvOptions(Urel)
10  );
11 
12  UrelEqn().relax();
13 
14  fvOptions.constrain(UrelEqn());
15 
16  solve(UrelEqn() == -fvc::grad(p));
17 
18  fvOptions.correct(Urel);
SRF
Info<< "Reading field p\n"<< endl;volScalarField p(IOobject("p", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Reading field Urel\n"<< endl;volVectorField Urel(IOobject("Urel", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Reading/calculating face flux field phi\n"<< endl;surfaceScalarField phi(IOobject("phi", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), linearInterpolate(Urel) &mesh.Sf());label pRefCell=0;scalar pRefValue=0.0;setRefCell(p, simple.dict(), pRefCell, pRefValue);mesh.setFluxRequired(p.name());Info<< "Creating SRF model\n"<< endl;autoPtr< SRF::SRFModel > SRF(SRF::SRFModel::New(Urel))
Foam::fvc::grad
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcGrad.C:52
p
p
Definition: pEqn.H:62
phi
surfaceScalarField & phi
Definition: setRegionFluidFields.H:8
UrelEqn
tmp< fvVectorMatrix > UrelEqn(fvm::div(phi, Urel)+turbulence->divDevReff(Urel)+SRF->Su()==fvOptions(Urel))
Foam::fvc::div
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcDiv.C:47
turbulence
autoPtr< compressible::turbulenceModel > turbulence
Definition: createFields.H:23
Urel
Urel
Definition: pEqn.H:45
solve
solve(UrelEqn()==-fvc::grad(p))
fvOptions
fv::options & fvOptions
Definition: setRegionFluidFields.H:16

Copyright © 2011-2018  OpenFOAM | OPENFOAM® is a registered trademark of OpenCFD Ltd.

Generated by   doxygen 1.8.17