MULES 是修正输运方程的正确选择吗？ [OpenFOAM]

Question

我想求解以下修改后的输运方程来计算灵敏度：

在 interFoam 中，我在求解 alphaEqn.H 之后立即实现了以下代码（连同初始字段 dAlpha 和 dU、dPhi 的计算和其他修改）：

{
    word dAlphaScheme("div(phi,dAlpha)");
    word dPhiScheme("div(dPhi,alpha)");

    #include "alphaSuSp.H"

    surfaceScalarField dAlphaPhi
    (
    fvc::flux(phi, dAlpha, dAlphaScheme)
    + fvc::flux(dPhi, alpha1, dPhiScheme)
    );

    MULES::explicitSolve
    (
    geometricOneField(),
    dAlpha,
    dAlphaPhi,
    Sp,
    Su
    );
}

有效，但我不确定 MULES::explicitSolve 是否是求解该方程的正确方法。经过几个步骤后，溶液会严重涂抹。你会推荐什么？

Answer 1

正如我在上面的评论中提到的，

cfd-online 可能是您最好的选择。看起来你正在尝试实现一个额外的标量传输方程，虽然当你说 'sensitivities'.

时不清楚你在计算什么物理

查看此处关于向 OpenFOAM 求解器添加额外的标量传输方程的示例：https://openfoamwiki.net/index.php/How_to_add_temperature_to_icoFoam 它对于 OpenFOAMv1.7 来说已经过时了，但是概念是 99.999% 相同的：

创建字段
向字典中添加适当的传输系数
向 fvSchemes 添加适当的方案
定义输运方程
对输运方程使用 .solve() 函数

MULES 是修正输运方程的正确选择吗？ [OpenFOAM]

Is MULES the right choice for modified transport equation? [OpenFOAM]

c++

openfoam

fluid-dynamics