[논문 리뷰] Unconditionally stable exponential time differencing schemes for the mass-conserving Allen-Cahn equation with nonlocal and local effects

Kun Jiang, Lili Ju외 2명|arXiv (Cornell University)|2021. 07. 10.

Solidification and crystal growth phenomena참고 문헌 39인용 수 9

한 줄 요약

이 논문은 국소적 및 비국소적 효과를 고려한 질량 보존 알렌-초앙 방정식에 대해 조건부 안정인 지수 시간 차분(ETD) 스킴을 제안한다. 선형 안정화 기법과 ETD 방법을 조합하여, 시간에 관계없이 최대 bound 원칙(MBP)과 질량 보존을 유지하는 일阶 및 이阶 시간 적분기법을 구축하였으며, 정규성 가정 하에 엄밀한 오차 추정과 2차원 및 3차원에서의 광범위한 수치 검증을 수행하였다.

ABSTRACT

It is well known that the classic Allen-Cahn equation satisfies the maximum bound principle (MBP), that is, the absolute value of its solution is uniformly bounded for all time by certain constant under suitable initial and boundary conditions. In this paper, we consider numerical solutions of the modified Allen-Cahn equation with a Lagrange multiplier of nonlocal and local effects, which not only shares the same MBP as the original Allen-Cahn equation but also conserves the mass exactly. We reformulate the model equation with a linear stabilizing technique, then construct first- and second-order exponential time differencing schemes for its time integration. We prove the unconditional MBP preservation and mass conservation of the proposed schemes in the time discrete sense and derive their error estimates under some regularity assumptions. Various numerical experiments in two and three dimensions are also conducted to verify the theoretical results.

연구 동기 및 목표

• 질량 보존 알렌-초앙 방정식에 대해 최대 bound 원칙(MBP)과 질량 보존을 모두 만족하는 조건부 안정 시간 적분 스킴을 개발하는 것.
• 상태 분리 모델의 장기 시뮬레이션에서 물리적 및 수학적 구조—특히 유계성과 질량 보존—유지를 도전 과제로 삼는 것.
• 정확한 질량 보존을 위해 필요로 하는 비국소 라그랑주 승수 항을 처리할 수 있도록 기존 ETD 스킴을 확장하면서도 안정성과 정확성을 유지하는 것.
• 정규성 가정 하에 정확한 해와 수치 해에 대해 엄밀한 오차 추정을 제공하고, 두 차원 및 세 차원에서 수치적으로 스킴을 검증하는 것.

제안 방법

• 수치 안정성을 향상시키기 위해 질량 보존 알렌-초앙 방정식을 선형 안정화 기법을 사용해 재구성한다.
• 시간 적분을 위한 일阶 및 이阶 지수 시간 차분 룬게쿠타(ETDRK1 및 ETDRK2) 스킴을 구축한다.
• 스티프한 비선형 항을 암시적으로 처리하기 위해 지수 시간 차분 방법을 적용하여 무조건적 안정성을 확보한다.
• 적분 제약 조건을 통해 비국소 라그랑주 승수 항을 통합하여 각 시간 단계에서 정확한 질량 보존을 확보한다.
• 에너지 기반 및 노름 기반 분석을 사용하여, 에너지 기반 및 노름 기반 분석을 사용하여 시간 이산화된 의미에서 최대 bound 원칙(MBP)과 질량 보존이 무조건적으로 유지됨을 증명한다.
• 정확한 해와 수치 해에 대해 적절한 정규성 가정 하에 오차 추정을 유도한다.

실험 결과

연구 질문

• RQ1 질량 보존 알렌-초앙 방정식에 대해 지수 시간 차분 스킴을 설계하여 최대 bound 원칙(MBP)을 조건부로 무조건적으로 유지할 수 있는가?
• RQ2 그러한 스킴은 시간 이산화된 설정에서 정확한 질량 보존을 보장할 수 있는가?
• RQ3 표준 정규성 가정 하에 제안된 ETD 스킴의 수렴 행동과 오차 추정은 어떠한가?
• RQ4 복잡한 역학을 보이는 상태 분리의 장기 시뮬레이션에서 스킴의 성능은 어떠한가?
• RQ5 선형 안정화 기법은 이 클래스의 비국소 비선형 PDE에 대해 정확성과 안정성을 유지하는 데 효과적으로 통합될 수 있는가?

주요 결과

• 제안된 일阶 및 이阶 ETD 스킴은 최대 bound 원칙(MBP)을 조건부로 무조건적으로 유지하여, 모든 시간에 걸쳐 해가 초기 값의 범위 내에 유계로 유지됨을 보장한다.
• 스킴은 질량 보존 형식의 요구 조건에 따라 각 시간 단계에서 정확한 질량 보존을 달성한다.
• 2D 및 3D에서의 수치 실험은 에너지가 단조롭게 감소하고, 질량은 일정하며, 해가 지정된 구간 내에서 유계로 유지됨을 확인한다.
• ETDRK2 스킴은 시간에 대해 이계 정확도를 달성하며, 매끄러운 해 조건 하에서 관측된 수렴률이 이론적 예측과 일치한다.
• 확장하는 버블 테스트에서 버블의 반지름은 시간에 따라 단조롭게 증가하며, t ≈ 18에 약 0.407의 기대 정적 상태 반지름에 도달한다.
• 스킴은 장기 시뮬레이션 동안에도 안정적이고 정확하며, 안정성 확보를 위한 시간 간격 제약 조건이 필요 없음을 확인하여 무조건적 안정성을 입증한다.