Structure-Preserving Hybrid Finite Element Methods

有限元方法是现代科学计算的基石之一，它将高性能与精度和稳定性的强大理论保证结合在一起。它们的应用包括联邦战略利益的几个领域：材料和制造、生物医学工程和生物技术、结构工程和民用基础设施、环境工程等。它们被学术界、工业界和国家实验室的科学家和工程师广泛用于模拟大型、复杂的物理系统。然而，为了使这些方法产生的模拟具有物理意义和可信度，它们最好能捕捉原始系统中存在的某些基本物理定律。某些系统满足“守恒定律”，即某些量(如质量、能量或电荷)可以在空间中移动，但不会自发地出现或消失，而另一些系统则满足“耗散定律”，要求这些量以一定的速度增长或衰减(例如，由于摩擦导致的能量衰减)。虽然经典的有限元方法很难表达这些物理定律，但PI已经证明了“混合”有限元方法提供了一种方法。在这个项目中，PI将对混合有限元方法进行进一步的研究，开发和分析保存这些重要物理结构的计算技术。该项目的成功将为各种高价值的科学应用带来新的计算方法和对现有方法的更好理解，对计算物理和工程具有重要的影响。许多偏微分方程(PDE)，特别是在物理应用中遇到的偏微分方程组，包含局部对称性、守恒定律和相关结构。乍一看，这些局部结构中的一些似乎很难被有限元方法捕捉到，因此，这方面的研究大多局限于矩形网格上的有限差分方法。然而，近年来，PI开发了一个成功的研究程序，表明杂交提供了一条绕过局部结构保留有限元的明显障碍的途径。这个项目将在这个方向上调查几个新的问题，围绕两个主要主题组织。(1)具有二次密度的局部守恒律在哈密顿偏微分方程组中是常见的，因为它们源于线性的点对称。然而，传统的有限元离散破坏了局部对称性，只产生了整体守恒律。PI的目的是通过使用混合方法来绕过这一障碍，并将其扩展到二次量被耗散而不是守恒的非保守系统。(2)PI和他的合作者最近成功地混合了有限元外微积分(FEEC)，它使用微分形式统一了几类拉普拉斯算子的方法。PI将研究混合FEEC方法在哈密顿偏微分方程组中的应用，特别关注多辛结构保持，并研究具有强(而不仅仅是弱)守恒通量的混合FEEC方法，包括非协调和可杂交不连续Galerkin(HDG)方法。该奖项反映了NSF的法定使命，并通过使用基金会的智力优势和更广泛的影响审查标准进行评估，被认为值得支持。

Finite element methods are among the cornerstones of modern scientific computing, combining high performance with strong theoretical guarantees of accuracy and stability. Their applications include several areas of Federal strategic interest: materials and manufacturing, biomedical engineering and biotechnology, structural engineering and civil infrastructure, environmental engineering, and more. They are widely used by scientists and engineers in academia, industry, and national laboratories to simulate large, complex physical systems. However, in order for the simulations produced by these methods to be physically meaningful and trustworthy, it is desirable that they capture certain fundamental physical laws present in the original systems. Certain systems satisfy "conservation laws," which state that some quantity (like mass, energy, or charge) can move through space but not appear or disappear spontaneously, while other systems satisfy "dissipation laws" that require these quantities to grow or decay at a certain rate (e.g., decay of energy due to friction). Although classical finite element methods make it difficult to express these physical laws, the PI has shown that "hybrid" finite element methods provide a way to do so. In this project, the PI will conduct further investigations into hybrid finite element methods, developing and analyzing computational techniques for preserving these important physical structures. The success of this project would lead to new computational methods and improved understanding of current methods for a wide variety of high-value scientific applications, with important ramifications for computational physics and engineering.Many partial differential equations (PDEs), particularly those encountered in physical applications, contain local symmetries, conservation laws, and related structures. At first glance, some of these local structures appear difficult for a finite element method to capture, and consequently, much of the research in this direction has been restricted to finite difference methods on rectangular grids. However, in recent years, the PI has developed a successful research program showing that hybridization provides a route around the apparent obstacles to local-structure-preserving finite elements. This project will investigate several new questions in this direction, organized around two main themes. (1) Local conservation laws with quadratic densities are common in Hamiltonian PDEs, because they arise from linear point symmetries. However, conventional finite element discretization breaks local symmetry, resulting only in global conservation laws. The PI aims to circumvent this obstacle by using hybrid methods, and to extend this to non-conservative systems where a quadratic quantity is dissipated rather than conserved. (2) The PI and his collaborators have recently succeeded in hybridizing finite element exterior calculus (FEEC), which uses differential forms to unify several families of methods for Laplace-type operators. The PI will investigate the application of hybridized FEEC methods to Hamiltonian PDEs, specifically focusing on multisymplectic structure preservation, and to investigate hybrid FEEC methods with strongly (rather than merely weakly) conservative fluxes, including nonconforming and hybridizable discontinuous Galerkin (HDG) methods.This award reflects NSF's statutory mission and has been deemed worthy of support through evaluation using the Foundation's intellectual merit and broader impacts review criteria.