观潮 · 模式预报与气象 AI 观察
熵稳定间断伽辽金方法：让数值模式同时守恒熵和总能量
数值模式的核心挑战之一是在离散化过程中保持物理守恒律。对于可压缩欧拉方程（描述大气和海洋运动的基本方程组），传统的数值方法往往只能守恒质量、动量和总能量，但无法同时守恒热力学熵——而熵守恒是保证数值解物理合理性的关键，特别是在长时间积分中。
这篇论文在间断伽辽金谱元方法（DGSEM）框架下，提出了一种新的数值通量构造方法，能够同时守恒热力学熵和总能量，并且对等温和等熵背景态是 well-balanced（均衡保持）的。
问题：为什么熵守恒这么难
可压缩欧拉方程在矢量不变形式下包含非守恒项（nonconservative terms），这些项在连续情况下不影响守恒律，但在离散化后如果处理不当会引入虚假的熵产生或耗散。
传统 DGSEM 方法使用 Lax-Friedrichs 或 Roe 等数值通量时，虽然能守恒总能量，但会引入数值耗散导致熵不守恒。这在短时间模拟中影响不大，但在气候尺度或长时间天气预报中，虚假的熵耗散会累积成系统性偏差。
作者要解决的问题是：如何构造一种数值通量，在离散化后同时严格守恒热力学熵和总能量？
方法：对称-反对称分解 + 两点通量
核心思路是将可压缩欧拉方程的非守恒项分解为对称和反对称两部分，然后分别处理：
- 对称部分：用中心通量处理，保持守恒性。
- 反对称部分：用迎风通量处理，引入必要的数值耗散以维持稳定性。
通过这种分解，作者推导出一种新的两点数值通量（two-point numerical flux），该通量满足：
- 熵守恒：在离散层面上严格守恒热力学熵。
- 总能量守恒：同时守恒总能量（动能 + 内能 + 势能）。
- Well-balanced：对等温（isothermal）和等熵（isentropic）背景态保持平衡，即如果初始状态是平衡态，数值解不会虚假地偏离。
此外，该方法还保证了密度和压力的正性（positivity），这是物理合理性的基本要求。
验证：从惯性重力波到斜压不稳定
作者通过一系列数值实验验证了方法的性能：
实验类型 · 验证内容 · 结果
熵和总能量守恒 · 长时间积分中熵和总能量的漂移 · 机器精度守恒
Well-balancedness · 等温/等熵背景态的保持 · 无数值漂移
惯性重力波 · 波动传播的准确性 · 与解析解一致
含地形结果 · 复杂地形下的稳定性 · 稳定且物理合理
斜压不稳定 · 天气尺度涡旋的发展 · 正确捕捉不稳定增长
特别值得注意的是斜压不稳定实验：这是大气中天气尺度涡旋（如温带气旋）形成的基本机制。新方法能够正确模拟这一过程，说明其不仅守恒律严格，而且动力学行为也是物理合理的。
与现有方法的关系
该工作属于 structure-preserving methods（结构保持方法）范畴，与以下研究方向密切相关：
- 熵守恒有限体积法：如 Tadmor (2003) 的熵守恒通量构造。
- Well-balanced DG 方法：如 Botta et al. (2004) 对浅水方程的工作。
- 矢量不变形式：如 Gassner et al. (2016) 对可压缩欧拉方程的矢量不变离散化。
本文的创新在于将熵守恒和 well-balancedness 统一在 DGSEM 框架下，并扩展到曲线网格（curvilinear meshes），使其适用于实际地球系统模式中的复杂几何。
作者承认的限制
该方法目前仅针对可压缩欧拉方程，尚未扩展到包含物理过程（如辐射、微物理、湍流）的完整地球系统模式。 作者指出，将熵守恒 DG 方法耦合到物理参数化方案是一个重要的未来方向。
Well-balancedness 仅对等温和等熵背景态验证。 实际大气中的背景态更加复杂，包含温度梯度和湿度变化，需要进一步研究。
计算成本高于传统 DGSEM 方法。 两点通量的计算比 Lax-Friedrichs 通量更昂贵，特别是在高分辨率模拟中。
一句话判断
对于需要长时间积分的气候模拟和天气预报，熵稳定 DG 方法提供了一种严格守恒热力学熵和总能量的数值框架，是结构保持方法在地球系统模式中的重要进展。
图 1：惯性重力波临界层附近垂直速度 w 在 4.5–12 小时的演变（80×180 网格，3 阶精度）。初始阶段（上排）波动结构清晰规整；随时间推移（下排）波动在临界层附近破碎，产生小尺度湍流结构——熵稳定格式在这一过程中保持了密度和压力的正性，没有出现虚假的数值振荡。