UEqns.H
Go to the documentation of this file.
1Info<< "Constructing face momentum equations" << endl;
2
3MRF.correctBoundaryVelocity(U1);
4MRF.correctBoundaryVelocity(U2);
5MRF.correctBoundaryVelocity(U);
6Info<< "Constructing face momentum equations" << endl;
7
8fvVectorMatrix U1Eqn(U1, rho1.dimensions()*U1.dimensions()*dimVolume/dimTime);
9fvVectorMatrix U2Eqn(U2, rho2.dimensions()*U2.dimensions()*dimVolume/dimTime);
10
11{
12 volScalarField Vm(fluid.Vm());
13
14 fvVectorMatrix UgradU1
15 (
16 fvm::div(phi1, U1) - fvm::Sp(fvc::div(phi1), U1)
17 + MRF.DDt(U1)
18 );
19
20 fvVectorMatrix UgradU2
21 (
22 fvm::div(phi2, U2) - fvm::Sp(fvc::div(phi2), U2)
23 + MRF.DDt(U2)
24 );
25
26 {
27 U1Eqn =
28 (
29 fvm::div(alphaRhoPhi1, U1) - fvm::Sp(fvc::div(alphaRhoPhi1), U1)
30 + MRF.DDt(alpha1*rho1, U1)
31 + phase1.turbulence().divDevRhoReff(U1)
32 + Vm*(UgradU1 - (UgradU2 & U2))
34 );
35 U1Eqn.relax();
36 fvOptions.constrain(U1Eqn);
37 U1.correctBoundaryConditions();
38 fvOptions.correct(U1);
39 }
40
41 {
42 U2Eqn =
43 (
44 fvm::div(alphaRhoPhi2, U2) - fvm::Sp(fvc::div(alphaRhoPhi2), U2)
45 + MRF.DDt(alpha2*rho2, U2)
46 + phase2.turbulence().divDevRhoReff(U2)
47 + Vm*(UgradU2 - (UgradU1 & U1))
49 );
50 U2Eqn.relax();
51 fvOptions.constrain(U2Eqn);
52 U2.correctBoundaryConditions();
53 fvOptions.correct(U2);
54 }
55}
fv::options & fvOptions
IOMRFZoneList & MRF
const volScalarField & alpha1
surfaceScalarField & phi2
twoPhaseSystem & fluid
phaseModel & phase1
volScalarField & rho2
const volScalarField & alpha2
surfaceScalarField & phi1
phaseModel & phase2
volVectorField & U1
volScalarField & rho1
volVectorField & U2
U
Definition: pEqn.H:72
surfaceScalarField & alphaRhoPhi2
surfaceScalarField & alphaRhoPhi1
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:82
messageStream Info
Information stream (stdout output on master, null elsewhere)
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:372
fvMatrix< vector > fvVectorMatrix
Definition: fvMatricesFwd.H:46