[论文解读] Augmented Lagrangian preconditioners for fictitious domain formulations of elliptic interface problems

Michele Benzi, Marco Feder+2 more|arXiv (Cornell University)|Mar 13, 2026

Advanced Numerical Methods in Computational Mathematics被引用 0

一句话总结

提出一种用于FD-DLM椭圆界面问题的增广拉格朗日预conditioner，以在大系数跃变下改善FGRES性能，并提出一个使用AMG对对角块进行折衷处理的修改型变体。

ABSTRACT

We present a novel augmented Lagrangian (AL) preconditioner for the solution of linear systems arising from finite element discretizations of elliptic interface problems with jump coefficients. The method is based on the Fictitious Domain with Distributed Lagrange Multipliers formulation and it is designed to improve the convergence of the Flexible Generalized Minimal Residual (FGMRES) method in the presence of large coefficient jumps. To reduce the computational cost, we also introduce a cheaper block-triangular variant of the preconditioner. We prove eigenvalue clustering for the ideal AL preconditioner and study the limiting behavior of the spectrum for the modified variant in terms of parameters and the size of the jumps. Numerical experiments on different immersed geometries confirm mesh-independent iteration counts and robustness over large coefficient jumps, with substantial reductions in wall-clock time for the modified approach.

研究动机与目标

• 通过分布式拉格朗格乘子将跳跃系数的椭圆界面问题转化为虚拟域方法来建模与求解。
• 开发并分析增广拉格朗日预conditioner，以提高由结果的鞍点系统的迭代求解器的收敛性。
• 引入计算成本更低的修改型AL预conditioner并评估其性能。
• 提供谱分析，展示特征值聚类性及对系数跃变的鲁棒性。
• 通过数值实验（包括3D弹性问题）验证所提出的预conditioners。

提出的方法

• 通过FD-DLM方法将椭圆界面问题定式为3×3块鞍点系统。
• 推导理想的增广拉格朗日(AL)预conditioner，使其达到2×2块结构以实现高效求解。
• 引入带块三角形形式的修改AL预conditioner，使对角块可使用标准AMG求解器。
• 选择加增算子W为M^2（M为乘子质量矩阵），以确保均匀的特征值界。
• 提供谱分析以确立预conditioned系统的特征值聚集性与鲁棒性。
• 通过在灵活的GMRES框架下描述理想预conditioner的近似实现来给出实际的实现路径。

实验结果

研究问题

• RQ1 如何为虚拟域FD-DLM形式的椭圆界面问题设计增广拉格朗日框架，以在大系数跃变下实现鲁棒性？
• RQ2 对于所得到的鞍点系统，理想和修改AL预conditioner的谱性质及其对网格参数与系数对比度的依赖性为何？
• RQ3 是否可以使用一种计算成本更低但仍有效的修改型AL预conditioner来保持鲁棒性并提供实际的求解时间优势？
• RQ4 在复杂3D问题及异质材料的FD-DLM设置中，所提出的预conditioners表现如何？

主要发现

• 理想的AL预conditioner在大系数跃变下表现出有利的特征值聚集性和鲁棒性。
• 修改型AL预conditioner通过利用标准AMG求解器处理对角块，在实现上显著简化，同时保持性能。
• 数值实验显示网格无关的迭代次数，以及在包括挑战性3D弹性测试在内的多种配置下的鲁棒性。
• 使用W = M^2可使增广块对称正定，并有助于使特征值与网格尺寸h及系数跃变无关地界定。
• 增广方法还消除了(2,2)块的奇异性，帮助求解器稳定性。
• 修改型方案在实际 wall-clock 时间方面相较理想形式有显著降低。