@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美元，就可以使用一年了。

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

@Samuel-Tu 比如你有一个大括号，包含很多行的内容，用一下就把你变成一行了。

@xpqiu 以前用sed，现在也爱上foamDictionary了，但是它会把注释弄没了，还强行帮我格式化一下，这点不太喜欢。

如果你的背景色也是白色的话

选用一个下限是白色的colorbar，就相当于隐藏了啊

@浪迹天大 在 openfoam添加湍流入口方法? 中说：

我把代码移植到OF-上了，并附上了我的算例以供参考：
https://github.com/ZhangYanTJU/boundaryConditions

少打了一个字。。。OF-7

我把代码移植到OF-7上了，并附上了我的算例以供参考：
https://github.com/ZhangYanTJU/boundaryConditions

我把代码移植到OF-上了，并附上了我的算例以供参考：
https://github.com/ZhangYanTJU/boundaryConditions

最近帮一个同事解决了类似问题，因为他对一个默认边界条件（calculated）的volScalarField量进行了fvc::grad操作，一个解决办法是fvc::grad之前先对这个量执行一下correctBoundaryConditions()。希望可以给你一个参考。

@东岳
我看到的=号也变边界啊
https://github.com/OpenFOAM/OpenFOAM-dev/blob/master/src/OpenFOAM/fields/GeometricFields/GeometricField/GeometricField.C#L1387
可以再详细谈一谈可能存在的坑吗？

苏老师培训讲义里有这个流程图，如果你能搞到 话。