[논문 리뷰] Stability analysis of RBF-FD and WLS based local strong form meshless methods on scattered nodes

Mitja Jančič, Gregor Kosec|arXiv (Cornell University)|2022. 02. 28.

Numerical methods in engineering참고 문헌 25인용 수 3

한 줄 요약

이 논문은 2D 및 3D 영역의 산점된 노드에서 두 가지 국소 강형 메쉬리스 방법—다항스플린을 보완한 RBF-FD(반경기저함수-유한차분)와 단지 단항식을 사용하는 WLS(가중 최소제곱)—의 안정성과 정확도를 비교한다. RBF-FD는 고차수 근사에서 훨씬 더 뛰어난 안정성과 높은 정확도를 보이며, WLS는 저차수 경우에서 더 안정적이고 효율적임을 발견한다.

ABSTRACT

The popularity of local meshless methods in the field of numerical simulations has increased greatly in recent years. This is mainly due to the fact that they can operate on scattered nodes and that they allow a direct control over the approximation order and basis functions. In this paper we analyse two popular variants of local strong form meshless methods, namely the radial basis function-generated finite differences (RBF-FD) using polyharmonic splines (PHS) augmented with monomials, and the weighted least squares (WLS) approach using only monomials. Our analysis focuses on the accuracy and stability of the numerical solution computed on scattered nodes in a two- and three-dimensional domain. We show that while the WLS variant is a better choice when lower order approximations are sufficient, the RBF-FD variant exhibits a more stable behavior and a higher accuracy of the numerical solution for higher order approximations, but at the cost of higher computational complexity.

연구 동기 및 목표

• 비정규적이고 산점된 노드 분포에서 두 가지 인기 있는 국소 강형 메쉬리스 방법—RBF-FD와 WLS—의 수치적 안정성과 정확도를 평가하기 위해.
• 안정성, 정확도, 계산 비용 측면에서 저차수와 고차수 근사에 최적의 방법을 결정하기 위해.
• 스티encil 크기와 노드 분포가 해의 신뢰성 및 수렴 행동에 미치는 영향을 조사하기 위해.
• 복잡한 기하구조를 포함하는 실제 시뮬레이션에서 RBF-FD와 WLS 사이에서 선택하는 데 실용적인 지침을 제공하기 위해.
• 다양한 메esh 수준과 고차원(2D 및 3D)에서 두 방법의 견고성을 평가하기 위해.

제안 방법

• 연구는 국소 강형 메쉬리스 프레임워크를 사용하며, 미분 연산자는 가중 최소제곱(WLS) 또는 반경기저함수를 이용한 유한차분(RBF-FD)을 통해 근사한다.
• WLS의 경우, 차수 m까지의 단항식이 기저 함수가 되며, 가중 최소제곱 오차 노름의 최소화를 통해 가중치가 계산된다.
• RBF-FD의 경우, 다항스플린(PHS)이 반경기저함수로 사용되며, 근사 안정성과 정확도 향상을 위해 단항식과 함께 보완된다.
• 근사는 이웃 노드로 구성된 국소 스티encil에서 구성되며, 노드 수가 기저 함수 수를 초과할 경우 선형 최소제곱 접근법으로 해를 구한다.
• 무한노름 오차를 계산하고 수렴성 및 안정성을 평가하기 위해 해석적 해가 존재하는 포아송 문제를 벤치마크로 사용한다.
• 안정성은 서로 다른 노드 분포를 가진 100회의 독립 실행에 걸친 오차의 정규화된 산란을 계산하여 정량화한다.

실험 결과

연구 질문

• RQ1 산점된 노드에서 저차수 및 고차수 근사에 대해 RBF-FD와 WLS 방법의 수치적 안정성은 어떻게 비교될 수 있는가?
• RQ2 기저 함수의 선택(단항식 대비 PHS+단항식)이 2D 및 3D 영역에서 수렴 속도와 오차 분포에 중대한 영향을 미치는가?
• RQ3 각 방법에 대해 노드 수와 근사 차수 증가에 따라 오차의 정규화된 산란은 어떻게 변화하는가?
• RQ4 어느 정도의 차수를 초과하면 한 방법이 다른 방법보다 훨씬 더 안정해지는 임계값이 존재하는가?
• RQ5 각 방법의 계산 복잡도는 근사 차수와 차원 증가에 따라 어떻게 변화하는가?

주요 결과

• RBF-FD 변형은 2D에서 모든 노드 수에 대해 약 1의 정규화된 산란을 보이며, 일관되고 안정된 성능을 나타내지만, WLS 변형은 산란이 약 두 개의 지수 차수 더 크고 더 예측 불가능하다.
• 고차수 근사(m > 2)의 경우, RBF-FD 방법은 훨씬 더 뛰어난 정확도와 안정성을 보이며, WLS 변형은 자주 수렴하지 못한다(예: 3D에서 m=6, N≈10^3일 때 무한노름 오차 약 10^1).
• 2D에서는 저차수 WLS 근사(m=2)가 RBF-FD보다 안정성과 정밀도 측면에서 뛰어나며, 더 낮은 정규화된 산란과 계산 비용을 보인다.
• 3D에서는 RBF-FD 방법이 노드 수 증가에 따라 안정적이고 감소하는 산란을 유지하지만, WLS 방법은 고차수, 특히 m=6일 때 산란이 증가하는 경향을 보인다.
• WLS 변형은 m=2일 때 저차수 근사에서 RBF-FD보다 계산 비용이 더 저렴하고 더 안정하지만, 고차수에서는 안정성이 떨어진다.
• 고차수 근사에서 RBF-FD 방법은 WLS보다 수십 개의 지수 차수 더 안정하여, 비균일한 노드 분포에 대한 견고성이 뛰어나고 고정밀 시뮬레이션에 더 적합하다.