Coupling Mesoscopic Boltzmann Transport Equation and Macroscopic Heat Diffusion Equation for Multiscale Phonon Heat Conduction

• Published in: Nanoscale and microscale thermophysical engineering Vol. 24; no. 3-4; pp. 150 - 167
• Main Authors: Cheng, W., Alkurdi, A., Chapuis, P.-O.
• Format: Journal Article
• Language: English
• Published: Abingdon Taylor & Francis 01-10-2020; Taylor & Francis Ltd
• Subjects: Boltzmann equation; Boltzmann transport equation; Cartesian coordinates; Computing time; Conduction heating; Conductive heat transfer; Convergence; discrete ordinate method; Domains; Engineering Sciences; Finite element method; Iterative methods; Mathematical analysis; Mechanics; Multiscale; phonon heat conduction; Phonons; Radiative transfer; Solvers; Thermics; Thermodynamics; Transport equations; Boltzmann equation; Phonon heat conduction; Multiscale; Discrete Ordinate Method
• ISSN: 1556-7265; 1556-7273
• DOI: 10.1080/15567265.2020.1836095

Phonon heat conduction has to be described by the Boltzmann transport equation (BTE) when sizes or sources are comparable to or smaller than the phonon mean free paths (MFPs). When domains much larger than MFPs are to be treated or when regions with large and small MFPs coexist, the computation time associated with full BTE treatment becomes large, calling for a multiscale strategy to describe the total domain and decreasing the computation time. Here, we describe an iterative method to couple the BTE, under the Equation of Phonon Radiative Transfer approximation solved by means of the deterministic Discrete Ordinate Method, to a Finite-Element Modeling commercial solver of the heat equation. Small-size elements are embedded in domains where the BTE is solved, and the BTE domains are connected to a domain where large-size elements are located and where the heat equation is applied. It is found that an overlapping zone between the two types of domains is required for convergence, and the accuracy is analyzed as a function of the size of the BTE domain. Conditions for fast convergence are discussed, leading to the computation time being divided by more than five on a study case in 2D Cartesian geometry. The simple method could be generalized to other types of solvers of the Boltzmann and heat equations.