[논문 리뷰] Multiscale Phonon Conduction in Nanostructured Materials Predicted by Bulk Thermal Conductivity Accumulation Function

Giuseppe Romano, Jeffrey C. Grossman|arXiv (Cornell University)|2013. 12. 30.

Thermal properties of materials참고 문헌 24인용 수 5

한 줄 요약

이 논문은 단일의 밀도 분포를 입력으로 사용하여 나노구조 재료의 열전도도를 효율적이고 정확하게 예측할 수 있는 다중스케일 열역학적 진동수 이동 모델인 MFP-BTE를 제안한다. 이 모델은 진동수 분포나 산란 시간 데이터를 요구하지 않으며, 구조적 복잡성에 관계없이 다양한 스케일에서의 열전도도 예측이 가능하다. 다공성 실리콘 막에서의 실험 측정치와 양호한 일치를 보이며, 메조스케일 및 나노스케일 영역에서의 적용 가능성을 입증한다.

ABSTRACT

We develop a computational framework, based on the Boltzmann transport equation, with the ability to compute the thermal transport in nanostructured materials of any geometry using as the only input the bulk thermal conductivity accumulation function. The main advantage of our method is twofold. First, while the scattering times and dispersion curves are unknown for most materials, the phonon mean free path distribution can be directly obtained by experiments. As a consequence, a wider range of materials can be simulated than with a frequency-dependent approach. Second, when phonon dispersions are available from first principles calculations, our approach allows one to include easily the whole Brillouen zone in the calculations without discretizing the phonon frequencies for all polarizations, reducing considerably the computational effort. Furthermore, after deriving the ballistic and diffusive limits of our model, we develop a multi-scale scheme that couples phonon transport across different scales, enabling efficient simulations of materials with wide phonon mean free path distributions. After we validate our model against the frequency-dependent approach, we apply the method to porous silicon membranes and find good agreement with experiments on mesoscale pores. By enabling the investigation of thermal transport in unexplored nanostructured materials, our method has potential to advance high-efficiency thermoelectric devices.

연구 동기 및 목표

• 나노구조 재료의 진동수 이동 열전도도를 예측하기 위한 계산적으로 효율적인 방법을 개발함으로써, 진동수 분포나 산란 시간 데이터에 의존하지 않도록 하는 것.
• 진동수 의존성 BTE(FD-BTE)의 한계를 극복함. FD-BTE는 세밀한 진동수 특성과 산란 시간 데이터를 필요로 하며, 특히 복잡한 재료에서는 높은 계산 비용을 유발함.
• 복잡한 격자 구조나 넓은 MFP 분포를 가진 재료의 시뮬레이션을 가능하게 하기 위해 실험적으로 측정 가능한 기초 MFP 분포를 활용함.
• MFP-BTE를 통해 구동되는 일관된 다중스케일 프레임워크를 수립하여, 구동, 산란, 중간 이동 영역을 Knudsen 수 기반으로 자연스럽게 연결함.
• 다양한 치수와 다공성 비율을 가진 다공성 실리콘 막에서의 실험 데이터와의 비교를 통해 모델의 예측 능력과 광범위한 적용 가능성을 검증함.

제안 방법

• 이 방법은 정적인 BTE를 진동수와 분극도 의존성 변수 대신 평균 자유로운 진동수(MFP) 기반 분포로 재구성함.
• 모델은 유일한 입력으로 기초 열전도도 누적 함수를 사용하며, 이는 실험적으로 직접 측정하거나 첫 원리 계산으로 산출할 수 있음.
• 모델은 구동 이동을 위한 구동 BTE 성분과 산란 이동을 위한 푸아송 법칙을 통합하며, Knudsen 수 기반 임계값에 따라 영역 간 전이를 제어하는 다중스케일 결합 기법을 적용함.
• FD-BTE와 비교해 진동수와 분극도의 이산화를 피함으로써, 특히 복잡한 브릴루앙 영역을 가진 재료에서 계산 비용을 크게 절감함.
• 넓은 MFP 분포를 다루기 위해 MFP 범위를 구동 영역, 전이 영역, 산란 영역으로 나누고 각 영역에 적절한 운동 모델을 적용하는 다중스케일 기법을 구현함.
• 다공성 실리콘 막에서의 다양한 기하구조와 기공 크기를 가진 시료에 대해 FD-BTE 예측값과 실험 데이터와의 비교를 통해 모델의 타당성을 검증함.

실험 결과

연구 질문

• RQ1 나노구조 재료에서의 진동수 이동 열전도도는 진동수 분포나 산란 시간 데이터 없이 기초 MFP 분포만으로 정확하게 예측할 수 있는가?
• RQ2 기존의 진동수 의존성 BTE(FD-BTE)에 비해 제안된 MFP-BTE 모델은 정확도와 계산 효율성 측면에서 어떻게 비교되는가?
• RQ3 다공성 재료에서의 고전적 치수 효과는 MFP-BTE 모델을 통해 어떤 정도로 나노스케일 및 메조스케일 영역에서 잘 포괄되는가?
• RQ4 Knudsen 수 기반 임계값을 기준으로 하여, 구동, 산란, 중간 이동 영역 간의 최적의 다중스케일 결합 전략은 무엇인가?
• RQ5 다공성 비율과 기하구조가 다양한 다공성 실리콘 막에서 MFP-BTE 모델은 실험적 열전도도 값을 얼마나 잘 예측하는가?

주요 결과

• MFP-BTE 모델은 메조스케일 다공성 실리콘 막에서 실험 데이터와 양호한 정량적 일치를 보이며, 예측한 효과적 열전도도는 56 W/mK이며 실험값 45 W/mK와 유사함.
• L=34 nm, φ=0.173 인 나노스케일 Si 나노메쉬의 경우 모델은 효과적 열전도도를 약 8 W/mK로 예측했으나, 실험값 2.85 W/mK보다 약간 높게 나옴. 이는 모델이 반응하지 않는 비정질 계면층 또는 파동 효과 기여 요소가 있을 수 있음을 시사함.
• 나노스케일 시스템의 경우 최적의 Knudsen 수 기준은 Kn_B = 10 (구동 영역) 및 Kn_D = 0.1 (산란 영역)이며, 메조스케일 시스템에서는 표면 산란 효과를 고려해 Kn_D = 0.025로 설정됨.
• 특히 복잡한 격자 단위세포를 가진 재료에서 FD-BTE보다 계산 효율성이 뛰어남. 이는 진동수 및 분극도 이산화를 피하고 MFP 샘플링 기반으로 확장되기 때문임.
• 모델은 다공성 실리콘에서의 구동에서 산란으로의 전이를 성공적으로 포착하며, 상하 표면에서의 산란으로 인해 열류가 중심 영역에 집중됨.
• 이 방법은 어떤 재료에도 일반적으로 적용 가능하며, 기초 MFP 분포가 알려져 있으면서, 열전도성 응용을 위한 새로운 나노구조 재료의 신속한 선별이 가능함.