动量方程有许多种形式。本文以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)