[논문 리뷰] Positivity-preserving, energy stable numerical schemes for the Cahn-Hilliard equation with logarithmic potential

Wenbin Chen, Cheng Wang외 2명|arXiv (Cornell University)|2017. 12. 08.

Solidification and crystal growth phenomena참고 문헌 32인용 수 17

한 줄 요약

이 논문은 로그 Flory-Huggins 잠재력과 함께 Cahn-Hilliard 방정식에 대한 일阶 및 이阶 유한차분 스킴을 제시하며, 양수성 보존과 무조건적인 에너지 안정성을 보장한다. 스킴은 비선형 로그 항과 표면 확산을 암시적으로 처리하며, 로그 잠재력의 특이성을 활용하여 단계 변수가 ±1에 도달하는 것을 방지하여 잘 정의된 문제와 시간 스텝 제약 없이 최적 수렴 속도를 보장한다.

ABSTRACT

We present and analyze finite difference numerical schemes for the Allen Cahn/Cahn-Hilliard equation with a logarithmic Flory Huggins energy potential. Both the first order and second order accurate temporal algorithms are considered. In the first order scheme, we treat the nonlinear logarithmic terms and the surface diffusion term implicitly, and update the linear expansive term and the mobility explicitly. We provide a theoretical justification that, this numerical algorithm has a unique solution such that the positivity is always preserved for the logarithmic arguments. In particular, our analysis reveals a subtle fact: the singular nature of the logarithmic term around the values of $-1$ and 1 prevents the numerical solution reaching these singular values, so that the numerical scheme is always well-defined as long as the numerical solution stays similarly bounded at the previous time step. Furthermore, an unconditional energy stability of the numerical scheme is derived, without any restriction for the time step size. The unique solvability and the positivity-preserving property for the second order scheme are proved using similar ideas, in which the singular nature of the logarithmic term plays an essential role. For both the first and second order accurate schemes, we are able to derive an optimal rate convergence analysis, which gives the full order error estimate. The case with a non-constant mobility is analyzed as well. We also describe a practical and efficient multigrid solver for the proposed numerical schemes, and present some numerical results, which demonstrate the robustness of the numerical schemes.

연구 동기 및 목표

• 로그 Flory-Huggins 에너지 잠재력이 있는 Cahn-Hilliard 방정식을 위한 강력한 수치 스킴을 개발하여 단계 변수의 물리적 양수성을 유지한다.
• 시간 스텝 제약 없이 일阶 및 이阶 시간 스킴에 대한 무조건적 에너지 안정성을 확립한다.
• φ = ±1 근처에서 로그 잠재력의 특이성에 기반하여 유일한 해 존재성과 양수성 보존을 증명한다.
• 모든 스킴에 대해 ℓ∞(0,T;H⁻¹ₕ) ∩ ℓ²(0,T;H¹ₕ 노름에서 최적 수렴 속도를 유도한다.
• 효율적인 다중격자 솔버를 설계하고 수치 실험을 통해 스킴을 검증한다.

제안 방법

• 일阶 시간 이산화에서 로그 항과 표면 확산 항을 암시적으로 처리하고, 확산 및 이동성 항은 명시적으로 갱신한다.
• 이阶 BDF 시간 스킴을 사용하며, 확산 항에 대해 이阶 Adams-Bashforth 외삽을 적용하고, 안정성을 향상시키기 위해 인위적 Douglas-Dupont 정규화 항을 추가한다.
• 균일한 직각격자에서 주기적 경계 조건을 갖는 유한차분 공간 이산화를 적용한다.
• φ = ±1 근처에서 ln(1±φ)의 특이성을 활용하여 수치 해가 이 값에 도달하는 것을 방지함으로써 로그 항의 정의역을 보장한다.
• 볼록-오목 분해 기법과 이산 에너지 소모 분석을 통해 에너지 안정성을 도출한다.
• V-사이클 스무딩을 갖는 FAS 다중격자 솔버를 구현하여 각 시간 단계에서 비선형 시스템을 효율적으로 해결한다.

실험 결과

연구 질문

• RQ1 일阶 유한차분 스킴을 구성할 수 있는가? 이 스킴은 1+φ와 1−φ의 양수성을 유지하면서도 무조건적 에너지 안정성을 확보할 수 있는가?
• RQ2 로그 잠재력의 특이성은 수치 해가 φ = ±1에 도달하는 것을 자연스럽게 방지하여 문제의 잘 정의됨을 보장하는가?
• RQ3 인위적 정규화를 갖는 이阶 BDF 스킴은 동일한 양수성 및 에너지 안정성 성질을 달성할 수 있는가?
• RQ4 제안된 스킴의 최적 수렴 속도는 관련 이산 노름에서 얼마인가?
• RQ5 다중격자 솔버는 스킴에서 발생하는 비선형 시스템을 다루는 데 얼마나 효율적이고 강력한가?

주요 결과

• 일阶 스킴은 로그 잠재력이 φ = ±1 근처에서 특이성을 띠므로, 로그 항의 양수성을 보장하며 무조건적 에너지 안정성을 확보한다.
• 인위적 정규화를 갖는 이阶 BDF 스킴은 유일한 해 존재성과 양수성 보존을 유지하며, 이론적 근거는 동일한 특이성에 기반한다.
• 모든 스킴에 대해 ℓ∞(0,T;H⁻¹ₕ) ∩ ℓ²(0,T;H¹ₕ 노름에서 O(Δt + h²)의 최적 수렴 속도가 증명된다.
• 수치 결과는 강력성과 효율성을 확인하였으며, BDF2 스킴은 뛰어난 정확도를 보였고, 다중격자 솔버는 시간 단계당 평균 약 3.5회의 V-사이클 반복을 달성하였다.
• 모든 스킴과 시뮬레이션에서 최대 단계 값 φ는 0.9967 이하로 유지되어 강력한 양수성 보존이 확인되었다.
• 이론적 증명 없이도 완전히 암시적인 스킴(예: 역 오일러)도 실질적으로 양수성을 유지함을 확인하였으며, 이는 기반 메커니즘이 더 넓은 적용 가능성을 지닌다는 것을 시사한다.