[논문 리뷰] Mesh-robustness of the variable steps BDF2 method for the Cahn-Hilliard model.

Hong-lin Liao, Bingquan Ji외 2명|arXiv (Cornell University)|2021. 02. 07.

Solidification and crystal growth phenomena참고 문헌 27인용 수 4

한 줄 요약

이 논문은 사용자 정의한 단계 비율 제한 $ r_{\text{user}} < 4.864 $ 하에서 에너지 소산과 메시 기반 안정성을 유지하는 가변 단계 BDF2 방법을 제안한다. 이는 이산 수직 컨볼루션 커널과 새로운 임bedding 부등식을 활용하여 시간 간격의 변화와 무관한 $ L^2 $-노름 오차 추정을 달성함으로써 메시 기반 수렴성을 보장하고, 적응형 시간 스텝을 통한 효율적인 장기 시뮬레이션을 가능하게 한다.

ABSTRACT

The two-step backward differential formula (BDF2) implicit method with unequal time-steps is investigated for the Cahn-Hilliard model by focusing on the numerical influences of time-step variations. The suggested method is proved to preserve a modified energy dissipation law at the discrete levels if the adjoint time-step ratios fulfill a new step-ratio restriction $0<r_k:= au_k/ au_{k-1}\le r_{\mathrm{user}}$ ($r_{\mathrm{user}}$ can be chosen by the user such that $r_{\mathrm{user}}<4.864$), such that it is mesh-robustly stable in an energy norm. We view the BDF2 formula as a convolution approximation of the first time derivative and perform the error analysis by using the recent suggested discrete orthogonal convolution kernels. By developing some novel convolution embedding inequalities with respect to the orthogonal convolution kernels, an $L^2$ norm error estimate is established at the first time under the updated step-ratio restriction $0<r_k\le r_{\mathrm{user}}$.The time-stepping scheme is mesh-robustly convergent in the sense that the convergence constant (prefactor) in the error estimate is independent of the adjoint time-step ratios. On the basis of ample tests on random time meshes, a useful adaptive time-stepping strategy is applied to efficiently capture the multi-scale behaviors and to accelerate the long-time simulation approaching the steady state.

연구 동기 및 목표

• 가변 단계 BDF2 스킴에서 Cahn-Hilliard 모델의 에너지 안정성과 수렴성을 유지하는 데 도전하는 문제를 해결하기 위해.
• 시간 간격의 변화가 장기 시뮬레이션에서 수치 안정성과 정확성에 미치는 영향을 분석하기 위해.
• 시간 간격 비율에 의존하지 않는 오차 한계를 갖는 메시 기반 수렴 이론을 수립하기 위해.
• 다양한 스케일의 역학을 효율적으로 포착할 수 있는 실용적인 적응형 시간 스텝 전략을 개발하기 위해.

제안 방법

• BDF2 방법은 가변 시간 간격을 사용하여 첫 번째 시간 도함수를 컨볼루션 근사로 설정한다.
• 이산 수직 컨볼루션 커널을 사용하여 시간 이산화 프레임워크 내에서 정밀한 오차 분석을 수행한다.
• 수직 컨볼루션 커널을 위한 새로운 컨볼루션 임베딩 부등식을 유도하여 시간 오차 성장을 제한한다.
• 에너지 소산을 이산 수준에서 유지하기 위해 $ 0 < r_k := \Delta t_k / \Delta t_{k-1} \leq r_{\text{user}} $ 이며 $ r_{\text{user}} < 4.864 $ 인 단계 비율 제한을 도입한다.
• 이 방법은 수정된 에너지 소산 법칙을 만족함으로써 에너지 노름에서 안정성을 보장한다.
• 임의의 시간 메시에서의 수치 실험 기반으로 적응형 시간 스텝 전략을 구현하여 안정 상태에 수렴하는 데 수렴 속도를 향상시킨다.

실험 결과

연구 질문

• RQ1 가변 단계 BDF2 방법이 임의의 시간 간격 변화에 대해 Cahn-Hilliard 모델의 에너지 소산을 유지할 수 있는가?
• RQ2 어떤 단계 비율 제한 조건이 에너지 및 $ L^2 $ 노름에서 메시 기반 안정성과 수렴성을 보장하는가?
• RQ3 이산 수직 컨볼루션 커널을 어떻게 활용하여 비균일 시간 메시에서 정밀한 오차 추정을 도출할 수 있는가?
• RQ4 수렴 상수는 시간 간격 비율에 얼마나 의존하는가? 그리고 이를 비율에 무관하게 만들 수 있는가?
• RQ5 적응형 시간 스텝 전략은 다중 스케일 역학을 해석하고 장기 시뮬레이션을 가속화하는 데 얼마나 효과적인가?

주요 결과

• 사용자 정의한 단계 비율 제한 $ r_k \leq r_{\text{user}} < 4.864 $ 하에서, 가변 시간 간격을 사용하는 BDF2 방법은 이산 수준에서 수정된 에너지 소산 법칙을 유지하여 에너지 안정성을 보장한다.
• 업데이트된 단계 비율 조건 하에서 첫 번째 시간 단계에서 새로운 $ L^2 $-노름 오차 추정이 확립되었으며, 수렴 상수는 시간 간격 비율에 영향을 받지 않는다.
• 수렴은 메시 기반 안정성이며, 오차 한계의 계수는 이웃 시간 간격 비율에 의존하지 않는다.
• 이중 이산 수직 컨볼루션 커널과 유도된 임베딩 부등식의 활용으로 $ L^2 $ 노름에서 최적의 수렴 속도를 달성한다.
• 임의의 시간 메시에서의 수치 실험은 다중 스케일 행동을 포착하는 데 있어 적응형 시간 스텝 전략의 강건성과 효율성을 확인한다.
• 이 방법은 안정성과 정확성을 유지하면서도 시간 간격을 동적으로 조정함으로써 장기 시뮬레이션을 가속화할 수 있다.