[논문 리뷰] GiftBTE: An efficient deterministic solver for non-gray phonon Boltzmann transport equation

Yue Hu, Ru Jia외 6명|arXiv (Cornell University)|2023. 06. 25.

Thermal properties of materialsMaterials Science인용 수 3

한 줄 요약

GiftBTE는 결정성 물질에서의 서미크론 온도 전도를 효율적이고 고정밀도로 시뮬레이션할 수 있도록 해주는 비회색 열파동 보른졸루션 방정식(BTE)에 대한 오픈소스이자 결정론적 해법기이다. 이는 정 steady-state 문제에 대해 이차 정확도를 갖는 음성 이산 체계 방법과 합성 반복 기법을 사용하고, 일시적 문제에 대해서는 명시적 기법을 사용하여 뛰어난 계산 효율성을 달성하며, 나노구조물, 트랜지스터, 레이저 가열 과정 등의 현실적인 3D 시뮬레이션을 가능하게 한다.

ABSTRACT

Advances in nanotechnology have facilitated the exploration of submicron thermal transport. At this scale, Fourier's law is no longer applicable, and the governing equation for thermal transport is the phonon Boltzmann transport equation (BTE). However, the availability of open-source solvers for the phonon BTE is limited, impeding progress in this field. This study introduces an open-source package, GiftBTE, for numerically solving the non-gray phonon BTE. GiftBTE employs deterministic solutions and provides both steady-state and transient solvers. For the steady-state solver, GiftBTE employs the implicit discrete ordinates method (DOM) with second-order spatial accuracy and the synthetic iterative scheme. For the transient solver, GiftBTE employs the explicit DOM with second-order spatial accuracy. This package demonstrates excellent computational efficiency, enabling realistic three-dimensional simulations of devices and materials. By interfacing with first-principles calculations, this solver enables parameter-free computation of submicron thermal transport. The application of GiftBTE includes, but is not limited to, computing the thermal conductivity of nanostructures, predicting temperature rises in transistors, and simulating laser heating processes.

연구 동기 및 목표

• 복잡한 기하학적 구조에서 비회색 열파동 BTE에 대한 일반 목적의 오픈소스 해법기가 부족한 문제를 해결한다.
• 기존 해법기가 밀도 또는 단순 기하학적 구조에 국한되어 있고, 일시적 해석이나 매개변수 없는 기능이 부족한 문제를 해결한다.
• 초기원리적 열파동 성질을 이용하여 실제 3D 재료 및 장치에서 서미크론 온도 전도를 정확하고 효율적이며 확장 가능한 방식으로 시뮬레이션할 수 있도록 한다.
• 나노전자소자 및 열전재료의 열전도 분석을 위한 CAE 유사 도구를 제공하여, 열전도도 예측 및 트랜지스터에서의 온도 상승 추정 등의 응용을 지원한다.
• 원자 척도 방법(MD, NEMD 등)과 매크로스코픽 열확산 간의 격차를 메우며, 서미크론 척도에서 결정론적이고 비푸아르르 방식의 시뮬레이션을 가능하게 한다.

제안 방법

• 비회색 열파동 BTE를 해소하기 위해 이산 체계 방법(DOM)에 기반한 결정론적 접근 방식을 사용하여 몬테카를로 방법에서 유래하는 통계적 노이즈를 피한다.
• 정상 상태 및 일시적 해법기 모두에서 이차 정확도의 공간 해석을 적용하여 메쉬 요구량을 줄이고 해석 정확도를 향상시킨다.
• 근접 확산 영역에서 표준 반복 기법이 느리게 수렴하는 문제를 해결하기 위해 정상 상태 해법기에 합성 반복 기법(SIS)을 적용한다.
• 일시적 시뮬레이션을 위해 이차 정확도를 갖는 명시적 DOM 설정을 사용하여 시간 해상도를 확보한 열역학적 동역학 분석을 가능하게 한다.
• phono3py 등과 같은 도구를 통해 초기원리적 열파동 계산에 연결하여 온도 및 波벡터 의존성 열파동 성질을 경험적 피팅 없이 확보한다.
• 임의의 결정성 재료 및 복잡한 3D 기하학적 형태(예: 나노다공성 매체, 박막, 장치 크기의 구조)를 지원한다.

실험 결과

연구 질문

• RQ1 결정론적이고 이차 정확도를 갖는 BTE 해법기가 현실적인 3D 서미크론 온도 전도 문제에 대해 높은 계산 효율성과 정확도를 달성할 수 있는가?
• RQ2 합성 반복 기법이 표준 반복 기법에 비해 근접 확산 영역에서 수렴 속도를 어떻게 향상시키는가?
• RQ3 비푸아르르 효과가 시간 영역 열반사도(TDTR) 실험과 같은 일시적 열응답에서 얼마나 뚜렷하게 나타나는가?
• RQ4 초기원리적 입력 자료를 사용하여 나노스케일 트랜지스터 및 레이저 가열 재료에서의 온도 상승을 고정밀도로 예측할 수 있는가?
• RQ5 GiftBTE 프레임워크가 단순 기하학적 형태를 초월한 복잡한 나노구조에서 매개변수 없는 열전도 시뮬레이션을 가능하게 하는가?

주요 결과

• 이차 정확도의 공간 해석과 합성 반복 기법을 갖춘 정상 상태 해법기는 일차 정확도 방법에 비해 메쉬 수요를 크게 줄여 효율성을 향상시킨다.
• 일시적 해법기는 비푸아르르 열역학적 동역학을 성공적으로 포착하여, TDTR 시뮬레이션에서 푸아르르 법보다 느린 온도 감쇠를 보이며, 뚜렷한 복사 운동 효과를 시사한다.
• TDTR 사례 연구에서 표면에서의 온도 상승은 열확산 방정식보다 열파동 BTE에서 더 천천히 감쇠되었으며, 64코어 클러스터에서 16시간이 소요된 시뮬레이션을 수행하였다.
• 초기원리적 열파동 데이터를 전부 활용하여 박막 및 다공성 재료와 같은 나노구조에서의 열전도도를 정확하게 계산할 수 있다.
• 모듈레이션된 펄스를 갖는 레이저 가열 과정을 성공적으로 시뮬레이션하여, 현실적인 빔 프로파일과 침투 깊이를 반영한 온도 반응에서 실험 유사 행동을 재현하였다.
• 복잡한 장치, 예를 들어 트랜지스터 및 열전 홀로구조에서의 3D 시뮬레이션에 대해 확장성과 견고성을 입증하였으며, 열파동 모드 해상도의 열류 및 온도 분포에 대한 전부 접근이 가능하다.