MULES 是修正输运方程的正确选择吗？ [开放泡沫]答案

【问题标题】：Is MULES the right choice for modified transport equation? [OpenFOAM]MULES 是修正输运方程的正确选择吗？ [开放泡沫]
【发布时间】：2019-06-27 08:23:54
【问题描述】：

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

在 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 是否是解决此等式的正确方法。经过几个时间步骤后，溶液会严重涂抹。你会推荐什么？

【问题讨论】：

这是 cfd-online.com 的 OpenFOAM 编程论坛的一个问题，因为它不是关于 C++ 编程，而是更多关于使用 OpenFOAM。
@EricBringley 我在那里尝试过，但直到现在我没有得到任何回复，也许你知道我可以问的类似论坛？

标签： c++ openfoam fluid-dynamics

【解决方案1】：

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

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

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

【讨论】：