0

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

https://latex.codecogs.com/gif.latex?%5Cfrac%7B%5Cpartial%7D%7B%5Cpartial%20t%7D%20%28%5Cdelta%20%5Calpha%29%20+%20%5Cnabla% 20%5Ccdot%20%28%5Cvec%7Bu%7D%20%5C%20%5Cdelta%20%5Calpha%29%20+%20%5Cnabla%20%5Ccdot%20%28%5Cdelta%20%5Cvec%7Bu %7D%20%5C%20%5Calpha%29

在 interFoam 中,我在求解后立即实现了以下代码alphaEqn.H(以及初始字段dAlphadU、计算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是解决这个等式的正确方法。经过几个时间步骤后,溶液会严重涂抹。你会推荐什么?

4

1 回答 1

1

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

看看这里关于向 OpenFOAM 求解器添加额外的标量传输方程的例子:https://openfoamwiki.net/index.php/How_to_add_temperature_to_icoFoam 它已经过时了,因为它对于 OpenFOAMv1.7,但这个概念是 99.999%相同的:

  • 创建字段
  • 在字典中添加适当的传输系数
  • 将适当的方案添加到 fvSchemes
  • 定义传输方程
  • 在传输方程上使用 .solve() 函数
于 2019-07-08T17:30:16.687 回答