vof的高手给点建议

开心果

   vof的高手给点建议：最近接触了一些vof的文章和程序，好像很多都是有限差分法为基础的，而我的程序是fvm的simple程序，各位知不知道有没有基于fvm的vof阿？我就怀疑是不是有限差分也无所谓，就是用现成的程序，求出C的控制方程，然后求出自由面的位置，然后再接过自己的程序施加边界条件就好了，是不是太天真了？

开心果

刘儒勋的《数值模拟方法和运动边界面追踪》

mhyuan

偶也发现不少VOF的文章中流场求解是用MAC或project一类的有限差分法。不仅VOF，其它基于固定网格的方法如marker粒子，level set 等的文章也不少是用了这些方法。我也有意用SIMPLE求流场，VOF求界面，不过我的问题界面上有辐射换热和相变，还不太确定如何做。

zcliuqd

我现在也在做这方面的工作，主要是蒸汽的冷凝，但是进展不是很大，希望大家可以交流。

开心果

我可找到知音了，我觉得好孤独，上面的几位朋友，我们大家能不能交流一下，我得qq：59842985

jcn

我用的是有限元和vof求解自由表面问题。应该说采用什么方法求解流场和采用什么方法求解自由表面没有什么关系。毕竟现在的vof法都是基于几何的方法求解的，也就是通过几何的方法求解下一时刻的单元的体积分数。
但是，从我现在的情况来看，用vof求解自由表面容易产生“气泡”现象，即单元过空或单元过满。从误差分析来看，这种现象主要是由于体积通亮不守恒引起的。但是，同样的程序，在没有自由表面的情况下（比如求解圆柱绕流问题），结果非常好，涡街的形状和脱落周期都还算完美。但是一同vof一起求解就会出现以上问题。我的vof程序也是经过验证的（在流畅速度已知，且严格守恒情况下），精度较高。现在还没有解决这个问题，希望高手能够指点一二。

anyone

我没做过VOF, 而且VOF有不同方法.
不过据我所知, 有限体积的VOF一般是加解一个F(体积分数)的输运方程, 所以体积通量是守恒的. 其实VOF没那么复杂, 其中一种做法是算单一流体, 但考虑物性随F变化.如果用压力修正方法:
1.解流体速度
2.压力修正
3.解F输运方程
4.根据F分布捕捉自由面
5.根据自由面加表面张力(体力)
6.根据F分布调整物性(密度,粘性..)
VOF部分可以考虑用显式, 也就是不在带回去修正速度压力.