위상장 모델(phase-field model)은 계면(interface) 문제를 해결하기 위해 사용되는 수학적 모델이다. 주로 응고(solidification) 역학 분야에 적용되어 왔으며, 그 외에도 점성 손가락화(viscous fingering), 파괴 역학, 수소 취성화, 그리고 소포(세포) 역학 등의 다양한 분야에 적용되었다.
정의와 기본 원리
위상장 모델은 계면에서의 경계 조건을 명시적으로 처리하는 대신, 보조 필드인 위상장(phase field)의 진화를 나타내는 편미분 방정식을 통해 문제를 해결하는 방법이다. 위상장은 질서 매개변수(order parameter) 역할을 하며, 각 상에서 두 가지 뚜렷한 값(예: +1과 −1)을 가진다. 계면 주변 영역에서는 두 값 사이의 부드러운 변화가 나타나, 유한한 너비를 가진 확산(diffuse) 상태를 이룬다. 계면의 이산적 위치는 위상장이 특정 값(예: 0)을 취하는 점들의 집합으로 정의될 수 있다.
위상장 모델은 일반적으로 계면 너비가 0으로 수렴하는 소위 "날카로운 계면 한계(sharp-interface limit)"에서 올바른 계면 역학이 복원되도록 구성된다. 이 접근법은 전체 시스템에 대한 편미분 방정식 세트를 통합함으로써 문제를 해결할 수 있게 하므로, 계면에서의 경계 조건을 별도로 처리할 필요가 없다는 장점이 있다.
역사
위상장 모델은 조지 픽스(George Fix)와 J. S. 랭거(Langer)에 의해 처음 소개되었다. 당시 랭거는 이 방법이 계면 두께가 일반적인 미세 구조의 크기에 비해 너무 작아서 실용적 용도가 없다고 판단하여 관련 연구를 출판하지 않았다. 그러나 이후 컴퓨터 성능의 향상과 함께 카르마(Karma)와 라펠(Rappel)이 얇은 계면 한계(thin-interface limit)를 도입하면서 위상장 모델을 사용한 실용적인 정량 시뮬레이션이 가능해졌다. 그 후 위상장 모델은 계면 문제의 수치 시뮬레이션을 위한 유용한 도구로 자리 잡았다.
수학적 구성
위상장 모델의 여러 공식은 질서 매개변수(위상장)와 확산 필드에 의존하는 열역학적 자유 에너지 함수(변분 공식)를 기반으로 한다. 모델의 방정식은 통계물리학의 일반적인 관계를 사용하여 얻어진다. 응고 문제의 간단한 경우, 자유 에너지 함수는 계면 에너지에 해당하는 기울기 항, 각 상의 벌크 자유 에너지 밀도를 나타내는 이중 우물 전위, 그리고 엔탈피 관련 항을 포함한다. 위상장의 진화 방정식과 확산 방정식은 이러한 자유 에너지 함수에 대한 변분 관계로부터 도출된다.
다른 공식(비변분 공식)은 열역학적 함수를 참조하지 않고 상장 방정식을 직접 작성하는 방식으로, 이 경우 유일한 참조는 날카로운 계면 모델이다.
주요 응용 분야
위상장 모델은 다음을 포함한 다양한 분야에 적용된다.
- 응고 역학: 나뭇가지 모양(dendrite) 성장 등 결정 성장 과정의 시뮬레이션
- 재료 공학: 고체상 변태, 결정립 성장, 재결정화, 오스테나이트에서 페라이트로의 상변태 등 미세 조직 진화의 예측
- 파괴 역학: 균열의 시작(nucleation), 성장, 분기(branching), 합체(coalescence) 과정을 모델링
- 수소 취성화: 수소 유도 균열(hydrogen-assisted cracking) 현상의 분석
- 세포 생물학: 집단 세포 이동(collective cell migration), 소포(vesicle) 역학 등
- 유체 역학: 점성 손가락화, 이상(non-miscible) 유체 유동, 이중상 나비에-스토크스 흐름의 모델링
다중 위상장 모델
다중 위상장(multiphase-field) 모델에서는 미세 구조가 질서 매개변수의 집합으로 설명되며, 각 매개변수는 특정 상 또는 결정학적 방향과 관련된다. 이 모델은 여러 결정립이 진화하는 과정, 예를 들어 결정립 성장이나 철 합금의 상변태와 같은 문제에 주로 사용된다. 여러 결정립을 설명할 수 있을 뿐만 아니라, 기술 합금 등급에서 발생하는 여러 열역학적 상을 고려할 수 있게 한다.
관련 소프트웨어
위상장 모델을 구현한 시뮬레이션 소프트웨어로는 PACE3D, MICRESS, MOOSE, OpenPhase, PRISMS-PF 등 다양한 패키지가 존재한다. 이들은 오픈 소스 또는 상업적으로 제공되며, 재료 과학 및 공학 분야의 미세 조직 시뮬레이션, 파괴 해석 등에 활용된다.
비판적 고려 사항
위상장 모델은 확산 계면(diffuse interface)을 사용하므로, 실제 계면 두께가 문제의 길이 스케일에 비해 매우 작은 경우 수치 계산에 상당한 계산 비용이 요구된다. 이러한 문제를 완화하기 위해 얇은 계면 한계 이론과 적응 메쉬(adaptive mesh) 기법 등이 개발되었다. 또한 자유 에너지 함수의 선택이 계면의 물리적 거동에 상당한 영향을 미칠 수 있어, 함수의 선정을 신중히 하여야 한다는 점이 지적된다.
참고 문헌
- Chen, Long-Qing (2002). "Phase-Field models for microstructure evolution", Annual Review of Materials Research 32, 113–140.
- Moelans, Nele; Blanpain, Bart; Wollants, Patrick (2008). "An introduction to phase-field modeling of microstructure evolution", Calphad 32, 268.
- Steinbach, Ingo (2009). "Phase-field models in materials science", Modelling and Simulation in Materials Science and Engineering 17.
- Provatas, Nikolas; Elder, Ken (2010). Phase-Field Methods in Materials Science and Engineering, Wiley-VCH.
