@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{\partial \rho k}{\partia 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美元，就可以使用一年了。

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

@starDust 如果最开始的时候，不输出为 raw 格式，而是 vtk 格式，就可以直接用 paraview 打开了。

@Henry-Zhou 我的意思很简单，就是需要用类似这样的操作：Y[0]，你需要找到你自己需要的那个下标啊。

这里的报错解读如下：
Y-Csat 里边的 Y 是一个 List场，用它减一个场，当然会出错。你需要对 Y 取下标，获取你需要的那个分量。

推荐一个可用的代码，支持 OF-4 和 OF-7
https://github.com/ZmengXu/sampledPlaneAverage

@hotairballon1997 我的意思是，找到能够获取 position的地方，然后通过函数参数，一层一层传到当前的函数。

@hotairballon1997 应该无法获取的，可以直接改接口，把 position 传进来。

@Zhy2022 大哥，可以麻烦找一下原始地址吗？这个网站是专门盗别人成果的。不要白白替他宣传了。

@同学博 是的，基于 freeCAD，专门给 OF 开发的前处理工具

@同学博 也许可以试试这个
https://github.com/jaheyns/CfdOF

renumberMesh 记得加overwrite啊