二维浅水波方程和欧拉方程数值激波不稳定性剖析.pdf

摘要 当使用能够精确分辨接触间断的数值格式来计算多维流体力学问题时，会发 现激波附近的扰动有时会出现剧烈增大的现象，这就是所谓的数值激波不稳定性 现象。许多文献对二维欧拉方程部分数值格式的不稳定性现象进行了分析和解 释，并给出了治愈激波不稳定性的方法。但几乎所有的文献均未讨论过浅水波方 程的激波不稳定性问题。本文对二维浅水波方程和欧拉方程数值激波不稳定性问 题进行了研究。 在对数值激波不稳定性研究现状进行回顾之后，论文的第一部分研究二维浅 水波方程的数值激波不稳定性问题。通过对部分数值格式的线性稳定性分析和相 关数值试验，我们发现临界稳定性与部分格式的激波不稳定现象存在密切联系。 部分格式在剪切粘性和非线性波粘性的不足，导致了二维浅水波方程的数值激波 不稳定性现象。利用临界稳定性分析，文中提出了一种能够消除数值不稳定性现 象的混合型数值格式。 论文的第二部分将二维浅水波方程数值激波不稳定性的研究方法推广N- 维欧拉方程。与浅水波方程研究不同的是，欧拉方程部分格式的稳定性除了需要 足够的剪切粘性，还需要对质量方程进行改进。与之类似的是，通过修改临界稳 定的特征值，构造了相应的混合型数值格式来消除二维欧拉方程的数值激波不稳 定性现象。 数值试验展示了混合格式对于消除二维浅水波方程和欧拉方程激波数值不 稳定现象的有效性和健壮性。 关键词：数值激波不稳定性，二维欧拉方程，二维浅水波方程，线性稳定性分析， 混合型通量格式。 万方数据 Abstract When the numericalschemeswhichCall using accurately to capturediscontinuity calculate multidimensionalfluid mechanics willfind problems，wethe perh】fbation nearthe shockwave increases iswhat dramatically．Thiswecall n_um鲥calShock instability．Severalhavebeenmadetounderstandandcure attempts me instabilitv on 2DEuler phenomenon few equations．Butliteraturesdiscussed t11enu】[11两cal shock of2D instabilityshallowwater