[논문 리뷰] Perturbative solution of a propagating interface in the phase field model

Mao Hiraizumi, Shin‐ichi Sasa|arXiv (Cornell University)|2021. 03. 25.

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

한 줄 요약

이 논문은 잠열 생성이 순서 매개변수와 온도 필드를 결합하는 초냉각 조건(∆ < 1)에서 전파되는 인터페이스에 대한 섭동 해를 제시한다. 작은 무차원 매개변수 η = ξṘ/D를 사용한 척도 가정을 통해, 1차 해를 유도하여 R(t) ∝ √t임을 보이고, 인터페이스 온도가 인터페이스 속도 Ṙ에 비례하여 평형 전이 온도 Tc에서 벗어남을 입증한다.

ABSTRACT

When a stable ordered phase and a metastable disordered phase are separated by a flat interface, the metastable state changes to the stable state through the propagation of the interface. For cases in which latent heat is generated, the interface displacement during some time interval is proportional to the square root of the time interval when the extent of supercooling is less than a certain value. We demonstrate this behavior by deriving a perturbative solution for a propagating interface in the phase field model. We calculate the leading-order contribution explicitly, and find that the interface temperature deviates from the equilibrium transition temperature in proportion to the interface velocity.

연구 동기 및 목표

• ∆ < 1 조건에서 단계 필드 모델에서 수치적으로 관측된 R(t) ∝ √t 의존성에 대한 이론적 유도.
• 초냉각 조건 하에서 전파되는 인터페이스에 대한 분석적 해의 부족을 해결.
• 잠열 존재 시 인터페이스 속도와 Tc로부터의 온도 편차 간의 관계 수립.
• 작은 매개변수 η = ξṘ/D를 사용한 척도 좌표와 함께 섭동 프레임워크 개발하여 1차 해에 대한 보정을 체계적으로 계산.

제안 방법

• 무차원 변수 도입: θ = cp(T - Tc)/L 및 ξ와 D를 사용하여 단계 필드 모델 방정식 (8)과 (9)를 재스케일링.
• 스케일드 공간 좌표 w = (x - R(t))/ξ 와 z = Ṙ(x - R(t))/D 정의로 인터페이스 및 열 길이 척도를 포괄.
• 정규 섭동 전개를 가능하게 하기 위해 작은 무차원 매개변수 η = ξṘ/D 도입.
• 해의 척도 형태 가정: φ(x,t) = Φ(w; η), θ(x,t) = Θ(z; η), R(t)는 1차적으로 R(t) = √(2DCt)로 진화.
• Φ0(w) 및 Θ0(z)에 대한 1차 방정식 (24)와 (28) 해석: Θ0(z) = -δ(z), Φ0(w)는 비선형 상미분방정식 만족.
• 다음 차수 보정 Φ1(w) 및 Θ1(z)에 대한 가용성 조건을 통합을 통해 유도하여 Θ1(0) 및 C1의 명시적 표현 도출.

실험 결과

연구 질문

• RQ1 단계 필드 모델에서 ∆ < 1 조건에서 수치적으로 관측된 R(t) ∝ √t 행동이 분석적으로 도출될 수 있는가?
• RQ2 잠열과 유한한 인터페이스 속도 존재 시 인터페이스 온도는 Tc에서 어떻게 벗어나는가?
• RQ3 무차원 초냉각 매개변수 ∆은 인터페이스 역학과 열 프로파일에 어떤 역할을 하는가?
• RQ4 η = ξṘ/D를 사용한 섭동 프레임워크는 1차 해에 대한 보정을 체계적으로 계산하는 데 어떻게 기여하는가?

주요 결과

• 1차 해는 R(t) = √(2DCt)를 도출하여 ∆ < 1 조건에서 인터페이스 위치의 √t 스케일링을 확인.
• 인터페이스 온도 T(R(t), t)는 인터페이스 속도 Ṙ에 비례하여 Tc에서 벗어나며, 이 편차는 T(R(t), t) = Tc + (LξṘ/cpD)Θ1(0)로 주어진다.
• 열 프로파일을 제어하는 매개변수 C는 ∆ → 1에서 아래로 발산하여 역학적 임계 전이를 나타낸다.
• 다음 차수 보정 Θ1(0)은 가용성 조건을 통해 명시적으로 계산되었으며, Θ1(0)은 자유 에너지 도함수를 포함하는 적분 I, J, K±에 의해 결정된다.
• D¨R/Ṙ³ 전개에서 계수 C1는 0으로 확인되어 섭동 전개의 일관성을 보장한다.
• Π = 0 (Ts = Tc) 경우, C(∆)는 계산되었으며 ∆ → 1에 수렴함에 따라 발산함을 보여, 임계 초냉각 임계점 근처에서 √t 스케일링의 붕괴를 확인한다.