[논문 리뷰] On Linear and unconditionally energy stable Algorithms for Variable Mobility Cahn-Hilliard Type Equation with Logarithmic Flory-Huggins Potential

Xiaofeng Yang, Jia Zhao|arXiv (Cornell University)|2017. 01. 25.

Solidification and crystal growth phenomena참고 문헌 20인용 수 10

한 줄 요약

이 논문은 로그형 플로리-후긴스 잠재력이 있는 변동 이동도 켈른-힐리아르 방정식에 대해, 불변 에너지 제곱화(IEQ) 접근법을 사용하여 선형이며 조건부 에너지 안정인 시간 적분 스킴을 제안한다. 비선형 잠재력을 변수 변환을 통해 이차형으로 재구성함으로써, 비선형 항을 반명시적으로 처리할 수 있게 되었으며, 이로 인해 선형, 대칭 양의정부형 시스템이 도출되어 시간에 관계없이 에너지 안정성이 보장되고 반복 해법을 통해 효율적으로 풀 수 있다. 이는 이 어려운 상전이 모델의 범주에 대해 선형, 고차, 증명된 조건부 에너지 안정성 스킴을 달성한 최초의 작업이다.

ABSTRACT

In this paper, we consider the numerical approximations for the fourth order Cahn-Hilliard equation with concentration dependent mobility, and the logarithmic Flory-Huggins potential. One challenge in solving such a diffusive system numerically is how to develop proper temporal discretization for nonlinear terms in order to preserve the energy stability at the time-discrete level. We resolve this issue by developing a set of the first and second order time marching schemes based on a novel, called "Invariant Energy Quadratization" approach. Its novelty is that the proposed scheme is linear and symmetric positive definite because all nonlinear terms are treated semi-explicitly. We further prove all proposed schemes are unconditionally energy stable rigorously. Various 2D and 3D numerical simulations are presented to demonstrate the stability, accuracy and efficiency of the proposed schemes thereafter.

연구 동기 및 목표

• 고무 및 상분리 시스템에서 널리 사용되는 로그형 플로리-후긴스 잠재력이 있는 변동 이동도 켈른-힐리아르 방정식에 대해 효율적이고 정확한 시간 적분 스킴을 개발하는 것.
• 기존 스킴의 한계를 극복하는 것 — 즉, 비선형(반복적 해법이 필요)이거나 조건부 안정성(작은 시간 간격이 필요)인 스킴은 두꺼운 인터페이스를 가진 강성 있는 시스템에 대해 특히 문제가 된다.
• 각 시간 단계에서 선형 시스템을 유지하면서도 시간에 관계없이 에너지 안정성을 확보함으로써 효율적이고 견고한 수치 시뮬레이션을 가능하게 하는 것.
• IEQ 접근법을 복잡한 경계가 있는 비선형 잠재력, 예를 들어 로그형 플로리-후긴스 잠재력과 같은 열역학적으로 일관된 모델에서 흔히 나타나는 것들로 확장하는 것.
• 복잡한 비선형 자유 에너지 항을 가진 그래디언트 플로우 문제에 대해 고차, 선형, 조건부 에너지 안정 스킴을 구성하는 일반적인 프레임워크를 제공하는 것.

제안 방법

• 비선형 잠재력을 보존하기 위해 새로운 보조 변수를 도입함으로써, 불변 에너지 제곱화(IEQ) 접근법을 사용해 로그형 플로리-후긴스 잠재력을 이차형으로 변환한다. 이는 비선형 항을 반명시적으로 처리할 수 있게 하며, 이를 통해 시스템의 에너지 소산 법칙이 유지되어 이산 수준에서도 에너지 안정성이 보존된다.
• 크랭크-니콜슨 시간 스텝을 사용한 1차 및 2차 시간 이산화 스킴(LS1-CN 및 LS2-CN)을 구성하며, 비선형 항은 IEQ 재구성에 의해 반명시적으로 처리된다.
• 결과로 도출된 선형 시스템은 대칭 양의정부형이므로, 질량 래핑 조절자와 함께 케리브-서브스페이스 방법을 사용해 효율적으로 해를 구할 수 있다.
• 표준 갈레르킨 유한요소 방법을 사용해 2D 및 3D 공간 이산화를 적용하였으며, 변분 형식이 일致성을 보장한다.
• 이론적 분석을 통해 제안된 모든 스킴이 시간 간격 크기와 관계없이 조건부 에너지 안정성을 보임을 증명하였다.

실험 결과

연구 질문

• RQ1 로그형 플로리-후긴스 잠재력이 있는 변동 이동도 켈른-힐리아르 방정식에 대해, 시간에 관계없이 에너지 안정성을 유지하는 선형, 고차 시간 이산화 스킴을 구성할 수 있는가?
• RQ2 IEQ 접근법이 기존의 볼록 분할 또는 안정화된 명시적 방법이 실패하거나 계산적으로 비효율적인 경우에도 안정적이고 효율적이며 정확한 수치 해법을 가능하게 하는가?
• RQ3 IEQ 기반 스킴이 비선형 해법을 피하면서도 조건부 에너지 안정성을 유지하고 대칭 양의정부형 선형 시스템을 보장할 수 있는가?
• RQ4 제안된 스킴은 2D 및 3D에서 상분리 역학의 장기 시뮬레이션에서 어떻게 성능을 발휘하는가? 특히 다양한 초깃값 조건에서의 성능은 어떻게 되는가?
• RQ5 IEQ 프레임워크는 복잡하고 경계가 있는 비선형 잠재력을 가진 다른 그래디언트 플로우 모델로 일반화 가능한가?

주요 결과

• 제안된 LS1-CN 및 LS2-CN 스킴은 변분 형식 기반 에너지 안정성 분석을 통해 엄밀히 조건부 에너지 안정성이 보장된다.
• 스킴은 선형이며 대칭 양의정부형이므로, 질량 래핑 조절자와 함께 케리브-서브스페이스 솔버를 사용해 효율적으로 해를 구할 수 있다.
• 2D 및 3D에서의 수치 시뮬레이션은 장기 안정성을 보여주며, 자유 에너지가 시간에 따라 단조 감소함을 확인하여 에너지 소산이 이루어짐을 입증한다.
• 3D에서의 상분리 역학은 φ₀ = 0.3 및 φ₀ = 0.7일 경우 구형 방울의 형성을 보이며, φ₀ = 0.5일 경우 이중층 구조를 형성한다. 이는 열역학적 기대와 일치한다.
• 2차 스킴인 LS2-CN은 δt = 0.001일 때 높은 정확도를 달성하였으며, 장기간의 통합에서도 불안정성이나 폭주 현상 없이 안정적인 결과를 보였다.
• 본 연구는 로그형 플로리-후긴스 잠재력이 있는 변동 이동도 켈른-힐리아르 시스템에 대해 선형, 고차, 조건부 에너지 안정 스킴을 제안한 최초의 작업으로, 문헌에서 중요한 격차를 메웠다.