邹东阳，林敬周，黄洁，刘君

1. 中国空气动力研究与发展中心 超高速空气动力研究所，绵阳 621000 2. 大连理工大学 航空航天学院，大连 116024

强扰动在空气中以超过声速的速度进行传播，表现为流体中出现一个参数的间断面，称为激波[1]。激波间断的存在给描述无黏流动的Euler方程的求解带来了很大挑战。过去半个多世纪里，为了发展出能够精确地描述激波的数值方法，科研人员做了大量的工作。如今这些工作是计算流体力学(Computational Fluid Dynamics, CFD)的重要组成部分。根据对激波处理方式的不同，CFD可以分为激波捕捉方法[2-8]和激波装配方法[9-11]2种。激波捕捉方法是以弱解理论为基础，全场采用统一的计算方法，无需特殊处理，具有通用性的特点。激波装配方法则是利用激波关系式(Rankine-Hugoniot，R-H)对激波进行描述，理论上是精确的，但是需要找到匹配激波间断解与其他光滑区域数值解的方法，存在兼容性的问题。对于复杂激波流动问题，特别是三维问题，激波装配方法存在应用上的困难。随着激波捕捉方法的成熟，特别是TVD(Total Variation Diminishing)格式的诞生，对于激波的计算逐步满足了一定的要求。过去30多年里，激波捕捉方法得到了迅猛发展，在航空航天、汽车建筑、天文地理等领域发挥了巨大的作用。

虽然激波装配方法的研究在相当长一段时间内寥若晨星，但是其在计算激波时表现出来的优势还是给许多激波捕捉方法研究提供了思路。像激波装配方法那样，将计算网格面与激波面相匹配来减少激波计算误差的处理手段在很多方法中都得到了应用。……