二维多介质可压缩流的RKDG有限元方法.pdfVIP

计 物 算 理 第23 咎第6 期 Vol.23 ,N..6 组踊年 11 月 CIDNESE JOURNAL OF COMPUTATIONAL PHYSICS Nov. ，组勘6 [ 文章编号1 lOOl.246X(2伽)06-0699-的 二维多介质可压缩流的 RKDG 有限元方法 陈荣三1 ， 蔚喜军2 (1.中西工程彷理研究院北京研究生部，北京 以削8; 2. 北京应用物理与计算数学研究所，计算物理实盘室，北京 阳刚的 [摘 要] 应用 RKDC (Runge- Kutta Oiscontinuous Galerkin) 有限元方法、Level Set 方法和 Ghost Fluid 方法数值模 拟二维多介质可压缩流，其中 Euler 方程组、Level Sel 方程和军新初始化方程的空间离散采用 DG (Di盹ontinuous Galerkin) 有限元方法，时间离散采用 Runge-Kuu8 方法 对二缘的气-气和气- 液两相流进行了数值计算，得到 了分辨率较高的计算结果 [关键词1 RKDG 有限元方法;多介质可压缩流; Gh耐F1uid 方法; Level Sel 方法 [中图分类号1 035 [文献标识码1 A 。 引言 间断有限元方法最早由 Reed 和 Hill 在 1973 年提出，用于求解中子输运方程.随后，一些计算工作者将 它推广到求解非线性双幽守恒律方程和方程组I 川J 另外，间断有限元方法被应用于求解椭圆方程I ， Hamilton-Jacobi 方程[的，对流扩散方程(1) 等. 间断有限元方法具有灵活处理间断和易于处理复杂的区域边界 和边值问题的能力. 在多介质问题的数值计算中，正确地追踪运动界面的位置相当重要.运动界丽的追踪有很多方法，直口格 子类方法、VOF 方法(8) 、Level Set 方法191 和波前追踪方法1101 等. Level Set 方法利用Level Set 函数 φ 的零等值 面表示两种流体的分界面，Level Set 函数值为该点到界面的符号距离，不同符号的Level Sel 函数值表示不同 lll 的流体.1999 年，Fedkiw ) 等在Level 如t 方法的基础上引人了 Ghost Fluid 方法.Ghost Fluid 方法的主要思想是 利用Level Set 函数将含有两种流体的流场分成只含有…种流体的两个流场，每个流场由真实流场和虚拟流 场组成，每个流场可以单独进行计算，从而避开了界面两边状态方程不同的问题. 本文将Level Set 方法和创ost Fl山d 方法相结合，处理多介质界面，用 RKDG 有限元求解Level Sel 方程、 重新初始化方程和流体力学方程组. 对二维的气 - 气和气-液的双介质可压缩流进行了数值计算，得到了分 辨率较高的计算结果. 1 方程 1.1 Euler 方程组 二维可压缩流体力学 Euler 方程组为 ρu I I ρu + p I I pUV I ( I ) + I I + I I = (ρ ) ll l m i pv I I pUV I I PV + P I E )， 飞 (E + p)U). 飞 (E + p) V), 这里ρ， p