高阶导数在DG方法中的数值通量.docVIP

高阶导数在DG方法中的数值通量

一、DG方法与数值通量的基础关联

不连续Galerkin（DiscontinuousGalerkin，DG）方法作为一种有限元方法，自20世纪70年代被提出以来，凭借其在处理复杂边界条件、自适应网格和间断问题上的独特优势，逐渐成为计算流体力学、计算电磁学等领域的重要数值工具。该方法的核心在于允许解在单元边界处不连续，这一特性使其能够高效捕捉激波、接触间断等物理现象，但同时也带来了如何在单元间传递信息的挑战。数值通量（NumericalFlux）正是解决这一问题的关键，它作为单元边界上的数值近似，承担着相邻单元之间信息交换的任务，其构造直接影响到DG方法的稳定性、精度和收敛性。

在传统的连续有限元方法中，解在单元边界处连续，因此可以通过直接的导数传递来实现单元间的信息交流。然而，DG方法中解的不连续性使得这种直接传递不再可行，必须引入数值通量来近似单元边界上的物理通量。数值通量的构造需要满足一定的条件，例如一致性、稳定性和守恒性等，以确保数值解能够准确反映物理问题的本质。高阶导数在DG方法中的引入，为数值通量的构造提供了新的思路和方法，使得DG方法能够在保持不连续性优势的同时，进一步提高数值解的精度。

二、高阶导数对DG方法精度的提升机制

（一）高阶导数与多项式基函数的关系

DG方法通常采用多项式基函数来逼近单元内的解，基函数的阶数直接决定了数值解