[논문 리뷰] Numerical Energy Dissipation for Time-Fractional Phase-Field Equations

Chaoyu Quan, Tao Tang외 1명|arXiv (Cornell University)|2020. 09. 14.

Solidification and crystal growth phenomena참고 문헌 31인용 수 11

한 줄 요약

이 논문은 비균일 시간 스텝을 사용할 때 시간 분수계 단상장 방정식에 대해 제1차 안정화된 반음향 L1 스킴을 사용하여 이산 에너지 법칙을 수립하고, 에너지 감소성과 유계성을 증명한다. 이는 (2−α)-차수의 L1-SAV 스킴으로 확장되어, 작은 α 값에서도 수정된 에너지가 여전히 유계이고 감소성을 유지함을 보여주며, 엄밀한 수치적 검증을 수반한다.

ABSTRACT

The numerical integration of phase-field equations is a delicate task which needs to recover at the discrete level intrinsic properties of the solution such as energy dissipation and maximum principle. Although the theory of energy dissipation for classical phase field models is well established, the corresponding theory for time-fractional phase-field models is still incomplete. In this article, we study certain nonlocal-in-time energies using the first-order stabilized semi-implicit L1 scheme. In particular, we will establish a discrete fractional energy law and a discrete weighted energy law. The extension for a $(2-α)$-order L1 scalar auxiliary variable scheme will be investigated. Moreover, we demonstrate that the energy bound is preserved for the L1 schemes with nonuniform time steps. Several numerical experiments are carried to verify our theoretical analysis.

연구 동기 및 목표

• 시간 분수계 단상장 방정식에 대해 이산 분수계 및 가중 에너지 법칙을 수립하기 위해.
• 비균일 시간 스텝을 사용하는 제1차 안정화된 반음향 L1 스킴의 에너지 안정성 분석을 위해.
• 분석을 더 높은 차수인 (2−α)-차수의 L1-SAV 스킴으로 확장하기 위해.
• 에너지의 유계성과 감소성이 수치적으로 유지되는지 확인하기 위해, 특히 작은 α 값에서의 경우를 중심으로.
• 시간 분수계 PDE에 적용 가능한 에너지 안정성 수치 스킴의 프레임워크를 제공하기 위해.

제안 방법

• 상태-분수계 시간 도함수를 추정하기 위해 제1차 안정화된 반음향 L1 스킴을 사용하여 단상장 방정식에 적용한다.
• 에너지 안정성 분석과 이산 에너지 법칙 유도를 위해 특수한 콜레스키 분해 기법을 적용한다.
• L1 스킴에 대해 이산 분수계 에너지 법칙과 이산 가중 에너지 법칙을 도출한다.
• 분석을 (2−α)-차수의 L1-SAV 스킴으로 확장하여 에너지의 유계성과 감소성을 보장한다.
• 시간에 따라 변화하는 메쉬(등급 기반)를 사용하여 t=0 근처의 해상도를 향상시킨다. 특히 작은 α 값에서 유리하다.
• Allen–Cahn, Cahn–Hilliard 및 MBE 방정식에 대한 수치 실험을 통해 이론적 결과를 검증한다.

실험 결과

연구 질문

• RQ1 시간 분수계 단상장 방정식에 적용된 제1차 L1 스킴에 대해 이산 분수계 에너지 법칙을 수립할 수 있는가?
• RQ2 비균일 시간 스텝을 사용하는 L1 스킴이 에너지의 유계성과 감소성을 유지하는가?
• RQ3 (2−α)-차수의 L1-SAV 스킴이 시간 분수계 단상장 모델의 에너지 안정성과 유계성을 유지할 수 있는가?
• RQ4 α 값을 감소시킬수록 L1 스킴에서 에너지 감소성과 수치적 해상도에 어떤 영향을 미치는가?
• RQ5 콜레스키 분해 기법이 시간 분수계 문제의 에너지 안정성 분석에 얼마나 효과적인가?

주요 결과

• 제1차 L1 스킴에 대해 이산 분수계 에너지 법칙이 수립되어 반연속 수준에서 에너지 감소성이 보장된다.
• 비균일 시간 스텝을 사용하는 L1 스킴이 작은 α 값에서도 에너지의 유계성을 유지함을 수치적으로 확인하였다.
• α=0.3일 때, 표준 L1 스킴에서는 해상도가 열 劣하여 에너지 감소성이 위반되지만, 여전히 에너지의 유계성은 유지된다.
• (2−α)-차수의 L1-SAV 스킴은 에너지의 유계성과 이산 분수계 에너지 법칙을 모두 보장한다.
• 수치 실험 결과, 모든 테스트된 α 값과 메쉬 유형에서 수정된 에너지가 초기 에너지 이하로 유계로 유지됨을 확인하였다.
• 콜레스키 분해 기법이 시간 분수계 단상장 문제의 에너지 안정성 분석에 효과적임을 입증하였다.