[논문 리뷰] Jacobi-Galerkin spectral method for eigenvalue problems of Riesz fractional differential equations

Lizhen Chen, Zhiping Mao외 1명|arXiv (Cornell University)|2018. 03. 09.

Fractional Differential Equations Solutions참고 문헌 17인용 수 3

한 줄 요약

이 논문은 동차 딜리클레 경계 조건을 가진 리프츠 분수 미분 방정식 고유값 문제를 해결하기 위한 자코비-갈레르킨 스펙트럴 방법을 제안한다. 특별히 정의된 소볼레프 공간 내에서 대칭적인 약한 형태로 문제를 재구성함으로써, 행렬의 대칭성과 양의 정부호성을 유지하여 엄밀한 오차 추정과 지수 수렴성을 보장한다. 주요 기여는 Weyl 유형의 점 游 법칙 $\lambda_n = \mathcal{O}(n^{2\alpha})$ 및 조건수 스케일링 $\chi_N = \mathcal{O}(N^{4\alpha})$의 수치적 검증이다. 이 방법은 광범위한 수치 실험을 통해 스펙트럴 정확도를 확보하고 이론적 수렴률을 확인한다.

ABSTRACT

An efficient Jacobi-Galerkin spectral method for calculating eigenvalues of Riesz fractional partial differential equations with homogeneous Dirichlet boundary values is proposed in this paper. In order to retain the symmetry and positive definiteness of the discrete linear system, we introduce some properly defined Sobolev spaces and approximate the eigenvalue problem in a standard Galerkin weak formulation instead of the Petrov-Galerkin one as in literature. Poincaré and inverse inequalities are proved for the proposed Galerkin formulation which finally help us establishing a sharp estimate on the algebraic system's condition number. Rigorous error estimates of the eigenvalues and eigenvectors are then readily obtained by using Babuška and Osborn's approximation theory on self-adjoint and positive-definite eigenvalue problems. Numerical results are presented to demonstrate the accuracy and efficiency, and to validate the asymptotically exponential oder of convergence. Moreover, the Weyl-type asymptotic law $ λ_n=\mathcal{O}(n^{2α})$ for the $n$-th eigenvalue $λ_n$ of the Riesz fractional differential operator of order $2α$, and the condition number $N^{4α}$ of its algebraic system with respect to the polynomial degree $N$ are observed.

연구 동기 및 목표

• 리프츠 분수 미분 방정식의 고유값 문제에 대해 고차수이고 안정적인 수치적 방법을 개발한다.
• 이산 시스템의 대칭성과 양의 정부호성을 유지하여 조건 수가 안정적이고 대수적 고유값 문제를 효율적으로 해결할 수 있도록 보장한다.
• 자기수반이고 양의 정부호인 고유값 문제에 대한 바부시카와 오스본의 근사 이론을 이용하여 고유값과 고유벡터에 대한 엄밀한 오차 추정을 수립한다.
• Weyl 유형의 점 游 법칙 $\lambda_n = \mathcal{O}(n^{2\alpha})$ 및 조건수 스케일링 $\chi_N = \mathcal{O}(N^{4\alpha})$의 수치적 검증을 수행한다.
• 다양한 분수 차수 $2\alpha \in (0,2)$ 및 최대 $5.6$까지의 고차수에 대해 방법의 지수 수렴 속도와 안정성을 입증한다.

제안 방법

• 경계에서 고유함수의 특이 행동을 수용할 수 있도록 특별히 구성된 소볼레프 공간 내에서 갈레르킨 약한 형태를 적용한다.
• 고유함수를 근사하기 위해 자코비 기저 함수 $\{(1-x^2)^\alpha J^{\alpha,\alpha}_n(x)\}$를 사용하여 $x = \pm 1$에서 $\sim (1-x^2)^\alpha$ 형태의 특이성을 정확히 반영한다.
• 이 형식은 얻어진 대수적 시스템의 대칭성과 양의 정부호성을 보장하여 일반화된 고유값 문제의 효율적이고 안정적인 해법을 가능하게 한다.
• 제안된 갈레르킨 공간에 대해 푸앵카레 및 역부등식을 유도하여 이는 이산 시스템의 조건수를 유계화하는 데 필수적이다.
• 자기수반이고 양의 정부호인 고유값 문제에 대한 바부시카와 오스본 이론을 이용하여 고유값과 고유벡터의 오차 추정을 유도한다.
• 다양한 다항식 차수 $N$과 분수 차수 $2\alpha$를 사용하여 고정밀도 해법기를 활용한 수치 실험을 수행하여 수렴성과 점근적 행동을 검증한다.

실험 결과

연구 질문

• RQ1 제안된 자코비-갈레르킨 스펙트럴 방법이 리프츠 분수 미분 방정식의 고유값에 대해 지수 수렴을 달성하는가?
• RQ2 이산 시스템의 조건수가 $N$이 증가함에 따라 $\mathcal{O}(N^{4\alpha})$로 정확히 예측되는가?
• RQ3 계산된 고유값이 큰 $n$에 대해 Weyl 유형의 점 游 법칙 $\lambda_n = \mathcal{O}(n^{2\alpha})$를 따르는가?
• RQ4 고차수 분수 차수 $2\alpha > 2$인 경우, 예를 들어 $3.6$ 및 $5.6$에서도 방법이 높은 정확도와 안정성을 유지하는가?
• RQ5 수치적 고유함수들이 끝점 근처에서 기대되는 특이 행동 $\sim (1-x^2)^\alpha$를 얼마나 잘 캐릭터라이즈하는가?

주요 결과

• 고유값에 대해 점차적으로 지수 수렴을 달성하며, 다항식 차수 $N$이 증가함에 따라 오차가 급격히 감소하는 수치 결과로 확인된다.
• 고유값은 $\lambda_n = \mathcal{O}(n^{2\alpha})$로 스케일링되며, $2\alpha \in (1.2, 2.0)$ 범위에서 $\lambda_n$과 $n^{2\alpha}$ 사이의 강한 선형 의존성이 수치 결과로 확인된다.
• 이산 시스템의 조건수는 $\chi_N = \mathcal{O}(N^{4\alpha})$로 스케일링되며, $\chi_N$과 $N$의 로그-로그 플롯이 예측된 거듭제곱 법칙과 일치함으로써 확인된다.
• $2\alpha = 2$일 경우, 첫 번째 고유값 $\lambda_1$은 정확한 값 $\left(\frac{\pi}{2}\right)^2$로 수렴하며, $n$이 증가함에 따라 $\rho_n = \lambda_n / \left(\frac{\pi n}{2}\right)^{2\alpha}$는 1로 수렴한다.
• 고차수 분수 차수인 $2\alpha = 3.6$ 및 $5.6$에서도 방법은 정확하고 효율적이며, 일관된 수렴 행동이 관찰된다.
• 수치 결과는 첫 번째 고유함수 $\psi_1(x)$가 경계 근처에서 $\sim (1-x^2)^\alpha$로 행동한다는 이론적 예측을 검증하며, 자코비 기저가 이를 효과적으로 캡처함을 보여준다.