UEqn.H
Go to the documentation of this file.
1 MRF.correctBoundaryVelocity(U);
2 
3 fvVectorMatrix UEqn
4 (
5  fvm::ddt(alphacRho, U)
6  + MRF.DDt(alphacRho, U)
7  - fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)
8  + fvm::div(rhoPhi, U)
9  + turbulence->divDevRhoReff(U)
10  ==
11  fvOptions(rho, U)
12  + cloudSU
13 );
14 
15 UEqn.relax();
16 fvOptions.constrain(UEqn);
17 
18 volScalarField rAUc(1.0/UEqn.A());
19 surfaceScalarField rAUcf(fvc::interpolate(rAUc));
20 
21 
22 surfaceScalarField phicForces
23 (
24  (fvc::interpolate(rAUc*cloudVolSUSu) & mesh.Sf())
25 );
26 
27 
28 if (pimple.momentumPredictor())
29 {
30  solve
31  (
32  UEqn
33  ==
34  fvc::reconstruct
35  (
36  phicForces/rAUcf
37  +
38  (
39  fvc::interpolate
40  (
41  mixture.sigmaK()
42  )*fvc::snGrad(alpha1)
43  - ghf*fvc::snGrad(rho)
44  - fvc::snGrad(p_rgh)
45  ) * mesh.magSf()
46  )
47  );
48 
49  fvOptions.correct(U);
50 }
phicForces
surfaceScalarField phicForces((fvc::interpolate(rAUc *cloudVolSUSu) &mesh.Sf()))
Foam::fvc::reconstruct
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> reconstruct(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcReconstruct.C:56
Foam::fvc::snGrad
tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > snGrad(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
Definition: fvcSnGrad.C:47
rhoPhi
rhoPhi
Definition: rhoEqn.H:10
ghf
const surfaceScalarField & ghf
Definition: setRegionFluidFields.H:18
turbulence
Info<< "Reading field U\n"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);volScalarField rho(IOobject("rho", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), thermo.rho());volVectorField rhoU(IOobject("rhoU", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *U);volScalarField rhoE(IOobject("rhoE", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *(e+0.5 *magSqr(U)));surfaceScalarField pos(IOobject("pos", runTime.timeName(), mesh), mesh, dimensionedScalar("pos", dimless, 1.0));surfaceScalarField neg(IOobject("neg", runTime.timeName(), mesh), mesh, dimensionedScalar("neg", dimless, -1.0));surfaceScalarField phi("phi", fvc::flux(rhoU));Info<< "Creating turbulence model\n"<< endl;autoPtr< compressible::turbulenceModel > turbulence(compressible::turbulenceModel::New(rho, U, phi, thermo))
Definition: createFields.H:94
Sp
zeroField Sp
Definition: alphaSuSp.H:2
alpha1
const volScalarField & alpha1
Definition: setRegionFluidFields.H:8
Foam::fac::div
tmp< GeometricField< Type, faPatchField, areaMesh > > div(const GeometricField< Type, faePatchField, edgeMesh > &ssf)
Definition: facDiv.C:50
rho
rho
Definition: readInitialConditions.H:88
MRF
IOMRFZoneList & MRF
Definition: setRegionFluidFields.H:22
pimple
pimpleControl & pimple
Definition: setRegionFluidFields.H:56
Foam::fvVectorMatrix
fvMatrix< vector > fvVectorMatrix
Definition: fvMatricesFwd.H:47
fvOptions
fv::options & fvOptions
Definition: setRegionFluidFields.H:23
p_rgh
volScalarField & p_rgh
Definition: setRegionFluidFields.H:15
Foam::volScalarField
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:57
mesh
dynamicFvMesh & mesh
Definition: createDynamicFvMesh.H:6
U
U
Definition: pEqn.H:72
Foam::surfaceScalarField
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
Definition: surfaceFieldsFwd.H:54
Foam::fac::ddt
tmp< GeometricField< Type, faPatchField, areaMesh > > ddt(const dimensioned< Type > dt, const faMesh &mesh)
Definition: facDdt.C:47
solve
solve(UEqn==-fvc::grad(p))
UEqn
fvVectorMatrix & UEqn
Definition: UEqn.H:13
rAUc
volScalarField rAUc(1.0/UEqn.A())
rAUcf
surfaceScalarField rAUcf(fvc::interpolate(rAUc))
mixture
Info<< "Creating temperaturePhaseChangeTwoPhaseMixture\n"<< endl;autoPtr< temperaturePhaseChangeTwoPhaseMixture > mixture
Definition: createFields.H:39
Foam::fac::interpolate
static tmp< GeometricField< Type, faePatchField, edgeMesh > > interpolate(const GeometricField< Type, faPatchField, areaMesh > &tvf, const edgeScalarField &faceFlux, Istream &schemeData)
Interpolate field onto faces using scheme given by Istream.