一类稳定混合型间断有限元方法：原理、特性与应用新析.docxVIP

一类稳定混合型间断有限元方法：原理、特性与应用新析

一、引言

1.1研究背景与意义

有限元方法作为一种强大的数值计算技术，在现代科学与工程领域中占据着举足轻重的地位。自1943年Courant提出三角形网格剖分的Dirichlet问题的分片线性逼近，以及我国计算数学家冯康先生独立发现该方法以来，有限元方法不断发展，已成为一门理论完善、应用广泛的数值计算方法。其基本原理是将原始问题转化为变分形式，在较弱的空间上求解，然后构造有限维空间逼近变分问题的求解空间。通过将求解区域剖分成许多小片，构造分片多项式逼近解，有限元方法能够有效处理各种复杂的数学物理问题，如二阶椭圆问题、抛物问题、双曲问题，以及流体中的Stokes问题等。

然而，传统有限元方法对单元间连续性要求较高，这在一定程度上限制了其应用范围。在处理具有奇性、振荡、边界层等特征的问题时，传统有限元方法往往难以精确逼近，导致计算精度下降。例如，在求解含有间断解的偏微分方程时，传统有限元方法会产生数值振荡，无法准确捕捉间断点的位置和性质。为了克服这些局限性，间断有限元方法应运而生。

间断有限元方法（DiscontinuousGalerkinMethod，简称DGM）在很大程度上继承了有限元方法的优点，同时极大地放松了对单元间连续性的要求。它允许有限元函数空间基函数完全间断，各个单元之间通过在单元边界上构造合适的数值流通量来实现通信。这种独特的特性使得间断有限元方法能够更精确地逼近具有奇性、振荡、边界层等特征的问题。例如，在计算流体力学中，间断有限元方法可以准确捕捉激波和接触间断，在光滑区域也能保证高精度，并且在间断区域保持数值无振荡，分辨率高。在处理具有复杂边界条件、波束转移、非线性等问题时，间断有限元方法也表现出明显的优势，其对物理界面的适应性更好、数值稳定性更高、计算效率更高且可并行化程度更高。

在众多间断有限元格式中，局部间断Galerkin有限元方法（LDG）具有独特的优势。它能够在形式上显式求解，容易实现并行计算，并且具有很好的稳定性。然而，不同的间断有限元格式在稳定性、误差估计和超收敛性等方面存在差异，对于一些复杂的问题，现有的间断有限元方法仍有待进一步完善和优化。

本文主要研究一类稳定混合型间断有限元方法，针对经典椭圆方程，给出一种由单元上的残量决定稳定项的离散格式。通过深入讨论该格式的有界性、稳定性及相容性，并给出在所定义范数下的最优误差估计，旨在为间断有限元方法的理论研究提供新的思路和方法，进一步拓展其在科学与工程计算中的应用范围。该研究对于提高数值计算的精度和效率，解决实际工程中的复杂问题具有重要的理论意义和应用价值。

1.2国内外研究现状

间断有限元方法的研究始于20世纪70年代，Reed和Hill在1973年为解决中子输运方程问题首次提出了间断有限元方法，这一开创性的工作为后续研究奠定了基础。此后，众多学者对间断有限元方法展开了深入研究，不断提出新的改进和发展。

在间断有限元方法的发展历程中，内惩罚(IP:InteriorPenalty)类方法在20世纪70年代被独立地提出来求解椭圆和抛物方程，后来被归为间断Galerkin方法的一种，记为内惩罚间断Galerkin(IPDG)方法。该方法的侧重点在于选择合适的惩罚项以保持格式的稳定性。

20世纪90年代以来，以Cockbum和舒其望为代表提出了Runge-Kutta间断Galerkin(RKDG)方法。该方法结合TVD(TVD:TotalVariationDiminishing)Runge-Kutta时间离散方法和间断有限元，用于求解一维双曲守恒律方程(组)乃至高维双曲守恒律方程(组)。RKDG方法能够适应复杂计算区域和边界条件，可以精确捕捉激波和接触间断，在光滑区域保证高精度，在间断区域保持数值无振荡且分辨率高，还可证明收敛到熵解，这些优点使其成为计算流体力学中流行的方法之一，并被广泛应用到气象学、海洋学、湍流、电磁学、石油勘探、水动力学、等离子物理和图像处理等众多领域。

1998年，Cockbum和舒其望提出了局部间断Galerkin有限元方法(LDG)。其思想来源于F.Bassi和S.Rebay求解可压缩Navier-Stokes方程的文章，通过引入辅助变量将含有高阶导数的微分方程写成只含有一阶导数的偏微分方程组，然后用DG方法进行空间离散。当LDG格式中基函数取为每个单元上的局部基函数时，辅助变量q是局部可解的，这也是该方法被称为“局部”间断Galerkin方法的原因。LDG方法能够在形式上显式求解，容易实现并行计算，并且具