[论文解读] Preconditioning of weighted H(div)-norm and applications to numerical simulation of highly heterogeneous media

Johannes Kraus, Raytcho Lazarov+3 more|arXiv (Cornell University)|Jun 17, 2014

Advanced Numerical Methods in Computational Mathematics参考文献 21被引用 4

一句话总结

本论文提出了一种针对高度异质介质中二阶椭圆问题的混合有限元离散化问题的鲁棒代数多重水平预条件子，采用基于加法Schur补近似的辅助空间多重网格（ASMG）方法。该方法在网格尺寸和材料对比度变化下均实现了与之无关的最优收敛性，数值实验验证了其在SPE10等基准问题上的有效性，问题规模最高达787,456自由度。

ABSTRACT

In this paper we propose and analyze a preconditioner for a system arising from a finite element approximation of second order elliptic problems describing processes in highly het- erogeneous media. Our approach uses the technique of multilevel methods and the recently proposed preconditioner based on additive Schur complement approximation by J. Kraus (see [8]). The main results are the design and a theoretical and numerical justification of an iterative method for such problems that is robust with respect to the contrast of the media, defined as the ratio between the maximum and minimum values of the coefficient (related to the permeability/conductivity).

研究动机与目标

• 解决由高度异质多孔介质中椭圆问题的混合有限元离散化所产生的大型不定鞍点系统所面临的挑战。
• 开发一种在材料系数极端变化（高对比度）下仍保持有效性和高效性的预条件子。
• 确保对网格细化（h）和系数对比度（κ）的鲁棒性，这对石油工程和环境模拟中的实际仿真至关重要。
• 在包含多尺度特征（如通道和夹杂物）的复杂、真实基准问题（如SPE10问题）上验证该方法。
• 展示在多个网格层次上的可扩展性与一致收敛性，确保总计算工作量与自由度呈线性关系。

提出的方法

• 使用Raviart-Thomas RT0元对速度进行离散化，压力使用分片常数元，构建混合有限元格式。
• 将问题表述为块结构的鞍点线性系统：[M_α -B_div^T; -B_div 0][u; p] = [0; f]，其中M_α表示加权H(div)内积。
• 应用基于加法Schur补近似的代数多重水平迭代（AMLI）循环辅助空间多重网格（ASMG）预条件子。
• 在ASMG求解器中采用W循环光滑器，包含一次前光滑和一次后光滑步骤，并使用MinRes作为全局系统的Krylov子空间求解器。
• 通过系数k(x)引入加权H(div)范数，并设计预条件子以对这一权重的变化保持鲁棒性。
• 利用离散Poincaré不等式和多重网格转移算子，确保在不同尺度间的一致收敛性。

实验结果

研究问题

• RQ1 能否设计出一种预条件子，使其在高度异质介质中对网格细化和极端材料对比度均保持收敛鲁棒性？
• RQ2 所提出的ASMG预条件子在包含多尺度特征的真实复杂基准问题（如SPE10问题）上的表现如何？
• RQ3 加法Schur补近似与AMLI循环多重网格的结合是否能对高对比度H(div)问题实现一致收敛？
• RQ4 非线性光滑化以及V循环与W循环配置对残差削减和迭代次数有何影响？
• RQ5 当右端项非零时，预条件子在多大程度上仍保持有效性，从而扩展理论假设的适用范围？

主要发现

• 所提出的ASMG预条件子在MinRes求解器中实现了鲁棒收敛，且MinRes与ASMG的迭代乘积在不同网格层次上几乎保持恒定。
• 对于SPE10基准问题（case [c]），当自由度从3,136增至787,456时，MinRes迭代次数仅从13次增至18次，表明具有良好的可扩展性。
• ASMG W循环配合一次光滑步骤在所有测试案例中（包括高对比度案例[a]和[b]）均实现了统一收敛，迭代次数在各层次间保持稳定。
• 即使不使用光滑步骤（m=0），非线性ASMG V循环在SPE10案例中仍表现出鲁棒收敛，证实了该方法的鲁棒性。
• 预条件系统条件数的数值界表明其独立于网格尺寸和对比度，即‖π_D̃‖²_Ã的上界为常数。
• 总计算工作量（以精细网格矩阵-向量乘法次数衡量）与自由度呈线性关系，证实了最优复杂度。