玻尔百科准确模拟物理世界——从机翼上的气流到穿越地球的波的传播——的探索，从根本上说，是求解偏微分方程（PDE）的挑战。几十年来，数值方法一直试图将这些连续的定律转化为计算机能够理解的语言，但许多传统方法都受到对连续性的严格要求的限制，在处理复杂几何形状和诸如激波之类的尖锐、间断解时力不从心。这一局限性造成了重大的知识鸿沟，阻碍了许多关键科学和工程领域的高保真模拟。
本文介绍了间断Galerkin (DG)方法，这是一个功能强大且灵活的框架，它颠覆了传统的约束。通过拥抱间断的自由，DG方法提供了前所未有的适应性和准确性。在接下来的章节中，您将踏上一段理解这项革命性技术的旅程。首先，在“原理与机制”中，我们将解构DG方法背后的核心理念，探索它如何利用破碎函数空间和数值通量来实现灵活性和物理一致性。随后，在“应用与跨学科联系”中，我们将领略该方法在不同科学和工程学科中的成功应用，见证其卓越的多功能性。
想象一下，你想用一些直边木块来建造一个光滑、弯曲雕像的精确复制品。传统方法，很像许多经典数值方法所采用的方法，会坚持要求每一块木块都必须与相邻木块完美贴合，不能有任何间隙或悬垂。边缘必须完美对齐。这听起来很合理，但却造成了极大的约束。为了捕捉一个细微的曲线，你需要大量无穷小且切割完美的木块。如果你想只给雕像的一小部分增加更多细节，你可能需要重新切割并重新安装所有周围的木块。
间断Galerkin (DG)方法始于一个极其解放、近乎颠覆性的想法：如果我们不强迫这些木块匹配呢？如果我们允许它们是不连续的呢？
DG方法的核心是一种“分而治之”的理念。我们将问题域——无论是一个固体、一体积的流体，还是一个电磁场——划分为一系列更小、更简单的形状，我们称之为单元。在每个单元内部，我们用一个简单的函数来近似我们的解（比如温度或速度），这个函数通常是某个次数为的多项式。
革命性的一步在于我们如何处理这些单元之间的边界。与强制严格整合的连续方法不同，DG方法允许相邻单元中的多项式近似在界面上不一致。这在解的值上产生了一个“跳跃”。我们所使用的函数集合，在每个单元内部是多项式，但可以跨越边界跳跃，这些函数据说存在于一个破碎函数空间中。
我们为什么要接纳这种看似混乱的状态？答案是灵活性。巨大的、实用的灵活性。
因为我们摆脱了整合的束缚，现在我们可以在不同区域使用不同大小的单元，这种技术称为-自适应。如果流体流动在我们的域的一个小角落里有一个复杂的涡流，我们可以在那里使用非常小的单元来捕捉细节，而在其他均匀流动的区域使用大的、粗糙的单元。更妙的是，DG自然地处理了一个大单元与几个小单元相邻的情况，这会产生所谓的悬挂节点。对于连续方法来说，这些悬挂节点是一场噩梦，需要复杂的约束来将解重新粘合在一起。而对于DG来说，这根本不是问题；跳跃就是跳跃，无论邻居的大小是否相同[@problem-id:3328245]。
这种自由也延伸到了我们构建块的复杂性上。我们可以在某些区域使用简单的低次多项式，而在其他区域使用复杂的高次多项式，以在最重要的地方获得精度。这被称为-自适应。两者的结合，即-自适应，为我们提供了一个强大的工具包，可以根据我们试图解决的问题的特征精确地定制我们的近似。
但这种自由是有代价的。最初的物理定律，即偏微分方程（PDE），描述了宇宙是如何连接的。它告诉我们，一个点的温度与其邻近点的温度相关。如果我们的数值构建块都是断开的，它们怎么可能遵守一个完全关乎连通性的定律呢？单元之间如何进行交流？
在DG世界里，交流的关键在于一个优美的概念，叫做数值通量。要理解它，让我们思考一下像热方程或守恒律这样的PDE在物理上代表了什么。它通常是一个衡算原则，一个平衡的陈述。对于任何一个小体积（我们的单元），其内部某个量（如热量或动量）的变化率等于流过其边界的总量。这种跨越边界的流动被称为通量。
当我们为一个单元内的多项式近似写下这个平衡定律时，一个称为分部积分的数学过程会自然地将一个涉及单元内部导数的项转化为一个涉及单元边界面上通量的项。这就是我们面临难题的地方。在两个单元之间的界面上，我们的间断近似有两个不同的值，一个来自左边，一个来自右边。物理通量取决于解的值，那么我们应该使用这两个值中的哪一个呢？
DG的优雅答案是：我们创造一个规则。我们定义一个特殊的配方，即数值通量，它接收来自左侧的状态（）和来自右侧的状态（），并将它们组合起来，为介导它们之间交换的通量产生一个单一、明确的值。这个数值通量就是通信协议，是让那些孤立的、不连续的单元能够相互对话，并共同近似一个连续物理世界行为的通用翻译器。
最简单的DG方法，使用分片常数近似（），提供了一个极好的洞见。在这种情况下，DG公式完全简化为众所周知的有限体积法（FVM），其中数值通量是一个熟悉的概念，用于计算单元之间的流动。因此，DG可以被看作是FVM的高阶推广，它将通量交换的核心思想扩展到更复杂的、基于多项式的近似。
这个数值通量不是一个随意的规则；它是方法的核心，其设计是一门精巧的艺术。一个好的数值通量必须满足两个关键属性。
首先，它必须是相容的。这是一个简单的合理性检查：如果碰巧解在界面处没有跳跃（），那么我们的数值通量必须退化为真实的物理通量。
其次，也是更深刻的一点，通量必须确保稳定性。让我们考虑一个从左向右传播的简单波。我们应该使用什么样的通量？一个看似公平的选择是中心通量，它简单地平均了左侧和右侧的物理通量。事实证明，这是一个糟糕的想法。中心通量导致了一个完全能量守恒的格式。虽然这听起来不错，但它意味着任何小的数值误差，任何微小的摆动或振荡，都无法被阻尼掉。它将永远存在，在域中四处反弹，污染整个解。模拟是不稳定的。
现在考虑一个更聪明的选择：迎风通量。这个规则观察波的移动方向，并简单地选择来自“迎风”一侧的值。对于从左向右移动的波，它只使用来自左侧的值，。事实证明，这个选择是神奇的。通过引入轻微的偏向，它就像一点点数值“摩擦”或耗散。这种耗散正是杀死非物理、高频摆动所需要的。能量分析表明，迎风通量保证了数值解的总能量永远不会增长；它只能保持不变或减少。这是一个稳定方法的标志。
这个原则是普适的。对于像结构力学这样的椭圆问题，稳定性是通过一个“惩罚”通量来实现的，该通量惩罚跳跃，实际上就像一个硬弹簧，将不连续的解拉到一起。对于像波传播这样的双曲问题，稳定性是通过像迎风通量或Lax-Friedrichs通量这样的耗散通量来实现的，它们在需要的地方，即界面处，添加了精确控制量的耗散。DG框架的美妙之处在于，这个关键的稳定性机制完全是局部的，包含在每个界面处数值通量的设计之中。
这个由破碎空间和数值通量构成的框架可能看起来抽象而复杂。为什么要费这么多功夫？回报是巨大的，使DG成为现代计算科学与工程中最强大和最流行的方法之一。
超凡的精度： 对于具有光滑解的问题，DG方法可以达到非常高的精度阶。通过简单地增加每个单元内的多项式次数，我们可以使误差以惊人的速度减小。对于固定数量的单元，对于非常光滑（解析）的解，误差可以随呈指数级下降。这通常意味着我们可以用比低阶方法少得多的自由度获得一个高精度的答案，从而节省巨大的计算成本。当然，这种能力也是有代价的：更高阶的方法需要更小的时间步长来保持稳定，这一约束由Courant-Friedrichs-Lewy (CFL) 条件描述，对于DG，它同时依赖于单元大小和多项式次数（通常为 ）。
几何灵活性和并行性： 正如我们所见，DG对间断的容忍使其成为处理具有非协调网格的复杂几何的理想选择。同样的局部性对于并行计算来说也是一个巨大的优势。当一个问题被分配到数千个计算机处理器上时，唯一需要的通信是在分区边界交换通量数据。需要发送的数据量与子域的表面积成正比，而不是其体积。对于使用次多项式的3D模拟，每个共享面上的通信量按比例缩放，而单元内部的计算工作量则按比例缩放。这种有利的表面积与体积之比意味着DG方法在最大的超级计算机上具有出色的可扩展性。
激波和间断的王者： DG最大的成就或许是在模拟带有尖锐移动锋面的问题，例如超音速喷气机尾气中的激波。全局近似方法在遇到间断时倾向于“振铃”，将虚假的振荡（吉布斯现象）扩散到整个域。DG的局部性，结合其固有的耗散数值通量，将这种现象控制在局部。间断被清晰地捕捉为跨越单元界面的尖锐但稳定的跳跃，而不会污染流动的光滑部分。现代DG方法变得更加智能，采用了问题单元指示器，其作用类似于局部传感器。这些指示器通过例如测量多项式系数衰减的速度来检测激波的存在。在被标记为“有问题”的单元中，该方法可以自动应用更鲁棒的、起稳定作用的限制器。在“良好”的单元中，它继续发挥其全部高阶能力。这使得该方法既能做到如剃刀般锋利，又能如磐石般稳固，根据解的性质动态地调整其特性[@problem-id:3425717]。
从一个简单、反直觉的想法——间断的自由——中，涌现出一个丰富而强大的框架。通过发明一种巧妙的通信协议——数值通量——并以对稳定性和耗散的深刻理解来设计它，我们创造了一台具有卓越灵活性、准确性和鲁棒性的机器。这就是间断Galerkin方法的故事：一段从局部混乱到全局和谐的旅程。
在经历了间断Galerkin (DG)方法的原理和机制之旅后，我们现在到达了探索中最激动人心的部分：见证这些思想的实际应用。一个物理或数学原理的真正美妙之处，不在于其抽象的公式，而在于其描述世界和解决实际问题的能力。DG方法，凭借其“分而治之”的理念和巧妙的解的“弱粘合”，是计算科学家的名副其实的瑞士军刀。它不仅仅是一个单一的工具，而是一个多功能的框架，能够以惊人的优雅应对那些让其他方法陷入困境的挑战。
现在，让我们巡览一下DG方法已变得不可或缺的广阔科学与工程领域。在此过程中，我们将看到我们所学的核心概念——分片多项式、数值通量和惩罚项——是如何在宏大的应用交响曲中反复出现的主旋律。
我们的世界很少是光滑和统一的。它充满了边界、界面和突变。穿过地球的地震波会遇到不同的岩层；海洋中的声波会从海床上反射；飞机机翼由不同的复合材料制成。许多数值方法都建立在连续性的假设之上，并被迫进行复杂的扭曲来处理这些自然的划分。相比之下，DG方法正是为这个世界而生。它从一开始就假设宇宙是由碎片组成的。
想象一下，将一束声波穿过一个复杂的地质构造，也许是为了寻找石油储备或了解地震传播。地面是由具有不同密度和刚度的材料拼接而成的。当波击中两种岩石类型之间的界面时，一部分反射，一部分透射。物理学规定，虽然压力必须在界面上是连续的，但粒子的运动更为复杂：由压力梯度和局部密度加权的法向通量才是守恒的。标准的连续Galerkin方法，它在任何地方都强制解的严格连续性，虽然可以处理变化的系数，但如果网格没有完美对齐，就会遇到困难。然而，DG方法以非凡的优雅处理这种情况。由于它允许不连续性，它利用界面上的数值通量，以一种弱的、积分的方式强制执行正确的物理跳跃条件——压力的连续性和法向通量的连续性。它不需要知道单元深处的细节；它只需要管理它们边界上的对话，并且它使用物理学的语言来做到这一点。
这种“分而治之”的理念从物理上的不连续性延伸到纯粹的几何不连续性。考虑模拟地下水流经一个包含地质断层的区域。在断层的两侧创建独立、不相关的计算网格通常是自然且高效的。由此产生的网格可能是“非匹配”的，节点和单元边缘在断层线上完全错位。对于依赖于全局连续函数空间的标准有限元方法来说，这是行不通的。但对于DG来说，这只是常规操作。数值通量是一个通用翻译器，它允许单元跨越非匹配界面的通信就像跨越匹配界面一样容易。这种从“协调网格的束缚”中解脱出来的自由，在模拟从断裂储层到复杂生物组织等复杂几何形状时，是一个革命性的优势。
有时，挑战不在于材料特性，而在于解本身的行为。例如，流体动力学方程可能产生远非光滑的解。
当流体以超音速运动时，无穷小的扰动会堆积起来并变陡，形成被称为激波的近乎不连续的锋面。这就是音爆的来源。模拟激波是一个出了名的难题；不谨慎的方法会在激波周围产生剧烈的、非物理的振荡。在这里，DG方法揭示了与基本物理学的深刻联系。对于这些被称为双曲守恒律的系统，必须满足一个额外的物理原理：热力学第二定律，它规定一个称为熵的量只能被创造，不能被毁灭。事实证明，人们可以设计特殊的“熵稳定”DG格式。通过在单元边界精心构造数值通量，可以使该方法满足这个熵不等式的离散版本。这不仅仅是防止振荡；它保证了该方法是非线性稳定的，并且即使在存在强激波的情况下也能收敛到物理上正确的解 ([@problem-id:3421650])。该方法不仅仅是在解一个方程；它是在尊重一个自然的根本法则。
拥抱非光滑性的同样原则也适用于一些看似简单却富有洞察力的问题。考虑近似一个带有尖锐“尖点”的函数，比如一个金融期权的支付函数，。如果你试图用一个单一、光滑的高次多项式在整个区间上近似这个函数，你不可避免地会在尖点附近得到波纹——即臭名昭著的吉布斯现象。DG的理念提供了一个绝妙的解决方案：“如果无法平滑它，就切开它。”通过在尖点处精确地放置一个单元边界，我们可以用一个分片多项式（一侧为常数零，另一侧为直线）来完美地表示该函数。如果尖点在时变问题（如美式期权的提前行权边界）中移动，这种完美可能不会持久，但它说明了核心策略：隔离困难点，并局部攻克 ([@problem-urlem:3367334])。
这种局部松弛的思想也解决了固体力学中令人烦恼的“体积锁定”问题。在模拟像橡胶这样的近不可压缩材料时，标准的低阶有限元可能会变得病态地刚硬。近乎零体积变化的数值约束对于简单的多项式来说太强了，无法满足，因此模拟的物体拒绝变形。DG方法可以通过几种方式规避这个问题。一种方法是引入压力作为一个独立的未知数，弱地强制不可压缩性。另一种更优雅的方式是设计一个原始DG方法，其中与体积变化相关的公式部分通过特殊的稳定项进行修改。这些项恰到好处地放松了不可压缩性约束，使单元能够正确变形，同时仍然防止虚假的体积变化。这是一个DG所提供的外科手术般精度的美丽例子：它可以对物理学的特定部分进行修复，而不干扰其余部分。
除了处理困难情况外，DG框架通常为那些在其他背景下需要专门工具的问题提供了一个更统一、更优雅的方法。
考虑梁或板弯曲的物理学。其控制方程是一个四阶偏微分方程，涉及位移的四阶导数。对于标准的有限元方法来说，这是一个大难题，因为它要求基函数具有连续的一阶导数（连续性），而这 notoriously 难以实现。DG方法完全回避了这个问题。它们在每个单元上重新构建问题，并使用其强大的界面机制来连接各个部分。在单元界面上添加惩罚项，以弱形式强制挠度和转角的连续性。这些惩罚项就像数字弹簧和扭转弹簧，确保不连续碎片的组合表现得像一个连续的物理梁。曾经的主要实现障碍变成了标准DG工具集的直接应用。
这种简化的主题延伸到计算电磁学。支配从无线电波到光的一切的麦克斯韦方程组，涉及具有复杂微分结构（旋度和散度）的矢量场。用于这些问题的协调有限元方法需要高度专业化的“边元”（如Nédélec元）来正确表示物理。DG提供了一个引人注目的替代方案。人们可以在每个单元上使用标准的、简单的多项式基，并让数值通量处理跨界面场切向和法向分量之间错综复杂的耦合。事实上，更深入的分析揭示了某些DG变体（如可杂交DG）与传统协调方法之间深刻的代数联系，表明DG不仅提供了一种替代方案，而且是一个可以从中推导出其他方法的更通用的框架。
DG框架的真正威力在其应用于现代科学和工程的重大挑战中达到顶峰，这些挑战几乎总是“多物理场”问题。
想一想模拟血液在柔性动脉中的流动，或飞机机翼在湍流中的颤振。这些都是流固耦合（FSI）问题，其中可变形固体与移动流体紧密耦合。这是DG方法大放异彩的领域。我们可以使用为流体量身定制的DG离散化和为固体量身定制的另一种离散化。它们之间的界面，不断移动和变形，可以被无缝处理，因为DG天生擅长处理非匹配和移动的网格。速度连续性和牵引力平衡的物理条件通过流固界面上的通量被弱强制。然后，人们可以选择在一个“整体”方法中一次性解决整个系统，这种方法鲁棒但计算量大；或者在一个“分区”方法中分别解决流体和固体部分并进行迭代，这种方法更模块化但可能存在稳定性问题。能够构建和比较这些复杂策略本身就证明了该方法的灵活性。
最后，DG方法的变分性质可以以微妙的方式带来更高的精度。在跟踪界面运动时，如在水平集方法中，简单笛卡尔网格上的离散化可能会引入取决于网格方向的误差。例如，一个在曲率流下演化的圆可能会因为网格的各向异性而想要变成一个正方形。因为DG方法是基于积分（一种弱形式）制定的，它们比基于点的有限差分方法更有效地平均掉这些方向性偏差。这导致了更忠实于几何的模拟，减少了虚假的数值效应并抓住了细节。
从地球物理波传播的宏大尺度，到流体与结构错综复杂的舞蹈，间断Galerkin方法已证明自己不仅仅是一种数值技术。它是计算科学的一个强大范式——一种拥抱复杂性、在离散层面尊重物理规律、并为解决重要问题提供统一、灵活、优雅框架的思维方式。