[OPFO]

For multiphase systems with 3 phases (gas, liquid, solid), how in the world can I access the phase fraction for each phase?

For interfoam, it's simple - alpha1 and alpha2. What is it for multiphase systems? alpha1, alpha2, alpha3 do not work.

I'm looking through the phaseModel and multiphaseSystem codes, but it's not as clear with my limited CPP knowledge.

[OPFO]

I haven't found anything on the forums and online so perhaps no one knows?

How can I explicitly reference the alpha.solid volume fraction?

From MulitphaseSystem.C: this prints out the various volume fractions for each phase.

for (phaseModel& phase : phases_)
{
volScalarField& alpha = phase;
alpha += alpha*sumCorr;

Info<< alpha.name() << " volume fraction = "
<< alpha.weightedAverage(mesh_.V()).value()
<< " Min(alpha) = " << min(alpha).value()
<< " Max(alpha) = " << max(alpha).value()
<< endl;
}

again, everything is always referenced in terms of the list of phases or phase pairs. Can anyone shed light on how I can determine the volume fraction of the solid phase only?

[geth03]

in multiphase euler foam you can access the alpha of the solid by

phases[index_Solid]. phases is the void fraction and because you have

multiple phases you need to know your index of the solid phase.

remember that in c++ you begin counting by zero. so if your solid is

in 3rd position you type phases[2].