InterFoam相方程代码疑惑

[小考拉]

各位朋友大家好！

本人最近在看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里面求解相方程的问题

[东岳]

去年我完成了MULES-FCT的相关研究，提交给Computers Fluids，审了半年只找到一个审稿人，被拒稿，后来就没有再投，也不想投了，我计划6月份之前把这部分内容直接挂在网上，MULES-FCT在网上现存一些讨论，你可以自己找找，包括一些sci，但我总觉得写的不详细不通透。这部分内容在论坛里面基本讨论不清楚，不是3 5个方程就能解释清楚的。如果可以等我后续会把详细推导发布。

这部分内容非常复杂，目前我建议你直接将其理解为对$\alpha$的显性离散后，求解即可。

[小考拉]

好的，谢谢东岳老师指点！