[mostanad]

Hi My Friends,
I am using icoReactingMultiphaseInterFoam in Openfoam v 2212, and somthing looks no clear to me. This solver is able to handle multicomponent multiphase reacting flows. From what I know, the species equation within each phase is defined as:



However, in the Mules solver for species equation in MultiComponentPhaseModel.C, we use:

MULES::explicitSolve
                    (
                        geometricOneField(),
                        Yi,
                        phi,
                        phiYiCorr,
                        Sp[i],
                        Su[i],
                        oneField(),
                        zeroField()
                    );

You see there is nothing about the of the phase related to that species, as a coefficient. Can anyone help me to undertand this?
Regards,
Mohammad

[mostanad]

No one has an opinion about this?

[mostanad]

OK. I realized it would be great to share my answer to the question above with you. Maybe in the future, someone will see this answer and find it useful for better clarification.

I do believe that OpenFOAM in icoReactingMultiphaseInterFoam solves instead of Y. This can be seen in MultiComponentPhaseModel.C, particularly in the following part:

fvScalarMatrix YiEqn
(
fv::EulerDdtScheme<scalar>(mesh).fvmDdt(Yi)
+ fv::gaussConvectionScheme<scalar>
(
mesh,
phi,
upwind<scalar>(mesh, phi)
).fvmDiv(phi, Yi)
==
Su[i] + fvm::Sp(Sp[i], Yi)
);

This part literally tells you that the governing equation is something like:

whereas, based on the equation in the first post, it should be:

And I don't know why! Even in test cases like evaporationMultiComponent in the tutorials, you can see that Y actually tries to follow Alpha values. This is because the explicit solution in MultiComponentPhaseModel.C is very similar to Alpha if we don't include Alpha in our equation: