[論文レビュー] 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

ひとこと要約

本稿は、2次元非局所拡散問題における、非局所ノイマン境界条件の新規な定式化を提案する。この方法は、ホライズンパラメータ $\delta \to 0$ の際に局所ノイマン問題へ $O(\delta^2)$ 収束を達成し、長年の問題であった最適でない収束性と不自然な表面効果を解消する。境界積分表現を導出し、適切なテイラー展開とカーネル重み付けを用いることで、局所極限を正確に捉える。この手法は、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$ の際に、非局所ノイマン境界条件が局所ノイマン問題へ最適収束するよう定式化すること。
• 不適切な非局所フラックス条件の適用によって生じる、境界付近での物理的に不自然な効果（例：境界軟化）を除去すること。
• 古典的ノイマン条件と数学的に整合性を持つ境界処理を構築し、漸近的整合性を維持すること。
• 角を持つ領域や不規則な形状の領域に対しても、2次収束性を保つように定式化を拡張すること。
• 混合境界値問題においても $O(\delta^2)$ 収束性を達成する、メッシュフリーで漸近的整合性を持つ離散化を構築すること。

提案手法

• 境界点における非局所ホライズン上での重み付き積分を用いて、古典的ノイマン条件を拡張することで、非局所ノイマン境界条件を導出。境界点における方向微分と法線成分を組み込む。
• 境界点周りの解のテイラー展開を用い、局所PDEと $O(\delta^2)$ まで一致するようにし、局所極限への一貫性を確保。
• 領域境界の幾何的性質、特に最近接境界点における法線成分と接線成分を考慮した修正されたカーネル積分を導入。
• 漸近的整合性理論を適用し、非局所解が $\delta \to 0$ の際に $O(\delta^2)$ の速度で局所解へ収束することを保証。
• 径数関数または粒子配置を用いたメッシュフリー離散化を構築。非局所作用素と境界条件は、数値積分とカーネル積分により近似。
• 角を持つ領域を用いた数値実験により、幾何的特異性に対しても2次収束性を確認。

実験結果

リサーチクエスチョン

• RQ1 ホライズンパラメータ $\delta \to 0$ の際に、非局所ノイマン境界条件が $O(\delta^2)$ の速度で局所ノイマン問題へ収束するように定式化可能か？
• RQ2 非局所境界条件は、ペリダイナミクスモデルにおける境界軟化などの不自然な表面効果を回避するようにどのように構築できるか？
• RQ3 角を持つような非滑らか境界を持つ領域においても、提案された境界処理は漸近的整合性を保っているか？
• RQ4 混合ディリクレ＝ノイマン問題へもこの手法を拡張可能か？ また、高次収束性を維持できるか？
• RQ5 得られたメッシュフリー離散化は安定的かつ最適な $O(\delta^2)$ 精度で収束するか？

主な発見

• 提案された非局所ノイマン境界条件は、$L^\infty(\Omega)$ ノルムにおいて、局所ノイマン問題へ $O(\delta^2)$ 収束を達成し、非局所方程式自体の最適収束速度と一致する。
• 従来の非局所境界処理に見られる物理的に不自然な表面効果（例：不自然な軟化）が、本手法により明確に除去された。
• 数値実験により、角を持つ領域に対しても2次収束性が確認され、幾何的不規則性に対して高いロバスト性を示した。
• 本手法の漸近的整合性は厳密に証明され、非局所解が $\delta \to 0$ の際に局所解へ収束することが保証された。
• 新しい境界条件に基づくメッシュフリー離散化は $O(\delta^2)$ 収束を達成し、理論的枠組みの有効性が検証された。