YEqns.H

Go to the documentation of this file.

{
    autoPtr<phaseSystem::massTransferTable>
        massTransferPtr(fluid.massTransfer());
 
    phaseSystem::massTransferTable&
        massTransfer(massTransferPtr());
 
    forAll(fluid.multiComponentPhases(), multiComponentPhasei)
    {
        phaseModel& phase = fluid.multiComponentPhases()[multiComponentPhasei];
 
        UPtrList<volScalarField>& Y = phase.YActiveRef();
        const volScalarField& alpha = phase;
        const tmp<volScalarField> trho = phase.rho();
 
        const auto& rho = trho();
 
        forAll(Y, i)
        {
            fvScalarMatrix YiEqn
            (
                phase.YiEqn(Y[i])
             ==
               *massTransfer[Y[i].name()]
              + fvOptions(alpha, rho, Y[i])
            );
 
            YiEqn.relax();
            YiEqn.solve("Yi");
        }
    }
}