好滴，谢谢啦~

@HydroJames 我暂时还查不到原因，测试过程中我还引入了其他差异，比如我进行过换硬盘的操作

forAll(walls.boundaryField(), patch) { walls.boundaryField()[patch] = L*hFlux[patch]/kappa; }

可以写如上代码计算，非常简单。上述代码没测试可能需要改动。

@李东岳

没，还是不知道那个公式咋来的，不过基本解决了flux limiter 和slope/gradient limiter的等价关系问题。

但是这玩意儿也是个大坑，对于多维网格带skew和non-uniform的情况，限制器设计和计算量都挺大。

OpenFOAM具体实现还没看太懂。有一个先有鸡还是先有蛋的问题，到底是interpolation先调用gradient呢？还是gradient调用interpolation？

我看的文章里是先leastSquare算gradient，不用Gauss的情况下是不用插值的。然后interpolation可以用一些limited gradient。

暂时只看懂了cellLimited, faceLimited没看懂。

@aufalpen 谢谢，公众号我早已经关注过了😄我去看看那个网页，多理解理解理论，还请多指教

边界条件不直接修改patchInternalField，而是通过改变离散后的系数矩阵起作用。

@李东岳 之前发过pdf，不过后来好像删了，你问问他。

单水滴破碎的倒是见过，多水滴碰撞的还没见过，楼主找到没有？

没用过gmsh...

网格体积：mesh.V()
网格面矢量 mesh.Sf()
网格面大小 mesh.magSf()

#问题1 自我解决，感觉像是bug.
0_1488141087886_streamline.png

@cfd-china 谢谢回复，我研究了源码已经搞懂了，顺便附上我总结的链接laplacian项离散解码

我这样写不对啊？ 不知道哪位可以给看看毛病在哪@administrators @李东岳

相对于计算发散，通常存在一个致命的因素可以解决。这种关于计算不准的东西，很难断定为什么，我经常计算不准。只能看sci，调节参数，调节网格，调节各种各种，毕竟人家算收敛了，收敛表明这个设置的结果，就是这样的。

Henry最近添加了部分描述：点击查看

确实，默认的rhoCentralDyMFoam没有使用correctPhi, 而且在其他可压流求解器中有使用correctPhi，不过需要注意的是，由于我在这里涉及到的是激波模拟，因此有压力间断，但是correctPhi似乎是通过压力进行通量修正，我不敢确定这是否会对压力间断产生影响，因此目前暂时还没有使用。我尝试过网格不运动，即直接使用rhoCentralFoam情况似乎是正确的。

只要原来的网格点都在就没问题

@100yearsalone
手头没有openfoam，但是在你使用dynamicMothion...代替solidBody....之后要定义solver，在你的字典里添加solver就可以了

solver xxx;

也有可能是：

solver { }

gedit打开文件很慢，可以换个编辑器，hoho

如果从0时间步开始计算就没问题？从非0时间步开始就有问题？是这个意思么