低密度低压力流动问题的高精度数值模拟方法研究-计算数学专业论文.docx

承诺书 本人声明所呈交的硕士学位论文是本人在导师指导下进 行的研究工作及取得的研究成果.除了文中特别加以标注和致 谢的地方外，论文中不包含其他人已经发表或撰写过的研究成 果，也不包含为获得南京航空航天大学或其他教育机构的学位 或证书而使用过的材料. 本人授权南京航空航天大学可以将学位论文的全部或部 分内容编入有关数据库进行检索，可以采用影印、缩印或扫描 等复制手段保存、汇编学位论文. （保密的学位论文在解密后适用本承诺书） 作者签名： 日 期： 南京航空航天大学硕士学位论文 南京航空航天大学硕士学位论文 低密度低压力流动问题的高精度数值模拟方法研究 低密度低压力流动问题的高精度数值模拟方法研究 I I PAGE PAGE IV 摘 要 在爆炸波等多介质流体力学问题的数值模拟中，会出现低密度低压力区域，由于数值误差 的影响，导致在计算过程中该区域的密度和压力出现负值，使得计算难以继续。X.X. Zhang 等 给出了求解理想气体欧拉方程的正保护 RKDG(Runge-Kutta Discontinuous Galerkin)方法。本文 在此基础上，在不改变正保护限制器性质的情况下，改进了正保护限制器，节省了计算量，给 出了多介质流体力学问题数值模拟的正保护 RKDG 方法，结合 RGFM(Real Ghost Fluid Method) 界面处理方法，对一维和二维多介质可压缩流体进行了数值试验，结果表明该方法能有效地避 免密度或压力出现负值，并能准确地捕捉各类间断。 关键词：RKDG 方法，低密度，低压力，正保护限制器，RGFM 方法 ABSTRACT In the numerical simulations of the multi-medium flow such as blast waves, the negative density or pressure may occur due to the numerical errors of the high-order scheme in the low density and pressure domain. Failure of preserving positivity of the density or pressure can cause blow-ups of the numerical algorithm. In this paper, based on the idea of Zhang for ideal gas, the limiter is modified and extended to the simulation of multi-medium flow. The obtained method is cost effective and keeps the property of the positivity-preserving limiter. The real Ghost Fluid method is coupled to treat the interface. Several one-dimensional cases and two-dimensional cases are given to test robustness and efficiency of the algorithm. Numerical results show that the obtained method can maintain the positivity of pressure and density and can capture the discontinuities accurately. Keywords: RKDG method, Low density, Low pressure, the Positivity-preserving Limiter, RGFM method III III PAGE PAGE IV 目 录 第一章 绪论 1 1.1 本文的选题背景 1 1.2 本文的主要工作 2 1.3 本文的内容安排 2 HYPERLINK \l _TOC_250003 第二章 RKDG 方法 4 2.1 控制方程 4 2.2 状态方程 4 HYPERLINK \l _TOC_250002 2.3 DG 有限元方法 5 HYPERLINK \l _TOC_250001 2.4 RK 时间离散方法 6 HYPERLINK \l _TOC_250000 2.5 TVB 限制器 6 第三章 正保护限制器 8 3.1 流体为理想气体 8 3.1.1 欧拉方程组无源项 9 3.1.2 欧