multiphase alpha.solid phase fraction
 

Hello,

I am using OpenFOAM-v1912 and the icoReactingMultiphaseInterFoam solver.

I am trying to convert a solver based on InterFoam where the term fvc::snGrad(alpha1) is used. However, for my icoReactingMultiphaseInterFoam solver, alpha1 won't work since I am using 3 phases and I believe the solver identifies phases based on pairs (e.g. ith phase).

Does anyone know how I can get the value of the solid phase fraction in volScalarField form for a multiphase solver?

Thanks in advance for any help!

Take a look into the source code in the surfaceTensionForce function and you will see a good example of referring to the phases. There, you will see that you have to loop over phases and via the iterators you will have access to alpha for each phase and also all the properties of that phase (rho, Cp, ...)