Research Papers: Radiative Heat Transfer

Solving Nongray Boltzmann Transport Equation in Gallium Nitride

[+] Author and Article Information
Ajit K. Vallabhaneni

G. W. Woodruff School of
Mechanical Engineering,
Georgia Institute of Technology,
Atlanta, GA 30332
e-mail: ajitkvallabhaneni@gmail.com

Liang Chen, Man P. Gupta, Satish Kumar

G. W. Woodruff School of
Mechanical Engineering,
Georgia Institute of Technology,
Atlanta, GA 30332

1Corresponding author.

Contributed by the Heat Transfer Division of ASME for publication in the JOURNAL OF HEAT TRANSFER. Manuscript received December 16, 2015; final manuscript received April 24, 2017; published online June 6, 2017. Assoc. Editor: Alan McGaughey.

J. Heat Transfer 139(10), 102701 (Jun 06, 2017) (8 pages) Paper No: HT-15-1796; doi: 10.1115/1.4036616 History: Received December 16, 2015; Revised April 24, 2017

Several studies have validated that diffusive Fourier model is inadequate to model thermal transport at submicron length scales. Hence, Boltzmann transport equation (BTE) is being utilized to improve thermal predictions in electronic devices, where ballistic effects dominate. In this work, we investigated the steady-state thermal transport in a gallium nitride (GaN) film using the BTE. The phonon properties of GaN for BTE simulations are calculated from first principles—density functional theory (DFT). Despite parallelization, solving the BTE is quite expensive and requires significant computational resources. Here, we propose two methods to accelerate the process of solving the BTE without significant loss of accuracy in temperature prediction. The first one is to use the Fourier model away from the hot-spot in the device where ballistic effects can be neglected and then couple it with a BTE model for the region close to hot-spot. The second method is to accelerate the BTE model itself by using an adaptive model which is faster to solve as BTE for phonon modes with low Knudsen number is replaced with a Fourier like equation. Both these methods involve choosing a cutoff parameter based on the phonon mean free path (mfp). For a GaN-based device considered in the present work, the first method decreases the computational time by about 70%, whereas the adaptive method reduces it by 60% compared to the case where full BTE is solved across the entire domain. Using both the methods together reduces the overall computational time by more than 85%. The methods proposed here are general and can be used for any material. These approaches are quite valuable for multiscale thermal modeling in solving device level problems at a faster pace without a significant loss of accuracy.

Copyright © 2017 by ASME
Your Session has timed out. Please sign back in to continue.


Chen, G. , 2005, Nanoscale Energy Transport and Conversion, Oxford University Press, New York.
Regner, K. T. , Sellan, D. P. , Su, Z. , Amon, C. H. , McGaughey, A. J. , and Malen, J. A. , 2013, “ Broadband Phonon Mean Free Path Contributions to Thermal Conductivity Measured Using Frequency Domain Thermoreflectance,” Nat. Commun., 4, p. 1640. [CrossRef] [PubMed]
Hua, C. , and Minnich, A. J. , 2014, “ Transport Regimes in Quasiballistic Heat Conduction,” Phys. Rev. B, 89(9), p. 094302. [CrossRef]
Hess, K. , 1988, “ Boltzmann Transport Equation,” The Physics of Submicron Semiconductor Devices, Springer, New York, pp. 33–43.
Ziman, J. M. , 1960, Electrons and Phonons: The Theory of Transport Phenomena in Solids, Oxford University Press, New York.
Maradudin, A. A. , 1971, Theory of Lattice Dynamics in the Harmonic Approximation, Academic Press, Cambridge, MA.
Narumanchi, S. V. , Murthy, J. Y. , and Amon, C. H. , 2005, “ Comparison of Different Phonon Transport Models for Predicting Heat Conduction in Silicon-on-Insulator Transistors,” ASME J. Heat Transfer, 127(7), pp. 713–723. [CrossRef]
Narumanchi, S. V. , Murthy, J. Y. , and Amon, C. H. , 2006, “ Boltzmann Transport Equation-Based Thermal Modeling Approaches for Hotspots in Microelectronics,” Heat Mass Transfer, 42(6), pp. 478–491. [CrossRef]
Ali, S. A. , Kollu, G. , Mazumder, S. , Sadayappan, P. , and Mittal, A. , 2014, “ Large-Scale Parallel Computation of the Phonon Boltzmann Transport Equation,” Int. J. Therm. Sci., 86, pp. 341–351. [CrossRef]
Ali, S. A. , and Mazumder, S. , 2015, “ Phonon Heat Conduction in Multidimensional Heterostructures: Predictions Using the Boltzmann Transport Equation,” ASME J. Heat Transfer, 137(10), p. 102401. [CrossRef]
Singh, D. , Murthy, J. Y. , and Fisher, T. S. , 2011, “ Effect of Phonon Dispersion on Thermal Conduction Across Si/Ge Interfaces,” ASME J. Heat Transfer, 133(12), p. 122401. [CrossRef]
Mishra, U. K. , Parikh, P. , and Wu, Y.-F. , 2002, “ AlGaN/GaN HEMTs—An Overview of Device Operation and Applications,” Proc. IEEE, 90(6), pp. 1022–1031. [CrossRef]
Kuzmík, J. , 2001, “ Power Electronics on InAlN/(In) GaN: Prospect for a Record Performance,” IEEE Electron. Device Lett., 22(11), pp. 510–512. [CrossRef]
Wu, Y.-F. , Kapolnek, D. , Ibbetson, J. P. , Parikh, P. , Keller, B. P. , and Mishra, U. K. , 2001, “ Very-High Power Density AlGaN/GaN HEMTs,” IEEE Trans. Electron. Devices, 48(3), pp. 586–590. [CrossRef]
Ni, C. , and Murthy, J. Y. , 2009, “ Parallel Computation of the Phonon Boltzmann Transport Equation,” Numer. Heat Transfer, Part B, 55(6), pp. 435–456. [CrossRef]
Sawetprawichkul, A. , Hsu, P.-F. , and Mitra, K. , 2002, “ Parallel Computing of Three-Dimensional Monte Carlo Simulation of Transient Radiative Transfer in Participating Media,” AIAA Paper No. 2002-2901.
Chui, E. , and Raithby, G. , 1992, “ Implicit Solution Scheme to Improve Convergence Rate in Radiative Transfer Problems,” Numer. Heat Transfer, Part B, 22(3), pp. 251–272. [CrossRef]
Fiveland, V. , and Jessee, J. , 1996, “ Acceleration Schemes for the Discrete Ordinates Method,” J. Thermophys. Heat Transfer, 10(3), pp. 445–451. [CrossRef]
Mazumder, S. , 2005, “ A New Numerical Procedure for Coupling Radiation in Participating Media With Other Modes of Heat Transfer,” ASME J. Heat Transfer, 127(9), pp. 1037–1045. [CrossRef]
Mathur, S. R. , and Murthy, J. Y. , 2009, “ An Acceleration Technique for the Computation of Participating Radiative Heat Transfer,” ASME Paper No. IMECE2009-12923.
Mathur, S. , and Murthy, J. , 1999, “ Coupled Ordinates Method for Multigrid Acceleration of Radiation Calculations,” J. Thermophys. Heat Transfer, 13(4), pp. 467–473. [CrossRef]
Adams, M. L. , and Larsen, E. W. , 2002, “ Fast Iterative Methods for Discrete-Ordinates Particle Transport Calculations,” Prog. Nucl. Energy, 40(1), pp. 3–159. [CrossRef]
Loy, J. M. , Mathur, S. R. , and Murthy, J. Y. , 2015, “ A Coupled Ordinates Method for Convergence Acceleration of the Phonon Boltzmann Transport Equation,” ASME J. Heat Transfer, 137(1), p. 012402. [CrossRef]
Loy, J. M. , Murthy, J. Y. , and Singh, D. , 2013, “ A Fast Hybrid Fourier–Boltzmann Transport Equation Solver for Nongray Phonon Transport,” ASME J. Heat Transfer, 135(1), p. 011008. [CrossRef]
Allu, P. , and Mazumder, S. , 2016, “ Hybrid Ballistic–Diffusive Solution to the Frequency-Dependent Phonon Boltzmann Transport Equation,” Int. J. Heat Mass Transfer, 100, pp. 165–177. [CrossRef]
Baloch, K. H. , Voskanian, N. , Bronsgeest, M. , and Cumings, J. , 2012, “ Remote Joule Heating by a Carbon Nanotube,” Nat. Nanotechnol., 7(5), pp. 316–319. [CrossRef] [PubMed]
Pollack, G. L. , 1969, “ Kapitza Resistance,” Rev. Mod. Phys., 41(1), p. 48. [CrossRef]
Zhang, Y. , Sun, M. , Liu, Z. , Piedra, D. , Lee, H.-S. , Gao, F. , Fujishima, T. , and Palacios, T. , 2013, “ Electrothermal Simulation and Thermal Performance Study of GaN Vertical and Lateral Power Transistors,” IEEE Trans. Electron Devices, 60(7), pp. 2224–2230. [CrossRef]
Parr, R. G. , and Yang, W. , 1989, Density-Functional Theory of Atoms and Molecules, Oxford University Press, New York.
Kresse, G. , and Furthmüller, J. , 1996, “ Software VASP, Vienna (1999),” Phys. Rev. B, 54(11), p. 169.
Lindsay, L. , Broido, D. , and Reinecke, T. , 2012, “ Thermal Conductivity and Large Isotope Effect in GaN From First Principles,” Phys. Rev. Lett., 109(9), p. 095901. [CrossRef] [PubMed]
Ruf, T. , Serrano, J. , Cardona, M. , Pavone, P. , Pabst, M. , Krisch, M. , D'astuto, M. , Suski, T. , Grzegory, I. , and Leszczynski, M. , 2001, “ Phonon Dispersion Curves in Wurtzite-Structure GaN Determined by Inelastic X-Ray Scattering,” Phys. Rev. Lett., 86(5), p. 906. [CrossRef] [PubMed]
Li, W. , Carrete, J. , Katcho, N. A. , and Mingo, N. , 2014, “ ShengBTE: A Solver of the Boltzmann Transport Equation for Phonons,” Comput. Phys. Commun., 185(6), pp. 1747–1758. [CrossRef]
Wu, X. , Lee, J. , Varshney, V. , Wohlwend, J. L. , Roy, A. K. , and Luo, T. , 2016, “ Thermal Conductivity of Wurtzite Zinc-Oxide From First-Principles Lattice Dynamics—A Comparative Study With Gallium Nitride,” Sci. Rep., 6(1), p. 22504. [CrossRef] [PubMed]


Grahic Jump Location
Fig. 4

Lattice temperature variation in the domain along the x-direction at the middle of the buffer (dotted line in the inset) aty = 1.9 × 10−6 μm. The Fourier model significantly under-predicts the hot-spot temperature.

Grahic Jump Location
Fig. 3

(a) Phonon dispersion of GaN from DFT calculations. The red dots shows the comparison with experiments from Ref. [32], and the blue lines show DFT results; (b) phonon relaxation time (τ) in picoseconds as a function of frequency at room temperature; (c) thermal conductivity of GaN as a function of temperature; (d) phonon mean free path of GaN as a function of frequency at room temperature (see color figure online).

Grahic Jump Location
Fig. 2

Flowchart showing the algorithm used for the adaptive BTE model

Grahic Jump Location
Fig. 1

(Left) block diagram of a GaN device with silicon substrate and AlGaN dielectric layer. (Right) schematic of the GaN buffer of size 50 μm × 2 μm divided into BTE and Fourier domain; the thick red line indicates the channel region where the energy is added. The exchange of temperature and heat flux information at the common interface is demonstrated in the bottom panel (see color figure online).

Grahic Jump Location
Fig. 5

Lattice temperature variation in the domain along the y-direction at the middle of the buffer (x = 1.0 × 10−6 μm). Inset shows the schematic of the domain with the GaN–Si interface along with boundary conditions. The thick red box indicates the region over which energy is added. The right-hand side panel shows the nonlinear temperature variation near the source on GaN side (see color figure online).

Grahic Jump Location
Fig. 6

Error in heat flux with respect to full BTE as a function of Kn cutoff

Grahic Jump Location
Fig. 7

(a) Spatial distribution of lattice temperature across the GaN Buffer layer using full BTE. Temperature distribution in the region highlighted in red from, (b) full BTE, (c) adaptive BTE, and (d) Fourier (the legend for temperature is same for the three cases (b), (c), and (d)) (see color figure online).

Grahic Jump Location
Fig. 8

Average time per iteration in seconds and error in maximum temperature as a function of Knudsen number cutoff for adaptive model (BTE is considered in only middle region of the domain, i.e., Fr + Adap. BTE + Fr model).

Grahic Jump Location
Fig. 9

Average time per iteration in seconds for the four methods considered




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In