VOF方程的疑问
-
各位老师同学大家好，最近遇到了点小问题，真心求帮助。
假设流体不可压缩。MULES算法中，添加了人工压缩项以保证界面尖锐：
按我的理解，其实就相当于在VOF方程中引入了一个源项S
虽然界面压缩项只在自由表面附近起作用，但是严格来讲，连续性方程应该修改成如下形式才对吧？(因为对于不可压缩流体，相体积分数的控制方程与连续性方程是一致的)
但是在interfoam里面并未发现类似的代码，请问这是怎么回事呢？是对于这个问题，这个源项可以忽略吗？
推广一下，如果想在VOF方程中添加其他的源项，是不是也要修改质量守恒方程？
此外，对于相体积分数为0的网格，上式中除以相体积分数这个操作是不是需要特殊处理一下？我试了一下，假设这样
程序运行倒是没什么问题，但是这个意义我没搞清楚，这么一除不是得到无限大的速度散度了吗？目前，我倾向于这么做：对于相体积分数为0的网格，我直接将源项S的值设为0，即源项不起作用。说的可能有点乱，但主要的意思就是这样。。。
谢谢大家！
-
额 好像发现个问题，质量守恒方程应该是这个
其中𝑆(𝜌)与 𝑆(𝛾)形式一样，只不过自变量由相体积分数变成了密度，不知道这样对不对？
感觉有点弄乱了，求大佬指点，不胜感激
