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



  • 好的,谢谢东岳老师指点!


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2