UEqns.H
Go to the documentation of this file.
1 Info<< "Constructing 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  autoPtr<phaseSystem::momentumTransferTable>
8  momentumTransferPtr(fluid.momentumTransfer());
9 
10  phaseSystem::momentumTransferTable&
12 
13  {
14  U1Eqn =
15  (
16  phase1.UEqn()
17  ==
18  *momentumTransfer[phase1.name()]
19  + fvOptions(alpha1, rho1, U1)
20  );
21  U1Eqn.relax();
22  fvOptions.constrain(U1Eqn);
23  fvOptions.correct(U1);
24  }
25 
26  {
27  U2Eqn =
28  (
29  phase2.UEqn()
30  ==
31  *momentumTransfer[phase2.name()]
32  + fvOptions(alpha2, rho2, U2)
33  );
34  U2Eqn.relax();
35  fvOptions.constrain(U2Eqn);
36  fvOptions.correct(U2);
37  }
38 }
phaseModel & phase1
Definition: createFields.H:12
phaseModel & phase2
Definition: createFields.H:13
volVectorField & U1
Definition: createFields.H:18
multiphaseSystem & fluid
Definition: createFields.H:10
Info<< "Constructing momentum equations"<< endl;fvVectorMatrix U1Eqn(U1, rho1.dimensions()*U1.dimensions()*dimVol/dimTime);fvVectorMatrix U2Eqn(U2, rho2.dimensions()*U2.dimensions()*dimVol/dimTime);{autoPtr< phaseSystem::momentumTransferTable > momentumTransferPtr(fluid.momentumTransfer())
volVectorField & U2
Definition: createFields.H:23
const dimensionSet dimTime(0, 0, 1, 0, 0, 0, 0)
Definition: dimensionSets.H:51
const dimensionSet dimVol(dimVolume)
Definition: dimensionSets.H:59
messageStream Info
fvMatrix< vector > fvVectorMatrix
Definition: fvMatricesFwd.H:45
fv::IOoptionList & fvOptions
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:251
rho1
Definition: pEqn.H:124
rho2
Definition: pEqn.H:125
volScalarField & alpha1
Definition: createFields.H:15
phaseSystem::momentumTransferTable & momentumTransfer(momentumTransferPtr())
alpha2
Definition: alphaEqn.H:112