@李东岳 至少论文里面是这样说的。The strain rate tensor S is treated implicitly in the modified flow solver.

请教大家一个问题。
如果将雷诺应力按如下方式分解
$$
\boldsymbol{\tau}^{m}=\nu_{t}^{L} \mathbf{S}+\left(\boldsymbol{\tau}-\boldsymbol{\tau}^{L}\right)+\operatorname{tr}(\boldsymbol{\tau})
$$
其中$\nu_{t}^{L}$可以理解为湍流涡黏系数（但有所不同）。$\mathbf{S}=\frac{1}{2}\left[\nabla \mathbf{U}+(\mathbf{\nabla} \mathbf{U})^{T}\right]$是应变率张量（strain-rate tensor）。
如果要求$\mathbf{S}$隐式离散，在OpenFOAM中应该如何实现？（The strain rate tensor $\mathbf{S}$ is treated implicitly in the modified flow solver.）

我个人的理解：
以不可压缩流动为例，OpenFOAM-7中对应的雷诺应力计算方式为：

template<class BasicTurbulenceModel>
Foam::tmp<Foam::fvVectorMatrix>
Foam::linearViscousStress<BasicTurbulenceModel>::divDevRhoReff
(
    volVectorField& U
) const
{
    return
    (
      - fvc::div((this->alpha_*this->rho_*this->nuEff())*dev2(T(fvc::grad(U))))
      - fvm::laplacian(this->alpha_*this->rho_*this->nuEff(), U)
    );
}

因此，只有$\nabla \mathbf{U}$是隐式离散，而$(\mathbf{\nabla} \mathbf{U})^{T}$其实是显式离散的。并且也没有直接使用$\mathbf{S}$，相当于是把 $2\mathbf{S}$ 做了分解。

如果要对$\mathbf{S}$进行隐式离散，是不是应使用fvm::div() ?
$\nu_{t}^{L} \mathbf{S}$ 应该如何表示 ？

好久没登录，想请教个问题，突然提示没有C币，只好来水一下了
水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水

@东岳 李老师，不是同一个人。

@溯光 icoDyMIbFoam

参考资料：
【1】 http://www.cfd-china.com/topic/231/openfoam并行计算如何手动划分区域

【2】 https://openfoam.com/releases/openfoam-v1712/parallel.php

问题来源：
foam-extend中使用IBM求解器时，并行总是出现问题，似乎需要对部分区域进行人工分区。

已解决：
参考【1】，使 setFields，对边界存在的区域分配指定的核。

存在的问题：

首先使用 decomposePar 对整个计算域进行分区，然后再使用 setFields 对特定区域进行分区，导致不同的核处理的网格量有很大差别，降低了计算效率。

待解决的问题：

能否对特定的区域指定分区方法，然后运行 decomposePar 同时分区，而不是必须先对整个计算域进行分区。如附图所示，对区域1使用人工分区，对区域2使用 scotch 分区，以保证每个核处理的网格量相近。我找到了一些参考资料【2】，但没有在 foam-extend4.0 中试验成功。