[논문 리뷰] A second-order accurate structure-preserving scheme for the Cahn-Hilliard equation with a dynamic boundary condition

Makoto Okumura, T. Fukao외 2명|arXiv (Cornell University)|2020. 07. 16.

Solidification and crystal growth phenomena참고 문헌 29인용 수 4

한 줄 요약

이 논문은 Cahn–Hilliard 방정식에 동적 경계 조건을 갖는 두 번째 차수 정확도를 가지며 구조를 유지하는 유한차분 스킴을 제안한다. 이는 이산 변분 도함수 방법(DVDM)을 사용하여 개발되었으며, 에너지 이산화를 재구성하고 합성-통분 공식을 적용함으로써, 공간적 두 번째 차수 정확도를 확보한다. 또한 질량 보존과 에너지 감소를 유지하며 안정성, 해의 존재성 및 유일성을 보장한다. 수치 예제를 통해 검증되었다.

ABSTRACT

We propose a structure-preserving finite difference scheme for the Cahn-Hilliard equation with a dynamic boundary condition using the discrete variational derivative method (DVDM). In this approach, it is important and essential how to discretize the energy which characterizes the equation. By modifying the conventional manner and using an appropriate summation-by-parts formula, we can use a standard central difference operator as an approximation of an outward normal derivative on the discrete boundary condition of the scheme. We show that our proposed scheme is second-order accurate in space, although the previous structure-preserving scheme by Fukao-Yoshikawa-Wada (Commun. Pure Appl. Anal. 16 (2017), 1915-1938) is first-order accurate in space. Also, we show the stability, the existence, and the uniqueness of the solution for the proposed scheme. Computation examples demonstrate the effectiveness of the proposed scheme. Especially through computation examples, we confirm that numerical solutions can be stably obtained by our proposed scheme.

연구 동기 및 목표

• 동적 경계 조건이 있는 Cahn–Hilliard 방정식을 위한 고차 정확도 수치 스킴을 개발하여 주요 물리적 구조를 유지하고자 한다.
• 이전의 구조 유지 스킴이 가지는 공간적 첫 번째 차수 정확도의 한계를 극복하고자 한다.
• 구조 유지 이산화를 통해 수치적 안정성, 해의 존재성 및 유일성을 확보하고자 한다.
• 장기 시뮬레이션 예제를 통해 스킴의 성능을 검증하고자 한다. 이는 안정적이고 물리적으로 일관된 행동을 보여준다.

제안 방법

• 스킴은 연속 문제의 에너지 감소 및 질량 보존 성질을 유지하기 위해 이산 변분 도함수 방법(DVDM)을 사용한다.
• 특정 합성-통분 공식과 결합된 수정된 에너지 이산화를 도입하여 공간적 두 번째 차수 정확도를 달성한다.
• 경계에서의 외부 법선 도함수는 이산 경계 조건을 철저히 다루어 표준 중심차분 연산자를 사용해 근사한다.
• 화학적 포텐셜 $ P $ 는 자유 에너지의 볼록 분할을 통해 반음성 형태로 갱신된다.
• 이산 경계 조건은 $ u $ 와 $ p $ 에 모두 동일한 노이만 조건을 강제하여 연속 문제와 일관성을 확보한다.
• 에너지 추정을 통해 스킴이 질량 보존과 에너지 감소를 유지하며, 안정성, 존재성 및 해의 유일성이 입증되었다.

실험 결과

연구 질문

• RQ1 동적 경계 조건이 있는 Cahn–Hilliard 방정식을 위한 구조 유지 유한차분 스킴이 공간적 두 번째 차수 정확도를 달성할 수 있는가?
• RQ2 고차 정확도를 달성하면서도 물리적 구조를 유지할 수 있도록 이산 경계 처리를 어떻게 설계할 수 있는가?
• RQ3 제안된 스킴이 장기 시뮬레이션에서 질량 보존과 에너지 감소를 유지하는가?
• RQ4 동적 경계 조건이 표준 노이만 조건과 비교하여 해의 장기적 행동에 어떤 영향을 미치는가?

주요 결과

• 제안된 스킴은 Fukao–Yoshikawa–Wada의 이전 구조 유지 스킴보다 공간적 첫 번째 차수 정확도를 초월하여 두 번째 차수 정확도를 달성한다.
• 수치 실험에서 질량 보존이 14자리 미만의 정밀도로 유지되어 스킴의 구조 유지 성격을 확인한다.
• 에너지 감소는 수치적으로 유지되며, 이산 에너지 함수가 시간에 따라 단조 감소한다.
• 장기 시뮬레이션 동안 해는 안정성을 유지하며 수치적 폭주나 진동이 관측되지 않는다.
• 수치 결과는 동적 경계 조건이 표준 노이만 경계 조건과 비교하여 다른 장기적 행동을 유도함을 확인한다.
• 이산 에너지 함수 $ A_{\rm d}^{(n)} $ 는 9자리 미만의 정밀도로 보존되어 스킴의 구조 유지 성격을 추가로 검증한다.