[论文解读] Modeling the cryosphere with FEniCS

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

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

一句话总结

本文提出了一种基于FEniCS的有限元框架，用于模拟冰冻圈过程，重点研究冰盖动力学和内部能量分布。该研究引入了CSLVR这一开源软件工具，通过稳定化有限元方法和PDE约束优化，实现对冰流、底部融化的高保真模拟，以及水含量分布的精确建模，并在南极和格陵兰冰盖几何结构上验证了结果的有效性。

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，简化冰冻圈建模中的网格生成、数据转换和可视化流程。

提出的方法

• 利用FEniCS有限元框架，通过Galerkin变分格式求解控制冰流和能量平衡的PDE。
• 采用Dolfin-Adjoint和IPOPT进行PDE约束优化，实现从观测数据反演速度场和能量场。
• 应用如泡函数增强和亚网格尺度模型（SSM）等稳定化技术，以处理对流主导和不稳定的数值问题。
• 为Stokes方程和全Stokes方程实现变分格式，边界条件包括无滑移和滑移-摩擦边界条件。
• 使用GMSH进行2D/3D网格生成，使用Paraview进行3D模拟结果的后处理与可视化。
• 提出一种新颖的优化流程，通过求解具有时间依赖水生成速率的能量平衡方程，实现对观测冰内水含量的匹配。

实验结果

研究问题

• RQ1 如何有效应用FEniCS中的有限元方法，对具有复杂几何形状的冰盖动力学进行建模？
• RQ2 对于数值不稳定的冰冻圈PDE（如对流-扩散-反应系统），哪些稳定化技术最为有效？
• RQ3 PDE约束优化框架能否准确反演由观测底部融化率推导出的内部能量和水含量分布？
• RQ4 速度数据的引入如何提升热-机械冰盖模拟的精度？
• RQ5 非连续热导率对暖冰中内部能量和水含量分布有何影响？

主要发现

• CSLVR软件成功实现了基于格陵兰和南极真实冰盖几何形状以及理想化几何形状的2D和3D冰流模拟。
• 稳定化有限元方法，特别是GLS和泡函数增强格式，在对流主导的能量与质量平衡问题中显著提升了数值稳定性。
• 经优化的能量平衡模型通过将水含量的时间导数与热通量和融化率相匹配，准确再现了观测到的冰内水含量。
• 当通过优化校准水生成速率时，由能量平衡模型推导出的底部融化率与观测数据高度一致。
• 对动量场和能量场的双重优化流程提升了模拟保真度，经L曲线分析和ISMIP-HOM基准测试验证。
• 基于质量跳跃关系和重心速度推导出的水通量边界条件，确保了模型中水含量、通量与底部融化的自洽性。