Barotropic-Baroclinic Splitting for Multilayer Shallow Water Models with Exchanges

技术摘要：具有交换机制的多层浅水模型的条纹-斜压分裂技术

问题陈述
全球海洋环流模型（OGCMs）面临着严重的计算约束：其时间步长受限于最快波速，即表面重力波（约 200 m/s），尽管主要研究现象（层结和洋流）是在更慢的时间尺度上（约 3 m/s）演化的。虽然此前已采用刚性盖模型（rigid lid approximation）或隐式表面高度法等方法，但这些方法存在局限性。条纹-斜压分解（barotropic-baroclinic decomposition）提供了一种解决方案，它将快速的表面重力波（条纹部分）与较慢的垂直调整（斜压部分）分离。然而，现有的分裂方法在处理具有质量交换的多层框架时，往往难以保持一致性、维持能量特性，并确保在不损失精度的情况下保证稳定性，特别是在地形跟随坐标系中。

方法论
作者对一种在具有层间质量交换的非线性多层框架内的条纹-斜压分裂进行了数值分析，该方法被表述为一种精确算子分裂。模型基于地形跟随（$\sigma$）坐标下的无粘欧拉方程，并假设密度恒定。

核心方法论涉及将层平均欧拉方程分解为两个不同的算子：

1. 条纹步（Barotropic Step）： 该步骤处理自由表面高度（$h$）和深度平均速度（$\bar{u}$）的演化。它被构建为一个包含压力项和地形源项的良好平衡（well-balanced）单层浅水模型，但不包含斜压强迫。至关重要的是，该步骤在平均速度 $\bar{u}$ 下处理各层速度（$u_\alpha$）和示踪物（$T_\alpha$）的输运。
2. 斜压步（Baroclinic Step）： 该步骤管理层间的垂直交换以及相对于平均速度的速度调整。它求解偏离平均速度的部分（$u_\alpha - \bar{u}$），并通过拉格朗日乘子（交换项 $G_{\alpha \pm 1/2}$）强制执行层高度与总高度成比例的约束（即 $h_\alpha = l_\alpha h$）。该系统是完全三维的，但缺乏压力和地形项，因此允许使用较大的时间步长。

数值方案采用了带有子循环策略的 Lie 分裂法：

• 斜压步在较大的时间步长 $\Delta t$ 上求解。
• 条纹步在相同的间隔内使用 $K$ 个受重力波速约束的小步长 $\delta t$ 进行求解。
• 一种修正机制确保了条纹步预测的平均速度与斜压步更新后的各层速度的平均值相一致。

主要贡献

• 精确算子分裂公式： 论文将分裂重新表述为精确算子和（$M_L = BT + BC$），明确了带有交换机制的多层系统的数学结构。
• 物理特性的保持： 作者证明了该分裂方法保持了总能量守恒（或耗散），并满足离散熵不等式。具体而言，条纹步满足能量等式，而斜压步由于交换项的耗散性质满足能量不等式。
• 离散最大值原理： 研究表明该方案满足速度和示踪物浓度的离散最大值原理，确保数值保持在物理合理的范围内。
• 良好平衡性（Well-Balancing）： 通过在条纹步中使用良好平衡求解器（静力学重建），整体方案保持了“静止湖泊”（lake-at-rest）平衡；当扩展应用最新的线性良好平衡策略时，该方案在粗网格上也能维持地转平衡。
• 双曲性分析： 论文研究了分裂系统的双曲性，证明了条纹系统是双曲的，并分析了斜压系统的特征值，指出在关于速度偏差的特定条件下，该系统保持严格双曲性。

结果
数值实验证明了所提方案的有效性：

• 计算效率： 在低弗劳德数（Froude number）模拟中（即流速与重力波速之比很小时），分裂法相比于非分裂方案大幅降低了计算成本。其增益与弗劳德数成正比。
• 精度： 尽管采用了分裂和子循环技术，该方法仍能保持与非分裂方案相当的精度。针对静力学欧拉方程解析解的收敛性研究显示其具有一阶收敛性。
• 扩散特性： 在某些测试案例（如风驱动腔体）中，发现分裂法比非分裂方案的扩散性更低；但在其他案例（如静态解）中，它表现出略高的扩散性。
• 地转平衡： 在条纹步中应用线性良好平衡策略，显著改善了长期积分中西边界流的解析能力，与标准方案相比，能在粗网格上捕捉到更窄、更快的流。

意义与主张
论文声称，这项工作提供了一种既具有计算效率又符合物理一致性的条纹-斜压分裂的严谨数值分析。其主要意义在于能够利用较大的时间步长（受限于缓慢的斜压波而非快速的重力波）来模拟长期海洋环流，且不会牺牲能量守恒或关键平衡态（静止湖泊和地转平衡）的保持。作者强调，用于条纹步的良好平衡策略被整体方案继承，从而允许在粗网格上对地转流进行精确模拟。该方法被呈现为一个稳健的框架，适用于以自由表面波并非主要关注点的全球及区域应用场景。