[논문 리뷰] A new class of efficient and robust energy stable schemes for gradient flows

Jie Shen, Jie Xu외 1명|arXiv (Cornell University)|2017. 10. 03.

Fluid Dynamics and Turbulent Flows참고 문헌 39인용 수 39

한 줄 요약

이 논문은 경사 하강 흐름을 위한 조건부 에너지 안정성, 선형성, 고차수 수치적 스킴을 구축하기 위해 새로운 스칼라 보조 변수(SAV) 방법을 제안한다. 비선형 항을 스칼라 보조 변수를 통해 분리함으로써, 일정 계수를 가진 선형 시스템을 효율적으로 풀 수 있게 되어, 조건부로 제2차 정확도를 달성하고, 높은 차수의 BDF 스킴을 안정적으로 적용할 수 있으며, 특히 적응형 시간 스텝과 조합될 경우 뛰어난 성능을 발휘한다.

ABSTRACT

We propose a new numerical technique to deal with nonlinear terms in gradient flows. By introducing a scalar auxiliary variable (SAV), we construct efficient and robust energy stable schemes for a large class of gradient flows. The SAV approach is not restricted to specific forms of the nonlinear part of the free energy, and only requires to solve {\it decoupled} linear equations with {\it constant coefficients}. We use this technique to deal with several challenging applications which can not be easily handled by existing approaches, and present convincing numerical results to show that our schemes are not only much more efficient and easy to implement, but can also better capture the physical properties in these models. Based on this SAV approach, we can construct unconditionally second-order energy stable schemes; and we can easily construct even third or fourth order BDF schemes, although not unconditionally stable, which are very robust in practice. In particular, when coupled with an adaptive time stepping strategy, the SAV approach can be extremely efficient and accurate.

연구 동기 및 목표

• 경사 하강 흐름을 해결하기 위한 강력하고 효율적인 수치적 프레임워크를 개발하여 에너지 안정성을 유지한다.
• 기존의 볼록 분할법과 안정화 방법의 한계를 극복하여, 매 시간 단계에서 비선형 해법이 필요하거나 제한적인 시간 스텝이 요구되는 문제를 해결한다.
• 에너지 안정성과 선형 시스템의 구조를 유지하면서 제2, 제3, 제4차 정확도를 달성할 수 있도록 한다.
• 비국소적 또는 다중 적분 자유 에너지 항을 포함한 복잡한 모델로의 적용 가능성을 확장한다.
• 고정 계수를 가진 분리된 선형 시스템을 통해 실용적인 구현을 용이하게 한다.

제안 방법

• 자유 에너지의 비선형 항의 제곱근으로 정의된 스칼라 보조 변수(SAV)를 도입하여 원래 시스템을 등가 시스템으로 변환한다.
• SAV와 원래 변수가 분리된 새로운 시스템으로 경사 하강 흐름을 재구성함으로써, 일정 계수를 가진 선형 방정식을 도출한다.
• SAV 방법을 이용해 조건부로 제2차 에너지 안정성 스킴을 구축하고, 이를 고차수 BDF 스킴으로 확장한다.
• 디지털 에너지 법칙이 연속적 법칙을 그대로 반영하도록 증명함으로써 에너지 안정성을 확보하며, SAV는 비선형 에너지 기여를 포괄한다.
• 매 시간 단계에서 비선형 해법을 피하고 고정 행렬을 가진 선형 시스템만을 풀어 계산 효율성을 높인다.
• SAV 프레임워크에 적응형 시간 스텝을 통합하여 장시간 시뮬레이션에서 효율성과 정확도를 향상시킨다.

실험 결과

연구 질문

• RQ1 비선형 시스템을 매 시간 단계에서 풀지 않고도 광범위한 경사 하강 흐름에 대해 조건부 에너지 안정성을 보장할 수 있는 새로운 수치적 방법을 개발할 수 있는가?
• RQ2 에너지 안정성 스킴에서 제3차 및 제4차 정확도를 달성하면서도 선형 구조를 유지할 수 있는가?
• RQ3 IEQ나 볼록 분할법과 비교해 볼 때, SAV 방법은 비국소적 또는 다중 적분 자유 에너지 항을 더 효과적으로 다룰 수 있는가?
• RQ4 적응형 시간 스텝은 장시간 시뮬레이션에서 SAV 기반 스킴의 효율성과 정확도에 어떤 영향을 미치는가?
• RQ5 전통적인 스킴인 볼록 분할법이나 안정화 방법과 비교해 SAV 방법은 성능와 안정성 면에서 어떻게 다른가?

주요 결과

• SAV 방법은 볼록 분할법보다 훨씬 효율적이고 구현이 용이하며, 조건부로 제2차 에너지 안정성 스킴을 제공한다.
• 이 방법은 앨런-카운과 케인-힐리아드 방정식을 포함한 다양한 경사 하강 흐름에 대해 조건부 안정성, 선형성, 제2차 정확도 스킴을 구축할 수 있다.
• SAV 프레임워크 기반의 제3차 및 제4차 BDF 스킴은 조건부 안정성이 아니지만 실용적으로 매우 안정적이다.
• 수치 결과는 SAV 방법이 특히 적응형 시간 스텝을 사용할 경우 제2차 BDF 스킴보다 더 높은 정확도를 달성함을 보여준다.
• SAV 방법은 비국소 자유 에너지와 다중 적분 항을 포함한 모델에도 적용 가능하여, IEQ 방법의 적용 범위를 초월한다.
• 적응형 시간 스텝과 결합될 경우 SAV 스킴은 매우 효율적이고 정확해지며, 장시간 시뮬레이션에서 표준 스킴을 능가한다.