[논문 리뷰] The generalized scalar auxiliary variable approach (G-SAV) for gradient flows

Qing Cheng|arXiv (Cornell University)|2020. 02. 01.

Solidification and crystal growth phenomena참고 문헌 35인용 수 7

한 줄 요약

이 논문은 그래디언트 플로우에 대한 조건부 에너지 안정 수치적 스킴을 구성하기 위한 새로운 프레임워크인 일반화된 스칼라 보조 변수(G-SAV) 방법을 제안한다. 기존 SAV 방법에서 보조 변수의 제한적인 제곱근 형태를 에너지 함수의 임의의 가역 함수로 대체함으로써, G-SAV는 각 시간 스텝에서 선형-일정 계수 해법을 효율적으로 수행하면서도 에너지 안정성을 유지할 수 있다. 이 방법은 에너지 유계성의 한계를 극복하고, 결합 Cahn-Hilliard 모델과 같은 강성 있는 시스템에 대해 더 높은 안정성과 탄력성을 제공한다.

ABSTRACT

We establish a general framework for developing, efficient energy stable numerical schemes for gradient flows and develop three classes of generalized scalar auxiliary variable approaches (G-SAV). Numerical schemes based on the G-SAV approaches are as efficient as the original SAV schemes \cite{SXY19,cheng2018multiple} for gradient flows, i.e., only require solving linear equations with constant coefficients at each time step, can be unconditionally energy stable. But G-SAV approaches remove the definition restriction that auxiliary variables can only be square root function. The definition form of auxiliary variable is applicable to any reversible function for G-SAV approaches . Ample numerical results for phase field models are presented to validate the effectiveness and accuracy of the proposed G-SAV numerical schemes.

연구 동기 및 목표

• 기존 SAV 및 IEQ 방법에서 보이는 보조 변수의 제한적인 제곱근 형태를 제거함으로써, 그래디언트 플로우의 에너지 안정적 시간 이산화를 위한 일반적 프레임워크를 개발하는 것.
• 자유 에너지가 아래로 유계가 아닌 시스템에 대해서도 에너지 안정적 스킴을 가능하게 하는 것, 이는 기존 SAV 접근 방식의 한계를 극복하기 위함이다.
• 각 시간 스텝에서 일정 계수를 가진 선형 해법만 요구함으로써, 원래 SAV 스킴과 유사한 계산 효율성을 유지하는 것.
• 원래 에너지 소산 법칙을 유지하는 것과 동시에 안정성, 구현의 단순성, 그리고 에너지 소산 법칙의 유지 간의 균형을 고려한 G-SAV의 여러 변형을 제공하는 것.
• Allen-Cahn 및 Cahn-Hilliard 상전이 모델, 특히 결합 BCP 모델에 대한 광범위한 수치 시뮬레이션을 통해 방법을 검증하는 것.

제안 방법

• G-SAV 프레임워크는 기존의 고전적 SAV 방법에서 사용하는 제곱근 형태가 아닌, 에너지 함수의 임의의 가역 함수로 정의된 보조 변수를 도입함으로써 제안된다.
• 이 방법은 일반화된 보조 변수를 도입함으로써 원래 그래디언트 플로우를 등가 시스템으로 변환하며, 시험 함수와 시험 함수 공간이 동일한 함수 공간에 속하는 변분 형식을 통해 에너지 안정성을 확보한다.
• 세 가지의 구체적 G-SAV 접근 방식이 제안된다: (i) 보조 변수를 암시적으로 처리하여 낮은 비용의 비선형 대수적 해법을 요구하는 방법; (ii) 보조 변수를 명시적으로 처리하여 완전히 명시적이고 반복이 없는 스킴을 만드는 방법; (iii) 원래 에너지 소산 법칙을 유지하는 변형(수정된 법칙이 아닌).
• 이 스킴들은 변분 형식과 시험 함수를 동일한 함수 공간에 놓은 이산 에너지 소산 분석을 통해 조건부 에너지 안정성이 증명된다.
• 시간 이산화에는 2차 BDF 스킴(BDF2)을 사용하며, 공간 이산화에는 유한요소법 또는 스펙트럼 방법과 같은 임의의 갈레르킨 유형 방법과 호환된다.
• 정확도와 효율성을 향상시키기 위해 장기 시뮬레이션에 적합한 안정화된 형식과 적응형 시간 스텝을 사용하여 방법을 구현한다.

실험 결과

연구 질문

• RQ1 보조 변수를 제곱근 함수로 제한하지 않고도 그래디언트 플로우에 대한 에너지 안정적 수치 스킴을 구성할 수 있는가?
• RQ2 보조 변수에 사용되는 가역 함수의 선택이 수치 스킴의 안정성, 정확도 및 탄력성에 어떤 영향을 미치는가?
• RQ3 일반화된 보조 변수 프레임워크는 IEQ 및 SAV 방법이 사용하는 수정된 에너지가 아닌 원래 에너지 소산 법칙을 유지할 수 있는가?
• RQ4 강성 있는 시스템인 결합 Cahn-Hilliard(BCP) 모델에 대해 G-SAV 접근 방식은 기존의 SAV 및 ETDRK2 스킴보다 어떻게 성능을 발휘하는가?
• RQ5 에너지 안정성과 계산 효율성을 유지하면서 G-SAV 프레임워크를 다성분 시스템으로 확장할 수 있는가?

주요 결과

• G-SAV 접근 방식은 보조 변수를 정의하는 데 어떤 가역 함수를 사용하든 간에 그래디언트 플로우에 대해 조건부 에너지 안정성을 달성하며, 에너지가 아래로 유계여야 한다는 요구 조건이 제거된다.
• G-SAV 기반 수치 스킴은 각 시간 스텝에서 일정 계수를 가진 선형 방정식만 풀면 되므로, 원래 SAV 스킴과 동일한 계산 효율성을 유지한다.
• 첫 번째 G-SAV 접근 방식(보조 변수를 암시적으로 처리)은 결합 Cahn-Hilliard 모델과 같은 강성 있는 시스템에 대해 두 번째 접근 방식(보조 변수를 명시적으로 처리)보다 더 높은 안정성을 보인다.
• 세 번째 G-SAV 접근 방식은 IEQ 및 SAV 방법이 사용하는 수정된 에너지가 아닌 원래 에너지 소산 법칙을 유지하며, [12]의 라그랑주 승수 방법과 동일한 결과를 낳는다.
• BCP 모델의 시뮬레이션 결과, $ G = e^{x/c} $ 및 $ G = \tilde{\sqrt{x+c}} $를 사용한 G-SAV 스킴은 구분할 수 없을 정도로 유사한 형태의 진화 패턴을 생성하며, 보조 함수의 선택에 대해 방법의 강건성을 확인한다.
• G-SAV 스킴은 ETDRK2 스킴과 원래 SAV 방법 수준의 정확도를 달성하며, 장기 시뮬레이션에서 더 높은 안정성과 효율성을 보인다.