Solution of the Boltzmann transport equation in two real-space dimensions using a spherical harmonic expansion in momentum space

Author(s):  
K. Rahmat ◽  
J. White ◽  
D.A. Antoniadis
VLSI Design ◽  
1998 ◽  
Vol 8 (1-4) ◽  
pp. 147-151 ◽  
Author(s):  
C.-H. Chang ◽  
C.-K. Lin ◽  
N. Goldsman ◽  
I. D. Mayergoyz

We perform a rigorous comparison between the Spherical Harmonic (SH) and Monte Carlo (MC) methods of solving the Boltzmann Transport Equation (BTE), on a 0.05 μm base BJT. We find the SH and the MC methods give very similar results for the energy distribution function, using an analytical band-structure, at all points within the tested devices. However, the SH method can be as much as seven thousand times faster than the MC approach for solving an identical problem. We explain the agreement by asymptotic analysis of the system of equations generated by the SH expansion of the BTE.


2017 ◽  
Vol 139 (10) ◽  
Author(s):  
Ajit K. Vallabhaneni ◽  
Liang Chen ◽  
Man P. Gupta ◽  
Satish Kumar

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.


Sign in / Sign up

Export Citation Format

Share Document