[논문 리뷰] Hermite spectral collocation methods for fractional PDEs in unbounded domains

Tao Tang, Huifang Yuan외 1명|arXiv (Cornell University)|2018. 01. 27.

Fractional Differential Equations Solutions인용 수 8

한 줄 요약

이 논문은 두 종류의 기저 함수인 과도 스케일링된 에르미트 함수 $e^{-x^2}H_n(x)$와 정규화된 에르미트 함수 $e^{-x^2/2}H_n(x)$를 사용하여 무한도메인 내 분수형 편미분방정식을 해결하기 위한 에르미트 스펙트럴 콜로케이션 방법을 제안한다. 재귀 관계를 통해 명시적 미분 행렬을 유도하고, 지수 감쇠를 보이는 해에 대해 스펙트럴 수렴을 보이며, 적응형 스케일링 인자로 정확도를 향상시킨다.

ABSTRACT

This work is concerned with spectral collocation methods for fractional PDEs in unbounded domains. The method consists of expanding the solution with proper global basis functions and imposing collocation conditions on the Gauss-Hermite points. In this work, two Hermite-type functions are employed to serve as basis functions. Our main task is to find corresponding differentiation matrices which are computed recursively. Two important issues relevant to condition numbers and scaling factors will be discussed. Applications of the spectral collocation methods to multi-term fractional PDEs are also presented. Several numerical examples are carried out to demonstrate the effectiveness of the proposed methods.

연구 동기 및 목표

• 기존 방법이 제한되어 있는 무한도메인 내 분수형 편미분방정식에 대해 효율적인 스펙트럴 콜로케이션 방법을 개발하기 위해.
• 무한공간 도메인 내 해의 감쇠 속도가 다양할 경우를 다루는 데 도전하기 위해.
• 에르미트 기반 스펙트럴 방법을 위한 명시적 미분 행렬을 구성하여 효율적인 계산을 가능하게 하기 위해.
• 스펙트럴 이산화에서 수렴 속도와 조건수에 미치는 스케일링 인자의 영향을 조사하기 위해.
• 뉴턴 반복을 통해 다항분수 라플라시안 방정식과 비선형 문제로의 방법 확장을 위해.

제안 방법

• 방법은 무한도메인 내 스펙트럴 근사에 사용하기 위한 글로벌 기저 함수로 에르미트 함수 $e^{-x^2}H_n(x)$와 $e^{-x^2/2}H_n(x)$를 사용한다.
• 분수형 편미분방정식을 약한 형태로 약화시키기 위해 가우스-에르미트 점에서 콜로케이션을 시행한다.
• 재귀 관계를 사용하여 명시적 미분 행렬을 유도함으로써 행렬-벡터 곱 연산을 효율적으로 수행할 수 있다.
• 해의 감쇠 속도와 일치시키기 위해 기저 함수에 스케일링 인자 $r$을 도입하여 수렴성을 향상시킨다.
• 뉴턴 반복을 통해 다항분수 라플라시안 방정식과 비선형 분수형 편미분방정식으로 방법을 확장한다.
• 비교를 위해 라그랑주 유형 기저 함수를 사용하는 등가 콜로케이션 체계도 논의된다.

실험 결과

연구 질문

• RQ1 어떻게 에르미트 함수를 사용하여 분수형 편미분방정식을 무한도메인에 적응시킬 수 있는가?
• RQ2 비표준 감쇠 속도를 가진 해에 대해 수렴성을 향상시키는 데 스케일링 인자가 수행하는 역할은 무엇인가?
• RQ3 도함수 행렬의 조건수는 차수 증가에 따라 어떻게 변화하는가?
• RQ4 무한도메인 내 에르미트 기반 스펙트럴 방법에 대해 명시적 미분 행렬을 재귀적으로 유도할 수 있는가?
• RQ5 이 방법은 $\mathbb{R}^d$ 내 비선형 및 다항분수형 편미분방정식에 대해 얼마나 효과적인가?

주요 결과

• 수치 예제에서 빠른 오차 감쇠를 통해 지수 감쇠를 보이는 해에 대해 스펙트럴 수렴을 달성한다.
• 다항분수 라플라시안 문제에서 스케일링 인자 $r = \sqrt{1.5}$ 선택이 가중치 및 최대 노름 모두에서 가장 빠른 수렴을 이끌어낸다.
• 정규화된 에르미트 함수 기저($e^{-x^2/2}H_n(x)$)는 비선형 분수형 편미분방정식 $(-\Delta)^{\alpha/2}u + u^2 = g$ 에서 스펙트럴 수렴을 달성한다.
• $(-\Delta)^{\alpha/2} + x^2$ 의 고유값 문제에서 오차는 고유값의 대수적 감쇠와 일치하는 대수적 감쇠를 보이며, 지수 감쇠가 아니다.
• 과도 스케일링된 기저($e^{-x^2}H_n(x)$)는 $N$ 증가에 따라 조건수가 급격히 증가함을 보이며, 이는 대규모 문제에서 조건부 행렬의 필요성을 시사한다.
• 명시적 미분 행렬을 사용한 제안된 방법은 비선형 문제에 대한 직접성과 적용 가능성 면에서 주파수-공간 콜로케이션보다 뛰어나다.