c++ - 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 是否是求解该方程的正确方法。经过几个步骤后，溶液会严重绘(repaint)抹。你会推荐什么？

Answer 1

正如我在上面的评论中提到的，cfd-online 可能是您最好的选择。看起来您正在尝试实现一个额外的标量传输方程，但不清楚您在说“灵敏度”时计算的是什么物理学。

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

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