openFOAM动量方程中张量的处理

动量方程有许多种形式。本文以rhoPimpleFoam为例子进行说明。
rhoPimpleFoam/UEqn.h内容如下：

文件内容

tmp<fvVectorMatrix> UEqn
(
    fvm::ddt(rho, U) 										(I)
       + fvm::div(phi, U)									(II)
       + MRF.DDt(rho, U)								(III)
       + turbulence -> divDevReff(U)		(IV)
       ==
       fvOptions(rho, U)									(V)
)

其中，(I)为非稳态项；(II)为对流项；(III)为MRF造成的附加修正；(IV)为应力张量；(V)用于处理fvOptions文件夹中的附加源项。

我们主要关注turbulence->divDevReff(U)函数。（如果是可压缩流体，那么这个函数就是turbulence->divDevRhoReff(U)）

turbulence->divDevReff(U)是什么

类的名称是turbulence，该类基于turbulenceModel类，并从imcompressibleTurbulenceModel类中衍生出来。turbulence类调用了函数divDevReff(U)。该函数内容如下：

tmp<fvVectorMatrix> laminar::divDevReff(volVectorField& U) const
{
    return
    (
      - fvm::laplacian(nuEff(), U)