气液两相流中存在气相空间压力不平衡的问题
-
简要描述：在计算两个容器中间的连通管内进行气液两相的重力流时，一开始设定的相同的压力0.0Pa表压，计算开始后，出现两个容器，和联通管上较大的压力差变化，导致两相流动不能顺利进行，不明白计算的压力变化从哪里来的，请各位大神指教
网格无关性基本验证了，网格划分加密5倍后，问题仍然存在。
以下是两个详细的记录文件，问题描述不清楚的，可以下载看看
某项目气液两相流模拟分析过程记录.pdf
网格划分生成报告.pdf
以下是较为详细的描述
项目背景：
在装置框架上，两台设备（KA103,VA118）高低布置，液体从高处设备(KA103)通过溢流堰，进入连通管，连通管接入低处设备(VA118)的顶部，VA118内部设导流管进入设备下部，液体向上返流后通过设备上部的环形堰，堰内液体通过管口流出，两台设备的气相区域设气相连通管，顶部均设排气口，通向外界。
项目目的：
从联通管进入VA118的液相，会夹带部分气体进入VA118，带入的气相在VA118内气液分离后会通过气相联通管回到KA103，核算液相通过连通管时的带气量，气相连通管能否满足气相联通要求。
模型概况：
模型建模采用GAMBIT，网格划分及计算采用ANSYS FLUENT 19.0。
设备基本规格：
KA103 DN3400mm H=16390mm
VA118 DN5000mm H=12338mm
为降低模型的复杂度，对于设备中不必要的管口进行简化，仅保留需要分析的接口及内部结构，具体情况需要对照设备图纸进行说明。
以下图片为设备的基本结构和内部结构的简图（忽略其中的不规格锯齿线，软件生成的，非结构线）
网格划分基本按照软件默认进行，四面积，未处理边界层，只在连通管等小直径区域进行了适当的加密处理，总体网格数量在161万左右，详细内容见软件生成的网格报告。
计算设定：
总体设定：采用基于压力的模型，瞬态，设置Z轴向重力-9.81m/s2;
计算模型：湍流模型采用标准k-ε模型，多相流采用VOF模型，气相为主相，液相为第二相，两相的交互系数考虑表面张力系数0.075n/m；其它系数保持默认；
介质：液相介质采用水，气相介质采用空气，介质性质来源于软件自带参数，未做修改；
单元区域条件：操作压力：101325pascal（绝压），其余保持默认；
边界条件：液相进口采用速度进口，流速1.7m/s，液体分率0.9；液体出口按自由流出口；两个气相出口按压力出口，设置参数为表压为0，与外界环境压力保持一致；
模型计算设定：
解算模型采用SIMPLE，压力采用PRESTO!，其余参数未更改，保持默认；初始化设置：
初始压力设定为表压0pascal（表压），对两台设备的两相界面进行patch操作，完成初始液面位置的设定
说明：界面不平滑，是因为靠近连通管位置网格粗糙
计算中发现的问题：
计算开始的0.005s内，每0.001s，20步迭代计算后未收敛，设定的两个气体出口存在气相回流，猜测的引发原因是两台设备内出现的低压，压力严重偏离初始设定的 0.000 pascal，设备内部出现负压的原因未知。之前尝试调整出口的边界条件，包括自由流，压力出口等，均出现以上问题。在后续计算的过程中，甚至出现了连通管内液体逆向回流的问题。
以下为初始阶段的几张软件截图
以下附图中为计算过程中截取的状态画面，便于开展过程状态的分析
![运行结果2.jpg]
