本人最近在看InterFoam里面alphaEqn.H代码（4.1版本），有一段代码不是很明白，希望大家给一些指点

if (alphaApplyPrevCorr && MULESCorr)
    {
        talphaPhiCorr0 = alphaPhi - talphaPhiCorr0;
    }

    if
    (
        word(mesh.ddtScheme("ddt(rho,U)"))
     == fv::EulerDdtScheme<vector>::typeName
    )
    {
        rhoPhi = alphaPhi*(rho1 - rho2) + phiCN*rho2;
    }
    else
    {
        if (ocCoeff > 0)
        {
            // Calculate the end-of-time-step alpha flux
            alphaPhi = (alphaPhi - (1.0 - cnCoeff)*alphaPhi.oldTime())/cnCoeff;
        }

        // Calculate the end-of-time-step mass flux
        rhoPhi = alphaPhi*(rho1 - rho2) + phi*rho2;
    }

上面这段代码 是利用MULES算法求解出修正的总通量alphaPhi 之后的 一段代码，

现在不太明白这段代码的作用，，同时希望大家推荐 有关这些代码的 一些参考公式推导 或者公式推导的思路

恳请好心人解答！同时也随时欢迎交流一下有关InterFoam里面求解相方程的问题