OpenFOAM: applications/solvers/incompressible/pimpleFoam/SRFPimpleFoam/UrelEqn.H Source File

OpenFOAM v3.0plus

The open source CFD toolbox

Go to the documentation of this file.

     // Relative momentum predictor
     tmp<fvVectorMatrix> UrelEqn
     (
         fvm::ddt(Urel)
       + fvm::div(phi, Urel)
       + turbulence->divDevReff(Urel)
       + SRF->Su()
      ==
         fvOptions(Urel)
     );
  
     UrelEqn().relax();
  
     fvOptions.constrain(UrelEqn());
  
     solve(UrelEqn() == -fvc::grad(p));
  
     fvOptions.correct(Urel);

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

Definition: pEqn.H:62

surfaceScalarField & phi

Definition: setRegionFluidFields.H:8

tmp< fvVectorMatrix > UrelEqn(fvm::div(phi, Urel)+turbulence->divDevReff(Urel)+SRF->Su()==fvOptions(Urel))

tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)

Definition: fvcDiv.C:47

autoPtr< compressible::turbulenceModel > turbulence

Definition: createFields.H:23

Urel

Definition: pEqn.H:45

solve(UrelEqn()==-fvc::grad(p))

tmp< GeometricField< Type, fvPatchField, volMesh > > ddt(const dimensioned< Type > dt, const fvMesh &mesh)

Definition: fvcDdt.C:45

fv::options & fvOptions

Definition: setRegionFluidFields.H:16