[논문 리뷰] Stability and Error estimates of the SAV Fourier-spectral method for the Phase Field Crystal Equation

Xiaoli Li, Jie Shen|arXiv (Cornell University)|2019. 07. 17.

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

한 줄 요약

이 논문은 단계 필드 크리스탈(PFC) 방정식에 대해 시간에 대해 SAV(스칼라 보조 변수) 접근법을 적용하고 공간에 대해 스펙트럼 이산화를 사용하여 조건부 에너지 안정성, 선형 스킴을 달성하는 안정화된 스칼라 보조 변수(SAV) 푸리에-스펙트럼 방법을 제안한다. 이는 완전 이산 선형 스킴에 대한 첫 번째 엄밀한 최적 오차 추정을 제공하며, 시간에 대해 2차 수렴성과 공간에 대해 스펙트럼 정확도를 보여주며, 안정화가 큰 시간 간격에서 정확도를 유지하는 데 필수적임을 밝힌다. 안정화가 안정성에는 필요하지 않지만 정확도를 확보하기 위해 필수적이다.

ABSTRACT

We consider fully discrete schemes based on the scalar auxiliary variable (SAV) approach and stabilized SAV approach in time and the Fourier-spectral method in space for the phase field crystal (PFC) equation. Unconditionally energy stability is established for both first- and second-order fully discrete schemes. In addition to the stability, we also provide a rigorous error estimate which shows that our second-order in time with Fourier-spectral method in space converges with order $O(Δt^2+N^{-m})$, where $Δt$, $N$ and $m$ are time step size, number of Fourier modes in each direction, and regularity index in space, respectively. We also present numerical experiments to verify our theoretical results and demonstrate the robustness and accuracy of the schemes.

연구 동기 및 목표

• 비선형 스킴의 복잡성 문제를 해결하기 위해 단계 필드 크리스탈(PFC) 방정식에 대해 무조건 에너지 안정성, 선형 수치 스킴을 개발한다.
• 기존 PFC 수치 해법에서 선형 스킴에 대한 수렴 분석과 오차 추정의 부족 문제를 해결한다.
• SAV 프레임워크와 스펙트럼 공간 이산화를 사용하여 완전 이산 스킴에 대한 엄밀한 오차 추정을 수립한다.
• 큰 시간 간격에서 정확도를 유지하면서 안정성에 영향을 주지 않는 안정화의 역할을 조사한다.
• 수치 실험을 통해 제안된 스킴의 강건성과 정확도를 입증한다.

제안 방법

• PFC 방정식을 선형 시스템으로 재구성하기 위해 스칼라 보조 변수(SAV) 접근법을 적용하여 에너지 안정성 확보와 시간 이산화의 단순화를 달성한다.
• 주기적 경계 조건을 활용하여 푸리에-스펙트럼 방법을 공간 이산화에 적용하여 공간에서 스펙트럼 정확도를 확보한다.
• 큰 시간 간격에서의 정확도 향상을 위해 안정화 항을 추가하여 안정화된 SAV(S-SAV) 공식을 도입한다.
• 크랭크-니콜슨 또는 후진 차분 공식을 시간에 사용하여 완전 이산 스킴을 구성하며, 무조건 에너지 안정성을 보장한다.
• 에너지 안정성에 기반한 이산 해에 대한 균일한 경계를 도출하여 엄밀한 오차 분석을 가능하게 한다.
• 이산 에너지 안정성과 정(regularity) 가정에 기반한 오차 분석을 수행하여 최적 수렴 속도를 도출한다.

실험 결과

연구 질문

• RQ1 SAV 접근법을 PFC 방정식에 확장하여 최적 수렴 속도를 갖는 선형, 무조건 에너지 안정성 스킴을 도출할 수 있는가?
• RQ2 S-SAV 공식에서 안정화 항의 역할은 큰 시간 간격에서 정확한 결과를 얻기 위해 어떤가?
• RQ3 고차수 및 비선형 특성을 지닌 문제에 대해 완전 이산 선형 스킴에 대한 엄밀한 오차 추정을 유도할 수 있는가?
• RQ4 시간 간격 크기의 선택이 안정화를 적용했을 때와 아닐 때 SAV-푸리에-스펙트럼 스킴의 정확도에 어떤 영향을 미치는가?
• RQ5 스펙트럼 공간 이산화가 스킴의 전체 수렴 속도와 수치적 효율성에 어느 정도 기여하는가?

주요 결과

• 제안된 2차 SAV-푸리에-스펙트럼 스킴은 각각 시간에 대해 $ O( au^2 + N^{-m}) $, 공간에 대해 $ O( au^2 + N^{-m}) $ 의 최적 수렴 속도를 달성한다. 여기서 $ au $ 는 시간 간격이고 $ N $ 은 푸리에 모드의 수이다.
• 첫 번째 및 두 번째 차수의 완전 이산 스킴에 대해 무조건 에너지 안정성이 증명되었으며, 이는 시간 간격 크기와 관계없이 강건함을 보장한다.
• S-SAV 공식에서의 안정화 항은 안정성이나 수렴성에는 필수적이지 않지만, 큰 시간 간격에서 합리적인 정확도를 확보하기 위해 필수적이다.
• 수치 실험을 통해 시간에 대해 2차 수렴성과 공간에 대해 스펙트럼 수렴성이 확인되었으며, 이는 이론적 오차 추정을 검증한다.
• 안정화가 적용된 경우, 조밀한 결정 성장 역학, 예를 들어 나선형 패턴, 불량 결함, 인터페이스 운동 등 복잡한 구조를 큰 시간 간격에서도 정확히 포착한다.
• 이 방법은 과냉각된 액체의 상전이 및 에피택셜 성장 시뮬레이션에서 낮은 수치적 확산을 보이며 강건성과 효율성을 입증한다.