基于THINC/QQ格式的两相界面流动数值模拟  

Numerical simulation of two-phase free-surface flow problems based on THINC/QQ scheme

作　　者：赵海洋 明平剑[1] 张文平[1] 齐文亮 ZHAO Haiyang;MING Pingjian;ZHANG Wenping;QI Wenliang(College of Power and Energy Engineering, Harbin Engineering University, Harbin 150001, China)

机构地区：[1]哈尔滨工程大学动力与能源工程学院,黑龙江哈尔滨150001

出　　处：《哈尔滨工程大学学报》2020年第12期1804-1810,共7页Journal of Harbin Engineering University

基　　金：国家自然科学基金项目(51479038);中央高校基本科研业务费重大项目(HEUCFP201711).

摘　　要：为了建立精确高效的两相自由面模拟程序,本文在实验室自主研发的N-S方程求解器基础上,开发了基于二次曲面和高斯积分的双曲正切函数界面捕捉模块。基于同位网格有限体积法对方程进行离散,采用SIMPLE算法处理速度压力耦合问题,应用连续表面力模型模拟两相流体间表面张力。对给定剪切速度场下液滴在非结构化网格中变形问题进行模拟,并对二维气泡上浮问题和三维双气泡融合问题进行仿真。实验表明:在经历较大变形后液滴仍能回到初始形状,证明了THINC/QQ格式的数值精度;数据结果与Level Set方法计算得到的结果吻合良好,证明了当前两相界面模拟程序的准确性。To establish an accurate and high-efficiency two-phase free-surface simulation program,the two-phase interface capturing scheme THINC/QQ(tangent of hyperbola interface capturing method with quadratic surface representation and Gaussian quadrature)was developed and employed based on the in-house software General Transport Equation Analyzer(GTEA)for solving N-S equations.The equations were discretized by the finite volume method on a collocated grid arrangement.The semi-implicit method for pressure-linked equations(SIMPLE)algorithm was used for the velocity-pressure coupling problem,and the surface tension force between the two-phase fluids was calculated using the continuum surface force(CSF)model.First,the deformation of a drop at a given shear velocity field was simulated.Then,the 2D bubble rising problem and the 3D coalescence problem of two rising bubbles were also simulated.The drop could return to its initial circle after great deformation,demonstrating the accuracy of the THINC/QQ scheme.The numerical results for bubble ring and coalescence cases agreed well with corresponding solutions calculated using the level set method in the literature,demonstrating the accuracy of the current two-phase simulating code.

关 键 词：两相流 非结构化网格 THINC/QQ 体积分数法 有限体积法 气泡上浮 气泡融合 数值仿真 

分 类 号：TK411[动力工程及工程热物理—动力机械及工程]