31 void Foam::solvers::multiphaseEuler::cellMomentumPredictor()
33 Info<<
"Constructing momentum equations" <<
endl;
37 autoPtr<phaseSystem::momentumTransferTable>
41 momentumTransfer(momentumTransferPtr());
58 *momentumTransfer[phase.name()]
63 UEqns[phase.index()].relax();
65 U.correctBoundaryConditions();
71 void Foam::solvers::multiphaseEuler::faceMomentumPredictor()
73 Info<<
"Constructing face momentum equations" <<
endl;
75 phaseSystem& fluid(fluid_);
77 autoPtr<phaseSystem::momentumTransferTable>
78 momentumTransferPtr(fluid.momentumTransferf());
81 momentumTransfer(momentumTransferPtr());
83 forAll(movingPhases, movingPhasei)
85 phaseModel& phase = movingPhases_[movingPhasei];
98 *momentumTransfer[phase.name()]
103 UEqns[phase.index()].relax();
105 U.correctBoundaryConditions();
117 UEqns.setSize(phases.size());
121 faceMomentumPredictor();
125 cellMomentumPredictor();
#define forAll(list, i)
Loop across all elements in list.
bool constrain(fvMatrix< Type > &eqn) const
Apply constraints to an equation.
virtual autoPtr< momentumTransferTable > momentumTransfer()=0
Return the momentum transfer matrices for the cell-based.
HashPtrTable< fvVectorMatrix > momentumTransferTable
Foam::fvModels & fvModels() const
Return the fvModels that are created on demand.
Foam::fvConstraints & fvConstraints() const
Return the fvConstraints that are created on demand.
PtrList< fvVectorMatrix > UEqns
Temporary phase momentum matrices.
const phaseSystem::phaseModelPartialList & movingPhases
Reference to the moving phases.
virtual void momentumPredictor()
Construct and optionally solve the momentum equation.
phaseSystem::phaseModelPartialList & movingPhases_
const phaseSystem & fluid
Reference to the multiphase fluid.
Foam::fvConstraints & fvConstraints(Foam::fvConstraints::New(mesh))
Foam::fvModels & fvModels(Foam::fvModels::New(mesh))
pimpleControl pimple(mesh)
Calculate the matrix for implicit and explicit sources.
volScalarField alpha(IOobject("alpha", runTime.name(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), lambda *max(Ua &U, zeroSensitivity))
VolField< vector > volVectorField
Ostream & endl(Ostream &os)
Add newline and flush stream.
VolField< scalar > volScalarField
fvMatrix< vector > fvVectorMatrix