CFD中文网

@李东岳 至少论文里面是这样说的。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}$ 应该如何表示 ？

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

@溯光 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 中试验成功。