一维多介质可压缩Euler方程的高精度RKDG有限元方法.pdfVIP

物 计 算 理 第 23 #筋1 ，目 Vol.23 . No .l CHINESE JOURNAL OF COMPUTATIONAL PHYSICS 氢陆华1 月 ~主垦 [文窜每) I∞1 .2.46X( 筑邸)01.ω43-07 一维多介质可压缩 Euler 方程的高精度 RKDG 有限元方法 陈荣三l ， 蔚喜军2 (1.中罔E抑物理研究院北京研究生部.北京 附XlI!8. 2. 北京应用物理与11 舞敛学研究所计算物穆实事盘空， Jt 京 100088) [第 要} 采用 RXOG 有限元方法、 l.evel !)t 方法初改进的帮 Isenω仰 修正的 CII05I Fluid 方法後拟了一维多 介l1t可压缩[nJ盯方程， 其中 Euler 方视、Level 1 方程和重新初始化方.fÈ~采用了三阶精度的 RKDG 有限元方法 进行.散，并对一维两种介质可压缩流体进行了数值实验，得到f较芮分辩凛的计算结果. [ 关键惆1 RKDC 有限元方法; l.evel Sel 方法; Gh四1 日山d 方法j Tsentropio 修正 [中图分自峰等) 035 (文献串串识码J 。 引言 间断有限元方法是20 世纪70 年代 Reed 和 HiU[ L) 提出的.并用于解快中子输运问题.20 世妃 80 年代. C侃kbu，和 Shu 构造了(k + 1 ) 阶的 RKDG 街限元方法求解一维标量守恒律方程IZI ， 并将该方法成功地应用 到了多维标量守恒律川和多维守恒律方程组;41 . RKDG 有限元方法的优点是具有灵活处理间断和易于处理 复朵的区域边界和边值问题的能为. 运动界丽的追踪有很多方法，如格子类方法r~J 、 VOF 方法川和Lcvcl Sel 方法川等.Level Set 方法是把随时 间运动的物质界面看作某个函数 φ(x ， t.) 的零等值函，并满足~定的方程.在任意时刻只要知道了 φ忡.t ) 的函撒值.就可以知道等值丽的位置I !!IJ 物质界面的位置.Fet:让iw[ ~J 等在LeveL Sel 方法的基础引入了 Ghost Fluid 方法. Chost Fluià 方法的主要思想是利用Le崎] Sel 函数将含有两种流体的梳场分成只含有一种流体的 两个流场，每个流场囱真实流场利虚拟流场组成，每个流场可以单独进行计算，从而避开了界面两边状态方 程不同的问题. 本文对 Ghost Fluid 方法作了适当的改进，与中般的 Gbωt fluid 方法不同的是，首先通过求解界面附近的 局部 Riemann 解得到界面的压力、速度以及界面左边和右边的密度，进而求出界面左边和右边的铺，再将界 丽的偏延拓到 Gh倒区域.为了防止飞vc由ellting~误差，本文采用LiU[91 提出的Isentropic修正， 多介质可压缩流的数值模拟在流体力学、材料科学、天体物理学以及受控核聚变等领域有重要的应用. 本文利用 nKDC 有限元方法和改进的带..L;entropic修正的 Ghost F1uid 方法敬值模拟一维多介质可压缩 ElJler 方程组，用Level Set 方法追踪运动界面，其中 RKDG 有限元方法均采用三阶的格式，并对一般密度比利大街 皮比的-维的气-气、气 - 液两种介质可压缩流体进行了数值实验，得到了具有较高分辨率的计算结果. 1 方程 1. 1 E叫方程组 一维可压缩流体力学副er 方程组为