[논문 리뷰] Phase Field Methods and Dislocations

David Rodney, A. Finel|arXiv (Cornell University)|2001. 01. 16.

Microstructure and mechanical properties인용 수 6

한 줄 요약

이 논문은 응력-free 변형이 그들의 버거스 벡터와 고리 법선과 동일한 등가인 판상 포함체로 디스locations를 모델링하는 단계장 방법을 제시한다. 이는 디스클로케이션 운동과 상전이 사이의 동적 결합을 가능하게 하며, 복잡한 디스클로케이션 역학, 예를 들어 프랭크-리드 소스 작동과 농도장 내에서의 고리 수축을 성공적으로 시뮬레이션한다. 이는 3차원 이방성 매질에서 탄성장과 미세구조 진화 사이의 견고한 결합을 보여준다.

ABSTRACT

We present a general formalism for incorporating dislocations in Phase Field methods. This formalism is based on the elastic equivalence between a dislocation loop and a platelet inclusion of specific stress-free strain related to the loop Burgers vector and normal. Dislocations are thus treated as platelet inclusions and may be coupled dynamically to any other field such as a concentration field. The method is illustrated through the simulation of a Frank-Read source and of the shrinkage of a loop in presence of a concentration field.

연구 동기 및 목표

• 디스클로케이션을 단계장 방법에 통합하기 위한 일반적인 수학적 형식을 개발하는 것.
• 특히 침전상과 디스클로케이션을 포함하는 복잡한 미세구조에 대해 동적 디스클로케이션 모델링의 부족을 해결하는 것.
• 실제 탄성 상호작용을 고려한 3차원 탄성 이방성 매질에서의 디스클로케이션 역학을 시뮬레이션할 수 있도록 하는 것.
• 디스클로케이션 소스와 변화하는 농도장과의 상호작용을 모델링한 방법의 능력을 입증하는 것.
• 메조스케일에서 결합된 디스클로케이션과 상전이 역학을 시뮬레이션하기 위한 통합된 프레임워크를 제공하는 것.

제안 방법

• 디스클로케이션은 응력-free 변형 텐서 $\epsilon^{oo}_{ij} = \frac{1}{2} \frac{b_i n_j + b_j n_i}{d}$ 를 가진 판상 포함체로 모델링되며, 여기서 $b$ 는 버거스 벡터, $n$ 은 고리 법선, $d$ 는 두께이다.
• 디스클로케이션 고리와 판상 포함체 사이의 탄성 등가성을 이용하여, 상전이와 동일한 수학적 형식을 통해 디스클로케이션을 단계장 프레임워크에 통합한다.
• 이 방법은 카하투리안의 단계장 미세탄성론을 사용하며, 탄성 에너지를 밀도장의 푸리에 변환의 기능으로 표현한다.
• 디스클로케이션 역학은 총 자유 에너지에서 유도된 운동 방정식에 의해 지배되며, 이는 화학적 및 탄성 기여를 포함한다.
• 시스템은 디스클로에이션 선 주변의 핵심 확산을 최소화하고 높은 응력장을 유지하기 위해 $d = 10b$ 로 설정된 3차원 정사각형 격자에서 수치적으로 해결된다.
• 농도에 의존하는 이동도를 가진 보존형 TDGL 방정식을 사용하여 농도장의 진화를 모델링하며, 이는 디스클로에이션과 상장에서 유래한 탄성 에너지와 결합된다.

실험 결과

연구 질문

• RQ1 3차원 이방성 탄성 매질에서 디스클로에이션을 단계장 시뮬레이션에 어떻게 동적으로 결합할 수 있는가?
• RQ2 디스클로에이션 고리와 판상 포함체 사이의 탄성 등가성을 단계장 프레임워크 내에서 디스클로에이션 역학을 모델링하는 데 사용할 수 있는가?
• RQ3 농도장이 존재할 경우 디스클로에이션은 상분리 및 침전핵 생성에 어떻게 영향을 미치는가?
• RQ4 상전이 중에 디스클로에이션의 변형장은 디스클로에이션 고리의 운동과 수축에 어떤 영향을 미치는가?
• RQ5 이 방법은 현실적인 역학을 갖춘 프랭크-리드 소스와 같은 복잡한 디스클로에이션 소스를 시뮬레이션할 수 있는가?

주요 결과

• 이 방법은 프랭크-리드 소스의 작동을 성공적으로 시뮬레이션하였으며, 비틀림 응력 하에서 디스클로에이션 세그먼트의 볼링 아웃 현상을 보여주었다.
• 디스클로에이션 핵심은 약 5개 격자 간격(두께 $d = 10b$ 일 경우 50b)으로 확산되며, 이는 일반적인 핵심 반경보다 크지만 메조스케일 시뮬레이션에는 수용 가능하다.
• 핵심 확산으로 인해 디스클로에이션 선 주변의 최대 응력은 감소하지만, 물리적 현실성을 확보하기 위해 충분한 단거리 상호작용 강도를 유지한다.
• 농도장이 존재할 경우, 탄성 상호작용으로 인해 디스클로에이션 고리가 수축하며, 모서리 성분은 잔류 침전물의 열을 남긴다.
• 디스클로에이션 변형장과 침전물 간의 상호작용은 고리를 锁정하는 경향이 있지만, 강한 내부 응력이 존재하면 고리가 이에 저항하고 면적을 감소시킬 수 있다.
• 이 방법은 공유된 탄성 에너지 항을 통해 디스클로에이션 역학과 상전이 동역학 사이의 자동적이고 자기 일관성 있는 결합을 가능하게 한다.