[논문 리뷰] An energy-stable parametric finite element method for simulating solid-state dewetting

Quan Zhao, Wei Jiang외 1명|arXiv (Cornell University)|2020. 03. 03.

Fluid Dynamics and Thin Films참고 문헌 44인용 수 5

한 줄 요약

이 논문은 표면 확산과 접선 이동을 고려한 2차원 고체상 융해를 시뮬레이션하기 위한 에너지 안정성 있는 매개변수 유한요소법(ES-PFEM)을 제안한다. 이론적 접선 각을 로빈형 경계조건으로 재구성하고, 조각별 선형 유한요소와 후진 오일러 시간 이산화를 사용함으로써, 다양체 거리에서 이차 수준의 공간 정확도, 조건부 에너지 안정성, 장시간에 걸친 균일한 메시 분포를 달성하였으며, 수치적 검증을 통해 이차 수렴성과 일차 수준의 접선 각 정확도를 확인하였다.

ABSTRACT

We propose an energy-stable parametric finite element method (ES-PFEM) for simulating solid-state dewetting of thin films in two dimensions via a sharp-interface model, which is governed by surface diffusion and contact line (point) migration together with proper boundary conditions. By reformulating the relaxed contact angle condition into a Robin-type boundary condition and then treating it as a natural boundary condition, we obtain a new variational formulation for the problem, in which the interface curve and its contact points are evolved simultaneously. Then, the variational problem is discretized in space by using piecewise linear elements. A full discretization is presented by adopting the backward Euler method in time, and the well-posedness and energy dissipation of the full discretization are established. The numerical method is semi-implicit (i.e., a linear system to be solved at each time step and thus efficient), unconditionally energy-stable with respect to the time step, and second-order in space measured by a manifold distance between two curves. In addition, it demonstrates equal mesh distribution when the solution reaches its equilibrium, i.e., long-time dynamics. Numerical results are reported to show accuracy and efficiency as well as some good properties of the proposed numerical method.

연구 동기 및 목표

• 이동하는 접선을 고려한 2차원 고체상 융해를 시뮬레이션하기 위한 강력하고 에너지 안정적인 수치적 방법을 개발한다.
• 열린 곡선의 안정적이고 정확한 진화를 가능하게 하는 변분형식에 이론적 접선 각 조건을 통합하는 데 도전한다.
• 정확한 평형형태를 확보하기 위해 수치적 평형형태와 정확한 평형형태 사이의 다양체 거리에서 이차 수준의 공간 수렴성을 확보한다.
• 장시간 시뮬레이션 동안 질량과 에너지 안정성을 유지하고, 평형 상태에서 균일한 메시 분포를 확보한다.
• 각 시간단계에서 선형 시스템을 한 번만 푸는, 반음성적이고 효율적인 스킴을 제공한다.

제안 방법

• 이론적 접선 각 조건을 로빈형 경계조건으로 재구성함으로써, 변분형식에 자연스럽게 통합할 수 있도록 한다.
• 경계면 곡선과 그 접점이 표면 확산과 접선 이동에 따라 동시에 진화하는 새로운 변분형식을 유도한다.
• 매개변수 메시 위에서 조각별 선형 유한요소를 사용하여 공간 이산화를 수행한다.
• 후진 오일러 방법을 시간에 적용하여 전 이산화의 조건부 에너지 안정성과 잘 정의된 성질을 확보한다.
• 구축 과정에서 에너지 감소를 보장하며, 각 시간단계에서 이산 에너지가 단조 감소하도록 한다.
• 각 시간단계에서 선형 시스템을 한 번만 푸는 반음성적 스킴을 구현하여 계산 효율성을 확보한다.

실험 결과

연구 질문

• RQ1 이동하는 접선이 있는 고체상 융해에서 열린 곡선 진화에 대해 에너지 안정성을 유지하는 매개변수 유한요소법을 구성할 수 있는가?
• RQ2 이론적 접선 각 조건을 변분 프레임워크에 자연스럽게 강제할 수 있는 방법은 무엇인가, 이를 통해 안정성과 정확도를 확보할 수 있는가?
• RQ3 수치적 평형형태와 정확한 평형형태 사이의 다양체 거리에서 이 방법의 수렴 속도는 얼마인가?
• RQ4 장시간 평형 상태에서 질량을 보존하고 균일한 메시 분포를 나타내는가?
• RQ5 이러한 열린 곡선 진화 문제 유형에서 이 방법이 이차 수준의 공간 정확도를 달성할 수 있는가?

주요 결과

• 제안된 ES-PFEM은 수치적 평형형태와 정확한 평형형태 사이의 다양체 거리 오차에서 이차 수렴성을 확보하여 높은 공간 정확도를 입증한다.
• 평형 상태에서의 접선 각 오차는 메시 크기에 대해 일차 수렴성을 보이며, 즉 |θₑ − θᵢ| = O(h)임을 확인하여 열역학적 평형을 정확히 캐릭터라이즈할 수 있음을 검증한다.
• 이 방법은 조건부 에너지 안정성을 확보하며, 이산 에너지가 시간에 따라 단조 감소함으로써 장시간에 걸친 안정적인 시뮬레이션을 보장한다.
• 평형 상태에서 수치해는 균일한 메시 분포를 보이며, 이는 방법이 정확한 장시간 역학을 올바르게 캐릭터라이즈하고 있음을 시사한다.
• 수치적 시뮬레이션 결과, 초기 고립형태가 초기 기하구조와는 무관하게 항상 기준면과 윤거각 θᵢ에서 만족하는 원호 형태로 수렴함을 확인하였다.
• 이 방법은 질량을 보존하며, 반음성적이고 선형 시스템 기반의 시간 스텝을 통해 뛰어난 계산 효율성을 보였다.