[논문 리뷰] An inverse theorem for compact Lipschitz regions in $R^d$ using localized kernel bases

Thomas Hangelbroek, Francis J. Narcowich외 2명|arXiv (Cornell University)|2015. 08. 12.

Numerical methods in engineering참고 문헌 26인용 수 4

한 줄 요약

이 논문은 $ vert\mathbb{R}^d$ 내의 컴acts Lipschitz 도메인에서 Matérn 및 표면 스플라인 RBFs로부터 유도된 국소화된 커널 기반을 사용하여, 반경 기반 함수(RBF) 근사에 대한 역추정(inverse estimates)을 수립한다. 이를 통해 이러한 유한차원 공간 내 함수의 소볼레프 노름이 채움 거리 $h$에 의해 스케일링된 $L_p$ 노름으로 유계임을 증명하며, 메쉬리스 방법에서 오랫동안 존재했던 격차를 해결하고, PDE 해법기에서 노름 안정성과 조건수 상한을 가능하게 한다.

ABSTRACT

While inverse estimates in the context of radial basis function approximation on boundary-free domains have been known for at least ten years, such theorems for the more important and difficult setting of bounded domains have been notably absent. This article develops inverse estimates for finite dimensional spaces arising in radial basis function approximation and meshless methods. The inverse estimates we consider control Sobolev norms of linear combinations of a localized basis by the $L_p$ norm over a bounded domain. The localized basis is generated by forming local Lagrange functions for certain types of RBFs (namely Matérn and surface spline RBFs). In this way it extends the boundary-free construction of Fuselier, Hangelbroek, Narcowich, Ward and Wright.

연구 동기 및 목표

• 메쉬리스 방법에서의 핵심적 격차를 메우기 위해, 유계이고 컴팩트한 Lipschitz 도메인에서 반경 기반 함수 근사에 대한 역추정을 개발하는 것.
• 이전에 경계가 없는 설정에서만 알려진 국소화된 RBF 기반 이론을 경계가 존재하는 도메인으로 확장하여 실용적인 PDE 응용에 더 부합하는 환경을 제공하는 것.
• RBF 공간 내 함수의 소볼레프 노름이 채움 거리 $h$에 의해 스케일링된 $L_p$ 노름으로 제어됨을 증명하여 안정성과 수렴성을 보장하는 것.
• 갈레르킨 방법에 필수적인 $L_p$ 안정성과 노름 등가 결과를 제공하여, 특히 혼합 유한요소 및 안장점 문제에 유리한 조건을 마련하는 것.
• 유계 도메인 설정에서의 추적 추정(trace estimates) 및 버르슈타인-니콜스키 부등식을 수립하여, 이전에 RBF 기반 방법에서는 가용하지 않았던 결과를 확보하는 것.

제안 방법

• 컴팩트한 Lipschitz 도메인에서 Matérn 및 표면 스플라인 커널에 대해 국소화된 라그랑주 함수를 이용하여 국소화된 RBF 기반을 구성한다.
• 단위 분할(partition of unity)과 반경 $\gamma q$의 구 $B_\xi$ 위에서 지지되는 국소화된 커널 $\chi_\xi$를 사용하며, $\gamma$ 는 $\rho$, $d$, $m$, $p$ 에 따라 달라진다.
• 멀리 떨어진 기반 함수 간의 상호작용을 제어하기 위해 준삼각부등식과 지수감쇠 추정을 적용한다.
• K-기능 방법과 $\ell_p$-노름 추정을 활용하여 $\|s\|_{W_p^\sigma(\Omega)}$ 를 $\|s\|_{L_p(\Omega)}$ 로 표현하는 bound를 유도한다.
• 기반의 $L_p$ 안정성을 증명하기 위해 $\|s\|_{L_p(\Omega)} \geq c_1 q^{d/p} \|\mathbf{a}\|_{\ell_p(\Xi)}$ 를 유도하며, 이때 $c_1 \geq C \rho^{\frac{d(d-2m)}{2\epsilon p}} (\log \rho)^{-\frac{d^2}{p\epsilon}}$ 이다.
• 기반 함수 간의 상호작용을 이중분할(dyadic decomposition)하고, RBF 커널의 지수감쇠 성질을 활용하여 비대칭항을 유계로 제한한다.

실험 결과

연구 질문

• RQ1 유계 Lipschitz 도메인에서 RBF 공간에 대해 $\|s\|_{W_p^\sigma(\Omega)} \leq C h^{-\sigma} \|s\|_{L_p(\Omega)}$ 형태의 역추정을 수립할 수 있는가?
• RQ2 갈레르킨 방법에서 $L_p$ 안정성과 강성행렬의 비대칭항 감쇠를 보장하기 위해 국소화된 RBF 기반을 어떻게 구성할 수 있는가?
• RQ3 역추정 상수의 정확한 의존성은 채움 거리 $h$, 차원 $d$, 그리고 부드러움 정도 $\sigma$ 에 따라 어떻게 달라지는가?
• RQ4 경계가 존재하는 도메인 설정에서의 RBF 기반 설정에서 추적 추정과 버르슈타인-니콜스키 부등식을 유도할 수 있는가? 이는 경계가 없는 경우에만 알려져 있었다.
• RQ5 경계가 없는 다양체에서의 국소화된 기반 프레임워크가 경계가 존재하는 도메인으로 얼마나 잘 일반화될 수 있는가?

주요 결과

• 논문은 컴팩트한 Lipschitz 도메인에서 RBF 공간 내 함수에 대해 $\|s\|_{W_p^\sigma(\Omega)} \leq C h^{-\sigma} \|s\|_{L_p(\Omega)}$ 의 역추정을 수립하였으며, 이때 $C$ 는 점 집합에 무관하지만 $d$, $p$, $\sigma$ 에 따라 달라진다.
• 국소화된 기반의 $L_p$ 안정성이 증명되었으며, $\|s\|_{L_p(\Omega)} \geq c_1 q^{d/p} \|\mathbf{a}\|_{\ell_p(\Xi)}$ 이고, 이때 $c_1 \geq C \rho^{\frac{d(d-2m)}{2\epsilon p}} (\log \rho)^{-\frac{d^2}{p\epsilon}}$ 이다. 여기서 $0 < \epsilon < m - d/2$ 이다.
• 역추정은 갈레르킨 방법에서 강성행렬의 조건수가 최대 $h^{-2\sigma}$ 의 속도로 증가함을 의미하며, 이는 수치적 안정성을 보장한다.
• 추적 추정과 버르슈타인-니콜스키 부등식은 역추정의 직접적 결과로서 도출되었으며, 이는 이들의 유효성을 유계 도메인으로 확장함을 의미한다.
• Matérn 및 표면 스플라인 RBFs에 대해 국소화된 라그랑주 함수를 이용한 국소화된 기반의 구성 방법이 경계가 존재하는 상황에서도 효과적이고 안정적임을 보였다.
• 이러한 결과는 메쉬리스 방법 문헌에서 오랫동안 부재했던, 특히 RBF-갈레르킨 및 콜로케이션 기반 방법에 대해 유계 도메인에서의 역추정이 해결되었음을 의미한다.