模拟相变问题的精细积分边界元法  

Simulation of phase change problems by precise integration boundary element method

作　　者：姚鸿骁 左冲 胡小飞[1] 姚伟岸[1] YAO Hong-xiao;ZUO Chong;HU Xiao-fei;YAO Wei-an(State Key Laboratory of Structural Analysis for Industrial Equipment,Dalian University of Technology,Dalian 116024,China;Sany Heavy Industry Co.,Ltd.,Wuhan 430100,China)

机构地区：[1]大连理工大学工业装备结构分析国家重点实验室,大连116024 [2]三一重工集团有限公司,武汉430100

出　　处：《计算力学学报》2022年第5期545-550,共6页Chinese Journal of Computational Mechanics

基　　金：国家自然科学基金(11672064)资助项目.

摘　　要：将精细积分边界元法和界面追踪法相结合求解相变问题。因为边界元法只需要将待求解空间域的边界离散,方便连续追踪移动界面位置和重构网格,所 以边界元法适合应用于移动边界问题的模拟。首先,利用精细积分边界元法在固相区域和液相区域分别求解相应的瞬态热传导控制方程,从而求得温度场 和边界热流密度。然后,根据固-液相变界面上的能量平衡方程,利用热流密度求得相变界面的移动速度,再采用界面追踪法预测移动相变界面的位置变 化。最后,给出了几个数值算例,并通过与参考解的对比验证本文方法的准确性。The precise integration boundary element method is combined with the front tracking method to simulate the phase change problems in this paper.One of the merits for the boundary element method is that only the boundary needs to be discretized,so it is easy for the implementation of grid deforming and remeshing.Firstly,the governing equations of heat conduction are solved by the precise integration boundary element method in the solid phase domain and the liquid phase domain respectively,and the moving velocity of the solid-liquid interface can be obtained by the heat flux on the interface and the energy balance equation.Then,an iterative algorithm of the front tracking method is adopted to track the moving interface.Finally,several numerical examples are presented to verify the validity and accuracy of the proposed method.

关 键 词：相变问题 精细积分法 径向积分边界元法 界面追踪法 

分 类 号：O241.8[理学—计算数学]