请问一下upwind<scalar>(mesh_, phi_).flux(alpha)是什么意思

[lrl3512]

请问一下这个代码是什么意思，alpha插值到面上，用迎风格式吗？我看求解相方程里面都有这个，加这一个标量有什么意义吗

[lrl3512]

还是就是Sf * Uf * alpha

[李东岳]

明确的指定一阶格式

[lrl3512]

@李东岳 东岳老师，那加这个迎风的通量目的是为了什么，在icoReactingMultiphaseFoam里面前面是一个alpha的通量flux::fvc(alpha,phi)的，后面又加了一个 ，原作者注释掉了upwind<scalar>(mesh_, phi).flux(alpha1)，实际上alpha1Eqn.flux的结果和upwind<scalar>(mesh_, phi).flux(alpha1)是一样的，同样在compressmultiphaseinterfoam里面也加上了这个通量，而在intermixingFoam里面是先减去在加上这个

[李东岳]

icoReactingMultiphaseFoam没详细看过。不过MULES需要使用一阶格式。需要一阶格式+高阶格式混合在一起。

[lrl3512]

@李东岳 东岳老师 ，我看了您关于MULES的解析，MULES的主要就是用两套通量去构建一个既保持高精度又不越界的通量，那具体要调用MULES::explicitSolve()要怎么输入参数呢，如果我要求解$\frac{\partial \alpha}{\partial t}+\nabla_{\bullet}\left(\alpha \vec{U}\right)+\nabla\cdot(\alpha_1\alpha_2c|\mathbf{U}|\frac{\nabla\alpha}{|\nabla\alpha|})=0$

template<class RhoType, class SpType, class SuType>
void Foam::MULES::explicitSolve
(
    const RhoType& rho,
    volScalarField& psi,
    const surfaceScalarField& phiPsi,
    const SpType& Sp,
    const SuType& Su
)

我在第三个参数phipsi要传入$\nabla_{\bullet}\left(\alpha \vec{U}\right)_{high}+\nabla\cdot(\alpha_1\alpha_2c|\mathbf{U}|\frac{\nabla\alpha}{|\nabla\alpha|})+\nabla_{\bullet}\left(\alpha \vec{U}\right)_{low}$这个总的通量吗

[lrl3512]

@lrl3512 自问自答一下，在MULES::explicitSolve()之前调用了

最后一个参数如果是true，

返回是通量的修正值，后面还是要加一个一阶迎风的通量，其实直接false也行，后面就不需要加了