高密度和高粘度比率下气液两相流动的数值模拟+研究.pdfVIP

巾同-】：程热物理学会 传热传质学 第十■届学术会议 编号：063193 高密度和高粘度比率下气液两相流动 的数值模拟 孙东亮 陶文铨 两安交通大学动力rJ：程多相流园家重点实验窒，陕西西安710049 联系电话：02982669106E-maih wqtao@mail．xjtu．edu．cn 摘要；本文在VOF方法的基础上，采用粗细两套网格对高密度和商粘度比率下的气液两相流动模拟进 彳亍了研究分析。在细网格中求解流体体积函数方程，在粗网格中采用交错列格求解动最方程和压力修 正方程，通过粗细网格间的数据传递获得求解动量方程时需要的准确的界面密度和粘度及控制体密度， 解决了高密度和高枯度比率下通过插值方法计算界面密度和枯度及控制体密度带来的较大误差，保证 ，质量和动量同时守恒。高密度和商粘度比率下气液两榴流动中气液交界面处密度、速度和压力急剧 变化，为了保证格式的有界性和稳定性，采用稳定的有界高阶组合格式STOIC。擐后模拟J，不同工况 下气泡在液体中的运动，并通过实验和模拟结果验证了方法的可行性及准确性。 关键词：VOF方法、粑细跚格、高密度和高粘度比率、气液两相流动 1．前言 得到了长足的发展，并在气液两相流模拟方面表现出其优越特性，但传统的VOF方法 在模拟高密度和高粘度比率下气液两相流动时会出现求解不稳定，．收敛困难，误差较大 等问题，所以在模拟分析中一般采用较低的密度和粘度比率，如文献【l】中密度和粘度比 率为1：80；文献[2]rp密度和粘度比率为1：100，如何模拟高密度和高粘度比率下气液两 相流动成为本文研究的重点。本文采用粗细两套网格对高密度和高粘度比率下的气液两 相流动模拟进行了研究分析，在细网格中求解流体体积函数方程，在粗网格中采用交错 网格求解动量方程和压力修正方程，通过粗细网格间的数据传递获得求解动量方程时需 要的准确的界面密度和粘度及控制体密度，为了保证格式的有界性和稳定性，采用稳定 的有界高阶组合格式STOIC，表面张力采用了CSF模犁，最后模拟了不同工况下气泡 在液体巾的运动，并通过实验结果验证了方法的可行性及准确性。 2．控制方程 求解气液两相流动的控制方程为： 流体体积函数方程 eC+V．(历)：0 (1) j‘ 矾 动量方程 掣+v．(户孑孑)：一vJp+v．r+pg+F (2) dr 连续性方程V∥=0 (3) 基金项FI：国家kl然科学基金 1085 划分为四个细网格。本文数值 ．I I ．．I I+I ，：3 ． ‘ ‘ ‘ 模拟方法的分析就在如图1 图1粗细网格及变量#4iz示意图 所示的网格系统巾进行。 3．1求解控制方程的算法 本文采用Projection算法求解描述气液两相流动的控制方程，计算分为四步： 第一步：求解流体体积函数方程c”1=c”一舀V．(c”材”) (4) (5) 第二步：求解动最方程：p川矿：p一万+国卜v．(p一驴歹)一VP一+v．f一+p一歹+可} 第三二步：求解压力修正方程 v．(vP，／∥+·)：(