好滴，谢谢啦~

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

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

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

@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时间步开始就有问题？是这个意思么

@zxyf1 是的，我的模型是个圆柱体，只需要让高度方向尺寸越来越小就可以，我已经实现了网格移动，就是把它用于interFoam时算的alpha场不对，目前 还在找原因，有兴趣可以一块探讨一下。