1 Info<<
"Constructing face momentum equations" <<
endl;
30 +
phase1.turbulence().divDevRhoReff(
U1)
31 + Vm*(UgradU1 - (UgradU2 &
U2))
36 U1.correctBoundaryConditions();
46 +
phase2.turbulence().divDevRhoReff(
U2)
47 + Vm*(UgradU2 - (UgradU1 &
U1))
52 U2.correctBoundaryConditions();
surfaceScalarField & phi2
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Ostream & endl(Ostream &os)
Add newline and flush stream.
const dimensionSet dimVol(dimVolume)
tmp< GeometricField< Type, fvPatchField, volMesh > > Sp(const volScalarField &sp, const GeometricField< Type, fvPatchField, volMesh > &vf)
surfaceScalarField & phi1
surfaceScalarField & alphaRhoPhi1
dimensionedScalar posPart(const dimensionedScalar &ds)
GeometricField< scalar, fvPatchField, volMesh > volScalarField
surfaceScalarField & alphaRhoPhi2
fvMatrix< vector > fvVectorMatrix
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
dimensionedScalar negPart(const dimensionedScalar &ds)