[논문 리뷰] $H^1$-norm stability and convergence of an L2-type method on nonuniform meshes for subdiffusion equation

Chaoyu Quan, Xu Wu|arXiv (Cornell University)|2022. 05. 12.

Differential Equations and Numerical Methods인용 수 4

한 줄 요약

이 논문은 일반적인 비균일 메esh에서 하위확산 방정식에 대한 L2형 방법에 대해 H¹-노름 안정성과 수렴성을 확립한다. 시간 간격 비율 제약 조건(0.457 ≤ ρₖ ≤ 3.562)이 만족될 경우, 핵심 이항형식의 양의 준정의성을 증명하여, r > 5/α − 1인 수정된 등급 메쉬에서 H¹-노름으로서 (5−α)/2의 수렴 차수를 갖는 장기적 H¹-안정성과 수렴성을 가능하게 한다.

ABSTRACT

This work establishes $H^1$-norm stability and convergence for an L2 method on general nonuniform meshes when applied to the subdiffusion equation. Under mild constraints on the time step ratio $ρ_k$, such as $0.4573328\leq ρ_k\leq 3.5615528$ for $k\geq 2$, the positive semidefiniteness of a crucial bilinear form associated with the L2 fractional-derivative operator is proved. This result enables us to derive long time $H^1$-stability of L2 schemes. These positive semidefiniteness and $H^1$-stability properties hold for standard graded meshes with grading parameter $15/α-1$. To the best of our knowledge, this study is the first work on $H^1$-norm stability and convergence of L2 methods on general nonuniform meshes for the subdiffusion equation.

연구 동기 및 목표

• 일반적인 비균일 시간 메쉬에서 하위확산 방정식에 대한 L2형 스킴의 H¹-노름 안정성과 수렴성을 확립하기 위해.
• 특히 낮은 정규성 해를 가진 경우에 비균일 메쉬에서 L2 방법의 H¹-안정성 분석 부족 문제를 다루기 위해.
• 시간 간격 비율 제약 조건 하에서 L2 분수도 도함수 연산자와 관련된 이항형식의 양의 준정의성을 증명하기 위해.
• 특정 등급 매개변수를 가진 등급 메쉬와 수정된 등급 메쉬로 안정성 및 수렴 결과를 확장하기 위해.
• 엄밀한 H¹-노름 오차 분석을 제공하고 수치적으로 수렴 속도를 확인하기 위해.

제안 방법

• k ≥ 2일 때 시간 간격 비율 제약 조건 ρₖ ∈ [0.4573328, 3.5615528] 하에서 이항형식 $\mathcal{B}_n(v,w) = \sum_{k=1}^n \langle L_k^\alpha v, \delta_k w \rangle$ 의 양의 준정의성을 증명한다.
• 장기적 H¹-안정성을 도출하기 위해 이산 분수도 그로워얼 불등식과 에너지 기법을 기반으로 한 새로운 분석 프레임워크를 사용한다.
• 표준 등급 메쉬에 분석을 적용하여 grading 매개변수 $r \leq 3.2016538$ 일 때 H¹-안정성을 확립한다.
• 일반적인 비균일 메쉬에 대해 H¹-노름 오차 추정을 유도하고, r > 5/α − 1일 때 수정된 등급 메쉬에서 수렴 차수 (5 − α)/2 를 증명한다.
• 고정밀 공간 해석을 위해 체비셰프–구아스–로바토 점들을 사용한 스펙트럴 콜로케이션 기법을 적용한다.
• 최대 H¹-오차를 측정하여 다양한 α, r, N 값에 대해 수정된 등급 메쉬에서 이론적 결과를 수치 실험으로 검증한다.

실험 결과

연구 질문

• RQ1 일반적인 비균일 메쉬에서 시간 간격 비율 제약 조건 하에서 L2 분수도 도함수 연산자와 관련된 이항형식이 양의 준정의인가?
• RQ2 하위확산 방정식에 대해 비균일 메쉬에서 암시적 L2 스킴의 장기적 H¹-안정성이 확립될 수 있는가?
• RQ3 수정된 등급 메쉬에서 L2 방법의 H¹-노름에서 최적 수렴 속도는 무엇이며, grading 매개변수 r에 어떻게 의존하는가?
• RQ4 수치 실험에서 이항형식의 양의 준정의성 결과가 여전히 날카로운가?
• RQ5 이론적 H¹-노름 수렴 속도 (5 − α)/2 는 실질적으로 관찰될 수 있으며, 관측된 수렴 차수와 일치하는가?

주요 결과

• k ≥ 2일 때 0.4573328 ≤ ρₖ ≤ 3.5615528 라는 온화한 제약 조건 하에서 이항형식 $\mathcal{B}_n$ 이 양의 준정의임을 증명하여 H¹-안정성 분석이 가능해진다.
• 일반적인 비균일 메쉬에서 동일한 경계 조건을 가진 하위확산 방정식에 대해 L2 스킴의 장기적 H¹-안정성이 확립된다.
• 표준 등급 메쉬의 경우 grading 매개변수 1 < r ≤ 3.2016538 를 만족할 때 H¹-안정성이 성립한다.
• r > 5/α − 1 일 때 수정된 등급 메쉬에서 H¹-노름으로서 (5 − α)/2 의 수렴 차수를 증명한다.
• 수치 결과는 H¹-노름에서 관측된 수렴 차수가 이론적 기대와 일치하는 $\min\{r\alpha, 3 - \alpha\}$ 약간의 값을 갖는 것으로 확인된다.
• λ_{\min} ≥ C N^{\alpha} 인 최소 고유값의 양의 정의성 경계는 수치적으로 최적임을 검증하였으며, λ_{\max} ≥ C N^{r\alpha} 이다.