湍流模型方程和代码问题

[李东岳]

@是嘿嘿牛啊 在 湍流模型方程和代码问题 中说：

图二中左边的代码和右边的公式感觉不匹配，求老师们指点一下，十分感谢！

我大略看了下，好像公式写错了，我明天详细看一下

图一中湍流模型方程左右都乘了一个相分数，那如果在两相流模型中，都是空气时，相分数为0，并且相分数的时间导数和空间导数也是0，这方程怎么解？

这个是著名的数值问题。每个版本的openfoam处理不一样。这个问题非常有意思，我这几天更新在ns方程笔记里面。

这就是phase intensive方法要处理的问题。有些文章叫做奇异。

[李东岳]

你第二个问题是哪个版本的openfoam，我这面openfoam-9跟你的代码不一样

[是嘿嘿牛啊]

@李东岳 李老师您好，我目前用的版本是v1812

[李东岳]

@是嘿嘿牛啊

1. openfoam9跟openfoam1812代码不一样

2. 我又细看了一下，虽然写的不一样，但是公式没问题，公式是哪里有问题呢？

[是嘿嘿牛啊]

@李东岳 公式没问题，只是左边红框内代码fvm::Sp(alpha()*rho()*epsilonByk(F1, tgradU()), k_)不知道是什么意思，是公式里面的哪一项？刚开始学of编程，求李老师指点一下，十分感谢！

[李东岳]

@是嘿嘿牛啊

fvm::Sp(alpha()*rho()*epsilonByk(F1, tgradU()), k_)

可以理解为alpha()*rho()*epsilonByk(F1, tgradU())*k
epsilonByk(F1, tgradU())的定义，要参考源代码的定义：

kOmegaSST<MomentumTransportModel, BasicMomentumTransportModel>::epsilonByk
(
    const volScalarField::Internal& F1,
    const volScalarField::Internal& F2
) const
{
    return betaStar_*omega_();
}

alpha()*rho()*epsilonByk(F1, tgradU())*k 就变成了alpha()*rho()*betaStar_*omega_()*k 。这个跟公式就可以对上了

[是嘿嘿牛啊]

@李东岳 李老师，下面那一行代码alpha()*rho()*betaStar_*omegaInf_*kInf_和alpha()*rho()*betaStar_*omega_()*k有什么区别吗？
我看上面代码中有对omegaInf_和kInf_定义，

    kInf_
    (
        dimensioned<scalar>::lookupOrAddToDict
        (
            "kInf",
            this->coeffDict_,
            k_.dimensions(),
            0
        )
    ),
    omegaInf_
    (
        dimensioned<scalar>::lookupOrAddToDict
        (
            "omegaInf",
            this->coeffDict_,
            omega_.dimensions(),
            0
        )
    )

那是否可以理解成omegaInf_和kInf_就是0？不明白alpha()*rho()*betaStar_*omegaInf_*kInf_这串代码的意义，请教一下老师，谢谢！

[李东岳]

是0 ，但是提供了一个用户自定义的接口

[是嘿嘿牛啊]

@李东岳 明白了，谢谢李老师！

[李东岳]

图一中湍流模型方程左右都乘了一个相分数，那如果在两相流模型中，都是空气时，相分数为0，并且相分数的时间导数和空间导数也是0，这方程怎么解？

http://dyfluid.com/reactingTwoPhaseEulerFoam.html 更新到这里了

[是嘿嘿牛啊]

@李东岳 正在学习，谢谢李老师！