[論文レビュー] A conservative, discontinuous Galerkin, tracer transport scheme using compatible finite elements

Timothy C. Andrews, Thomas M. Bendall|arXiv (Cornell University)|Mar 19, 2026

Advanced Numerical Methods in Computational Mathematics被引用数 0

ひとこと要約

論文は、互換有限要素フレームワーク内で質量を守る保守的トレーサー輸送スキームを提示し、DG離散化とトレーサー密度を輸送するための特別な空間写像を用いて一貫性と非負性を保持します。視点は輸送と比較して、保守性の向上を実証するテスト設定で示されます。

ABSTRACT

This paper outlines a conservative transport scheme for scalar tracers within a compatible finite element model for geophysical fluid equations. Instead of using the advective transport equation for a mixing ratio, a conservative transport equation is solved for the tracer density of the mixing ratio multiplied by the dry density. This ensures mass conservation in the continuous equations, which can be preserved in the discrete equations with a discontinuous Galerkin transport scheme. Our method is designed to work for two placements of the mixing ratio in a Charney-Phillips vertical staggering: either co-located with the dry density or vertically staggered from it. The new scheme is designed to conserve the tracer density and ensure consistency by maintaining a constant mixing ratio. Additionally, a mass-conserving limiter is developed to ensure non-negativity in the co-located configuration. Tests with terminator toy chemistry and a moist rising bubble show the use of the new transport scheme with physics terms and its ability to accurately model mass conservation of moisture species in a dynamical core setup.

研究の動機と目的

• 地球物理流体モデルにおける質量を守るトレーサー輸送の必要性を動機付ける。
• グローバル質量を保存するトレーサー密度を用いた保守的輸送形式を提案する。
• 関数空間を跨いで密度と混合比を結ぶ空間写像演算子（保守的Galerkin射影、回復、埋め込みDG）を開発する。
• Charney-Phillips staggered フレームワーク内で一貫性と非負性を確保する。

提案手法

• トレーサーを密度rho_X = rho mとして表し、rho mの保守的輸送方程式を導出する。
• 乾燥密度とトレーサー密度輸送の両方に対して不連続Galerkin弱形式を用い、局所およびグローバルな保守性を保証する。
• 密度と混合比を異なる関数空間に配置できるようにし、保守的Galerkin射影、回復、埋め込みDG写像を適用して共通空間で輸送する。
• 次に低次の空間で2次精度を実現するための回復スキームを導入し、回復空間J = I R + I(1 - P~R)を用いる。
• 同時配置設定で非負性を強制する平均混合比リミッターを導入する。

実験結果

リサーチクエスチョン

• RQ1 rho mの保守的トレーサー輸送方程式を解くことで輸送中にトレーサー質量がグローバルに保存されるか？
• RQ2 空間的に一定の混合比が輸送下で変化せず一貫性を保てるか？
• RQ3 共存配置と鉛直的な staggered配置のような異なる関数空間配置が質量保守性と一貫性にどう影響するか？
• RQ4 回復および埋め込みDG写像が最も低次の場の精度に及ぼす影響はどの程度か？
• RQ5 質量保全を崩さず非負性リミッターを実装できるか？

主な発見

• 保守的トレーサー輸送スキームはテストでトレーサー密度（rho m）を機械精度近く保存する。
• 回復スキームを介して次に低い階の空間で輸送する場合、k=0要素の2次精度を達成する。
• k=1空間では元の空間での輸送が2次精度を超える収束を示す。
• 教示的輸送スキームは保守的トレーサー scheme に比べてテストで質量のずれが大きいことを示す。
• 同時配置設定で輸送下の空間的に一定の混合比を維持することで一貫性を示す。