YEqns.H
Go to the documentation of this file.
1 {
2  autoPtr<phaseSystem::massTransferTable>
3  massTransferPtr(fluid.massTransfer());
4 
5  phaseSystem::massTransferTable&
6  massTransfer(massTransferPtr());
7 
9  {
10  phaseModel& phase = phases[phasei];
11 
12  PtrList<volScalarField>& Y = phase.Y();
13  const volScalarField& alpha = phase;
14  const volScalarField& rho = phase.rho();
15 
16  forAll(Y, i)
17  {
18  tmp<fvScalarMatrix> YiEqn(phase.YiEqn(Y[i]));
19 
20  if (YiEqn.valid())
21  {
22  YiEqn =
23  (
24  YiEqn
25  ==
26  *massTransfer[Y[i].name()]
27  + fvOptions(alpha, rho, Y[i])
28  );
29 
30  YiEqn->relax();
31  YiEqn->solve(mesh.solver("Yi"));
32  }
33  }
34  }
35 
36  fluid.massTransfer(); // updates interfacial mass flow rates
37 }
multiphaseSystem & fluid
Definition: createFields.H:10
label phasei
Definition: pEqn.H:37
dynamicFvMesh & mesh
fv::IOoptionList & fvOptions
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:52
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
PtrList< volScalarField > & Y
Definition: createFields.H:36
phaseSystem::massTransferTable & massTransfer(massTransferPtr())
forAll(phases, phasei)
Definition: YEqns.H:8
multiphaseSystem::phaseModelList & phases
Definition: createFields.H:11