UrelEqn.H
Go to the documentation of this file.
1 // Relative momentum predictor
2
3 tmp<fvVectorMatrix> tUrelEqn
4 (
5 fvm::div(phi, Urel)
6 + turbulence->divDevReff(Urel)
7 + SRF->Su()
8 ==
10 );
11 fvVectorMatrix& UrelEqn = tUrelEqn.ref();
12
13 UrelEqn.relax();
14
15 fvOptions.constrain(UrelEqn);
16
17 if (simple.momentumPredictor())
18 {
19 solve(UrelEqn == -fvc::grad(p));
20
21 fvOptions.correct(Urel);
22 }
fv::options & fvOptions
surfaceScalarField & phi
volScalarField & p
compressible::turbulenceModel & turbulence
Urel
Definition: pEqn.H:56
tmp< fvVectorMatrix > tUrelEqn(fvm::ddt(Urel)+fvm::div(phi, Urel)+turbulence->divDevReff(Urel)+SRF->Su()==fvOptions(Urel))
fvVectorMatrix & UrelEqn
Definition: UrelEqn.H:11
const dictionary & simple
CEqn solve()
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, pimple.dict(), pRefCell, pRefValue);mesh.setFluxRequired(p.name());Info<< "Creating SRF model\n"<< endl;autoPtr< SRF::SRFModel > SRF(SRF::SRFModel::New(Urel))