All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
YEqns.H
Go to the documentation of this file.
1 {
2  autoPtr<phaseSystem::specieTransferTable>
3  specieTransferPtr(fluid.specieTransfer());
4 
5  phaseSystem::specieTransferTable&
6  specieTransfer(specieTransferPtr());
7 
8  fluid.correctReactions();
9 
10  forAll(fluid.multiComponentPhases(), multiComponentPhasei)
11  {
12  phaseModel& phase = fluid.multiComponentPhases()[multiComponentPhasei];
13 
14  UPtrList<volScalarField>& Y = phase.YActiveRef();
15  const volScalarField& alpha = phase;
16  const volScalarField& rho = phase.rho();
17 
18  forAll(Y, i)
19  {
20  fvScalarMatrix YiEqn
21  (
22  phase.YiEqn(Y[i])
23  ==
24  *specieTransfer[Y[i].name()]
25  + fvModels.source(alpha, rho, Y[i])
26  );
27 
28  YiEqn.relax();
29  fvConstraints.constrain(YiEqn);
30  YiEqn.solve("Yi");
32  }
33  }
34 
35  fluid.correctSpecies();
36 }
fvMatrix< scalar > fvScalarMatrix
Definition: fvMatricesFwd.H:42
volScalarField alpha(IOobject("alpha", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE), lambda *max(Ua &U, zeroSensitivity))
forAll(fluid.multiComponentPhases(), multiComponentPhasei)
Definition: YEqns.H:10
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:58
Foam::fvConstraints & fvConstraints
phaseSystem & fluid
Definition: createFields.H:11
phaseSystem::specieTransferTable & specieTransfer(specieTransferPtr())
bool constrain(fvMatrix< Type > &eqn) const
Apply constraints to an equation.
Foam::fvModels & fvModels
PtrList< volScalarField > & Y