二维浅水波方程的数值激波不稳定性.pdf

第29卷第l期 计 算 物 理 V01．29．No．1 PHYSICS 2012年1月 CHINESE OFCOMPUTATIONAL Jan．．2012 JOURNAL 文章编号：1001．246x(2012)0l_0025—11 二维浅水波方程的数值激波不稳定性 沈智军1’2， 胡立军3， 闫 伟1 (1．北京应用物理与计算数学研究所计算物理实验室，北京100088； 2．北京大学应用物理与技术研究中心，北京10087l；3．中国工程物理研究院研究生部，北京100088) 摘 要：研究二维浅水波方程的数值激波不稳定性问题．线性稳定性分析和数值实验表明，格式的临界稳定性与 数值激波的不稳定现象有重要的联系．基于扰动量的增长矩阵分析，本文将高分辨率的数值格式和HLL格式进 行特定的加权，设计一类新的混合型数值格式．其中可以调节非线性波速的HLLc与HLL的混合格式，数值试 验展示了消除浅水波方程激波不稳定现象的有效性和鲁棒性． 关键词：激波不稳定性；浅水波方程；线性稳定性分析；混合通量格式 中图分类号：0241 文献标识码：A O 前言 在近几十年中，双睦型守恒律方程的计算格式得到了蓬勃的发展．以通量分裂格式和Godunov型格式为 代表的流体力学数值方法取得了长足的进步．这些格式在大多数情况下都表现出很好的健壮性．然而，当计 算二维问题时，它们还是会出现一些非物理现象．数值激波不稳定性就是其中的现象之一． 3在二维绕柱钝体流动的计算中首先发现，使用Roe格式会产生所谓的“瘢痕现象” Peery和Imlay[1 似Riemann解，总是会产生一种“奇偶失联”(odd．even 与“瘢痕现象”总是相互伴随，产生前一种现象的数值格式，一定会产生后一种现象，反之亦然．由于这些非 物理现象常常出现在激波的波阵面附近，所以也常被称为激波不稳定现象．在双马赫反射计算和环绕斜体的 半圆锥激波的计算中口。4]，激波不稳定性以另外的形式出现．这些现象的产生与计算格式的内在性质相关， 因此引起了人们的高度关注．探索这种现象产生的原因和治愈这种缺陷的方法一直是流体力学数值计算中 的一个热点m18]． 上述研究激波不稳定性问题的文章，均未涉及浅水波方程的激波不稳定性问题．研究浅水波方程，除了 本身的重要性以外，还可以加深对不稳定现象的理解．与欧拉方程相比，浅水波方程仅包括剪切间断而不包 括接触间断．数值计算中仍然出现了激波不稳定现象，将表明激波不稳定现象的产生，确实与接触间断没有 本质联系，而可能是源于其它波系的不稳定性． 本文讨论二维浅水波方程的数值激波不稳定性问题．在第1节简介了二维浅水波方程和数值格式后，第 2节利用von 增长矩阵的特征值．相比谱半径小于l的严格稳定，类似文[3]引出临界稳定性概念．第3节根据扰动量的 增长矩阵分析，确定产生数值激波不稳定性的网格界面位置，并给出了激波不稳定现象的一种解释．第4节 针对二维浅水波方程提出了一种新的数值格式，其主要思想是将较高接触波分辨率的数值格式和HLL格式 进行加权求和得到混合型通量格式，将临界稳定格式改为严格稳定格式．在第5节，进一步的数值实验表明， 非线性波上耗散不足，也是导致数值激波不稳定的重要因素．可以调节非线性波速的HLLc与HLL的组合 收稿日期：2011一03—18；修回日期：2011一05—27 验室基金资助项目 作者简介：沈智军(1966一)，男，辽宁，博士，研究员，从事计算流体力学方法研究 26 计 算 物 理 第29卷 格式，展示了消除浅水波方程激波不稳定现象的有效性和鲁棒性． 1 二维浅水波方程和数值方法 掣+掣旦+鼍旦：o， (1)、 。