[论文解读] An Asymptotically Compatible Approach For Neumann-Type Boundary Condition On Nonlocal Problems

Huaiqian You, Xin Lü+2 more|arXiv (Cornell University)|Aug 11, 2019

Numerical methods in engineering参考文献 69被引用 10

一句话总结

该论文提出了一种用于二维非局部扩散问题的新型非局部诺伊曼边界条件，当非局部参数 $\delta \to 0$ 时，实现了与局部诺伊曼问题的 $O(\delta^2)$ 收敛率，解决了长期存在的次优收敛率和人为表面效应问题。该方法通过推导边界积分形式，结合细致的泰勒展开与核函数加权，精确捕捉局部极限，确保渐近相容性，经无网格离散化验证，达到二阶精度。

ABSTRACT

In this paper we consider 2D nonlocal diffusion models with a finite nonlocal horizon parameter $δ$ characterizing the range of nonlocal interactions, and consider the treatment of Neumann-like boundary conditions that have proven challenging for discretizations of nonlocal models. While existing 2D nonlocal flux boundary conditions have been shown to exhibit at most first order convergence to the local counter part as $δ ightarrow 0$, we present a new generalization of classical local Neumann conditions that recovers the local case as $O(δ^2)$ in the $L^{\infty}(Ω)$ norm. This convergence rate is optimal considering the $O(δ^2)$ convergence of the nonlocal equation to its local limit away from the boundary. We analyze the application of this new boundary treatment to the nonlocal diffusion problem, and present conditions under which the solution of the nonlocal boundary value problem converges to the solution of the corresponding local Neumann problem as the horizon is reduced. To demonstrate the applicability of this nonlocal flux boundary condition to more complicated scenarios, we extend the approach to less regular domains, numerically verifying that we preserve second-order convergence for domains with corners. Based on the new formulation for nonlocal boundary condition, we develop an asymptotically compatible meshfree discretization, obtaining a solution to the nonlocal diffusion equation with mixed boundary conditions that converges with $O(δ^2)$ convergence.

研究动机与目标

• 解决在 $\delta \to 0$ 时，非局部诺伊曼边界条件向局部诺伊曼问题最优收敛的挑战。
• 消除因非局部通量条件不当施加而引起的非物理表面效应（例如边界附近的软化现象）。
• 开发一种在数学上与经典诺伊曼条件一致且保持渐近相容性的边界处理方法。
• 在保持二阶收敛率的前提下，将该公式推广至具有角点和非规则几何形状的区域。
• 构建一种无网格、渐近相容的离散化方法，实现混合边界值问题的 $O(\delta^2)$ 收敛率。

提出的方法

• 通过在非局部作用域上对经典诺伊曼条件进行加权积分扩展，推导出非局部诺伊曼边界条件，引入边界点处的方向导数和法向分量。
• 利用边界点附近解的泰勒展开，将局部极限匹配至 $O(\delta^2)$ 阶，确保与局部PDE的一致性。
• 引入一种改进的核积分形式，考虑区域边界几何特性，特别是最近边界点处的法向与切向分量。
• 应用渐近相容性理论，确保当 $\delta \to 0$ 时，非局部解以 $O(\delta^2)$ 速率收敛至局部解。
• 采用径向基函数或粒子配点法实现无网格离散化，通过数值积分与核函数积分近似非局部算子及边界条件。
• 在具有角点的区域上进行数值验证，结果表明尽管存在几何奇点，仍保持二阶收敛率。

实验结果

研究问题

• RQ1 能否构造一种非局部诺伊曼边界条件，使其在 $\delta \to 0$ 时以 $O(\delta^2)$ 速率收敛至局部诺伊曼问题？
• RQ2 如何构建非局部边界条件以避免在非局部模型中出现如边界软化等人为表面效应？
• RQ3 所提出的边界处理方法在具有非光滑边界（如角点）的区域中是否仍能保持渐近相容性？
• RQ4 该方法能否推广至混合狄利克雷-诺伊曼问题，同时保持高阶收敛率？
• RQ5 所得到的无网格离散化是否稳定且收敛，且具备最优的 $O(\delta^2)$ 精度？

主要发现

• 所提出的非局部诺伊曼边界条件在 $L^\infty(\Omega)$ 范数下实现了与局部诺伊曼问题的 $O(\delta^2)$ 收敛率，与非局部方程本身的最优收敛率一致。
• 该方法成功消除了早期非局部边界处理中常见的非物理表面效应（如人为软化）。
• 数值实验表明，即使在具有角点的区域中，该方法仍保持二阶收敛率，展现出对几何不规则性的鲁棒性。
• 该方案的渐近相容性得到严格证明，确保当 $\delta \to 0$ 时，非局部解收敛至局部解。
• 基于新边界条件的无网格离散化实现了 $O(\delta^2)$ 收敛率，验证了理论框架的有效性。