CFD中文网

浪

@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 格式。

0_1467122753910_无标题.png
另外， Davidson提到，SIMPLE算法也是允许的：