[논문 리뷰] A symmetric fractional-order reduction method for direct nonuniform approximations of semilinear diffusion-wave equations

Pin Lyu, Seakweng Vong|arXiv (Cornell University)|2021. 01. 24.

Fractional Differential Equations Solutions참고 문헌 35인용 수 4

한 줄 요약

이 논문은 일반적인 비균일 시간 메esh에서 반선형 시간 분수형 확산-파동 방정식에 대해 고차 정확도 수치 해법을 구축하기 위해 대칭 분수형 순환(SFOR) 방법을 제안한다. 원래 문제를 순서 $\alpha/2$의 분수형 도함수를 갖는 결합된 시스템으로 줄임으로써, 저자들은 약간의 메쉬 및 정규성 조건 하에서도 최적의 $\mathcal{O}(\tau^2)$ 시간 수렴성을 달성하는 선형화된 L1 및 아리카노 스킴을 개발하였다. 이는 적응형 시간 스텝 및 수치 실험을 통해 검증되었다.

ABSTRACT

We introduce a symmetric fractional-order reduction (SFOR) method to construct numerical algorithms on general nonuniform temporal meshes for semilinear fractional diffusion-wave equations. By using the novel order reduction method, the governing problem is transformed to an equivalent coupled system, where the explicit orders of time-fractional derivatives involved are all $α/2$ $(1

연구 동기 및 목표

• 비균일 시간 메쉬에서 약한 초기 특이성을 갖는 반선형 확산-파동 방정식에 대해 최적 수렴성을 달성하는 데 도전하는 문제를 해결한다.
• 초기 특이성으로 인해 해의 정규성이 낮더라도 고차 정확도를 유지하는 강력한 수치 해법을 개발한다.
• L1 및 아리카노 스킴 기반의 적응형 시간 스텝 전략을 설계하여 효율성과 정확도를 향상시킨다.
• 최소한의 정규성 및 메쉬 제약 조건 하에서 $H^2$ 에너지 방법을 사용하여 엄밀한 수렴 분석을 수립한다.
• 이중 분수형 도함수를 갖는 균형 잡힌 시스템으로 원래의 미분 방정식을 단순화하는 대칭 분수형 순환 기법을 제공한다.

제안 방법

• 원래의 $\alpha$-순서 시간 분수형 확산-파동 방정식을 $\alpha/2$-순서 도함수를 갖는 등가 결합 시스템으로 변환하기 위해 대칭 분수형 순환(SFOR) 방법을 제안한다.
• 일반적인 비균일 시간 메쉬에서 감소된 시스템을 이산화하기 위해 선형화된 L1 스킴과 아리카노 스킴을 적용한다.
• 약한 메쉬 등급 및 정규성 가정 하에서 최적의 $\mathcal{O}(\tau^2)$ 수렴 속도를 유도하기 위해 $H^2$ 에너지 방법을 사용한다.
• 시간 스텝을 동적으로 조정하기 위해 (빠른) 선형화된 L1 및 아리카노 스킴 기반의 적응형 시간 스텝 전략을 도입한다.
• L1 및 아리카노 스킴에 대해 전역 일致성 오차 추정을 사용하여 시간 이산화 오차를 정의하고 분석한다.
• 이산화 컨볼루션 커널 $P_{n-j}^{(n)}$ 및 $A_{j-k}^{(j)}$를 포함하는 합계 부등식을 통해 오차 한계를 수립한다.

실험 결과

연구 질문

• RQ1 비균일 메쉬에서 반선형 확산-파동 방정식에 대해 정확도와 안정성을 향상시키기 위해 대칭 분수형 순환 방법을 개발할 수 있는가?
• RQ2 약한 정규성 및 메쉬 제약 조건 하에서 선형화된 L1 및 아리카노 스킴이 이러한 문제에 적용되었을 때 달성 가능한 최적 수렴 속도는 무엇인가?
• RQ3 계산 비용을 줄이면서도 높은 정확도를 유지하기 위해 어떻게 적응형 시간 스텝 전략을 설계할 수 있는가?
• RQ4 이격 오차 추정치가 메쉬 등급 파라미터 $\gamma$ 및 해의 정규성 파라미터 $\sigma_1, \sigma_2$에 얼마나 의존하는가?
• RQ5 $H^2$ 에너지 방법이 비균일하고 비균일적으로 등급이 매겨진 메쉬에서 분수형 초확산 문제에 대해 최적 수렴 결과를 도출하는 데 효과적으로 적용될 수 있는가?

주요 결과

• 선형화된 L1 스킴은 조건 $\gamma\sigma_1 \geq 2$ 하에서 $\mathcal{O}(\tau^2)$ 수렴 속도를 달성하며, 오차 한계는 $\sum_{j=1}^{n}P_{n-j}^{(n)}\|\Delta_h({\cal T}_f)^{n-\theta}\| \leq C\tau^{\min\{2\gamma\sigma_1,2\}}$ 로 주어진다.
• 아리카노 스킴은 $\gamma\sigma_1 \geq 3 - \beta$ 조건 하에서 $\mathcal{O}(\tau^2)$ 수렴 속도를 달성하며, 오차 한계는 $\sum_{j=1}^{n}P_{n-j}^{(n)}\|\Delta_h({\cal T}_u)^{j-\theta}\| \leq C\tau^{\min\{3-\beta,\gamma\sigma_1\}}$ 로 주어진다.
• 공간 오차는 $\|{\cal S}^n\| \leq C_u h^2$ 로 유계이며, 이는 이차 공간 수렴을 확인한다.
• L1 스킴의 경우, 정규성 가정 $|g''(t)| \leq C_g(1 + t^{\sigma-2})$ 하에서 시간 이격 오차는 $C\tau^{\min\{2-\beta,\gamma\sigma_2\}}$ 로 유계이다.
• 아리카노 스킴의 경우, 가정 $|g^{\prime\prime\prime}(t)| \leq C_g(1 + t^{\sigma-3})$ 하에서 이격 오차는 $C\tau^{\min\{3-\beta,\gamma\sigma_1\}}$ 로 유계이며, 추가로 $\tau_1^{\sigma}$ 및 $\tau_k^3$ 항이 포함된다.
• L1 및 아리카노 스킴 기반의 적응형 시간 스텝 전략은 수치 예제를 통해 정확도와 효율성을 유지하는 데 성공적으로 기여하였다.