[balkrishna]

I wish to modify bubble Foam of 2 phases to n phase and add heat transfer and species mixing to it .. One stumbling block i came across was the way the alpha equation was programmed .
it consists of the following line in the header file alphaEqn.H

 word scheme("div(phi,alpha)");

    surfaceScalarField phir = phia - phib;

This means they perform the phase calculation as stated in the wiki here . How does one extend the relative flux to n phase ...