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里面求解相方程的问题