[논문 리뷰] Low artificial anisotropy cellular automaton model and its applications to the cell-to-dendrite transition in directional solidification

Lei Wei, Xin Lin외 2명|arXiv (Cornell University)|2014. 10. 15.

Solidification and crystal growth phenomena참고 문헌 1인용 수 3

한 줄 요약

이 논문은 방향 응고 중 미세구조 진화를 시뮬레이션하기 위해 제한된 이웃 고체 분율(LNSF) 캡처 규칙을 적용한 낮은 인위적 이방성 셀룰러 오토마타(CA) 모델을 제시한다. 시뮬레이션과 SCN-ace 및 SCN-카무포르 시스템의 실험 데이터를 통해 검증된 결과, 용질 분할 계수 k₀가 임계 임계값 0.125를 초과할 경우 세포에서 나선형으로의 전이 속도(Vcd)가 k₀에 비례하는 scaling(Vcd ∝ Vc/k₀)에서 일정한 값(Vcd = 8Vc)으로 전이됨을 밝혀냈다.

ABSTRACT

A low artificial anisotropy cellular automaton (CA) model is developed for the simulation of microstructure evolution in directional solidification. The CA model's capture rule was modified by a limited neighbor solid fraction (LNSF) method. Various interface curvature calculation methods have been compared. The simulated equilibrium shapes agree with the theoretical shapes, when the interface energy anisotropy coefficient is ε=0.01, ε=0.03 and ε=0.05, respectively.The low artificial anisotropy CA model is used in the numerical simulation of the cell-to-dendrite transition (CDT) in directional solidification. The influence of physical parameters (Γ, Dl, k0, ml) on CDT has been investigated. The main finding in this paper is the discovery of the changing behavior of the Vcd when the solute partition coefficient k0 is larger than a critical value. When k0 is less than 0.125, the Vcd follows the Kurz and Fisher criterion Vc/k0; while when k0>0.125, the Vcd equals to 8Vc. The experimental data of succinonitrile-acetone (SCN-ace, k0=0.1) and SCN-camphor (k0=0.33) support the conclusion from CA simulations.

연구 동기 및 목표

• 방향 응고에서 정확한 미세구조 시뮬레이션을 위해 최소한의 인위적 이방성을 갖춘 셀룰러 오토마타 모델을 개발하기 위해.
• 세포에서 나선형으로의 전이(CDT)에 영향을 미치는 핵심 물리적 매개변수(Γ, Dl, k₀, ml)의 영향을 조사하기 위해.
• 용질 분할 계수 k₀에 따른 세포에서 나선형으로의 전이 속도(Vcd)의 전이 행동을 명확히 하기 위해.
• 다양한 k₀ 값을 갖는 시클로헥산 기반 시스템의 실험 데이터와 모델의 시뮬레이션 결과를 비교 검증하기 위해.

제안 방법

• 모델의 인위적 이방성을 감소시키기 위해, 제한된 이웃 고체 분율(LNSF) 방법을 도입하여 CA 캡처 규칙을 수정하였다.
• 정확한 나선형 성장 시뮬레이션을 확보하기 위해 여러 경계 곡률 계산 방법을 평가하고 비교하였다.
• 낮은 이방성 계수(ε = 0.01, 0.03, 0.05)에서 이론적 평형 형상과의 일치를 확인하기 위해 이론적 평형 형상에 기반한 모델 校정을 수행하였으며, 이론과 양호한 일치를 보였다.
• k₀, Dl, ml, 및 Γ와 같은 다양한 물리적 매개변수를 포함한 세포에서 나선형으로의 전이에 대한 수치 시뮬레이션을 수행하였다.
• k₀ = 0.1 (SCN-ace) 및 k₀ = 0.33 (SCN-카무포르)인 시스템에 대해 모델을 적용하여 실험 관측 결과와 일치하는 결과를 도출하였다.

실험 결과

연구 질문

• RQ1 용질 분할 계수 k₀는 방향 응고에서 세포에서 나선형으로의 전이 속도(Vcd)에 어떤 영향을 미치는가?
• RQ2 Vcd 행동이 Vc/k₀ scaling에서 일정한 값으로 전이되는 k₀의 임계 값은 얼마인가?
• RQ3 낮은 인위적 이방성 CA 모델은 낮은 이방성에서 이론적 평형 형상과 얼마나 정확하게 재현되는가?
• RQ4 SCN-ace 및 SCN-카무포르 시스템에 대한 시뮬레이션된 Vcd 값은 실험 데이터와 얼마나 잘 일치하는가?
• RQ5 경계 에너지 이방성(ε)은 성장하는 나선형의 안정성과 형태에 어떤 역할을 하는가?

주요 결과

• 낮은 이방성 계수(ε = 0.01, 0.03, 0.05)에서 이론적 평형 형상을 성공적으로 재현하여 인위적 이방성이 최소화되었음을 확인하였다.
• k₀ < 0.125일 경우, 세포에서 나선형으로의 전이 속도(Vcd)는 Kurz와 Fisher 기준에 따라 Vc/k₀ 비례 scaling을 따른다.
• k₀ > 0.125일 경우, Vcd는 일정한 값 8Vc로 안정화되며, 성장 역학에서 명백한 전이가 일어남을 시사한다.
• SCN-ace(k₀ = 0.1) 및 SCN-카무포르(k₀ = 0.33)에 대한 시뮬레이션된 Vcd 값은 실험 관측 결과와 일치하여 모델의 타당성을 뒷받침한다.
• k₀ = 0.125에서의 전이는 용질 확산 제한에서 열역학적 제어로의 성장 메커니즘 전이를 나타내며, 중요한 전환점이다.