Barotropic-Baroclinic Splitting for Multilayer Shallow Water Models with Exchanges

技術要約：交換を伴う多層浅水モデルのためのバロトロピック・バロクリニック分割

問題提起
全球海洋循環モデル（OGCM）は、深刻な計算上の制約に直面している。すなわち、タイムステップが最も速い波の速度、具体的には表面重力波（約200 m/s）によって制限されることである。これに対し、主要な関心対象である現象（成層や流速）は、はるかに遅い時間スケール（約3 m/s）で進化する。剛体蓋近似や、暗黙的な表面高度法などの従来のアプローチも用いられてきたが、それらには限界がある。バロトロピック・バロクリニック分解は、速い表面重力波（バロトロピック）と遅い垂直調整（バロクリニック）を分離することで、この問題に対する解決策を提供する。しかし、既存の分割手法は、質量交換を伴う多層フレームワークにおいて、一貫性の維持、エネルギー特性の保存、および精度を損なうことなく安定性を確保することに苦慮する場合が多い。特に地形追従座標系においてはその傾向が顕著である。

手法
著者らは、質量交換を伴う非線形多層フレームワーク内でのバロトロピック・バロクリニック分割の数値解析を提案する。これは、定数密度を仮定した、地形追従（$\sigma$）座標における非粘性オイラー方程式として定式化された厳密な演算子分割である。

核となる手法は、層平均オイラー方程式を2つの異なる演算子に分解することである。

1. バロトロピック・ステップ： このステップは、自由表面高度（$h$）および層平均速度（$\bar{u}$）の進化を扱う。これは、圧力項と地形ソース項を含む、ウェル・バランストな1層浅水モデルとして定式化されるが、バロクリニックな強制力は含まれない。決定的なのは、このステップが個々の層の速度（$u_\alpha$）およびトレーサー（$T_\alpha$）の輸送を、平均速度（$\bar{u}$）において扱う点である。
2. バロクリニック・ステップ： このステップは、層間の垂直方向の交換、および平均速度に対する速度の偏差の調整を管理する。これは、平均速度からの偏差（$u_\alpha - \bar{u}$）を解き、層の高さが全高に比例するという制約（$h_\alpha = l_\alpha h$）をラグランジュ乗数（交換項 $G_{\alpha \pm 1/2}$）を介して強制する。このシステムは完全な3次元であるが、圧力項と地形項を欠いているため、より大きなタイムステップが可能となる。

数値スキームは、サブサイクリング戦略を用いたLie分割を採用している。

• バロクリニック・ステップは、大きなタイムステップ $\Delta t$ で解かれる。
• バロトロピック・ステップは、重力波の速度によって制約される $K$ 個の小さなサブステップ $\delta t$ を用いて、同じ間隔内で解かれる。
• 補正メカニズムにより、バロトロピック・ステップによって予測された平均速度が、バロクリニック・ステップによって更新された個々の速度の平均と一致することを保証する。

主な貢献

• 厳密な演算子分割の定式化： 本論文は、交換を伴う多層システムの数学的構造を明確にするために、分割を厳密な演算子の和（$M_L = BT + BC$）として再定式化した。
• 物理的特性の保存： 著者らは、この分割が全エネルギーの保存（または散逸）を保持し、離散エントロピー不等式を満たすことを証明した。具体的には、バロトロピック・ステップはエネルギー等式を満たし、バロクリニック・ステップは交換項の散逸的な性質によりエネルギー不等式を満たす。
• 離散最大原理： 本スキームは、速度およびトレーサー濃度に関する離散最大原理を満たすことが示されており、値が物理的な範囲内に留まることを保証している。
• ウェル・バランシング： バロトロピック・ステップにウェル・バランスト・ソルバー（静水圧再構成）を利用することで、全体的なスキームは「湖の静止状態（lake-at-rest）」の平衡を保持し、また、最新の線形ウェル・バランシング戦略を適用することで、粗いメッシュ上でも地衡平衡を維持する。
• 双曲性解析： 本論文は分割された系の双曲性を調査し、バロトロピック系が双曲的であることを証明し、速度偏差に関する特定の条件下でバロクリニック系の固有値が厳密に双曲性を維持することを分析している。

結果
数値実験は、提案されたスキームの有効性を実証している。

• 計算効率： 低フルード数シミュレーション（流速と重力波速度の比が小さい場合）において、分割法は分割なしのスキームと比較して計算コストを劇的に削減する。この利得はフルード数に比例する。
  تعتبر。
• 精度： 分割およびサブサイクリングにもかかわらず、本手法は分割なしのスキームと同等の精度を維持する。静水圧オイラー方程式の解析解に対する収束研究は、1次収束を示している。
• 拡散特性： 特定のテストケース（風駆動のキャビティなど）では、分割法は分割なしのスキームよりも拡散が少ないことが判明したが、他のケース（定常解）では、わずかに高い拡散を示すこともあった。
• 地衡平衡： バロトロピック・ステップにおける線形ウェル・バランシング戦略の適用は、長期間の積分において西岸境界流の解像度を大幅に向上させ、標準的なスキームと比較して、粗いメッシュ上でより狭く速い流を捉えることができる。

意義と主張
本論文は、計算効率が高く、かつ物理的に一貫した、バロトロピック・バロクリニック分割の厳密な数値解析を提供することを主張している。主な意義は、重力波（速い波）ではなく、遅いバロクリニック波によって制約される大きなタイムステップを用いても、エネルギーの保存や重要な平衡（湖の静止状態および地衡平衡）の維持を犠牲にすることなく、長期的な海洋循環をシミュレートできる能力にある。著者らは、バロトロピック・ステップで使用されるウェル・バランシング戦略が全体のスキームに継承されるため、粗い格子上で地衡流を正確にシミュレートできることを強調している。本手法は、自由表面波が主要な関心の対象ではない、全球および地域的な応用における堅牢なフレームワークとして提示されている。