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
8 forAll(fluid.multiComponentPhases(), multiComponentPhasei)
9 {
10 phaseModel& phase = fluid.multiComponentPhases()[multiComponentPhasei];
11
12 UPtrList<volScalarField>& Y = phase.YActiveRef();
13 const volScalarField& alpha = phase;
14 const volScalarField& rho = phase.rho();
15
16 forAll(Y, i)
17 {
18 fvScalarMatrix YiEqn
19 (
20 phase.YiEqn(Y[i])
21 ==
22 *massTransfer[Y[i].name()]
23 + fvOptions(alpha, rho, Y[i])
24 );
25
26 YiEqn.relax();
27 YiEqn.solve(mesh.solver("Yi"));
28 }
29 }
30}
fv::options & fvOptions
twoPhaseSystem & fluid
PtrList< volScalarField > & Y
dynamicFvMesh & mesh
phaseSystem::massTransferTable & massTransfer(massTransferPtr())
volScalarField & alpha
#define forAll(list, i)
Loop across all elements in list.
Definition: stdFoam.H:333