UEqns.H
Go to the documentation of this file.
1 Info<< "Constructing momentum equations" << endl;
2 
3 PtrList<fvVectorMatrix> UEqns(phases.size());
4 
5 {
6  autoPtr<phaseSystem::momentumTransferTable>
7  momentumTransferPtr(fluid.momentumTransfer());
8 
9  phaseSystem::momentumTransferTable&
11 
13  {
14  phaseModel& phase = phases[phasei];
15 
16  const volScalarField& alpha = phase;
17  const volScalarField& rho = phase.rho();
18  volVectorField& U = phase.U();
19 
20  UEqns.set
21  (
22  phasei,
23  new fvVectorMatrix
24  (
25  phase.UEqn()
26  ==
27  *momentumTransfer[phase.name()]
28  + fvOptions(alpha, rho, U)
29  )
30  );
31 
32  UEqns[phasei].relax();
33  fvOptions.constrain(UEqns[phasei]);
34  fvOptions.correct(U);
35  }
36 }
PtrList< fvVectorMatrix > UEqns(fluid.phases().size())
multiphaseSystem & fluid
Definition: createFields.H:10
forAll(phases, phasei)
Definition: UEqns.H:12
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())
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
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:52
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
phaseSystem::momentumTransferTable & momentumTransfer(momentumTransferPtr())
GeometricField< vector, fvPatchField, volMesh > volVectorField
Definition: volFieldsFwd.H:55
multiphaseSystem::phaseModelList & phases
Definition: createFields.H:11
U
Definition: pEqn.H:82
int phasei
Definition: UEqns.H:6