[논문 리뷰] Asymptotically compatible energy of variable-step fractional BDF2 formula for time-fractional Cahn-Hilliard model

Hong-lin Liao, Nan Liu외 1명|arXiv (Cornell University)|2022. 10. 22.

Fractional Differential Equations Solutions인용 수 6

한 줄 요약

이 논문은 Caputo 도수 $\alpha \in (0,1)$의 시간분수 Cahn-Hilliard 모델을 위한 가변단계 분수 BDF2 스킴을 제안한다. 분수 BDF2 공식의 국소-비국소 분할을 통한 새로운 이산 기울기 구조를 도입하여, 약한 단계비 제약 하에 이산 에너지 소산 법칙을 수립하고, $\alpha \to 1^{-}$일 때 고전적 Cahn-Hilliard 에너지 법칙과 점점 가까워지는 점근적 호환성을 입증한다. 이는 적응형 시간 스텝을 통해 장시간 안정성과 정확성을 확보한다.

ABSTRACT

A new discrete energy dissipation law of the variable-step fractional BDF2 (second-order backward differentiation formula) scheme is established for time-fractional Cahn-Hilliard model with the Caputo's fractional derivative of order $α\in(0,1)$, under a weak step-ratio constraint $0.4753\le τ_k/τ_{k-1}

연구 동기 및 목표

• 일반적인 비균일 시간 메esh를 사용하는 시간분수 Cahn-Hilliard 모델에 대해 안정적이고 두 번째 차수의 시간 이산화를 개발한다.
• 가변 시간 스텝 하에서 물리적 에너지 감쇠 성질을 유지하는 이산 에너지 소산 법칙을 수립한다.
• $\alpha \to 1^{-}$일 때 고전적 Cahn-Hilliard 에너지 법칙과 점근적 호환성을 확보한다.
• 초기 시점에서의 급격한 변화와 이후의 응집 단계를 잘 포착하는 적응형 시간 스텝 전략을 통해 효율적인 장시간 시뮬레이션을 가능하게 한다.

제안 방법

• 분수 BDF2 공식을 국소-비국소 분할하여, 고전적 이중단계 후진미분과 유사한 국소 부분과 Caputo 도수의 L1 공식과 유사한 비국소 부분으로 나눈다.
• 이 분할 기반으로 새로운 이산 기울기 구조를 구성하여 일致성 있고 안정적인 에너지 소산 법칙을 유도한다.
• 안정성과 에너지 감쇠를 확보하기 위해 약한 단계비 제약 조건인 $0.4753 \leq \tau_k / \tau_{k-1} < r^*(\alpha)$를 유도하며, $r^*(\alpha) \geq 4.660$이다.
• 초기 특이성과 동적 변화를 해결하기 위해 $\tau_{\text{ada}} = \max\{\tau_{\min}, \tau_{\max}/\Pi_\phi\}$를 사용하는 적응형 시간 스텝 전략을 구현하며, $\Pi_\phi = \sqrt{1 + \eta \|\partial_\tau \phi^n\|^2}$이다.
• 공간 이산화에 푸리에 가짜스펙트럴 방법을 사용하고, $128 \times 128$ 그리드를 적용하며, $\varepsilon = 10^{-12}$로 분수 BDF2 공식의 빠른 평가를 위해 SOE 기법을 활용한다.

실험 결과

연구 질문

• RQ1 시간분수 Cahn-Hilliard 모델에 대해 이산 에너지 소산 법칙을 유지하는 두 번째 차수의 가변단계 분수 BDF2 스킴을 구성할 수 있는가?
• RQ2 이산 에너지 법칙이 $\alpha \to 1^{-}$일 때 고전적 Cahn-Hilliard 에너지 소산 법칙과 점점 가까워지는가?
• RQ3 이 스킴은 계단형 메쉬와 적응형 메쉬에서 최적의 시간 정확도 $O(\tau^{\min\{\gamma, 3 - \alpha\}})$를 유지하는가?
• RQ4 적응형 시간 스텝 전략은 초기 특이성과 장시간 응집 동역학을 얼마나 효과적으로 포착하는가?

주요 결과

• 제안된 가변단계 분수 BDF2 스킴은 계단형 및 적응형 메쉬에서 $O(\tau^{\min\{\gamma, 3 - \alpha\}})$의 두 번째 차수 시간 정확도를 확보하며, $\alpha = 0.4$와 $\alpha = 0.7$에 대해 수치적으로 확인되었다.
• 약한 단계비 제약 조건 $0.4753 \leq \tau_k / \tau_{k-1} < r^*(\alpha)$ 하에 새로운 이산 에너지 소산 법칙이 수립되었으며, $\alpha \in (0,1)$에 대해 $r^*(\alpha) \geq 4.660$이다.
• 이산 에너지 및 에너지 소산 법칙은 $\alpha \to 1^{-}$일 때 고전적 Cahn-Hilliard 모델과 점점 가까워지며, 정수계수 극한에서 일致성을 보장한다.
• 적응형 시간 스텝 전략을 $\eta = 10^3$로 적용하여 초기 단계에서의 급격한 에너지 감쇠와 장시간에 걸친 응집 동역학을 성공적으로 포착하였으며, 시간 스텝이 해의 기울기에 따라 동적으로 조정됨을 확인하였다.
• 수치 결과는 응집 속도와 상호작용 진화가 $\alpha$에 의해 크게 영향을 받으며, $\alpha \to 1^{-}$일 때 고전적 $O(t^{-1/3})$ 속도에 점점 수렴함을 보여주었다.