[논문 리뷰] Analysis of the L1 scheme for fractional wave equations with nonsmooth data

Binjie Li, Tao Wang외 1명|arXiv (Cornell University)|2019. 08. 24.

Fractional Differential Equations Solutions참고 문헌 36인용 수 5

한 줄 요약

이 논문은 초기 데이터가 부드럽지 않은 시간 분수형 파동 방정식에 대한 L1 스킴을 분석하여, 새로운 안정성 추정을 유도하고 표준 L1 스킴의 시간 수렴 차수 $\mathcal{O}(\tau^{3-\alpha})$ 를 증명한다. 수정된 L1 스킴은 $\mathcal{O}(\tau^2)$ 정확도를 달성하며, 균일한 시간 격자와 유계인 $\tau^\alpha / h_{\min}^2$ 조건 하에서 $L^2(\Omega)$-노름에서 두 스킴에 대해 수렴성을 확립한다. 결과는 수치적으로 검증된다.

ABSTRACT

This paper analyzes the well-known L1 scheme for fractional wave equations with nonsmooth data. A new stability estimate is obtained, and the temporal accuracy $ \mathcal O(τ^{3-α}) $ is derived for the nonsmooth initial data. In addition, a modified L1 scheme is proposed, and stability and temporal accuracy $ \mathcal O(τ^2) $ are derived for this scheme with nonsmooth initial data. The convergence of the two schemes in the inhomogeneous case is also established. Finally, numerical experiments are performed to verify the theoretical results.

연구 동기 및 목표

• 비부드러운 초기 자료를 가진 분수형 파동 방정식에 L1 스킴를 적용할 때 $L^2(\Omega)$-노름에서의 수렴 분석 부족 문제를 해결하기 위해.
• 해의 최소한의 정규성 가정 하에서 L1 스킴의 향상된 시간 정확도를 도출하기 위해.
• 비부드러운 자료에 대해 $\mathcal{O}(\tau^2)$ 수렴을 달성하는 수정된 L1 스킴을 제안하고 분석하기 위해.
• 비동차 케이스($f \not\equiv 0$)에서 두 스킴에 대한 수렴성을 확립하기 위해.
• 비율 $\tau^\alpha / h_{\min}^2$ 가 시간 정확도를 결정하는 데 핵심적인 역할을 하는지 규명하기 위해.

제안 방법

• 비부드러운 자료를 가진 분수형 파동 방정식에 적용된 L1 스킴에 대해 새로운 안정성 추정을 유도한다.
• 분석은 리만-리우빌 분수도함수 $\operatorname{D}_{0+}^{\alpha-1}$ 을 사용하며, 방정식 $\operatorname{D}_{0+}^{\alpha-1}(u' - u_1)(t) - \Delta u(t) = f(t)$ 에 적용한다.
• 분수도함수의 수정된 구적법을 사용하여 시간 정확도를 향상시킨 수정된 L1 스킴을 제안한다.
• 완전한 이산화에서는 시간에 대한 L1 또는 수정된 L1 스킴과 공간에 대한 표준 $P_1$-유한요소 방법을 조합한다.
• 균일한 시간 격자와 유계인 $\tau^\alpha / h_{\min}^2$ 조건 하에서 오차 추정을 도출한다. 여기서 $h_{\min}$ 은 최소 요소 직경이다.
• 다양한 $\alpha$ 와 메쉬 크기에서 두 스킴의 이론적 수렴 속도를 검증하기 위해 수치 실험을 수행한다.

실험 결과

연구 질문

• RQ1 비부드러운 초기 자료를 가진 분수형 파동 방정식에 대해 $L^2(\Omega)$-노름에서 표준 L1 스킴의 시간 수렴 차수는 무엇인가?
• RQ2 수정된 L1 스킴은 비부드러운 자료에 대해 $\mathcal{O}(\tau^2)$ 시간 정확도를 달성할 수 있는가?
• RQ3 비율 $\tau^\alpha / h_{\min}^2$ 는 L1 스킴의 시간 정확도에 어떤 영향을 미치는가?
• RQ4 비동차 케이스($f \not\equiv 0$)에서 L1 및 수정된 L1 스킴에 대한 수렴성이 확립되었는가?
• RQ5 비균일한 시간 격자 또는 변수 계수의 영향을 받아 L1 스킴의 오차 추정에 어떤 영향을 미치는가?

주요 결과

• 균일한 시간 격자와 유계인 $\tau^\alpha / h_{\min}^2$ 조건 하에서 비부드러운 초기 자료에 대해 표준 L1 스킴은 $L^2(\Omega)$-노름에서 시간 수렴 차수 $\mathcal{O}(\tau^{3-\alpha})$ 를 달성한다.
• 새로운 안정성 추정이 도출되었으며, 이는 $\mathcal{O}(\tau^{3-\alpha})$ 수렴 속도를 증명하는 데 핵심적이다.
• 수정된 L1 스킴은 비부드러운 자료에 대해 $\mathcal{O}(\tau^2)$ 시간 정확도를 달성하여 표준 L1 스킴의 $\mathcal{O}(\tau^{3-\alpha})$ 수렴 속도를 향상시킨다.
• 두 스킴의 수렴성이 비동차 케이스에서 확립되었으며, 오차 추정은 소스 항 $f$ 의 정규성에 따라 달라진다.
• 수치 실험을 통해 예측된 수렴 차수를 확인하였다: 표준 L1 스킴의 경우 $\alpha \leq 3/2$ 에서 $\mathcal{O}(h^2)$, $\alpha > 3/2$ 에서는 $\mathcal{O}(h^{6/\alpha - 2})$ 이며, 수정된 L1 스킴의 경우 모든 $\alpha \in (1,2)$ 에 대해 $\mathcal{O}(h^2)$ 이다.
• 실험 결과, $\tau^\alpha / h_{\min}^2$ 가 균일하게 유계가 아닐 경우 두 스킴의 시간 정확도가 떨어지며, 이는 이론적 조건이 확인됨을 시사한다.