@李东岳

没，还是不知道那个公式咋来的，不过基本解决了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时间步开始就有问题？是这个意思么

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

@金石为开
不太懂CFx，因此也无法对比。不过速度入口，压力出口，我觉得你应该会知道的吧？inlet: fixedValue, outlet: zeroGradient.

@队长别开枪  好的，谢谢你的回复，我再看看