[mostanad]

Hi foamers,
I am working on a case with three phases. Because of dealing with multiple phases, I found multiphaseInterFoam a proper option for my aim.
But in my simulation case, just for one existing phase, and not between two of them, I need to add a source term to that specific phase alpha equation.
We know alpha equation for interFoam is solved in alphaEqn.H. However, in multiphaseInterFoam, the calculation of this equation using MULES solver has been included in multiphaseMixture.C.
https://github.com/OpenFOAM/OpenFOAM...phaseMixture.C
When we look at its line 657, it seems that Su and Sp for source terms is set zeroFields in the explicit solver.

MULES::explicitSolve
        (
            geometricOneField(),
            alpha,
            alphaPhi,
            zeroField(),
            zeroField()
           );

and it is not possible to add a source term for multiphase solver, even for a single phase.
So can anyone guide my on this matter whether it is feasible to have source term in a multiple phase solver using multiphaseMixture?

[mostanad]

No comment from OF experts?!