[논문 리뷰] A fourth-order maximum principle preserving operator splitting scheme for three-dimensional fractional Allen-Cahn equations

Dongdong He, Kejia Pan외 1명|arXiv (Cornell University)|2018. 04. 19.

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

한 줄 요약

이 논문은 3차원 공간 분수형 Allen-Cahn 방정식을 위한 네째계수 최대원칙 유지 연산자 분할 스킴을 제안한다. Strang 분할, 분수형 확산 단계에서의 Crank-Nicolson ADI, 네째계수 유한차분, Richardson 보간을 조합함으로써, 시간 스텝 제약 조건 하에 무조건적 안정성과 이산 최대원칙 유지가 가능해지며, 2차원 및 3차원 테스트에서 네째계수 정확도와 최대원칙 만족도를 수치적으로 검증하였다.

ABSTRACT

In this paper, by using Strang's second-order splitting method, the numerical procedure for the three-dimensional (3D) space fractional Allen-Cahn equation can be divided into three steps. The first and third steps involve an ordinary differential equation, which can be solved analytically. The intermediate step involves a 3D linear fractional diffusion equation, which is solved by the Crank-Nicolson alternating directional implicit (ADI) method. The ADI technique can convert the multidimensional problem into a series of one-dimensional problems, which greatly reduces the computational cost. A fourth-order difference scheme is adopted for discretization of the space fractional derivatives. Finally, Richardson extrapolation is exploited to increase the temporal accuracy. The proposed method is shown to be unconditionally stable by Fourier analysis. Another contribution of this paper is to show that the numerical solutions satisfy the discrete maximum principle under reasonable time step constraint. For fabricated smooth solutions, numerical results show that the proposed method is unconditionally stable and fourth-order accurate in both time and space variables. In addition, the discrete maximum principle is also numerically verified.

연구 동기 및 목표

• 3차원 공간 분수형 Allen-Cahn 방정식을 위한 고차수 수치 스킴을 개발하여 이산 최대원칙을 유지하고자 한다.
• 무조건적 안정성을 유지하면서 시간 및 공간에서 네째계수 정확도를 달성하고자 한다.
• 합리적인 시간 스텝 제약 조건 하에 수치 해가 [−1, 1] 범위 내에서 유한하게 유지되도록 하고자 한다.
• 비선형 해법기법을 피하거나 에너지 안정성과의 상충을 피하는 선형적이고 효율적이며 무조건적 안정성 보장 스킴을 제공하고자 한다.
• 이산 최대원칙 유지 및 고차수 수렴 특성을 수치적으로 검증하고자 한다.

제안 방법

• Strang의 두 번째계수 연산자 분할을 적용하여 분수형 Allen-Cahn 방정식을 세 개의 부분문제로 분해한다: 두 개는 해석적 상미분방정식을 포함하고, 나머지 하나는 3차원 선형 분수형 확산 방정식을 포함한다.
• 3차원 분수형 확산 단계는 Crank-Nicolson 교차방향의료(ADI) 방법을 사용하여 다차원 문제를 일련의 1차원 문제로 줄인다.
• 공간 분수도를 고차수 유한차분 스킴으로 이산화하여 고차수 공간 정확도를 확보한다.
• 시간 적분의 정밀도를 네째계수로 높이기 위해 Richardson 보간을 적용한다.
• 푸리에 분석을 통해 매끄러운 해에 대해 전체 스킴의 무조건적 안정성을 증명한다.
• 이산 최대원칙을 수치적으로 유지하기 위해 시간 스텝 제약 조건을 유도한다.

실험 결과

연구 질문

• RQ1 3차원 공간 분수형 Allen-Cahn 방정식에 대해 이산 최대원칙을 유지하는 고차수(네째계수) 수치 스킴을 구성할 수 있는가?
• RQ2 분수도가 존재하는 상황에서 고차수 시간 및 공간 정확도를 확보하면서도 무조건적 안정성을 유지할 수 있는가?
• RQ3 제안된 선형 연산자 분할 스킴이 합리적인 시간 스텝 제약 조건 하에 해의 범위 [−1, 1]을 유지하는가?
• RQ4 분수도 α가 계면 두께와 상분리 시뮬레이션에서의 응집 동역학에 어떤 영향을 미치는가?
• RQ5 Richardson 보간이 안정성이나 최대원칙 유지에 영향을 주지 않으면서 분할 스킴의 시간 정밀도를 효과적으로 향상시킬 수 있는가?

주요 결과

• 푸리에 분석을 통해 3차원 분수형 Allen-Cahn 방정식의 매끄러운 해에 대해 제안된 스킴이 무조건적 안정성을 확보함을 증명하였다.
• 제작된 매끄러운 해를 사용한 수렴 테스트를 통해 시간 및 공간 변수에서 네째계수 정확도를 확인하였다.
• 이산 최대원칙은 수치적으로 검증되었으며, 시간 스텝이 유도된 제약 조건을 만족할 경우 해가 [−1, 1] 내에 유한하게 유지된다(예: α=1.7일 때 Δt ≤ 0.4945).
• 시간 스텝이 안정성 제약 조건을 초과할 경우(예: Δt = 2), 해의 최대값이 1을 초과하여 이산 최대원칙을 위반한다.
• 분수도 α를 감소시키면 계면 두께가 얇아지고 상분리 응집 속도가 느려지며, α = 1.2일 때 가장 비균일하고 세밀한 상 구조를 생성한다.
• 이 스킴은 높은 정밀도와 유한한 해를 유지하면서 복잡한 3차원 상분리 동역학, 즉 계면 진화 및 응집을 성공적으로 캐릭터라이즈하였다.