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 }
U
Definition: pEqn.H:83
multiphaseSystem & fluid
Definition: createFields.H:11
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())
multiphaseSystem::phaseModelList & phases
Definition: createFields.H:12
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:253
GeometricField< vector, fvPatchField, volMesh > volVectorField
Definition: volFieldsFwd.H:55
fv::options & fvOptions
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:52
int phasei
Definition: UEqns.H:6
phaseSystem::momentumTransferTable & momentumTransfer(momentumTransferPtr())
forAll(phases, phasei)
Definition: UEqns.H:12
PtrList< fvVectorMatrix > UEqns(fluid.phases().size())
fvMatrix< vector > fvVectorMatrix
Definition: fvMatricesFwd.H:45
messageStream Info
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].