v3/a05002_source.html

 {
     fvScalarMatrix rhoEqn
     (
         fvm::ddt(rho)
       + fvm::div(phi, rho)
     );

     rhoEqn.solve();

     rhoPhi = rhoEqn.flux();

     Info<< "max-min rho: " << max(rho).value()
         << " " << min(rho).value() << endl;

     rho == max(rho, rhoMin);
 }
Foam::Info
messageStream Info

Foam::fvc::ddt
tmp< GeometricField< Type, fvPatchField, volMesh > > ddt(const dimensioned< Type > dt, const fvMesh &mesh)
Definition: fvcDdt.C:45

rhoPhi
rhoPhi
Definition: rhoEqn.H:10

rho
rho
Definition: readInitialConditions.H:96

Foam::fvScalarMatrix
fvMatrix< scalar > fvScalarMatrix
Definition: fvMatricesFwd.H:42

Foam::endl
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:251

Foam::max
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)

phi
surfaceScalarField & phi
Definition: setRegionFluidFields.H:8

rhoMin
PtrList< dimensionedScalar > rhoMin(fluidRegions.size())

Foam::fvc::div
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcDiv.C:47

Foam::min
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)