UEqns.H
Go to the documentation of this file.
1 Info<< "Constructing face momentum equations" << endl;
2 
3 fvVectorMatrix U1Eqn(U1, rho1.dimensions()*U1.dimensions()*dimVol/dimTime);
4 fvVectorMatrix U2Eqn(U2, rho2.dimensions()*U2.dimensions()*dimVol/dimTime);
5 
6 {
7  volScalarField Vm(fluid.Vm());
8 
9  fvVectorMatrix UgradU1
10  (
12  + MRF.DDt(U1)
13  );
14 
15  fvVectorMatrix UgradU2
16  (
18  + MRF.DDt(U2)
19  );
20 
21  const volScalarField dmdt12(posPart(fluid.dmdt()));
22  const volScalarField dmdt21(negPart(fluid.dmdt()));
23 
24  {
25  U1Eqn =
26  (
28  + fvm::Sp(dmdt12, U1) + dmdt21*U2
29  + MRF.DDt(alpha1*rho1, U1)
30  + phase1.turbulence().divDevRhoReff(U1)
31  + Vm*(UgradU1 - (UgradU2 & U2))
32  );
33  U1Eqn.relax();
34  fvOptions.constrain(U1Eqn);
35  U1.correctBoundaryConditions();
36  fvOptions.correct(U1);
37  }
38 
39  {
40  U2Eqn =
41  (
43  - fvm::Sp(dmdt21, U2) - dmdt12*U1
44  + MRF.DDt(alpha2*rho2, U2)
45  + phase2.turbulence().divDevRhoReff(U2)
46  + Vm*(UgradU2 - (UgradU1 & U1))
47  );
48  U2Eqn.relax();
49  fvOptions.constrain(U2Eqn);
50  U2.correctBoundaryConditions();
51  fvOptions.correct(U2);
52  }
53 }
phaseModel & phase1
Definition: createFields.H:12
surfaceScalarField & alphaRhoPhi2
Definition: createFields.H:26
phaseModel & phase2
Definition: createFields.H:13
volVectorField & U1
Definition: createFields.H:18
multiphaseSystem & fluid
Definition: createFields.H:10
volVectorField & U2
Definition: createFields.H:23
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
Definition: dimensionSets.H:51
dimensionedScalar posPart(const dimensionedScalar &ds)
const dimensionSet dimVol(dimVolume)
Definition: dimensionSets.H:59
tmp< GeometricField< Type, fvPatchField, volMesh > > Sp(const volScalarField &sp, const GeometricField< Type, fvPatchField, volMesh > &vf)
Definition: fvcSup.C:67
messageStream Info
surfaceScalarField & phi2
Definition: createFields.H:24
fvMatrix< vector > fvVectorMatrix
Definition: fvMatricesFwd.H:45
fv::IOoptionList & fvOptions
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:251
IOMRFZoneList & MRF
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:52
dimensionedScalar negPart(const dimensionedScalar &ds)
volScalarField & alpha1
Definition: createFields.H:15
surfaceScalarField & phi1
Definition: createFields.H:19
surfaceScalarField & alphaRhoPhi1
Definition: createFields.H:21
alpha2
Definition: alphaEqn.H:112
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcDiv.C:47