[논문 리뷰] Linking Phase-Field and Atomistic Simulations to Model Dendritic Solidification in Highly Undercooled Melts

J. Bragard, Alain Karma외 2명|arXiv (Cornell University)|2001. 12. 10.

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

한 줄 요약

이 연구는 분자역학에서 유도된 비등방성 계면 특성을 사용하여 상전이 필드와 원자적 시뮬레이션을 융합하여 고도로 초냉각된 니켈 액체에서의 수상성 응고를 모델링한다. 이는 운동학적 비등방성이 실험적 수상성 끝단 속도를 정량적으로 일치시키는 데 핵심적임을 보여주며, 표준 확산 모델이 관측된 속도의 비선형 변화를 재현하지 못함을 시사하여 잡질이나 유동 효과가 원인일 수 있음을 시사한다.

ABSTRACT

Even though our theoretical understanding of dendritic solidification is relatively well developed, our current ability to model this process quantitatively remains extremely limited. This is due to the fact that the morphological development of dendrites depends sensitively on the degree of anisotropy of capillary and/or kinetic properties of the solid-liquid interface, which is not precisely known for materials of metallurgical interest. Here we simulate the crystallization of highly undercooled nickel melts using a computationally efficient phase-field model together with anisotropic properties recently predicted by molecular dynamics simulations. The results are compared to experimental data and to the predictions of a linearized solvability theory that includes both capillary and kinetic effects at the interface.

연구 동기 및 목표

• 모르쇠 계면 비등방성이 알려져 있지 않은 상황에서 초냉각된 액체에서의 수상성 응고를 정량적으로 모델링하는 오랜 도전 과제를 해결하기 위해.
• 원자적 기반 유도된 비등방성 계면 특성을 상전이 필드 시뮬레이션에 통합하여 이론 모델과 실험 관측 간 격차를 메우기 위해.
• 니켈에서 널리 관측되는 V–ΔT 곡선의 속도 비선형 변화가 표준 확산 운반 모델로 설명될 수 있는지, 아니면 잡질이나 유동과 같은 추가 효과가 필요하는지 테스트하기 위해.
• 분자역학에서 유도된 정밀한 계면 파arameter를 사용하여 실험 데이터와 안정성 이론 예측과 상전이 필드 모델을 검증하기 위해.

제안 방법

• 상전이 필드 모델은 분자역학 시뮬레이션에서 유도된 비등방성 표면장력 및 운동학 계수를 사용하여 고도로 초냉각된 니켈 액체에서의 수상성 응고를 시뮬레이션한다.
• 모델은 곡률과 운동학 비등방성을 포함한 지브스-톰슨 조건을 사용하며, 원자적 연구에서 유도된 계면 에너지 비등방성 εc = 0.018 및 운동학 비등방성 εk = 0.13을 적용한다.
• 모의 실험은 초냉각도 범위(ΔT/(L/cp) = 0.2에서 0.8)와 계면 두께(λ = 2에서 30)를 다양하게 설정하여 안정 상태에서 속도와 형태를 분석한다.
• 결과는 룸 등과 윌네커 등이 발표한 실험 데이터 및 표면장력과 운동학 효과를 모두 포함한 선형화된 안정성 이론 예측과 비교된다.
• 운동학 비등방성을 변화시켜 수상성과 조밀한 가지가지 형태 간의 전이를 분석하여 안정된 수상성 성장에 대한 임계 조건을 규명한다.
• 3D 시뮬레이션에서 계면 진화와 속도를 추적하며, 스냅샷과 속도 프로파일을 사용하여 형태 안정성과 선택 메커니즘을 분석한다.

실험 결과

연구 질문

• RQ1 원자적 기반 유도된 비등방성 계면 특성을 갖춘 상전이 필드 시뮬레이션은 초냉각된 니켈에서 실험적 수상성 끝단 속도를 정량적으로 재현할 수 있는가?
• RQ2 왜 실험적으로 관측된 V–ΔT 곡선의 속도 비선형 변화가 순수하게 확산 운반 모델을 사용한 표준 상전이 필드 모델에서는 나타나지 않는가?
• RQ3 운동학 비등방성은 특히 조밀한 가지가지 형태로의 전이 근처에서 수상성 성장의 안정성과 형태에 어떤 영향을 미치는가?
• RQ4 표면장력과 운동학 비등방성을 모두 포함할 경우, 안정성 이론 예측이 상전이 필드 시뮬레이션 결과와 어느 정도 일치하는가?
• RQ5 니켈 응고에서 관측된 속도 비선형 변화는 내재된 확산 물리에 기인한 것인가, 아니면 잡질이나 유체역학적 효과와 같은 외부 요인에 기인한 것인가?

주요 결과

• 원자적 기반 유도된 비등방성 파arameter(εc = 0.018, εk = 0.13)를 사용한 상전이 필드 시뮬레이션은 다양한 초냉각도에서 실험적 수상성 끝단 속도를 정량적으로 일치시킨다.
• 실험적으로 관측된 V–ΔT 곡선의 속도 비선형 변화는 표준 확산 모델에서는 재현되지 않지만, 현재의 상전이 필드 모델에서는 재현되지 않아 누락된 물리적 메커니즘이 존재함을 시사한다.
• 실험 데이터와 일치시키기 위해 운동학 비등방성 εk는 크기가 커야 하며(≥0.13) 이는 분자역학 시뮬레이션의 예측을 확인한다.
• 낮은 운동학 비등방성(εk = 0)에서 수상성에서 조밀한 가지가지 형태로의 전이가 발생하며, 임계 초냉각도에서 안정성 해면의 분기 현상으로 나타난다.
• 표면장력과 운동학 비등방성을 모두 포함할 경우, 안정성 이론 예측은 상전이 필드 시뮬레이션 결과와 정량적으로 일치한다.
• 시뮬레이션에서 속도 비선형 변화가 나타나지 않는 것은, 실험적으로 관측된 행동이 내재된 확산 물리보다는 잡질이나 유동 효과에 기인할 수 있음을 시사한다.