1Info<<
"Constructing momentum equations" <<
endl;
3MRF.correctBoundaryVelocity(
U1);
4MRF.correctBoundaryVelocity(
U2);
5MRF.correctBoundaryVelocity(
U);
21 +
phase1.turbulence().divDevRhoReff(
U1)
27 - fvm::Sp(fvc::div(
phi1),
U1)
33 U1Eqn += fvm::Sp(
Kd,
U1);
35 U1.correctBoundaryConditions();
45 +
phase2.turbulence().divDevRhoReff(
U2)
51 - fvm::Sp(fvc::div(
phi2),
U2)
57 U2Eqn += fvm::Sp(
Kd,
U2);
59 U2.correctBoundaryConditions();
const volScalarField & alpha1
surfaceScalarField & phi2
const volScalarField & alpha2
surfaceScalarField & phi1
const volScalarField Kd(fluid.Kd())
surfaceScalarField & alphaRhoPhi2
surfaceScalarField & alphaRhoPhi1
GeometricField< scalar, fvPatchField, volMesh > volScalarField
messageStream Info
Information stream (stdout output on master, null elsewhere)
Ostream & endl(Ostream &os)
Add newline and flush stream.
fvMatrix< vector > fvVectorMatrix