求解双曲型守恒律方程的高性能数值方法的分析.pdf

The Research of High Performance Numerical Methods for Hyperbolic Conservation Laws A Dissertation Submitted for the Master Degree of Science Candidate：Zhu Xioajun Supervisor：Prof. Feng Jianhu Chang’an University, Xi’an, China 摘 要 数值模拟方法的研究对计算流体力学有十分重要的意义，其中非线性双曲守恒律方 程由于其应用广泛和求解困难，一直以来受到广泛关注。以往的数值方法往往忽略了物 理系统的一个重要因素：热力学第二定律，即熵稳定条件，导致一些非物理现象的产生。 本文从物理概念出发，研究满足熵稳定条件的一类高精度高分辨率的熵相容格式，应用 于各种数值试验，主要创新点如下： (1) 熵相容格式比熵稳定格式进一步地控制了激波处的熵增量，可有效消除膨胀激波 及间断处的伪振荡等现象。通过在单元交界面处进行高阶 WENO/CWENO 重构，得到 一类高精度的熵相容格式。通过数值算例表明，改进后的格式具有可靠、高精度、高分 辨率、基本无振荡等特点。 (2) 提出了一种求解双曲型守恒律方程的自适应网格方法，即在数值解变化平缓的 区域网格点自动变得稀疏，而在数值解变化剧烈的区域网格点自动变得稠密。数值算例 表明该方法有效地消除了 “膨胀激波”、“负密度”和 “负压力”等非物理现象，说明 了该方法的鲁棒性、可靠性和无振荡性。 (3) 提出了求解二维双曲型守恒律方程的熵相容格式的构造方法和相关理论。首先 介绍了求解二维标量方程的构造方法，并用其求解Burgers 方程，取得了比较好的结果。 其次介绍了求解二维Euler 方程组的构造方法。最后给出相应的数值算例，算例结果表 明该方法能很好地求解二维Euler 方程组，有效地消除了 “红斑”现象。 关键词：双曲守恒律，WENO/CWENO 重构，比较原则，自适应方法，熵守恒/熵 稳定/熵相容格式, 旋转黎曼求解器 I Abstract Numerical simulation methods are of great significance for Computational Fluid Dynamics. Much attention has been paid to nonlinear hyperbolic conservation laws due to its wide range of applications and being difficult to solve. The previous numerical methods tend to ignore an important factor of physical systems—the second law of thermodynamics, which is entropy stable condition, and have generated some non-physical phenomena. Based on physical concepts, a new method called entropy consistent scheme which was entropy stable, accurate and high resolution was studied and applied to various equations. The main achievements of this work were as follows: (1) Compared with entropy stable scheme, entropy c