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