关于两相流的数值计算

hilllike

本人自己编了一个三维两相流的数值计算程序。
poison方程求解总不能收敛，有没有高手给看下。。

1  投影法求解的ns方程；
2  vof法追踪液面；验证过vof正常。
3  模型是基于two phase model ，但是在计算的时候，我没加表面张力，（老师说不用加，不影响）

我一直怀疑是缺少表面张力相才导致压力方程无法收敛，但是缺少证据来反驳老师。有没有高手给指点下，两相流的资料太杂，一直没有看到系统的书籍 。

zhmeww

1 表面张力影响的确很小，可以不考虑；
2 采用投影法求解不可压问题是不稳定的，特别是非定常问题，建议采用simple或piso算法，若采用投影法建议采用petsc工具库，还需将压力收敛值设为10E-7左右。
3 建议先不考虑两相流，测试投影法。
水平有限，仅供参考

hilllike

谢谢回帖！

1，从前我是用的piso做的，老师说他们都在用mac方法，所以我最好改成投影法，因为他们懂。
2， 我的模型是一体模型，计算时候空气部分是和液体部分一样计算到的，密度用vof插值得到的。没有自由表面边界条件能行吗？

zouchu

找个VOF2D,或VOF3D代码来看看，
现在网上MAC+VOF和SIMPLE+VOF都可以找到代码