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$的显性离散后，求解即可。

小考拉

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

CFD中文网

InterFoam相方程代码疑惑