[논문 리뷰] A Variational Lagrangian Scheme for a Phase Field Model: A Discrete Energetic Variational Approach

Chun Liu, Yiwei Wang|arXiv (Cornell University)|2020. 03. 23.

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

한 줄 요약

이 논문은 에너지-소산 구조를 유지하기 위해 이산 에너지적 변분 방법을 활용한 수정된 Allen-Cahn 상전이 모델을 위한 변분 라그랑주 방법을 제안한다. 라그랑주 사상과 정규화 항을 통해 동역학을 기술함으로써, 이 방법은 에너지 안정성을 확보하고 최소한의 메쉬 포인트로 얇은 분산 경계면을 정확히 포착하여 평형 상태의 안정적인 계산을 가능하게 한다.

ABSTRACT

In this paper, we propose a variational Lagrangian scheme for a modified phase-field model, which can compute the equilibrium states for the original Allen-Cahn type model. Our discretization is based on a prescribed energy-dissipation law in terms of the flow map. By employing a discrete energetic variational approach, this scheme preserves the variational structure of the original energy-dissipation law and is energy stable. Plentiful numerical tests show that, by choosing the initial value properly, our methods can produce the desired equilibrium and capture the thin diffuse interface with a small number of mesh points.

연구 동기 및 목표

• Allen-Cahn 유형의 상전이 모델에서 평형 상태를 계산하기 위한 구조를 유지하는 수치적 스킴을 개발하기 위해.
• 에울러리안 방법이 최소한의 계산 비용으로 얇은 분산 경계면을 해상도하는 데에 한계를 가진다는 점을 극복하기 위해.
• 라그랑주 사상에 자연스러운 변분 구조가 없는 $L^2$-강하류에 대해 라그랑주 방법을 확장하기 위해.
• 에너지-소산 법칙에 정규화 항을 도입하여 에너지 안정성과 메쉬 품질을 확보하기 위해.
• 최소한의 메쉬 해상도로 날카운 경계면을 포착하고 정확한 평형 상태로 수렴하는 스킴의 효과성을 입증하기 위해.

제안 방법

• 스킴은 흐름 사상에 기반한 정의된 에너지-소산 법칙을 바탕으로 하여 연속적인 변분 구조의 유지 보장을 보장한다.
• 이산 에너지적 변분 방법을 적용하여 시간 이산화 스킴을 유도하며, 이는 이산 수준에서 연속적인 에너지-소산 법칙을 유지한다.
• 메쉬 뒤틀림을 방지하고 안정성을 확보하기 위해 소산 항에 정규화 항 $\nu|\nabla\mathbf{u}|^2$ 를 통합한다.
• 시스템은 라그랑주 사상 $\mathbf{X}(\mathbf{X}_0, t)$ 를 사용하여 기술되며, 상 함수 $\varphi$ 는 $\varphi_t + \nabla\varphi \cdot \mathbf{u} = 0$ 를 통해 진화한다.
• 이산 시스템은 일관된 질량 행렬 $\mathbf{M}^e$ 와 강성 행렬 $\mathbf{역행렬}^e$ 를 사용한 유한요소 방법으로 해결되며, $\nu > 0$ 일 때 양의 정부호성을 보장한다.
• 이 방법은 자가 적응형이며, 메쉬 포인트를 적응 메쉬 세분화가 필요 없이 분산 경계면에 집중시킨다.

실험 결과

연구 질문

• RQ1 변분 라그랑주 스킴은 이산 수준에서 Allen-Cahn 유형의 상전이 모델의 에너지-소산 구조를 유지할 수 있는가?
• RQ2 에너지-소산 법칙에 자연스러운 변분 구조가 없는 $L^2$-강하류에 대해 라그랑주 방법은 어떻게 적응시킬 수 있는가?
• RQ3 제안된 스킴은 최소한의 메쉬 해상도로 평형 상태를 정확히 계산하고 메쉬 뒤틀림을 피할 수 있는가?
• RQ4 정규화 항 $\nu|\nabla\mathbf{u}|^2$ 는 수치적 안정성과 메쉬 품질 확보에 어떤 역할을 하는가?
• RQ5 초기 조건 $\varphi_0$ 의 선택이 라그랑주 시뮬레이션에서 정확한 평형 상태로의 수렴에 어떻게 영향을 미치는가?

주요 결과

• 제안된 스킴은 이산 에너지적 변분 형식을 통해 연속적인 에너지-소산 법칙의 변분 구조를 유지하며 에너지 안정성을 확보한다.
• 적절하게 선택된 초기 조건을 사용할 경우, 이 방법은 Allen-Cahn 모델의 정확한 평형 상태를 성공적으로 계산한다.
• 스킴은 적은 수의 메쉬 포인트로 얇은 분산 경계면을 포착하여 높은 해상도 효율성을 보여준다.
• 수치 실험 결과, 이 방법은 단순한 라그랑주 접근 방식과 달리 메쉬 뒤틀림을 피하고 늦은 시점까지 양호한 메쉬 품질을 유지함을 보였다.
• 정규화 항 $\nu|\nabla\mathbf{u}|^2$ 는 이산 소산 행렬 $\mathbf{D}$ 의 양의 정부호성을 보장하여 수치적 붕괴를 방지하는 데 필수적이다.
• 부적절한 초기 조건 $\varphi_0$ 를 사용한 실패한 사례는 $\varphi_0$ 가 $[-1,1]$ 범위에 속하지 않을 경우 라그랑주 방법이 정확한 평형 상태에 도달하지 못할 수 있음을 보여주며, 초기 자료 선택의 중요성을 강조한다.