[논문 리뷰] Explicit temperature coupling in phase-field crystal models of solidification

Maik Punke, Steven M. Wise외 2명|arXiv (Cornell University)|2022. 05. 30.

Solidification and crystal growth phenomena참고 문헌 54인용 수 10

한 줄 요약

이 논문은 온도에 따라 변화하는 격자 상수를 통해 온도 변화와 열팽창을 명시적으로 결합한 단상장 결정(PFC) 모델을 제시한다. 이는 응고 중 열응력과 형태 변화를 시뮬레이션할 수 있도록 한다. 전통적인 PFC 모델을 발전시켜 열역학적 유속과 온도에 의존하는 주기성을 통합하였으며, Eshelby 포함체 및 나선형 성장 시뮬레이션을 통해 검증되었으며, 정확한 열응력 예측과 온도 변화에 따른 격자 간격 조절을 입증하였다.

ABSTRACT

We present a phase-field crystal (PFC) model for solidification that accounts for thermal transport and a temperature-dependent lattice parameter. Elasticity effects are characterized through the continuous elastic field computed from the microscopic density field. We showcase the model capabilities via selected numerical investigations which focus on the prototypical growth of two-dimensional crystals from the melt, resulting in faceted shapes and dendrites. This work sets the grounds for a comprehensive mesoscale model of solidification including thermal expansion.

연구 동기 및 목표

• 메조스케일 PFC 모델을 개발하여 열역학적 전도를 결정 성장 역학과 명시적으로 결합한다.
• PFC 자유 에너지 기능에 온도에 따라 변화하는 격자 상수를 통합하여 열팽창 효과를 포함한다.
• 확산 시간 스케일에서 응고하는 결정에서 열응력과 탄성 회복을 시뮬레이션할 수 있도록 한다.
• 기존의 해석적 해(예: Eshelby 포함체)와의 비교를 통해 모델을 검증하고, 나선형 및 면체 성장에 대한 예측 능력을 입증한다.

제안 방법

• 파동수 q0가 β(T, T0)q0로 대체된 수정된 PFC 자유 에너지 기능을 사용하며, β는 εth(T, T0) = αth(T − T0)를 통해 열팽창을 표현한다.
• 결합된 진화 방정식은 밀도 장 ψ에 대한 보존적 H−1 기울기 유동과 온도 장 T에 대한 확산-반응 방정식을 포함하며, 열유속은 열확산도 M를 통해 모델링된다.
• 온도 장 T는 무차원이며 융해점에서 T = 1로 정규화되며, ψ와 T 사이에 γ∂tT − ϑ∂tψ 항을 통해 명시적인 결합이 이루어진다.
• 수치적 해는 스펙트럼 공간 이산화와 후진 오일러 시간 적분을 사용한 푸리에 가짜 스펙트럼 방법으로 계산된다.
• 연속적인 탄성 장은 밀도 장 ψ로부터 응력 텐서를 유도하여 기계적 응력 성분 σij를 계산할 수 있도록 한다.
• 시뮬레이션은 공간적으로 균일한 T = T0와 액체(A) 또는 고체(A, φ±) 상으로 설정된 초기 조건을 사용하며, 정오각형 및 삼각형 결정 구조를 재현하도록 매개변수를 조정한다.

실험 결과

연구 질문

• RQ1 온도에 따라 변화하는 격자 간격은 액상에서 고체로의 2차원 결정 성장 중 형태 및 응력 진화에 어떻게 영향을 미치는가?
• RQ2 PFC 모델은 고체 상에서 열팽창으로 인한 열응력을 어느 정도 정확하게 재현할 수 있는가?
• RQ3 열유속과 온도 변화를 포함함으로써, 비등온 PFC 모델과 비교했을 때 나선형 및 면체 성장 패턴은 어떻게 변화하는가?
• RQ4 열팽창과 온도 기울기가 존재하는 조건에서 모델은 Eshelby 포함체 문제를 정확하게 포착할 수 있는가?

주요 결과

• 모델은 Eshelby 포함체 문제를 성공적으로 재현하였으며, 응력장에 대한 해석적 및 수치적 해가 뛰어난 일치를 보였다.
• 열응력 σα_ij는 정확히 예측되었으며, T > T0일 경우 양수, T < T0일 경우 음수로 나타나, A1(T < T0)과 A2(T > T0)로 초기화된 결정에서 반대 부호를 가지며 물리적 일관성을 확인하였다.
• 격자 간격은 온도에 따라 변화하며, 크기는 σα_ij ≈ 10⁻⁵ 수준이었으며, 이는 열팽창이 고체-액체 인터페이스 간격과 결정 융합 역학을 조절함을 보여주었다.
• 열기울기 조건 하에서 나선형 및 면체 성장 형태가 자연스럽게 발생하였으며, 고체-액체 인터페이스에서 특징적인 응력 패tern이 관찰되었다.
• 액체 상에서는 응력장 σij가 0이 되며, 고체 상에서는 강한 공간적 변화를 보였다. σxx(A1)와 σxx(A2)는 각각 대칭성과 반대칭성을 가지는 σ0_ij 및 σα_ij 성분으로 분리되었으며, 이는 대칭성과 반대칭성을 반영한다.
• 모델은 탄성 회복과 확산 역학을 분리할 수 있도록 하였으며, 표준 PFC와 일치하는 느린 회복을 보였으며, 주요 방법론의 변경 없이 응력 기여도를 해석할 수 있도록 하였다.