小白对于openfoam k-epsilon湍流模型程序编写的困惑

lllwonderliquid

大家好，我最近在阅读openfoam4.1中的k-epsilon湍流模型程序。。在读到epsilon方程和k方程时有些小困惑，以epsilon方程为例，对于扩散项和耗散项，程序中将其分为三项来写，如下所示：

tmp<fvScalarMatrix> epsEqn
    (
        fvm::ddt(alpha, rho, epsilon_)
      + fvm::div(alphaRhoPhi, epsilon_)
      - fvm::laplacian(alpha*rho*DepsilonEff(), epsilon_)
     ==
        C1_*alpha()*rho()*G*epsilon_()/k_()
      - fvm::SuSp(((2.0/3.0)*C1_ + C3_)*alpha()*rho()*divU, epsilon_)
      - fvm::Sp(C2_*alpha()*rho()*epsilon_()/k_(), epsilon_)
      + epsilonSource()
      + fvOptions(alpha, rho, epsilon_)
    );

中的这三项：

 C1_*alpha()*rho()*G*epsilon_()/k_()
      - fvm::SuSp(((2.0/3.0)*C1_ + C3_)*alpha()*rho()*divU, epsilon_)
      - fvm::Sp(C2_*alpha()*rho()*epsilon_()/k_(), epsilon_)

为什么这三项需要使用不同的差分格式？为什么第一项不需要fvm::，第二项使用fvm::SuSp而第三项却使用fvm::Sp？

李东岳

fvm表示隐性离散，不带fvm的表示显性离散。具体显性还是隐性和你的方程有关。第一项如果隐性离散的话，会使你的矩阵对角系数降低，因此显性离散。第二项由于不确定正负号因此使用fvm::SuSp自动调节显性/隐性离散，第三项因为为负号可以提高矩阵的对角占优因此使用fvm::Sp。

非常好的问题。

lllwonderliquid

@李东岳那显性离散的不需要加fvc::吗:surprised: ，经常看到有些量加有些量却不加:crying:

李东岳

如果你的函数不需要操作，比如就是-U，那就不需要加。如果你的函数需要进行张量操作，如果显性离散，就需要加fvc

lllwonderliquid

@李东岳谢谢东岳前辈！:cheeky:

Calf.Z-DNS

我还想问下epsilonSource()跟 fvOptions(alpha, rho, epsilon_) 在代码里面的作用是什么呢？

strive_DUT

@Calf-Z-DNS 请问您了解这两个函数的意义了吗？

李东岳

@strive_DUT 分别是epsilon的源项以及fvOptions源项，如果没有源项可以不管，主要是方便普适性编程

strive_DUT

@东岳谢谢李老师！请问里面的代码含义是不是表示量纲的乘除？比如说dimVolume代表体积量纲，dimTime代表时间量纲，k_.dimensions()代表湍动能的量纲？

dimVolume*this->rho_.dimensions()*k_.dimensions()
            /dimTime

李东岳

是的

strive_DUT

@东岳谢李老师！

CFD中文网

小白对于openfoam k-epsilon湍流模型程序编写的困惑