[论文解读] A rigorous Multiscale Method for semi-linear elliptic problems

Patrick Henning, Axel Målqvist+1 more|arXiv (Cornell University)|Nov 27, 2012

Advanced Mathematical Modeling in Engineering参考文献 28被引用 3

一句话总结

本文提出了一种严格的多尺度方法，用于求解具有高度振荡、非均匀系数的半线性椭圆问题。该方法通过在大小为 H log(H⁻¹) 的局部区域中进行细尺度计算，构建广义有限元基函数，实现了相对于粗网格尺寸 H 的线性 H¹ 收敛，并采用阻尼牛顿法求解由此产生的多尺度系统。

ABSTRACT

In this paper we propose and analyze a new Multiscale Method for solving semi-linear elliptic problemswith heterogeneous andhighly variable coefficient functions. For this purposewe construct a generalized finiteelement basisthat spansalow dimensionalmultiscale space. The basis is assembled by performing localized linear fine-scale computations in small patches that have a diameter of order Hlog(H −1) where H is the coarse mesh size. Without any assumptions on the type of the oscillations in the coefficients, we give a rigorous proof for a linear convergence of the H 1-error with respect to the coarse mesh size. To solve the arising equations, we propose an algorithm that is based on a damped Newton scheme in the multiscale space.

研究动机与目标

• 开发一种针对具有高度变化和非均匀系数的半线性椭圆问题的鲁棒多尺度方法。
• 构建一个低维多尺度空间，以准确捕捉细尺度特征，且无需对系数振荡施加限制性假设。
• 在最小正则性假设下，确保方法的收敛性并进行严格的误差分析。
• 设计一种高效的求解算法，用于求解由此产生的多尺度系统，采用阻尼牛顿法。

提出的方法

• 通过在直径为 H log(H⁻¹) 的局部区域中求解局部细尺度问题，构建广义有限元基函数，其中 H 为粗网格尺寸。
• 通过组合能够反映非均匀系数细尺度行为的局部解，组装多尺度空间。
• 在降维的多尺度空间中建立全局多尺度问题，保留原始问题的关键特征。
• 应用阻尼牛顿法求解在多尺度空间中通过 Galerkin 投影得到的非线性系统。
• 通过控制步长确保迭代求解器的收敛性，以保持全局收敛性质。
• 进行严格的误差分析，建立 H¹ 误差相对于粗网格尺寸 H 的线性收敛性。

实验结果

研究问题

• RQ1 是否可以构建一种多尺度方法，用于求解具有高度振荡系数的半线性椭圆问题，而无需假设周期性或尺度分离？
• RQ2 在多尺度基函数构造中，局部计算的最优区域大小是什么，以在精度和计算效率之间取得平衡？
• RQ3 所提出的方法是否在 H¹ 范数下相对于粗网格尺寸 H 实现线性收敛？
• RQ4 如何为由此产生的非线性多尺度系统设计一种鲁棒且高效的迭代求解器？
• RQ5 是否可以在不施加对系数可变性限制性假设的前提下，严格建立收敛性分析？

主要发现

• 该方法实现了 H¹ 误差相对于粗网格尺寸 H 的线性收敛，且不依赖于系数振荡的类型。
• 多尺度基函数通过大小为 H log(H⁻¹) 的局部区域计算构建，确保了计算效率和可扩展性。
• 该方法无需对系数振荡的性质作任何假设，例如周期性或尺度分离。
• 阻尼牛顿法确保了在多尺度公式中产生的非线性系统的鲁棒收敛。
• 由此产生的多尺度空间为低维空间，能够捕捉准确逼近解所必需的细尺度关键特征。
• 收敛性分析是严格的，并且在对系数函数的假设最少的情况下依然成立。