[논문 리뷰] An Asymptotically Compatible Approach For Neumann-Type Boundary Condition On Nonlocal Problems

Huaiqian You, Xin Lü외 2명|arXiv (Cornell University)|2019. 08. 11.

Numerical methods in engineering참고 문헌 69인용 수 10

한 줄 요약

이 논문은 2D 비국소 확산 문제에 대해 $O(\delta^2)$ 수렴성을 확보하는 새로운 비국소 노이만 경계 조건을 제안한다. 이는 경계 조건의 수렴 속도가 $\delta \to 0$일 때 최적의 수렴 속도를 달성하며, 오랫동안 지속된 부분 수렴성 문제와 인위적인 표면 효과를 해결한다. 이 방법은 정확한 테일러 전개와 커널 가중치를 통해 국소 극한을 정확히 포착하는 경계 적분 형식을 유도함으로써 점차적 일致성(Asymptotic Compatibility)을 보장한다. 메쉬프리 이산화를 통해 두 번째 차수 정확도를 확보한 결과가 검증되었다.

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)$까지 일致시키며, 비국소 편미분방정식과의 일致성을 확보한다.
• 영역 경계의 기하학적 특성, 특히 가장 가까운 경계점에서의 법선 및 격자 성분을 고려한 수정된 커널 적분을 도입한다.
• 점차적 일치성 이론을 적용하여 비국소 해가 $\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)$ 수렴성을 확보하였으며, 이론적 프레임워크의 타당성을 검증하였다.