气液欧拉双流体模型动量方程奇异问题

[刘宗库]

各位老师好，我再用气液双流体模型做管道充水过程，在对空管道充水过程中，水柱通过后就全成液相，而未充水的管道还是气相，这个时候一相体积分数几乎为0，这个时候求解双流体模型就会出现问题。另外我是利用SIMPLE算法求解守恒型控制方程。请问如何从数值计算上解决这个问题呢？麻烦大家了

[李东岳]

不是OpenFOAM？可以借鉴OpenFOAM的处理方式

http://dyfluid.com/reactingTwoPhaseEulerFoam.html

[刘宗库]

@李东岳 谢谢李老师回复，我自己用Matlab编程，那我参考一下OpenFOAM里面的处理方式，我也在学习李老师的相关论文

[刘宗库]

@李东岳 李老师您好，我看了您的论文，想请教您问题。在您的论文《Simulation of bubbly flows with special numerical treatments of the semi-conservative and fully conservative two-fluid model》中，您对于双流体模型动量方程奇异的问题，提到了两种处理方法并指出他的不足，即第一种限制相分数，第二种指定限定速度。他们的不足分别是不满足全域守恒和计算结果具有随机性。在接下来，您提到了给动量方程除以各自体积分数将方程转换为非守恒形式，这样不用指定速度。通过学习OpenFOAM的处理方式，我有两个问题，1，OpenFOAM中我看要指定速度，是不是采取了您提到的第二种方法？2，您提到的采用非守恒形式的动量方程是不是就可以自动求解，而不用限制相分数或者指定某一相的速度，也就是说您提到的动量方程的非守恒处理方法是不是不同于您提到的前两种方法？

[李东岳]

1，OpenFOAM中我看要指定速度，是不是采取了您提到的第二种方法？

新版本的OpenFOAM是守恒形式，就是这里写的 http://dyfluid.com/reactingTwoPhaseEulerFoam.html

2，您提到的采用非守恒形式的动量方程是不是就可以自动求解，而不用限制相分数或者指定某一相的速度，也就是说您提到的动量方程的非守恒处理方法是不是不同于您提到的前两种方法？

非守恒形式就是你的图中的方程10

[刘宗库]

@李东岳 李老师，请问如果求解非守恒型的动量方程，那么对于其中一相为0的情况，是不是就直接求解就好了，不用再做别的处理了，如让相为0的相的速度等于按照模型去求出的速度。