2212/multiphase_2MPPICInterFoam_2UEqn_8H_source.html

 MRF.correctBoundaryVelocity(U);

 fvVectorMatrix UEqn
 (
     fvm::ddt(alphacRho, U)
   + MRF.DDt(alphacRho, U)
   - fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)
   + fvm::div(rhoPhi, U)
   + turbulence->divDevRhoReff(U)
   ==
     fvOptions(rho, U)
   + cloudSU
 );

 UEqn.relax();
 fvOptions.constrain(UEqn);

 volScalarField rAUc(1.0/UEqn.A());
 surfaceScalarField rAUcf(fvc::interpolate(rAUc));


 surfaceScalarField phicForces
 (
     (fvc::interpolate(rAUc*cloudVolSUSu) & mesh.Sf())
 );


 if (pimple.momentumPredictor())
 {
     solve
     (
         UEqn
      ==
         fvc::reconstruct
         (
             phicForces/rAUcf
           +
             (
                 fvc::interpolate
                 (
                     mixture.sigmaK()
                 )*fvc::snGrad(alpha1)
                 - ghf*fvc::snGrad(rho)
                 - fvc::snGrad(p_rgh)
             ) * mesh.magSf()
         )
     );

     fvOptions.correct(U);
 }
solve
solve(UEqn==-fvc::grad(p))

rho
rho
Definition: readInitialConditions.H:88

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.

Foam::fac::div
tmp< GeometricField< Type, faPatchField, areaMesh > > div(const GeometricField< Type, faePatchField, edgeMesh > &ssf)
Definition: facDiv.C:43

mixture
Info<< "Creating temperaturePhaseChangeTwoPhaseMixture\"<< endl;autoPtr< temperaturePhaseChangeTwoPhaseMixture > mixture
Definition: createFields.H:39

Foam::fac::ddt
tmp< GeometricField< Type, faPatchField, areaMesh > > ddt(const dimensioned< Type > dt, const faMesh &mesh)
Definition: facDdt.C:40

fvOptions
fv::options & fvOptions
Definition: setRegionFluidFields.H:21

ghf
const surfaceScalarField & ghf
Definition: setRegionFluidFields.H:16

Foam::volScalarField
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:84

phicForces
surfaceScalarField phicForces((fvc::interpolate(rAUc *cloudVolSUSu) &mesh.Sf()))

mesh
dynamicFvMesh & mesh
Definition: createDynamicFvMesh.H:6

MRF
IOMRFZoneList & MRF
Definition: setRegionFluidFields.H:20

rAUcf
surfaceScalarField rAUcf(fvc::interpolate(rAUc))

pimple
pimpleControl & pimple
Definition: setRegionFluidFields.H:54

turbulence
Info<< "Reading field U\"<< 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\"<< endl;autoPtr< compressible::turbulenceModel > turbulence(compressible::turbulenceModel::New(rho, U, phi, thermo))
Definition: createFields.H:94

rAUc
volScalarField rAUc(1.0/UEqn.A())

Foam::fvc::reconstruct
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> reconstruct(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcReconstruct.C:49

U
U
Definition: pEqn.H:72

UEqn
fvVectorMatrix & UEqn
Definition: UEqn.H:13

rhoPhi
rhoPhi
Definition: rhoEqn.H:10

Foam::fvVectorMatrix
fvMatrix< vector > fvVectorMatrix
Definition: fvMatricesFwd.H:40

Foam::surfaceScalarField
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
Definition: surfaceFieldsFwd.H:64

p_rgh
volScalarField & p_rgh
Definition: setRegionFluidFields.H:13

Foam::fvc::snGrad
tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > snGrad(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
Definition: fvcSnGrad.C:40

Sp
zeroField Sp
Definition: alphaSuSp.H:2

alpha1
const volScalarField & alpha1
Definition: setRegionFluidFields.H:6