융합용 텅스텐의 방사선‑구동 재결정 메커니즘을 예측하는 다중스케일 모델

초록

**
본 논문은 결정소성, 확률적 클러스터 동역학, 그리고 정점 기반 입계 역학을 결합한 물리 기반 다중스케일 프레임워크를 개발한다. 이를 통해 고온 가공·방사선 손상·입계 이동이 동시에 일어나는 실제 텅스텐 다결정 미세구조의 변형·재결정·입계 성장 과정을 3‑D(2‑D) 수준에서 직접 시뮬레이션한다. 온도가 입계 이동성을 지배하고, 변형 전 이력은 약하게, 중성자에 의한 Re 원소 침착은 재결정을 크게 억제한다는 핵심 결과를 제시한다.

**

상세 분석

**
이 연구는 융합 플라즈마와 직접 접촉하는 텅스텐(PFC) 재료가 직면하는 세 가지 복합 현상을 하나의 통합 모델로 풀어낸 점에서 혁신적이다. 첫째, 유한요소 기반 결정소성(FE‑CP) 모듈은 각 입계마다 서로 다른 결정학적 방향성을 고려해 변형률·전위·전위 경화(격자 마찰, 포레스트 경화, 입계 제한 경화, 방사선 결함 경화)를 정량화한다. 특히 방사선에 의해 생성된 점·클러스터 결함을 ‘분산 장벽 경화’ 모델(τ_irr)로 구현함으로써 전위 흐름이 손상 축적에 따라 급격히 감소하는 메커니즘을 재현한다.

둘째, 확률적 클러스터 동역학(SCD) 모듈은 중성자 플럭스, 온도, 전위 상태를 입력으로 받아 점결함·다중클러스터(Vo, SI, 크루시블 등)의 생성·성장·소멸을 시뮬레이션한다. 이때 전위에 의한 재결합률(k₂)과 온도 의존성(Arrhenius 형태)이 결합되어, 고온·고전위 조건에서 결함이 빠르게 소멸하고 저장에너지(전위·점결함) 축적이 억제되는 현상을 포착한다.

셋째, 정점 동역학(Vertex Dynamics, VD) 모듈은 2‑D 입계 네트워크를 삼각형 정점과 삼중점(TJ)으로 이산화하고, 화학적 구동력(입계 에너지 차, Re·W 전이 원소 농도 구배)과 탄성 구동력(전위·점결함에 의한 응력 차이)을 기반으로 입계 이동 속도 v = M·f 를 계산한다. 이동도 M은 온도 의존적 Arrhenius 식으로 정의되며, Re 침착에 의한 Zener 고정 효과를 추가해 입계 이동성을 실질적으로 감소시킨다.

모듈 간 상호작용은 순환적으로 수행된다. FE‑CP가 제공하는 전위·전위 밀도장은 SCD에 입력되어 결함 생성률을 조정하고, SCD가 산출한 결함 농도·크기는 FE‑CP의 경화 항 τ_irr에 반영된다. 동시에, 전위·결함에 의해 발생한 저장에너지와 화학적 구배는 VD에 전달되어 입계 구동력을 산출한다. 이러한 피드백 루프는 시간·공간적으로 이질적인 재결정 핵생성·성장 과정을 자연스럽게 재현한다.

시뮬레이션 결과는 세 가지 주요 인사이트를 제공한다. (1) 온도는 입계 이동성 M의 지배 변수이며, 1200 K 이상에서 재결정 속도가 급격히 증가한다. (2) 사전 변형(전위 밀도) 은 재결정 초기 핵생성률을 약간 높이지만, 온도에 비해 영향이 미미하다. (3) 중성자에 의해 생성된 Re 원소가 입계에 축적되면, 화학적 구동력 ψ가 감소하고 입계 이동도가 1‑2 order of magnitude 낮아져 재결정이 현저히 지연된다. 특히, 방사선 환경에서 비방사선 상태 대비 재결정 온도가 200 °C 이하로 낮아지는 반면, Re 축적이 이를 상쇄하거나 오히려 100 °C 이상 상승시킨다.

이 모델은 기존 평균장(mean‑field) 접근법이 가정하던 고정된 입계 미스오리엔테이션·핵생성률을 탈피하고, 실제 입계 곡률·다중입계 상호작용·불균일 전위·결함 분포를 그대로 반영한다. 따라서 실험적 Avrami 지수 n≈1 (핵생성·성장이 동시에 진행되는 방사선‑보조 재결정)과도 정량적으로 일치한다. 또한, 모델 파라미터(τ_l, k₁, M₀, Q_m 등)를 실험 데이터와 역학적으로 보정함으로써, 향후 설계 단계에서 합금 조성·열처리·작동 온도 창을 최적화하는 예측 도구로 활용 가능하다.

**