[논문 리뷰] An Energy Stable Finite-Difference Scheme for Functionalized Cahn-Hilliard Equation and its Convergence Analysis

Wenqiang Feng, Zhen Guan외 3명|arXiv (Cornell University)|2016. 10. 08.

Solidification and crystal growth phenomena참고 문헌 6인용 수 3

한 줄 요약

이 논문은 페르소나형 칸힐리아르(FCH) 방정식, 즉 복합계에서의 상분리 현상을 모델링하는 6차 비선형 포아르비아 방정식에 대해 조건부 에너지 안정성과 수렴성을 보장하는 유한차분 스킴을 제안한다. 비볼록 에너지 함수를 보정하기 위해 보조 항을 도입함으로써, 비볼록성 문제를 해결하고 볼록 분할 접근법을 적용하여 유일한 해 존재성과 조건부 에너지 안정성을 확보하였다. 이는 전역적 $H_{\rm per}^{2}$ 안정성과 효율적인 조절된 최강하강법을 사용한 $H^{-1}$-기울기 흐름 시스템의 해법을 통해 수렴성이 증명되었다.

ABSTRACT

We present and analyze an unconditionally energy stable and convergent finite difference scheme for the Functionalized Cahn-Hilliard equation. One key difficulty associated with the energy stability is based on the fact that one nonlinear energy functional term in the expansion appears as non-convex, non-concave. To overcome this subtle difficulty, we add two auxiliary terms to make the combined term convex, which in turns yields a convex-concave decomposition of the physical energy. As a result, an application of the convex splitting methodology assures both the unique solvability and the unconditional energy stability of the proposed numerical scheme. To deal with a 4-Laplacian solver in an $H^{-1}$ gradient flow at each time step, we apply an efficient preconditioned steepest descent algorithm to solve the corresponding nonlinear systems. In addition, a global in time $H_{ m per}^2$ stability of the numerical scheme is established at a theoretical level, which in turn ensures the full order convergence analysis of the scheme. A few numerical results are presented, which confirm the stability and accuracy of the proposed numerical scheme.

연구 동기 및 목표

• 복합계 상분리에서 발생하는 매우 비선형적이며 6차의 기능화된 칸힐리아르(FCH) 방정식을 위한 수치적 안정성과 수렴성을 보장하는 스킴을 개발하기 위해.
• 에너지 함수의 비볼록성으로 인해 에너지 안정성 분석과 수치적 구현이 복잡해지는 문제를 해결하기 위해.
• 시간 스텝 제약 없이도 조건부 에너지 안정성과 유일한 해 존재성을 보장하기 위해.
• 각 시간 단계에서 발생하는 $H^{-1}$-기울기 흐름 시스템을 효율적으로 해결하기 위한 비선형 해법을 설계하기 위해.
• 수치적 스킴의 전역적 $H_{\rm per}^{2}$ 안정성에 기반한 엄밀한 수렴성 분석을 수립하기 위해.

제안 방법

• 에너지 함수에 두 개의 보조 항을 도입하여 볼록-비볼록 분해를 만들어내어 볼록 분할 기법의 적용을 가능하게 하였다.
• 시간에 대해 완전히 암시적이고, 공간에 대해 중심 차분을 사용한 유한차분 스킴을 설정하여 조건부 에너지 안정성을 확보하였다.
• 각 시간 단계에서 발생하는 $H^{-1}$-기울기 흐름으로부터 유도된 비선형 시스템을 해결하기 위해 조절된 최강하강(PSD) 알고리즘을 적용하였다.
• 에너지 함수로부터 화학적 포텐셜 $\mu$를 유도하기 위해 변분 형식을 사용하였으며, 이는 4라플라시안 항과 고차 비선형성을 포함한다.
• 수치적 해의 전역적 $H_{\rm per}^{2}$ 안정성을 확립하여, 전체 차수의 수렴성 분석의 기초를 마련하였다.
• 주기적 경계 조건을 구현하고, 질량 보존과 에너지 감소를 유지하기 위해 반암시적 시간 이산화를 사용하였다.

실험 결과

연구 질문

• RQ1 에너지 함수의 비볼록성에도 불구하고 FCH 방정식에 대해 조건부 에너지 안정성을 확보하는 수치적 스킴을 구성할 수 있는가?
• RQ2 비볼록 에너지 항은 어떻게 재구성하여 볼록 분할이 가능하게 하고, 해의 유일성 보장을 확보할 수 있는가?
• RQ3 제안된 스킴의 수렴 차수는 무엇이며, 안정성 추정치를 기반으로 엄밀하게 증명할 수 있는가?
• RQ4 실제로 $H^{-1}$-기울기 흐름으로부터 유도된 고차 비선형 시스템은 얼마나 효율적으로 해결될 수 있는가?
• RQ5 장기간의 시뮬레이션 동안 수치적 스킴이 에너지 감소와 질량 보존과 같은 핵심 물리적 성질을 유지하는가?

주요 결과

• 보조 항을 통한 에너지 함수의 볼록-비볼록 분해를 통해 제안된 스킴은 조건부 에너지 안정성과 유일한 해 존재성을 확보하였다.
• 스킴은 시간에 대해 1차 수렴성을 보였으며, 공간에 대해서는 약 2.0–2.5의 수렴률을 관측하여 이론적 예측을 확인하였다.
• 조절된 최강하강 해법은 강력한 수렴성을 보였으며, 반복마다 잔차가 거의 일정한 요소만큼 감소하였다.
• 수치적 시뮬레이션은 장기간의 진화 동안 에너지 감소와 질량 보존을 확인하였으며, 물리적 기대와 일치하였다.
• 스킴은 밀레트 네트워크, 이중층-실린더 접합부, Y형 접합부와 같은 복잡한 형태의 구조를 성공적으로 포착하였으며, 기존 실험 및 시뮬레이션 결과와 일치하였다.
• 수치적 스킴의 전역적 $H_{\rm per}^{2}$ 안정성이 이론적으로 확립되었으며, 이는 전체 차수의 수렴성 분석을 가능하게 하였다.