在两相流代码中添加空化模型的一些问题
-
原两相流代码是基于有限差分方法，
速度压力解耦是使用的投影法，
第一步：求中间速度V*
第二步：
》》》》
求压力
第三步：求解速度
VOF 使用的是THINC/WLIC方法，（https://doi.org/10.1016/j.jcp.2007.06.020）通量更新如下
现在想把SchnerrSauer模型加进去，求解的方程是如下四个
我的思路如下：
首先在某一点，利用时间步n上的值，r,p利用方程（4）计算得到m，把带入方程（1）（2），利于前面提到的投影法的思路计算中间速度，用方程1，3,得到压力泊松方程（方程里面带有m），使用隐式迭代方法计算压力，用压力计算得到速度，然后使用前面提到的THINC 方法更新相分数，这样整体思路有问题吗，
如果没有，有几个涉及的细节问题，
1，因为m里面带有体积分数α，压力P，在求解压力泊松方程的时候如果不考虑空化离散的方程如下
，使用迭代法求解，考虑空化的情况会有m相关的项，m里面涉及的压力p需要如何考虑，就直接用算出来的m放到压力泊松方程里面，然后 迭代求解吗，
2.同样相分数更新的时候涉及m的项里面也有α，这个如何考虑，也是直接用算出来的m做一个源项计算吗，
3，在上面提到的THINC 方法中对相分数的跟新使用的是x,y方向交替更新的，加入空化后涉及到源项m，我在交替更新的时候x,y方向，每一个方向考虑0.5*m可以吗，
问题有点多，如果有对相关问题熟悉或者了解的前辈，希望不吝赐教，非常感谢。
或者有涉及到求解附带空化模型的N-S方程算法的文章可以分享吗，可能是我找文献方式有问题，没找到关于如何求解附带空化模型的N-S方程算法的文章
