UEqn.H
Go to the documentation of this file.
1  fvVectorMatrix UEqn
2  (
3  betav*fvm::ddt(rho, U) + fvm::div(phi, U)
4  + MRF.DDt(rho, U)
5  + turbulence->divDevTau(U)
6  ==
7  betav*rho*g
8  + betav*fvModels.source(rho, U)
9  );
10 
11  fvConstraints.constrain(UEqn);
12 
13  volSymmTensorField invA(inv(I*UEqn.A() + drag->Dcu()));
14 
15  if (pimple.momentumPredictor())
16  {
17  U = invA & (UEqn.H() - betav*fvc::grad(p));
18  U.correctBoundaryConditions();
19  fvConstraints.constrain(U);
20  K = 0.5*magSqr(U);
21  }
Foam::fvConstraints::constrain
bool constrain(fvMatrix< Type > &eqn) const
Apply constraints to an equation.
Definition: fvConstraintsTemplates.C:29
invA
volSymmTensorField invA(inv(I *UEqn.A()+drag->Dcu()))
UEqn
fvVectorMatrix & UEqn
Definition: UEqn.H:11
fvConstraints
Foam::fvConstraints & fvConstraints(Foam::fvConstraints::New(mesh))
fvModels
Foam::fvModels & fvModels(Foam::fvModels::New(mesh))
MRF
IOMRFZoneList MRF(mesh)
pimple
pimpleControl pimple(mesh)
drag
Info<< "Reading strained laminar flame speed field Su\n"<< endl;volScalarField Su(IOobject("Su", runTime.name(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);Info<< "Reading field betav\n"<< endl;volScalarField betav(IOobject("betav", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE), mesh);Info<< "Reading field Lobs\n"<< endl;volScalarField Lobs(IOobject("Lobs", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE), mesh);Info<< "Reading field CT\n"<< endl;volSymmTensorField CT(IOobject("CT", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE), mesh);Info<< "Reading field Nv\n"<< endl;volScalarField Nv(IOobject("Nv", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE), mesh);Info<< "Reading field nsv\n"<< endl;volSymmTensorField nsv(IOobject("nsv", mesh.facesInstance(), mesh, IOobject::MUST_READ, IOobject::NO_WRITE), mesh);IOdictionary PDRProperties(IOobject("PDRProperties", runTime.constant(), mesh, IOobject::MUST_READ_IF_MODIFIED, IOobject::NO_WRITE));autoPtr< PDRDragModel > drag
Definition: createFields.H:194
K
K
Definition: pEqn.H:75
U
U
Definition: pEqn.H:72
Foam::fvc::ddt
tmp< VolField< Type > > ddt(const dimensioned< Type > dt, const fvMesh &mesh)
Definition: fvcDdt.C:45
Foam::fvc::grad
tmp< VolField< typename outerProduct< vector, Type >::type > > grad(const SurfaceField< Type > &ssf)
Definition: fvcGrad.C:46
Foam::fvc::div
tmp< VolField< Type > > div(const SurfaceField< Type > &ssf)
Definition: fvcDiv.C:47
Foam::I
static const Identity< scalar > I
Definition: Identity.H:93
Foam::fvVectorMatrix
fvMatrix< vector > fvVectorMatrix
Definition: fvMatricesFwd.H:45
Foam::inv
dimensionedSphericalTensor inv(const dimensionedSphericalTensor &dt)
Definition: dimensionedSphericalTensor.C:71
Foam::volSymmTensorField
VolField< symmTensor > volSymmTensorField
Definition: volFieldsFwd.H:64
Foam::magSqr
dimensioned< scalar > magSqr(const dimensioned< Type > &)
turbulence
autoPtr< incompressible::momentumTransportModel > turbulence(incompressible::momentumTransportModel::New(U, phi, viscosity))
p
volScalarField & p
Definition: createFieldRefs.H:5