非线性守恒方程两种稳定化谱元方法pdf.pdfVIP

总 38 卷　第 4 期 数　学　研　究 V ol. 38　N o. 4 2005 年 12 月 Journal of M athem atical Study D ec. 2005 非线性守恒方程的两种稳定化谱元方法 杨炜铖　　许传炬 ( 厦门大学数学科学学院, 福建 厦门 361005) 摘　要　考虑非线性守恒方程的高阶数值解法, 介绍了基于谱元法的两种稳定性方法, 一种是 谱粘性消去法(SVV ) , 另一种是过滤法. 在 SVV 方法中, 我们推广并分析了传统的基于单区域的 SVV 算子的定义. 在过滤法中, 我们分析了 SVV - H elm holtz 过滤算子的性质. 文中从分析和计算 两方面对两种方法进行了比较, 建立了两者之间的关系. 最后通过一系列数值试验说明方法的有效 性. 关键词　谱元法(SEM ) ; 谱粘性消去法(SVV ) ; 过滤法 中图分类号　 174. 52　　　　文献标识码　 O A 1　引　　言 谱方法是一种求解偏微分方程的高阶数值方法. 谱方法最显著的优势在于它的“无穷阶” 收敛性, 即如果问题的解无限光滑, 那么谱方法求得的数值解将以 1 的任意次幂速度收敛 N ( ) 于精确解, 这个优点是其它有限阶方法 如有限差分法和有限元法 所无法比拟的. 由于上述 特点, 谱方法已经被广泛应用于椭圆、抛物型问题的求解中. 但传统的谱方法的局限在于无法 灵活处理复杂的计算区域, 80 年代末发展的谱元法很好地弥补了这个不足. 谱元法在快速收 敛的谱方法基础上结合了有限元法的区域剖分技术, 其基本思想是先将复杂区域剖分成若干 个较小的区域, 然后在每个区域上应用谱方法. 谱元法的灵活特点及高效性在计算流体力学 ( ) 中得到了很好的体现 [ 4 ] . 然而, 谱元法在处理守恒方程等双曲型方程时仍面临很大的困 难: 首先, 与椭圆、抛物型不同的是, 双曲型问题本身不具备物理耗散性, 即使很小的误差都可 ( ) 能引起数值解的不稳定 . . [ 1 ] . 更为主要的是, 当解出现断点时, 谱元解会产生 现 e g Gibbs 象, 导致无法得到收敛解. 而且, 即使解能够保持长时间的稳定, 计算结果也只有一阶精度, 这 使得谱元法作为高阶方法的有效性受到质疑. 为了克服这些困难, 人们提出许多改进的方法, 其中包括谱粘性消去法和过滤法. 谱粘法消去法(SVV ) 是通过引入一个人工粘性参与数值计算来控制数值解的不稳定. 该 方法在单区域情形下有比较完备的理论基础, 并得到了数值试验的验证. SVV 法有多种实现 方式, 例如 . . [ 3 ]提出的 法和郭本喻等人[ 2 ]提出的改进的 法等等, 均能 S M Kaber L SVV SV 　收稿日期: 2004- 04- 08