@Aeronastro 我觉得是这样的：
where |S|= sqrt(2SijSij)，这句话可能来自LES，这只是一个定义，它不是正常意义上的magnitude.
正常的magnitude就是|S|= sqrt(SijSij)

并行计算时，利用 gMax 获取了某个 volScalarField 的全局最大值，然后我想找到这个网格的 cell_ID，因为我想输出这个网格的其它的一些 volScalarField 值，我写的代码只能在串行下实现这个功能，请问如何在并行下实现这个想法呢？
我的串行代码如下：

scalar theMax = gMax(oneVolScalarField);
forAll (oneVolScalarField,cellI)
{
    if ( oneVolScalarField[cellI] >= theMax )  
    {
        Info<<" When oneVolScalarField peaks, anotherVolScalarField ======"<<anotherVolScalarField[cellI]<<endl;  
    } 
}

@cfd-china 我是直接把图片拖进来

@wwzhao
谢谢您更新了我的认知！
现在我消化一下：
假如我们要求解方程 $\frac{\ \rho k}{\p t} = -30 + 50$。
左边这项当然变成了：fvm::ddt(rho, k)
右边第一项可以变成：- fvm::Sp(30/k, k)，移动这一项到左边，即把负号去掉了，然后把 Sp 的定义代入，即变成了 fvm.diag += V*30/k。这项加到对角元素上去了，和我们要求解的方程是一致的，并且使对角更占优了。
右边第二项可以变成：fvm::Su(50, k)，移动这一项到左边，并且把 Su 的定义代入，即变成了 fvm.source += V*50，与我们要求解的方程是一致的。

如果右边第一项写成这样：- fvm::Su(30, k)，那么它移动到左边后，将会使 fvm.source -= V*30，与我们要求解的方程是一致的。
如果右边第二项写成这样：fvm::Sp(50/k, k)，那么它移到左边后，将会使 fvm.diag -= V*50/k。这不利于对角占优，因此不推荐这样做。
结论：对于方程右边的项：
负数，写成 Sp——推荐这样做，因为可以使对角更占优。
负数，写成 Su——可以这样做。
正数，写成 Sp——不推荐这样做，因为会损害对角占优。
正数，写成 Su——应该这样做。
还有一种没讨论的情况：正负未知就写成 SuSp。
@wwzhao 请问我这里的理解正确吗？

在Chalmers大学的Lars Davidson教授的讲义《Fluid mechanics, turbulent flow and turbulence modeling》
里边提到，RANS 代码改成 LES，需要做的就是，1.加入SGS模型，2.修改对流项离散格式，3.修改时间离散格式。
修改这两个格式是因为不能引入额外的 numerical dissipation。

他还给了一个例子，把对流项离散格式改成中心差分，时间离散改成Crank-Nicolson 格式。

另外， Davidson提到，SIMPLE算法也是允许的：

@alex 有关注某个帖子的按钮啊，并且还可以关注某条回答，加入书签即可。这些都可以在你的个人主页查看。

@浪迹天大 在 最快速的上google的方式是什么？ 中说：

购买国外服务器，推荐digitalOcean啊，5美元1个月。另外github学生优惠可以获得50美元。首次使用必须充值5美元，加上github的50美元，就可以使用一年了。

这是我昨晚折腾一个加速软件之后的网速测试结果。

@shiyu 一个有点冗余的办法，定义一个vector场，每个网格都等于你那个gradP。这样就可以在别处lookup到了。

@星星星星晴 在 DecomposePar 的问题 中说：

uncollated

既然你提到了 uncollated，那么猜测你的算例应该是用了 collated 模式，但是有些文件却使用的 uncollated 的方式读取。
如果是这个问题，解决办法可以参考：
https://github.com/ZhangYanTJU/boundaryConditions/blob/master/LEMOSinflowGeneratorMod/decayingTurbulenceFvPatchVectorField.C#L140

这也有家公司：
https://www.nextsilicon.com/
https://github.com/nextsilicon/nextOpenFOAM

从 ESI 把这个功能偷过来了，现在在 OF-7下能用了。但是限制比较大，需要完全均一的网格。
https://github.com/ZhangYanTJU/functionObjects
这是我做的结果：

探索过一点，总结在这里了：field
应该是和本帖问题相关。

https://zhuanlan.zhihu.com/p/136258767

@Micro 我有个朋友写过相关的博文，感兴趣可以看看
https://zmeng.tech/2020/10/DynamicLibraries/

@马乔 都只是对官方的进行同级别的扩充。确实没有自己额外创造过结构啊。

@Micro 哈哈，到现在也没有写过自己的RTS

@freeman writeFormat改为binary呢？