[논문 리뷰] A fast synthetic iterative scheme for the stationary phonon Boltzmann transport equation

Chuang Zhang, Songze Chen외 2명|arXiv (Cornell University)|2018. 12. 16.

Thermal properties of materials참고 문헌 53인용 수 3

한 줄 요약

이 논문은 정적 진동자 BTE(보울츠만 운동방정식)를 해결하기 위한 빠른 합성 반복 기법을 제안한다. 이 기법은 BTE의 제0차 및 제1차 모멘트 방정식에서 유도된 밀도적으로 결합된 매크로스코픽 확산 방정식을 사용하며, 위상공간 전반에 걸친 진동자 결합을 향상시켜 수렴 속도를 높인다. 이로 인해 근접 확산 영역에서 기존의 암시적 이산 순서 방법보다 1~3개 정도의 지수 차수 빠르게 100회 이내에 수렴을 달성한다.

ABSTRACT

In this paper, a fast synthetic iterative scheme is developed to accelerate convergence for the implicit DOM based on the stationary phonon BTE. The key innovative point of the present scheme is the introduction of the macroscopic synthetic diffusion equation for the temperature, which is obtained from the zero- and first-order moment equations of the phonon BTE. The synthetic diffusion equation, which is asymptomatically preserving to the Fourier's heat conduction equation in the diffusive regime, contains a term related to the Fourier's law and a term determined by the second-order moment of the distribution function that reflects the non-Fourier heat transfer. The mesoscopic kinetic equation and macroscopic diffusion equations are tightly coupled together, because the diffusion equation provides the temperature for the BTE, while the BTE provides the high-order moment to the diffusion equation to describe the non-Fourier heat transfer. This synthetic iterative scheme strengthens the coupling of all phonons in the phase space to facilitate the fast convergence from the diffusive to ballistic regimes. Typical numerical tests in one-, two-, and three-dimensional problems demonstrate that our scheme can describe the multiscale heat transfer problems accurately and efficiently. For all test cases convergence is reached within one hundred iteration steps, which is one to three orders of magnitude faster than the traditional implicit DOM in the near-diffusive regime.

연구 동기 및 목표

• 정적 진동자 BTE에 대한 암시적 이산 순서 방법(DOM)이 근접 확산 영역에서 수렴 속도가 느린 문제를 해결한다.
• 다양한 波벡터와 공간 위치 간의 진동자를 결합하는 데에 있어 전통적인 반복 기법의 비효율성을 해결한다.
• 미세역학적 운동론과 매크로스코픽 확산 물리학을 통합한 합성 기법을 개발하여 수렴성과 정확도를 향상시킨다.
• 1D, 2D, 3D 문제에서 확산 영역에서 구동 영역까지의 다스케일 열전달 시뮬레이션을 효율적이고 정확하게 수행할 수 있도록 한다.

제안 방법

• 진동자 BTE의 제0차 및 제1차 모멘트 방정식에서 유도된 합성 매크로스코픽 확산 방정식을 도출하여, 확산 극한에서 푸아송의 법칙이 점점 유지되도록 보장한다.
• 확산 방정식의 열류를 온도 기울기에 의존하는 푸아송 항과 분포 함수의 제2차 모멘트에 의존하는 비푸아송 항으로 분해한다.
• 미세역학적 BTE와 매크로스코픽 확산 방정식을 밀접하게 결합한다: BTE는 확산 방정식에 제2차 모멘트를 제공하고, 확산 방정식은 BTE에 온도를 공급한다.
• 반복 알고리즘을 구현하여 각 이산 파동벡터에 대해 BTE를 풀고, 모멘트 정보를 사용해 확산 방정식을 갱신함으로써 빠른 전역 수렴을 가능하게 한다.
• 고해상도 3D 시뮬레이션을 처리하기 위해 고체 각도 공간에 걸쳐 MPI 병렬 처리를 구현한다.
• 각 반복 단계에서 모멘트 기반 갱신을 강제하여 에너지 보존과 수치적 안정성을 확보한다.

실험 결과

연구 질문

• RQ1 근접 확산 영역에서 정적 진동자 BTE에 대한 암시적 DOM에 대해 합성 반복 기법이 수렴 속도를 크게 향상시킬 수 있는가?
• RQ2 모멘트 방정식에서 유도된 매크로스코픽 확산 방정식과 미세역학적 BTE를 결합함으로써 수렴 속도와 정확도는 어떻게 향상되는가?
• RQ3 합성 기법이 수렴 횟수를 줄이면서도 전체 운반 영역(확산에서 구동까지)에서 정확도를 얼마나 잘 유지하는가?
• RQ4 이 기법은 최소한의 계산 오버헤드로 대규모 3D 다스케일 열전달 문제를 효율적으로 처리할 수 있는가?
• RQ5 실제 3D 구성에서 기존의 암시적 DOM에 비해 수렴 속도 향상 정도는 정량적으로 얼마인가?

주요 결과

• 모든 3D 테스트 케이스에서 합성 반복 기법은 최대 12×12×6 μm 영역을 포함한 대규모 문제까지 72~99회의 반복 단계 내에 수렴을 달성한다.
• 1D, 2D, 3D 문제에 걸쳐 다양한 격자 크기와 각도 이산화 수준을 포함한 모든 테스트 구성에서 수렴이 100회 이내에 일관되게 달성된다.
• 근접 확산 영역에서 기존의 암시적 DOM에 비해 수렴 속도가 1~3개 지수 차수 빠르게 향상된다.
• 기법은 온도가 높은 소스 근처에서의 온도 점프와 같은 비푸아송 열전달 효과를 정확히 포착하여, 이 지역에서 푸아송의 법칙이 실패함을 시사한다.
• 합성 확산 방정식은 정확히 푸아송 극한을 복원하며, 확산 및 비확산 열전달 물리 현상을 모두 정확하게 모델링한다.
• MPI 병렬 처리를 활용하여 192개 코어에서 효율적인 계산이 가능해져 고해상도로 대규모 3D 시뮬레이션을 실현할 수 있게 되었다.