[논문 리뷰] A Second-Order Energy Stable Backward Differentiation Formula Method for the Epitaxial Thin Film Equation with Slope Selection

Wenqiang Feng, Cheng Wang외 2명|arXiv (Cornell University)|2017. 06. 06.

Fluid Dynamics and Thin Films참고 문헌 17인용 수 4

한 줄 요약

이 논문은 경사 선택을 고려한 에피택셜 박막 방정식에 대해 두 번째 차수의 에너지 안정적인 역차분 공식(BDF) 유한차분 스킴을 제안한다. 이 스킴은 일정 계수를 가진 안정화 항을 도입하여 무조건적인 에너지 안정성을 확보한다. 이 방법은 최적 수렴성을 달성하며, 비선형 4라플라시안 시스템을 효율적으로 해결하기 위한 전처리된 반복 해법을 사용한다. 수치 시뮬레이션을 통해 장기간에 걸친 안정성과 정확성이 입증되었다.

ABSTRACT

In this paper, we study a novel second-order energy stable Backward Differentiation Formula (BDF) finite difference scheme for the epitaxial thin film equation with slope selection (SS). One major challenge for the higher oder in time temporal discretization is how to ensure an unconditional energy stability and an efficient numerical implementation. We propose a general framework for designing the higher order in time numerical scheme with unconditional energy stability by using the BDF method with constant coefficient stabilized terms. Based on the unconditional energy stability property, we derive an $L^\infty_h (0,T; H_{h}^2)$ stability for the numerical solution and provide an optimal the convergence analysis. To deal with the 4-Laplacian solver in an $L^{2}$ gradient flow at each time step, we apply an efficient preconditioned steepest descent algorithm and preconditioned nonlinear conjugate gradient algorithm to solve the corresponding nonlinear system. Various numerical simulations are present to demonstrate the stability and efficiency of the proposed schemes and slovers.

연구 동기 및 목표

• 경사 선택을 고려한 에피택셜 박막 방정식에 대해 두 번째 차수의 시간 이산화를 개발하여 무조건적인 에너지 안정성을 보장한다.
• BDF 프레임워크에 안정화된 일정 계수 항을 도입하여 고차수 시간 스킴에서 에너지 안정성을 유지하는 데 도전 과제를 해결한다.
• 수치 해의 $L^{ty}_{h}(0,T;H_{h}^{2})$ 안정성을 보장하는 엄밀한 수렴 분석을 제공한다.
• 매 시간 단계에서 발생하는 비선형 4라플라시안 시스템을 효율적으로 해결하기 위한 반복 해법을 설계한다. 전처리된 최강하강법과 비선형 공액 기울기 방법을 활용한다.
• 스펙트럼 방법과 빠른 푸리에 변환을 사용한 광범위한 수치 시뮬레이션을 통해 스킴의 장기 안정성과 정확성을 입증한다.

제안 방법

• 경사 선택을 고려한 에피택셜 박막 방정식에 대해 일정 계수 안정화 항을 포함한 두 번째 차수 BDF 스킴을 수립하여 무조건적인 에너지 안정성을 확보한다.
• 비선형 항을 명시적인 시간 분할이나 외삽 없이 처리하기 위해 BDF 방법과 안정화를 융합한 새로운 에너지 안정성 프레임워크를 적용한다.
• 주기적 경계 조건을 가진 유한차분 공간 이산화를 사용하고, 고차수 정확도를 확보하기 위해 푸리에 스펙트럼 방법을 적용한다.
• 매 시간 단계에서 발생하는 4라플라시안 항에 기인한 비선형 시스템을 해결하기 위해 전처리된 최강하강 알고리즘과 전처리된 비선형 공액 기울기 방법을 구현한다.
• 반복 해법에서 선형 시스템을 빠르게 해결하기 위해 빠른 푸리에 변환(FFT)을 활용하여 계산 효율성을 향상시킨다.
• 무조건적인 에너지 안정성과 스킴의 일致성에 기반하여 $L^{ty}_{h}(0,T;H_{h}^{2})$ 노름에서 최적 수렴 속도를 도출한다.

실험 결과

연구 질문

• RQ1 경사 선택을 고려한 에피택셜 박막 방정식에 대해 무조건적인 에너지 안정성을 보장하는 두 번째 차수의 BDF 스킴을 구성할 수 있는가?
• RQ2 명시적인 시간 분할이나 큰 안정화 계수에 의존하지 않고, BDF 프레임워크에 안정화된 일정 계수 항을 어떻게 통합하여 에너지 안정성을 유지할 수 있는가?
• RQ3 에너지 안정성 조건 하에서 제안된 스킴의 $L^{ty}_{h}(0,T;H_{h}^{2})$ 노름에서 최적 수렴 속도는 얼마인가?
• RQ4 매 시간 단계에서 비선형 4라플라시안 시스템을 다루는 데 있어 전처리된 반복 해법의 효율성은 어떠한가?
• RQ5 제안된 스킴은 박막 에피택시의 느린 역학을 시뮬레이션하기 위해 요구되는 장기 안정성과 정확성을 유지하는가?

주요 결과

• 제안된 두 번째 차수 BDF 스킴은 일정 계수 항을 포함한 새로운 안정화 프레임워크를 통해 무조건적인 에너지 안정성을 달성한다.
• 무조건적인 에너지 안정성과 일致성에 기반하여 $L^{ty}_{h}(0,T;H_{h}^{2})$ 노름에서 최적 수렴 속도가 증명되었으며, 이는 스킴의 정확성과 강건성을 확인한다.
• 스킴은 무조건적으로 안정하고 장기간의 시뮬레이션 동안 에너지 감쇠를 유지하며, 물리적 예측과 일치하는 $t^{-1/3}$ 에너지 감쇠 경향을 반영한다.
• 전처리된 최강하강법과 비선형 공액 기울기 해법은 매 시간 단계에서 4라플라시안 시스템을 해결하는 데 있어 효율성을 크게 향상시킨다.
• 수치 시뮬레이션을 통해 스킴의 안정성, 정확성, 효율성이 확인되었으며, 장시간에 걸쳐도 최소한의 정확도 저하가 관찰되었다.
• 특히 수정된 크랭크-니콜슨 접근법에서 관찰된 열악한 성능을 피하기 위해, 현존하는 두 번째 차수 스킴들에 비해 에너지 안정성과 장기적 행동에서 뛰어난 성능을 보였다.