用于求解不可压的pseudo-compressible method, 被认为有更好的收敛特性，见A. J. Chorin, JCP, 1967, pp. 12-26

http://www.cfd-china.com/topic/2744

我重新推出来一个带偏导数的压力方程..

ef44221d-062f-4315-b9c1-50e96688d88e-image.png

一个公式真是折腾死了。

我看了好几个人讲关于von Neumann稳定性分析。

12步入门Python这个女老师讲的。 MIT qiqi wang讲的PDE课程

还是有些地方很模糊。

@G8S7

究竟能不能用欧拉-欧拉模型来说明这个体系下固体颗粒的运动情况，

如果有实验图片的话更好判断。取决于固体颗粒粒径，相分数，粘度等。如果用你所说的VOF+欧拉-欧拉才能描述这一个过程，也是可以的。另一种方式是三相欧拉，在皮层熔体中对界面进行尖锐挤压处理。

另外，用什么模型取决你你关注的重点，如果关注内部的运动，皮层溶体不是重点用欧拉无妨，VOF也是有耗散。

At last, 几何很简单，是验证算法的利器，不用OpenFOAM可惜了。

@余正东 我都没听过 。。。

好的，谢谢

@肖艺 只是看过相关的影片.. :zoule:

我没说清楚。P23和P25的页码写在了下面。而前面所有的奇数页都在上面。没统一。

好像不能发动图，我试试截几个图吧（前两张是计算“正常”时，后面两张是出现了震荡和最后崩溃时的情况）
diver1.jpg diver2.jpg diver3.jpg diver4.jpg

@wwzhao 谢谢两位老师，文献DOI: 10.1115/1.1471361中有提到离散求解的

@netbuger :huahua: :huahua:

咦，这个工作我应该是2015年做的，怎么2017年才发帖..

整个域中的摩擦速度用的是当地的还是壁面的呢？

壁面的，我刚刚验证过，确定

替代文字

@wwzhao 再次感谢:xinxin:

@金石为开 Any updates?

看起来是的，找了下资料没看出来理由，看看大家什么看法?

倒是发现一篇0阶精度的文章..一会发个CFD界

@东岳 我也在做测吹方面的研究，但是E-E有点算不起（因为我的吹入速度很大)东岳先生的结晶器图片我看过，很漂亮。

是的，我觉得，本质上，还是应该走看书 - 手算 - 编程这个过程，如果看书直接去编程，可能理不清思路。

@东岳 过两天我测试一下6142

又研究了一下明白了：如果手动添加各项异性的全隐式源项，
需要拷贝一份A()并且构造成为张量，然后添加隐式源项，之后再构造HbyA，这样就能将各向异性的源线用全隐的方式添加到面通量上面了。