求教并讨论VOF界面重构

天空之城

[这个贴子最后由天空之城在 2008/04/24 09:55pm 第 1 次编辑]

各位大侠们，因为论文需要，自己编程用VOF算了一些模型，虽然算得还可以，但在编程过程中有一些疑问，望大侠们赐教。
1、刘儒勋老师那本自由液面追踪书后的代码我试过（剪切流流那个），为什么体积函数C的初始值是各点到圆心的距离？那样得出的C是一个连续分布的，而界面应该是间断，C应取是0或者1。
2、关于平移（同心框）、旋转（马蹄型）、剪切（圆）的那三个典型模型的计算结果界面为什么是一条线？而不是云图？按照书上的代码得出的文件，用TECPLOT画出来的应该是连续分布的云图。我们输出的文件中只能输出C的值，并不能输出具体在一个单元中界面是如何分布的（也没有必要），这样在显示软件中（TECPLOT）得出的也应该是云图，应该通过云图的边界来表达界面，在流体内部（C=1）和外部（C=0）中自然形成一条边界。
3、当然如果要说通过设置等值线来表达边界，比如设置C=0.5，那么我觉得这应该算是level set的思想。另外，C的分布是不连续的，我怎么就知道C=0.5的网格就一定能够形成一天边界呢，可能在边界上一些单元是0.5,令一些单元是0.9,0.1，即使设两条等值线也不一定能够保证能够完全把边界包括起来。所以我觉得仍应用云图，通过原图的边界来显示。
4、除了FLAIR，YOUNG这几个施主-受主模型中，界面重构直接影响传输的流体体积量，其他模型，如平均积分法等所谓的高阶格式中，界面重构与传输量有关系吗？他们的体积流通量中体现出了界面的形状和方向吗？既然体积输运的方程中没有体现出界面形状和方向，那么重构又有什么目的。