[논문 리뷰] A universal solution scheme for fractional and classical PDEs

Yixuan Wu, Yanzhi Zhang|arXiv (Cornell University)|2021. 01. 30.

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

한 줄 요약

이 논문은 일반화된 역다중함수(GIMQ) 라디얼 기저 함수를 기반으로 한 통합된 메쉬리스 편미분법을 제안한다. 이 방법은 $\alpha \in (0,2]$ 에 대해 $(-\Delta)^{\frac{\alpha}{2}}$ 연산자를 갖는 고전적 및 분수형 PDE를 해결하기 위한 것이다. GIMQ 함수의 라플라스 연산을 가우스 제타함수를 통해 해석적으로 표현함으로써, 국소적 및 비국소적 라플라스 연산의 이산화를 통합하고, 초특이적 적분의 수치 근사를 피하며, 차원에 관계없이 적은 수의 미지수와 낮은 계산 비용으로 고정밀도를 달성할 수 있다.

ABSTRACT

We propose a unified meshless method to solve classical and fractional PDE problems with $(-Δ)^{\fracα{2}}$ for $α\in (0, 2]$. The classical ($α= 2$) and fractional ($α< 2$) Laplacians, one local and the other nonlocal, have distinct properties. Therefore, their numerical methods and computer implementations are usually incompatible. We notice that for any $α\ge 0$, the Laplacian $(-Δ)^{\fracα{2}}$ of generalized inverse multiquadric (GIMQ) functions can be analytically written by the Gauss hypergeometric function, and thus propose a GIMQ-based method. Our method unifies the discretization of classical and fractional Laplacians and also bypasses numerical approximation to the hypersingular integral of fractional Laplacian. These two merits distinguish our method from other existing methods for the fractional Laplacian. Extensive numerical experiments are carried out to test the performance of our method. Compared to other methods, our method can achieve high accuracy with fewer number of unknowns, which effectively reduces the storage and computational requirements in simulations of fractional PDEs. Moreover, the meshfree nature makes it free of geometric constraints and enables simple implementation for any dimension $d \ge 1$. Additionally, two approaches of selecting shape parameters, including condition number-indicated method and random-perturbed method, are studied to avoid the ill-conditioning issues when large number of points.

연구 동기 및 목표

• 같은 계산적 기반으로 고전적 ($\alpha=2$) 및 분수형 ($\alpha<2$) PDE를 모두 해결할 수 있는 통합 수치 프레임워크를 개발하는 것.
• 기존 국소적 및 비국소적 라플라스 연산자에 대한 수치 방법 간의 부적합성 문제를 해결하여 별도의 구현이 필요 없도록 하는 것.
• 분수형 라플라스 연산 계산에서 초특이적 적분의 수치 근사가 필요 없도록 하는 것.
• 저장소 및 계산 비용을 줄이며 효율적이고 고정밀도로 분수형 PDE를 시뮬레이션할 수 있도록 하는 것.

제안 방법

• 차원 $d \geq 1$ 에 대해 거듭제곱 $\beta = -(d+1)/2$ 를 갖는 일반화된 역다중함수(GIMQ) 라디얼 기저 함수를 사용한다.
• 가우스 제타함수를 통한 분석적 표현을 활용하여, GIMQ 함수에 $(-\Delta)^{\frac{\alpha}{2}}$ 를 적용한 결과를 정확하게 계산할 수 있으며, 이는 수치 적분이 필요 없음을 의미한다.
• 해법은 GIMQ 기저 함수를 산점된 노드 중심으로 선형 조합으로 근사하는 편미분 스펙트럴 콜로케이션 방법을 사용한다.
• 디리클레 경계 조건은 $\alpha=2$ 인 경우 도메인 경계($\partial\Omega$)에 직접 적용하고, $\alpha<2$ 인 경우 보완 영역($\Omega^c$)에 적용한다. 이는 분수형 라플라스 연산자의 점별 정의에 기반한다.
• 불안정성 방지를 위해 두 가지 형태의 형태 매개변수 선택 전략을 도입: 조건수 기반 선택 및 무작위 변형 방법. 이는 큰 점 집합에서의 불안정성 문제를 완화한다.
• 이 방법은 메쉬 기반 방식이 아니므로, 기하학적 제약 없이 임의의 차원 $d \geq 1$ 에서 쉽게 구현할 수 있다.

실험 결과

연구 질문

• RQ1 $(-\Delta)^{\frac{\alpha}{2}}$ 를 따르는 고전적 및 분수형 PDE를 동시에 해결할 수 있는 단일 수치적 기반을 개발할 수 있는가?
• RQ2 초특이적 적분의 수치적 통합 없이도 분수형 라플라스 연산자의 비국소적 성질을 정확히 포착할 수 있는가?
• RQ3 경계가 도메인과 접촉하지 않는 경우, 특히 분수형 확산에서 경계 조건이 해의 도메인에 미치는 영향은 어떠한가?
• RQ4 다양한 형태 매개변수 선택 전략이 대규모 문제에 대한 GIMQ 기반 방법의 조건수 및 정밀도에 미치는 영향은 어떠한가?
• RQ5 기존 RBF 또는 가우시안 기반 방법에 비해 제안된 방법이 정밀도 및 계산 효율성 측면에서 얼마나 뛰어나게 성능을 발휘하는가?

주요 결과

• 제안된 방법은 최근의 가우시안 기반 방법과 유사한 스펙트럴 정밀도를 달성하지만, 더 적은 수의 미지수와 낮은 계산 비용을 제공한다.
• GIMQ 기반 방법은 가우스 제타함수를 통한 해석적 표현을 활용하여 분수형 라플라스 연산에서 초특이적 적분의 수치 평가를 피한다.
• 균일하게 분포된 중심점에 대해, 가우시안 기반 방법의 최적 형태 매개변수는 점의 수에 더 민감하게 반응하는 반면, GIMQ 기반 방법은 그 민감도가 낮다.
• 조건수 기반 형태 매개변수 선택 전략은 큰 점 집합에서도 불안정성을 효과적으로 억제하고 정밀도를 유지하지만, 추가적인 계산 시간이 소요된다.
• 무작위 변형 형태 매개변수 방법은 계산 시간을 줄이면서도 만족스러운 정밀도와 안정성을 유지한다.
• 수치 결과는 분수형 라플라스 연산의 비국소적 성질로 인해 경계 영역이 도메인과 분리되어 있더라도 경계 조건이 해에 상당한 영향을 미치며, 이는 $\alpha$ 가 작을수록 또는 거리가 짧을수록 더욱 뚜렷하게 나타난다.