[论文解读] 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）进行了设计。它在 1D 和 2D 基准测试中表现出稳健性能，包括溃坝事件和真实世界现场数据，通过开放源码和版本控制实现了验证过的准确性与可复现性。

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）实现对降雨、入渗和变化的底地形等复杂水文过程的精确模拟。
• 通过在开放许可证下发布完整源代码并配合版本控制与回归测试，确保科学研究的可复现性。
• 为水文学、洪水模拟及环境工程领域的研究人员和实践者提供一个用户友好且可扩展的框架。

提出的方法

• 基于守恒形式的浅水方程实现有限体积法，采用 Godunov 型格式与 Riemann 求解器。
• 集成满足平衡的重构技术，以保持稳态并确保在干湿界面处水深的正性。
• 将降雨、入渗和底摩擦（如 Manning 公式）的源项直接整合到数值通量中。
• 支持基于 CFL 条件的结构化笛卡尔网格上的 1D 和 2D 模拟，以及自适应时间推进。
• 采用开源软件工程实践，包括版本控制（Git）、错误追踪与持续集成，以保障代码质量与可扩展性。
• 使用 MPI 实现并行化的 2D 版本，以加速大规模模拟，例如 Malpasset 溃坝案例。

实验结果

研究问题

• RQ1 一个满足平衡的开源有限体积格式能否准确模拟具有干湿边界转换和小水深的浅水流动？
• RQ2 FullSWOF 在模拟塞内加尔蒂耶斯地区田间地块的地表径流等真实水文事件时表现如何？
• RQ3 该软件在大规模溃坝情景（如 Malpasset 事件）中，能否准确再现复杂的流动动力学？
• RQ4 即使在存在陡峭梯度和干态的情况下，该软件是否能保持数值稳定性和准确性，而无需数据滤波？
• RQ5 该开源模型在促进可复现性与社区驱动开发方面，在水文建模中的有效性如何？

主要发现

• FullSWOF_2D 在塞内加尔蒂耶斯地区一块 400 m² 的沙质地块上成功模拟了 70 mm/h 降雨强度下的地表径流，结果在定性上与先前研究一致，且无需数据滤波。
• 该软件准确捕捉了 Malpasset 溃坝基准测试中九个测站的最大水位高度，与缩尺物理模型的实验数据吻合。
• 并行化的 FullSWOF_2D 版本将大区域（17.3 km × 8.4 km，共 486,000 个单元）的模拟时间缩短至数小时，实现了高效的大规模模拟。
• FullSWOF_1D 和 FullSWOF_2D 在多种测试案例中均表现出稳定与鲁棒性，包括解析解、现场数据与复杂地形。
• 该软件在无需校准的情况下正确处理了入渗（θs−θi=0.12）与地表摩擦（n=0.033 m⁻¹/³s），展示了良好的物理一致性。
• 通过版本控制与回归测试的开源发布确保了可复现性，并支持社区贡献，从而提升了长期可维护性与科学透明度。