[論文レビュー] FullSWOF: A free software package for the simulation of shallow water flows

Olivier Delestre, Frédéric Darboux+4 more|arXiv (Cornell University)|Jan 16, 2014

Meteorological Phenomena and Simulations参考文献 22被引用数 13

ひとこと要約

FullSWOF は、浅水域の流れをシミュレートするためのオープンソース C++ ソフトウェアパッケージであり、良好にバランスの取れた有限体積法を用いて、小さな水深、湿乾遷移、降雨、浸透、デジタル高程モデル（DEMs）を特に扱うように設計されている。1次元および2次元のベンチマーク、ならびに実際の現場データにおいて、オープンソースコードとバージョン管理を用いた検証済みの正確性と再現可能性を示しており、強固な性能を発揮している。

ABSTRACT

Numerical simulations of flows are required for numerous applications, and are usually carried out using shallow water equations. We describe the FullSWOF software which is based on up-to-date finite volume methods and well-balanced schemes to solve this kind of equations. It consists of a set of open source C++ codes, freely available to the community, easy to use, and open for further development. Several features make FullSWOF particularly suitable for applications in hydrology: small water heights and wet-dry transitions are robustly handled, rainfall and infiltration are incorporated, and data from grid-based digital topographies can be used directly. A detailed mathematical description is given here, and the capabilities of FullSWOF are illustrated based on analytic solutions and datasets of real cases. The codes, available in 1D and 2D versions, have been validated on a large set of benchmark cases, which are available together with the download information and documentation at http://www.univ-orleans.fr/mapmo/soft/FullSWOF/.

研究の動機と目的

• 水文学および地球物理モデリング分野において、再現可能でオープンソースのツールが不足しているという問題に対処すること。
• 安定状態を保持し、湿乾境界における水深の正負を維持する、強固で良好にバランスの取れた有限体積ソルバーを開発すること。
• デジタル高程モデル（DEMs）を用いて、降雨、浸透、変化する底面地形などの複雑な水文的プロセスを正確にシミュレートすること。
• バージョン管理とレグRESSIONテストを併用したオープンライセンスでの完全なソースコード公開により、科学的再現性を確保すること。
• 水文学、洪水モデリング、環境工学分野の研究者および実務家が使いやすく、拡張可能なフレームワークを提供すること。

提案手法

• 保存形の浅水域方程式に基づく有限体積法を実装し、リーマンソルバを用いたゴドノフ型スキームを採用する。
• 安定状態を保持し、湿乾界面における水深の正負を保証するため、良好にバランスの取れた再構成技術を統合する。
• 降雨、浸透、底面摩擦（例：マンニンの法則）の源項を数値フラックスに直接統合する。
• 構造的直交グリッドを用いた1次元および2次元のシミュレーションを実施し、CFL 条件に基づく適応的時間刻みを採用する。
• コード品質と拡張性を確保するため、バージョン管理（Git）、バグトラッキング、継続的インテグレーションを含むオープンソースソフトウェア開発手法を採用する。
• MPI を用いた並列化された2次元バージョンを開発し、マルドゥー・ダム破壊事例のような大規模なシミュレーションを高速化する。

実験結果

リサーチクエスチョン

• RQ1 良好にバランスの取れたオープンソース有限体積スキームは、湿乾遷移および小さな水深を伴う浅水域流れを正確にシミュレートできるか？
• RQ2 FullSWOF は、セネガルのティエスで実施されたフィールドプロットにおける地上流出を、実際の水文的イベントとして正確にシミュレートできるか？
• RQ3 FullSWOF は、マルドゥー事例のような大規模なダム破壊シナリオにおいて、複雑な流れのダイナミクスをどの程度正確に再現できるか？
• RQ4 急激な勾配や乾燥状態が存在する状況でも、データフィルタリングを一切行わずに、数値的安定性と正確性を維持できるか？
• RQ5 オープンソースモデルは、水文モデリング分野における再現可能性とコミュニティ主導の開発をどの程度効果的に促進できるか？

主な発見

• FullSWOF_2D は、セネガルのティエスで実施された 400 m² の砂質土壌フィールドに 70 mm/h の降雨強度を適用した地上流出を、データフィルタリングなしで、先行研究と定性的に整合する結果を示している。
• ソフトウェアは、マルドゥー・ダム破壊ベンチマークの9か所の水位計で観測された最大水位を、スケールド物理モデルの実験データと一致させる形で正確に捉えている。
• 並列化された FullSWOF_2D バージョンは、17.3 km × 8.4 km の大規模領域（486,000 セル）に対して数時間のシミュレーション時間を短縮し、大規模なシミュレーションを効率的に行えるようにしている。
• FullSWOF_1D および FullSWOF_2D は、解析解、現場データ、複雑な地形を含む多様なテストケースにおいて、安定的かつ強固に動作している。
• ソフトウェアは、校正なしに浸透（θs−θi=0.12）および表面摩擦（n=0.033 m⁻¹/³s）を正しく扱っており、物理的整合性を示している。
• バージョン管理とレグRESSIONテストを併用したオープンソースリリースにより、再現性が確保され、コミュニティの貢献が可能となり、長期的な保守性と科学的透明性が向上している。