欧拉方程roe格式与高精度半拉氏方法研究.pdfVIP

分类号……………… 密级……………… UDC 编号……………… 中国工程物理研究院 学位论文 欧揎赢程．胝．格式点赢精度半撞氏方法砑究 齐进 指导教师姓各……王双虎…研究员…………………………… 申请学位级别………顾±…………专业名称……让算数学…… 论文提交日期……．20．1．虬5…………论文答辩日期…～2Q1Q：§… 学位授予单位……史．国王程物理砑巍院……………………… 答辩委员会主席…江…松一研究员 万方数据 I删煅 独创性声明 本人声明所呈交的学位论文是本人在导师指导下进行的研究工作及取 得的研究成果。据我所知，除了文中特另,J30n以标注和致谢的地方外，论文中 不包含其他人已经发表或撰写过的研究成果，也不包含为获得中国工程物理 研究院或其他教育机构的学位或证书使用过的材料。与我--MI作的同志对 本研究所做的任何贡献均已在论文中作了明确的说明并表示谢意。 学位论文作者签名：首边 签字日期．≯f哞／月彤日 学位论文版权使用授权书 本学位论文作者完全了解并接受中国工程物理研究院研究生部有关保 存、使用学位论文的规定，允许论文被查阅、借阅和送交国家有关部门或机 构，同时授权中国工程物理研究院研究生部可以将学位论文全部或部分内容 编入有关数据库进行检索，可以影印、缩印或扫描等复制手段保存、汇编学 位论文。 0 、 学位论文作者签名： 前地 J 靳虢钞父了阮 签字日期：肋19年易月谚日 签字日期：力步年占月吃阳 万方数据 摘要 多尺度多介质复杂流动问题在惯性约束聚变、天体物理等科学与工程领域有 着广泛的应用背景，由于理论分析和实验的困难，数值模拟已经成为研究这类问 题的主要手段。在以有限体积方法为基础的流体力学数值方法中，其中一个关键 问题是如何求解网格边界的数值通量，Godunov方法是目前最常用的方法之一。 由此本文分析了几种适用于一般状态方程的Godunov型格式的数值特性，在此 基础上给出了一种修正的Roe格式，并发展了一种新的有限体积半拉氏方法。 目前国内外已发展了许多基于近似黎曼解的高阶数值方法，广泛应用于各类 实际问题的计算，但其数值性能仍有待深入系统的分析比较。Roe．Pike、HLL、 HLLC方法是目前较为常用的可以用于一般状态方程求解的近似黎曼解方法。 MUSCL和PPM方法是两种常用的能够抑制振荡的高精度数值方法。本文利用 这两种方法的高阶无振荡重构技术，基于上述三种近似黎曼解方法，构造了一系 列的高精度Godunov型方法。随后通过有代表性的数值算例对这几种方法进行 了分析比较，考察了它们的数值特性和优缺点，从而将有助于在实际应用中选择 适宜的方法，并可为新方法的设计提供了参考。 针对Roe矩阵构造过程的复杂性，本文给出了一种修正的Roe方法，在保 证守恒性的同时，简化了Roe矩阵的限制条件并构造了新的Roe矩阵，而且能 够方便地推广到一般状态方程和多维情形。另外，为了有效模拟低密度流动问题， 本文给出了一种简单有