니코어싱을 억제하는 새로운 상 필드 시뮬레이션 접근법

초록

본 연구는 실험적으로 얻은 SOFC 양극 미세구조를 초기 조건으로 삼아 상‑필드 모델을 통해 니켈‑이산화이트륨(YSZ) 양극의 코어싱(Ki coarsening) 동역학을 정량적으로 예측한다. 샘플 크기와 접촉각 측정 오차를 검토하고, 활성·폐쇄·고립된 상 클러스터별 TPB 길이와 니켈 표면적, 기공 토르성(Tortuosity)의 변화를 분석한다. 결과는 큰 클러스터가 느리게 성장하고, 코어싱이 전반적인 전해질 전도성에는 미미한 긍정적 영향을 주지만, TPB는 국부 형태에 민감하게 감소한다는 점을 보여준다.

상세 분석

이 논문은 고체산화물 연료전지(SOFC) 양극에서 니켈(니) 입자들의 코어싱이 전극 성능 저하의 주요 원인임을 전제로, 상‑필드(Phase‑Field) 방법을 이용해 코어싱 메커니즘을 정량적으로 규명한다. 초기 구조는 실제 제조된 양극의 단면을 X‑ray CT 등으로 획득한 3차원 이미지에서 추출했으며, 이를 격자 기반 상‑필드 모델에 직접 매핑함으로써 실험과 시뮬레이션 간의 격차를 최소화했다. 모델은 Cahn‑Hilliard 방정식과 Allen‑Cahn 방정식을 결합해 물질 이동과 계면 이동을 동시에 기술하고, 니와 YSZ 사이의 계면 에너지와 접촉각을 파라미터화하였다. 특히, 접촉각 측정 시 발생할 수 있는 오차를 통계적으로 분석해 모델 민감도에 미치는 영향을 정량화하였다. 샘플 크기 효과에 대해서는 여러 부피(10 µm³~100 µm³)에서 시뮬레이션을 수행해 경계 조건과 통계적 대표성 사이의 trade‑off를 평가했으며, 충분히 큰 부피에서는 TPB 길이와 니 표면적이 수렴함을 확인했다.

코어싱 과정에서 클러스터를 ‘활성(active)’, ‘폐쇄(dead‑end)’, ‘고립(isolated)’ 세 종류로 구분하고, 각각의 부피, 표면적, TPB 기여도를 추적했다. 결과는 특성 길이(클러스터 평균 직경)가 큰 클러스터일수록 곡률 구동 코어싱 속도가 낮아, 전체 니 표면적 감소율이 예상보다 작다는 것을 보여준다. 반면, 작은 클러스터는 빠르게 사라지면서 TPB 길이는 비선형적으로 감소한다. TPB 감소는 전기화학 반응 면적 손실로 직결되지만, 전체 니 표면적 감소는 상대적으로 작아 전극 전도성에 미치는 영향은 제한적이다.

기공 구조 변화는 토르성(tortuosity) 분석을 통해 정량화했으며, 코어싱이 진행될수록 기공 연속성이 향상되고 토르성이 약 5~8 % 감소한다. 이는 가스 확산 저항 감소와 연관되어 전반적인 전극 전도성에 미세한 긍정적 효과를 가져온다. 그러나 TPB 감소가 전기화학 반응 속도에 미치는 부정적 영향이 더 크게 작용한다는 결론을 도출했다.

이러한 결과는 양극 설계 시 ‘큰 클러스터 유지 + 작은 클러스터 최소화’ 전략이 코어싱에 대한 내성을 높이고 TPB 손실을 억제할 수 있음을 시사한다. 또한, 상‑필드 모델에 실험 기반 초기 구조와 접촉각 보정값을 적용함으로써, 실제 제조 공정에서 발생하는 미세구조 변동성을 반영한 예측이 가능함을 입증했다. 향후 연구에서는 고온·고전압 조건에서의 전기화학적 반응과 연계된 다중 물리 모델을 결합해, 코어싱과 전극 성능 저하 사이의 인과관계를 더욱 정밀하게 규명할 필요가 있다.