[论文解读] FreeFem++, a tool to solve PDEs numerically

Georges Sadaka|arXiv (Cornell University)|May 7, 2012

Advanced Numerical Methods in Computational Mathematics参考文献 1被引用 6

一句话总结

FreeFem++ 是一款开源、用户友好的有限元软件，通过变分格式数值求解偏微分方程（PDEs）。它支持二维和三维问题，具备基于各向异性适应的自动网格生成，并支持多种有限元类型（包括不连续和混合元），可通过高级脚本语言和内置求解器（如 UMFPACK 和 GMRES）高效求解线性和非线性椭圆、抛物线及双曲型 PDEs。

ABSTRACT

FreeFem++ is an open source platform to solve partial differential equations numerically, based on finite element methods. It was developed at the Laboratoire Jacques-Louis Lions, Universit\\'e Pierre et Marie Curie, Paris by Fr\\'ed\\'eric Hecht in collaboration with Olivier Pironneau, Jacques Morice, Antoine Le Hyaric and Kohji Ohtsuka. The FreeFem++ platform has been developed to facilitate teaching and basic research through prototyping. FreeFem++ has an advanced automatic mesh generator, capable of a posteriori mesh adaptation; it has a general purpose elliptic solver interfaced with fast algorithms such as the multi-frontal method UMFPACK, SuperLU . Hyperbolic and parabolic problems are solved by iterative algorithms prescribed by the user with the high level language of FreeFem++. It has several triangular finite elements, including discontinuous elements. For the moment this platform is restricted to the numerical simulations of problems which admit a variational formulation. We will give in the sequel an introduction to FreeFem++ which include the basic of this software. You may find more information throw this link http://www.freefem.org/ff++.

研究动机与目标

• 为科学计算中的教学和原型设计提供一种易于访问的开源数值 PDE 求解工具。
• 使研究人员和学生能够无需低级编程即可快速实现和测试有限元方法。
• 支持通过变分格式求解复杂 PDE 系统，包括非线性和时间依赖问题。
• 集成先进的数值特性，如自适应网格、质量矩阵对角化及高效的稀疏求解器。
• 促进 PDE 领域的数值实验，用于模型验证、数学分析和算法开发。

提出的方法

• 使用高级脚本语言通过其变分形式定义 PDE，与数学符号高度相似。
• 采用基于 Delaunay-Voronoi 的自动网格生成器，并结合后验各向异性适应以提高精度。
• 支持多种有限元类型，包括连续和不连续的拉格朗日元（P1, P2, P1dc, P2dc）以及 Raviart-Thomas 元。
• 集成快速稀疏求解器（如 UMFPACK、SuperLU、CG 和 GMRES），以高效求解线性系统。
• 通过用户定义的迭代格式实现时间积分，包括对抛物型问题具有稳定性分析的 θ-格式。
• 提供网格到网格插值、误差估计及可视化工具（如 gnuplot 等外部工具）

实验结果

研究问题

• RQ1 一个高级、开源的平台在多大程度上可以简化研究人员和教育工作者对 PDE 数值求解的流程？
• RQ2 当在 FreeFem++ 中使用 θ-格式求解热方程时，其在空间和时间上的收敛行为如何？
• RQ3 自动网格适应与质量矩阵对角化在提升时间依赖型 PDE 求解的稳定性和效率方面效果如何？
• RQ4 FreeFem++ 在通过变分格式求解包含非线性和耦合项的复杂 PDE 系统方面，其支持程度如何？
• RQ5 FreeFem++ 在求解泊松方程和非线性椭圆方程等基准问题时，其精度和鲁棒性如何？

主要发现

• 当 θ = 0（显式）或 θ = 1（后向欧拉）时，θ-格式在空间上呈现二阶收敛，在时间上为一阶收敛；当 θ = 0.5（Crank-Nicolson）时，空间和时间均达到二阶收敛。
• 当 θ ∈ [0, 0.5) 时，θ-格式在满足 CFL 条件 Δt ≤ CFL·(Δx)²/(4μ·(1−2θ)) 下保持稳定，数值测试结果确认时间收敛率为约 1.0，空间收敛率为约 2.0（当 Δx = Δy 时）。
• 通过 qf1pTlump 实现质量矩阵对角化后，质量矩阵变为对角矩阵，易于求逆，从而提升计算效率，并扩展显式格式的 CFL 条件。
• 热方程的数值结果表明，当 θ = 0 时，空间收敛率为约 1.99992，时间收敛率为约 0.99996，与理论预期一致。
• 当 θ = 1/2 时，空间收敛率为 2.00195，时间收敛率也为 2.00195，表明在两个维度上均达到二阶精度。
• 软件能够以高精度重现精确解，表现为 L² 误差随网格细化和时间步长减小而呈比例下降。