turbulence-&gtdivDevReff(U)代码解读

目前，PIMPLE算法所采用的动量方程（不可压缩流）代码如下：

fvVectorMatrix UEqn
(
    fvm::ddt(U) + fvm::div(phi, U)
  + MRF.DDt(U)    //MRF是多参考坐标系，在动网格中会用到（像旋转机械等），一般情况下可忽略
  + turbulence->divDevReff(U)
 ==
    fvOptions(U)//源项或约束
);

本文旨在分析解决turbulence->divDevReff(U)这一行代码。

由Navier-Stokes方程，动量方程表达式如下：
$$\frac{\partial (\rho \mathbf{U})}{\partial t}=-\nabla \cdot (\rho \mathbf{U}\mathbf{U})-\nabla \cdot \tau -\nabla p+\rho \mathbf{f}$$

• ρ为密度，u为速度，p为压力，τ为剪切应力，f为体积力。

对于不可压缩流，有：
$$\frac{\partial \mathbf{U}}{\partial t}+\nabla \cdot (\mathbf{U}\mathbf{U})+\frac{1}{\rho }\nabla \cdot \tau =-\nabla \frac{p}{\rho }+\mathbf{f}$$

$$\frac{1}{\rho }\nabla \cdot \tau =\nabla \cdot \frac{\tau }{\rho }$$

令 $p=\frac{p}{\rho }$（运动压力），
$$\frac{\partial \mathbf{U}}{\partial t}+\nabla \cdot (\mathbf{U}\mathbf{U})+\nabla \cdot \frac{\tau }{\rho }=-\nabla p+\mathbf{f}$$

我们这里对剪切力项进行分析。将剪切应力用新符号表示：

$$\frac{\tau }{\rho }={\mathbf{R}}^{eff}$$

该项包含两部分：

1. 分子扩散项
2. 雷诺应力的偏分量的散度。 Boussinesq假设允许我们将雷诺应力加到剪切应力上。该应力张量（包含剪切应力和雷诺应力）表示为 $R^{eff}$，即有效应力张量。RAS和LES湍流模型都是基于Boussinesq假设的。雷诺应力主分量的散度归结到了压力项中，这是大多数雷诺时均RAS和大涡模型LES实现的一贯做法。因此压力比层流模型中多了一项,那就是雷诺应力的主分量，通常被称为湍动压力。

对于不可压缩流体，τ 的迹满足连续性方程：

由此，根据矩阵的偏置

代码实施：

注：

• 上述张量公式的运算规则可在东岳流体内查看
• fvOptions—OpenFOAM 中引入的一个新的框架，允许用户选择任何可以表示为控制方程的源或约束的物理量，例如多孔介质、磁流变液体和体积力等。具体内容可参考fvOptions代码分析。

参考网址：

1. http://www.jku. at/pfm/
2. http://blog.sina.com.cn/s/blog_5fdfa7e60100d8ah.html
3. https://openfoam.org/release/2-2-0/fv-options/