[논문 리뷰] A stabilized second order exponential time differencing multistep method for thin film growth model without slope selection

Wenbin Chen, Weijia Li외 3명|arXiv (Cornell University)|2019. 07. 04.

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

한 줄 요약

이 논문은 인위적 안정화를 통합한 푸리에 가짜스펙트럴 공간 이산화와 ETD 기반 시간 적분을 사용하는 안정화된 두 번째 차수의 지수 시간 차분 다단계(sETDMs2) 방법을 제안한다. 이 방법은 무조건적인 장기 에너지 안정성을 확보하며, $\ell^{\infty}(0,T;\ell^2)$ 노름에서 $O(\tau^2)$ 수렴성을 보이며, 수치 실험을 통해 스케일링 법칙에서 $\lesssim 4\%$ 이내의 상대 오차로 정확한 장기 응고 동역학을 재현한다.

ABSTRACT

In this paper, a stabilized second order in time accurate linear exponential time differencing (ETD) scheme for the no-slope-selection thin film growth model is presented. An artificial stabilizing term $Aτ^2\frac{\partialΔ^2 u}{\partial t}$ is added to the physical model to achieve energy stability, with ETD-based multi-step approximations and Fourier collocation spectral method applied in the time integral and spatial discretization of the evolution equation, respectively. Long time energy stability and detailed $\ell^{\infty}(0,T; \ell^2)$ error analysis are provided based on the energy method, with a careful estimate of the aliasing error. In addition, numerical experiments are presented to demonstrate the energy decay and convergence rate.

연구 동기 및 목표

• 무한한 산기 경사가 존재하는 장기 응고를 보이는 '기울기 선택 없음' 박막 성장 모델에 대해 고차수이자 무조건적인 에너지 안정성 시간 적분 방법을 개발한다.
• 비선형 오목한 항에 대해 에너지 안정성 제2차 방법을 구성하는 데 도전하는 문제를 해결하기 위해 인위적 안정화 항을 도입한다.
• 에너지 방법을 사용하여 엄밀한 장기 에너지 안정성과 $\ell^{\infty}(0,T;\ell^2)$ 노름에서 $O(\tau^2)$ 수렴성을 확립한다.
• 에너지 감소, 수렴 속도, 표면 거칠기 및 기울기의 장기 스케일링 법칙을 정확히 포착하는 것을 통해 수치적으로 방법을 검증한다.

제안 방법

• 이산 설정에서 에너지 안정성을 확보하기 위해 물리적 PDE에 인위적 안정화 항 $A\tau^2 \frac{\partial \Delta^2 u}{\partial t}$ 를 도입한다.
• 시간 이산화에 두 번째 차수의 지수 시간 차분(ETD) 다단계 방법을 적용하여 선형 항은 정확히 통합하고 비선형 항은 명시적으로 처리한다.
• 공간 이산화에 푸리에 가짜스펙트럴 방법을 사용하여 높은 정확도와 스펙트럴 수렴성을 확보한다.
• 연속적인 에너지 감소 성질 $\frac{dE}{dt} = -\|u_t\|_{L^2}^2$ 를 모방하는 이산 에너지 기능을 구성한다.
• 에너지 방법을 사용한 세밀한 오차 분석을 수행하며, 비선형 항에서의 앨리아싱 오차를 신중하게 추정한다.
• 시간 간격을 적응적으로 조절하여 $t < 1000$ 에서 $\tau = 0.01$, $1000 \leq t < 2000$ 에서 $\tau = 0.02$, $t \geq 2000$ 에서 $\tau = 0.04$ 로 설정하고, 시간 간격 변경 후 이전 단계를 재사용한다.

실험 결과

연구 질문

• RQ1 기울기 선택 없음 박막 성장 모델에 대해 장기 에너지 안정성을 보장하는 두 번째 차수의 시간 정확한 방법을 구성할 수 있는가?
• RQ2 인위적 안정화 항을 추가함으로써 비선형 항을 명시적으로 처리하는 두 번째 차수의 ETD 방법에서 에너지 안정성이 달성되는가?
• RQ3 제안된 방법의 $\ell^{\infty}(0,T;\ell^2)$ 노름에서의 수렴 속도는 얼마이며, 어떤 시간 간격 제약 조건 하에서 성립하는가?
• RQ4 이 방법은 표면 거칠이와 기울기의 장기 응고 동역학, 특히 $h_{\mathcal{N}} \sim t^{1/2}$ 와 $m_{\mathcal{N}} \sim t^{1/4}$ 스케일링 법칙을 정확히 포착할 수 있는가?
• RQ5 장기 영역에서 인위적 안정화 파rameter $A$ 의 다양한 값에 대해 방법의 안정성은 어떻게 되는가?

주요 결과

• sETDMs2 방법은 기울기 선택 없음 박막 성장 모델에 대해 무조건적인 장기 에너지 안정성을 확보하며, 이산 에너지가 시간이 지남에 따라 단조롭게 감소한다.
• 에너지 방법과 앨리아싱 오차 통제를 통해 $\tau \leq \frac{1}{8}$ 일 때 $\ell^{\infty}(0,T;\ell^2)$ 노름에서 $O(\tau^2)$ 수렴성이 수학적으로 증명된다.
• 수치 실험 결과는 이론적 예측과 일치하는 에너지 감소와 수렴 속도를 확인하며, $\ell^2$ 오차 노름이 두 번째 차수 속도로 감소함을 보여준다.
• 평균 표면 거칠이 $h_{\mathcal{N}}$ 은 $t^{0.4865 \sim 0.4879}$ 로 스케일링되고, 평균 기울기 $m_{\mathcal{N}}$ 은 $t^{0.2609 \sim 0.2614}$ 로 스케일링되며, 이는 이론적 $t^{1/4}$ 와 $t^{1/2}$ 법칙과 매우 유사하다.
• 에너지 $E_{\mathcal{N}}$ 는 $-39.26 \sim -39.38 \ln t - 54.36$ 로 스케일링되며, 큰 $t$ 에서 이론적 $O(-\ln t)$ 감소와 일치한다.
• 다양한 안정화 파rameter $A$ 값(예: $A = 1/8$, $10^{-2}$, $\frac{2+\sqrt{3}}{6}$)을 사용한 결과 장기 스케일링 행동이 거의 동일하게 나타나, 방법의 강건성을 입증한다.