[논문 리뷰] Arbitrarily High-order Unconditionally Energy Stable Schemes for Gradient Flow Models Using the Scalar Auxiliary Variable Approach

Yuezheng Gong, Jia Zhao외 1명|arXiv (Cornell University)|2019. 07. 09.

Numerical methods for differential equations참고 문헌 46인용 수 6

한 줄 요약

이 논문은 에너지 안정성을 조건 없이 유지하는 고차원 시간정확도 스킴을 위한 새로운 가족인 고차원 스칼라 보조변수(HSAV) 방법을 제안한다. SAV 접근법과 고차원 시간적분기법을 조합함으로써, 기존 SAV 스킴보다 훨씬 큰 시간스텝을 허용하면서도 최적의 수렴률을 달성하고, 다양한 벤치마크 모델에서 조건 없이 에너지 안정성과 높은 정확도를 유지한다.

ABSTRACT

In this paper, we propose a novel family of high-order numerical schemes for the gradient flow models based on the scalar auxiliary variable (SAV) approach, which is named the high-order scalar auxiliary variable (HSAV) method. The newly proposed schemes could be shown to reach arbitrarily high order in time while preserving the energy dissipation law without any restriction on the time step size (i.e., unconditionally energy stable). The HSAV strategy is rather general that it does not depend on the specific expression of the effective free energy, such that it applies to a class of thermodynamically consistent gradient flow models arriving at semi-discrete high-order energy-stable schemes. We then employ the Fourier pseudospectral method for spatial discretization. The fully discrete schemes are also shown to be unconditionally energy stable. Furthermore, we present several numerical experiments on several widely-used gradient flow models, to demonstrate the accuracy, efficiency and unconditionally energy stability of the HSAV schemes. The numerical results verify that the HSAV schemes can reach the expected order of accuracy, and it allows a much larger time step size to reach the same accuracy than the standard SAV schemes.

연구 동기 및 목표

• 에너지 안정성을 보존하는 임의의 고차원 시간정확도 수치 스킴을 구성하기 위한 일반적인 프레임워크를 개발하는 것.
• 기존 SAV 스킴이 일반적으로 2차 정확도에 국한되어 있는 한계를 극복하여 더 높은 시간 수렴 차수를 가능하게 하는 것.
• 시간스텝 크기와 관계없이 조건 없이 에너지 안정성을 보장함으로써 효율적인 장시간 시뮬레이션을 가능하게 하는 것.
• 자유 에너지 함수의 구체적 형태에 대한 가정 없이도 열역학적으로 일관된 다수의 기울기 흐름 모델에 대해 SAV 접근법의 적용 범위를 넓히는 것.
• Allen-Cahn, Cahn-Hilliard, MBE 방정식과 같은 널리 사용되는 모델들에 대한 수치 실험을 통해 방법의 효율성과 정확도를 입증하는 것.

제안 방법

• HSAV 방법은 원래 기울기 흐름 모델을 스칼라 보조변수(SAV) 접근법을 통해 재구성하여, 자유 에너지가 새로운 변수들에 대해 이차형이 되도록 보조변수를 도입한다.
• SAV 재구성 시스템에 고차원 선형적임의 룬게쿠타 또는 다단계 스킴(예: 가우스-레지오드 또는 역차분 공식)을 적용하여 고차원 시간적분을 달성한다.
• 결과적으로 유도된 반연속 스킴은 수정된 자유 에너지 기능에 대해 조건 없이 에너지 안정성을 보장하며, 이는 이산 수준에서 에너지 감소를 보장한다.
• 공간 이산화에는 푸리에 가짜스펙트럴 방법을 사용하여 공간에서 스펙트럴 정밀도를 확보하고 FFT를 통한 효율적 계산을 가능하게 한다.
• 완전히 이산화된 스킴은 시간적분기법과 SAV 재구성의 구조보존 성질 덕분에 조건 없는 에너지 안정성 성질을 그대로 이어받는다.
• 이 방법은 자유 에너지 함수의 구체적 형태에 대한 가정이 필요 없어 다양한 기울기 흐름 모델에 적용 가능하다.

실험 결과

연구 질문

• RQ1 SAV 프레임워크와 고차원 시간적분을 조합하여 기울기 흐름 시뮬레이션에서 높은 정확도와 조건 없는 에너지 안정성을 동시에 달성할 수 있는가?
• RQ2 제안된 HSAV 방법이 Allen-Cahn, Cahn-Hilliard, MBE 모델을 포함한 다양한 기울기 흐름 모델에서 큰 시간스텝에서도 에너지 안정성을 유지하는가?
• RQ3 동일한 정확도 수준과 시간스텝 크기에서 HSAV 스킴의 수치 오차와 계산 비용은 기존 표준 SAV 스킴보다 어떻게 비교되는가?
• RQ4 HSAV 방법은 에너지 소산과 기대 수렴 속도를 유지하면서 얼마나 큰 시간스텝을 허용할 수 있는가?
• RQ5 HSAV 프레임워크는 시험된 예제를 넘어서 다른 열역학적으로 일관된 모델로 일반화될 수 있는가?

주요 결과

• HSAV 스킴은 시간 정밀도 실험에서 기대하는 정확도 순서를 달성하며, 4차 및 6차 스킴이 시간 정밀도 분석에서 최적의 수렴 속도를 보였다.
• 유사한 시간스텝에서 표준 SAV-CN 스킴에 비해 HSAV 스킴은 수치 오차가 수 개의 주자리만큼 작아서 극적으로 낮은 오차를 보였다.
• 경사 선택이 있는 MBE 모델의 경우, 6차 HSAV 스킴은 시간스텝 δt = 0.05에서도 에너지 감소를 정확히 포착했지만, SAV-CN 스킴은 유사한 정확도를 확보하기 위해 δt = 10⁻⁴ 가 필요했다.
• 시간스텝 δt = 0.025일 때 4차 HSAV 스킴의 총 CPU 시간은 79.45초였고, 시간스텝 δt = 0.0001인 SAV-CN 스킴은 672.39초였으며, 이는 5~8배의 속도 향상을 의미한다.
• 시간스텝 δt = 0.025인 4차 HSAV 스킴을 사용한 MBE 모델의 수치 해는 알려진 문헌 결과와 일치하는 형태의 일관성 있는 패턴을 생성하여 방법의 물리적 정확성을 검증했다.
• 모든 시험된 모델에서 HSAV 방법은 조건 없는 에너지 안정성을 유지하였으며, 시간스텝 크기와 관계없이 에너지가 단조롭게 감소하는 경향을 보였다.