[論文レビュー] Modeling the cryosphere with FEniCS

Evan Cummings|arXiv (Cornell University)|Sep 6, 2016

Cryospheric studies and observations参考文献 21被引用数 4

ひとこと要約

本稿では、FEniCSを用いた有限要素フレームワークを提示し、氷河力学および内部エネルギー分布のモデル化に焦点を当てている。本研究は、安定化された有限要素法とPDE制約最適化を用いて、氷の流れ、底面融解、水分含量分布の高精度なシミュレーションを可能にするオープンソースソフトウェアツールCSLVRを紹介している。結果は、南極およびグリーンランド氷床の幾何形状で妥当性が検証されている。

ABSTRACT

This manuscript is a collection of problems and solutions related to modeling the cryosphere using the finite element software FEniCS. Included is an introduction to the finite element method; solutions to a variety of problems in one, two, and three dimensions; an overview of popular stabilization techniques for numerically-unstable problems; and an introduction to the governing equations of ice-sheet dynamics with associated FEniCS implementations. The software developed for this project, Cryospheric Problem Solver (CSLVR), is fully open-source and has been designed with the goal of simplifying many common tasks associated with modeling the cryosphere. CSLVR possesses the ability to download popular geological and geographical data, easily convert between geographical projections, develop sophisticated two- or three-dimensional finite-element meshes, convert data between many popular formats, and produce production-quality images of data. Scripts are presented which model the flow of ice using geometry defined by mathematical functions and observed Antarctic and Greenland ice-sheets data. A new way of solving the internal energy distribution of ice to match observed intra-ice water contents within temperate regions is thoroughly explained.

研究の動機と目的

• 複雑な氷圏システム、特に氷床力学と熱力学的結合をモデル化するための堅牢でオープンソースの計算フレームワークの開発。
• サブスケールモデルやGLS法を含む高度な安定化技術を用いて、氷圏シミュレーションにおける数値的不安定性を解消すること。
• エネルギー収支方程式の最適化により、温暖帯氷領域における内部エネルギー分布および水分含量の正確なシミュレーションを可能にすること。
• 観測データ（例：グリーンランドおよび南極からのデータ）を有限要素メッシュに統合し、現実的なシミュレーションを実現すること。
• メッシュ生成、データ変換、可視化を単一のツールであるCSLVRで簡素化する統合型ツールの提供。

提案手法

• ガラーキン変分式を用いて、氷の流れとエネルギー収支を記述するPDEを解くためにFEniCSの有限要素フレームワークを活用する。
• PDE制約最適化のためDolfin-AdjointとIPOPTを用い、観測データから速度場およびエネルギー場を逆算可能にする。
• 顕著な対流優勢および不安定な問題に対処するため、バブル拡張およびサブスケールモデル（SSM）などの安定化技術を適用する。
• スティーブス方程式およびフル・スティーブス方程式の変分形式を、滑らかでない境界条件および滑り摩擦境界条件を含めて実装する。
• 2次元および3次元メッシュ生成にはGMSHを、3次元シミュレーション出力の後処理および可視化にはParaviewを用いる。
• 時間依存の水分生成率を伴うエネルギー収支方程式を解くことで、観測された内部氷内水分含量に一致する新しい最適化手順を導入する。

実験結果

リサーチクエスチョン

• RQ1 FEniCSにおける有限要素法は、複雑な幾何形状を有する氷床の動態を効果的にモデル化するためにどのように適用可能か？
• RQ2 対流優勢で不安定な氷圏PDE（例：対流拡散反応系）に対して、最も効果的な安定化技術は何か？
• RQ3 PDE制約最適化フレームワークは、観測された底面融解率から内部エネルギーおよび水分含量分布を正確に推定できるか？
• RQ4 速度データの組み込みは、熱力学的結合を伴う氷床シミュレーションの精度をどのように向上させるか？
• RQ5 不連続な熱伝導率が、温暖帯氷内における内部エネルギーおよび水分含量の分布に与える影響は何か？

主な発見

• CSLVRソフトウェアは、グリーンランドおよび南極の実際の氷床幾何形状に加え、理想化された幾何形状を用いた2次元および3次元の氷の流れを成功裏にシミュレートした。
• 特にGLS法およびバブル拡張形式を用いた安定化された有限要素法は、対流優勢なエネルギー収支および質量収支問題における数値的安定性を顕著に向上させた。
• 最適化されたエネルギー収支モデルは、水分含量の時間微分を熱フラックスおよび融解率に一致させることで、観測された内部氷内水分含量を正確に再現した。
• 最適化により水生成率をキャリブレーションした場合、エネルギー収支モデルから導かれた底面融解率は、観測データと強く一致した。
• 運動量場およびエネルギー場の二重最適化手順により、シミュレーションの忠実性が向上し、L曲線解析およびISMIP-HOMベンチマーク比較によって確認された。
• 質量ジャンプ関係および重心速度に基づく導出された水分フラックス境界条件は、モデル内での水分含量、フラックス、底面融解の整合性を保証した。