A Novel Approach for Lattice Boltzmann Modeling of Energy Transport

Author(s):  
Dadong Wang ◽  
Yanbao Ma

It is well known that Fourier law breaks down for the prediction of heat conduction in nano-scale, where the length scale is comparable to the mean free path of energy carriers. Over the past decade, Boltzmann transport equation (BTE) has been used to predict thermal transport in dielectrics and semiconductors at micro-scale and nano-scale. In this work, a new modified gray model is obtained from BTE. The implicit lattice Boltzmann method (LBM) is developed to simulate the thermal transport process. Based on the new model, we can derive Guyer-Krumhansl equation. Transient heat conduction through a thin nano-film and hotspot self-heating in sub-micron transistors are examined. The numerical results are compared with those provided by Fourier, Cattaneo, and Guyer-Krumhansl equation.

Author(s):  
Ankur Chattopadhyay ◽  
Arvind Pattamatta

Heat transport at nanoscales departs substantially from the well established classical laws governing the physical processes at continuum level. The Fourier Law of heat conduction cannot be applied at sub-continuum level due to its inability in modeling non-equilibrium energy transport. Therefore one must resort to a rigorous solution to the Boltzmann Transport Equation (BTE) in the realm of nanoscale transport regime. Some recent studies show that a relatively inexpensive and accurate way to predict the behavior of sub continuum energy transport in solids is via the discrete representation of the BTE referred to as the Lattice Boltzmann method (LBM). Although quite a few numerical simulations involving LBM have been exercised in the literature, there has been no clear demonstration of the accuracy of LBM over BTE; also there exists an ambiguity over employing the right lattice configurations describing phonon transport. In the present study, the Lattice Boltzmann Method has been implemented to study phonon transport in miniaturized devices. The initial part of the study focuses upon a detailed comparison of the LBM model with that of BTE for one dimensional heat transfer involving multiple length and time scales. The second objective of the present investigation is to evaluate different lattice structures such as D1Q2, D1Q3, D2Q5, D2Q8, D2Q9 etc. for 1-D and 2-D heat conduction. In order to reduce the modeling complexity, gray model assumption based on Debye approximation is adopted throughout the analysis. Results unveil that the accuracy of solution increases as the number of lattice directions taken into account are incremented from D2Q5 to D2Q9. A substantial increase in solution time with finer directional resolutions necessitates an optimum lattice. A novel lattice dimension ‘Mod D2Q5’ has been suggested and its performance is also compared with its compatriots. It is also demonstrated that the inclusion of the center point within a particular lattice structure can play a significant role in the prediction of thermal conductivity in the continuum level. However, as the size of the device comes down to allow high Knudsen numbers, in the limiting case of ballistic phonon transport, the choice of lattice seems to have negligible effect on thermal conductivity.


Author(s):  
Rodrigo A. Escobar ◽  
Cristina H. Amon ◽  
Amador M. Guzma´n

Numerical simulations of time-dependent thermal energy transport in semiconductor thin films are performed using the Lattice Boltzmann Method applied to phonon transport. The discrete Lattice Boltzmann Method is derived from the continuous Boltzmann transport equation assuming nonlinear, frequency-dependent phonon dispersion for acoustic and optical phonons. Results indicate that the heat conduction in silicon thin films displays a transition from diffusive to ballistic energy transport as the characteristic length of the system becomes comparable to the phonon mean free path, and that the thermal energy transport process is characterized by the propagation of multiple, superimposed phonon waves. The methodology is used to characterize the time-dependent temperature profiles inside films of decreasing thickness. Thickness-dependent thermal conductivity values are computed based on steady-state temperature distributions obtained from the numerical models. It is found that reducing feature size into the subcontinuum regime decreases the thermal conductivity when compared to bulk values, at a higher rate than what was displayed by the Debye-based gray Lattice Boltzmann Method.


Author(s):  
Sartaj S. Ghai ◽  
Rodrigo A. Escobar ◽  
Myung S. Jhon ◽  
Cristina H. Amon

The lattice Boltzmann method (LBM) is used to examine multi-length scale, confined heat conduction problems in one dimension for which sub-continuum effects are important. This paper describes the implementation of the method and its application to electronic devices. A silicon-on-insulator device with internal heat generation is used as a case study to illustrate the advantages of the LBM. We compare our results with various hierarchical equations of heat transfer such as Fourier, Cattaneo, and Boltzmann transport equations, as well as with experimental and numerical data from the literature. Our results provide excellent agreement with other methodologies, at a far less computational effort.


2008 ◽  
Vol 130 (9) ◽  
Author(s):  
Rodrigo A. Escobar ◽  
Cristina H. Amon

Numerical simulations of time-dependent thermal energy transport in semiconductor thin films are performed using the lattice Boltzmann method applied to phonon transport. The discrete lattice Boltzmann Method is derived from the continuous Boltzmann transport equation assuming nonlinear, frequency-dependent phonon dispersion for acoustic and optical phonons. Results indicate that the heat conduction in silicon thin films displays a transition from diffusive to ballistic energy transport as the characteristic length of the system becomes comparable to the phonon mean free path and that the thermal energy transport process is characterized by the propagation of multiple superimposed phonon waves. The methodology is used to characterize the time-dependent temperature profiles inside films of decreasing thickness. Thickness-dependent thermal conductivity values are computed based on steady-state temperature distributions obtained from the numerical models. It is found that reducing feature size into the subcontinuum regime decreases thermal conductivity when compared to bulk values, at a higher rate than what was displayed by the Debye-based gray lattice Boltzmann method.


2016 ◽  
Vol 138 (4) ◽  
Author(s):  
Cheng Chen ◽  
James Geer ◽  
Bahgat Sammakia

In this paper, a lattice Boltzmann method (LBM)-based model is developed to simulate the subcontinuum behavior of multidimensional heat conduction in solids. Based on a previous study (Chen et al., 2014, “Sub-Continuum Thermal Modeling Using Diffusion in the Lattice Boltzmann Transport Equation,” Int. J. Heat Mass Transfer, 79, pp. 666–675), phonon energy transport is separated to a ballistic part and a diffusive part, with phonon equilibrium assumed at boundaries. Steady-state temperature/total energy density solutions from continuum scales to ballistic scales are considered. A refined LBM-based numerical approach is applied to a two-dimensional simplified transistor model proposed by (Sinha et al. 2006, “Non-Equilibrium Phonon Distributions in Sub-100 nm Silicon Transistors,” ASME J. Heat Transfer, 128(7), pp. 638–647), and the results are compared with the Fourier-based heat conduction model. The three-dimensional (3D) LBM model is also developed and verified at both the ballistic and continuous limits. The impact of film thickness on the cross-plane and in-plane thermal conductivities is analyzed, and a new model of the supplementary diffusion term is proposed. Predictions based on the finalized model are compared with the existing in-plane thermal conductivity measurements and cross-plane thermal conductivity molecular dynamics (MD) results.


2020 ◽  
Vol 24 (6 Part A) ◽  
pp. 3749-3756
Author(s):  
Ya Han ◽  
Shuai Li ◽  
Hai-Dong Liu ◽  
Weipeng Cui

In order to deeply investigate the gas heat conduction of nanoporous aerogel, a model of gas heat conduction was established based on microstructure of aerogel. Lattice Boltzmann method was used to simulate the temperature distribution and gas thermal conductivity at different size, and the size effects of gas heat conduction have had been obtained under micro-scale conditions. It can be concluded that the temperature jump on the boundary was not obvious and the thermal conductivity remained basically constant when the value of Knudsen number was less than 0.01; as the value of Knudsen number increased from 0.01 to 0.1, there was a clear temperature jump on the boundary and the thermal conductivity tended to decrease and the effect of boundary scattering increased drastically, as the value of Knudsen number was more than 0.1, the temperature jump increased significantly on the boundary, furtherly, the thermal conductivity decreased dramatically, and the size effects were significantly.


2016 ◽  
Vol 64 (3) ◽  
pp. 599-606 ◽  
Author(s):  
A. Piasecka Belkhayat ◽  
A. Korczak

Abstract In the paper a description of heat transfer in one-dimensional crystalline solids is presented. The lattice Boltzmann method based on Boltzmann transport equation is used to simulate the nanoscale heat transport in thin metal films. The coupled lattice Boltzmann equations for electrons and phonons are applied to analyze the heating process of thin metal films via laser pulse. Such approach in which the parameters appearing in the problem analyzed are treated as constant values is widely used, but in the paper the interval values of relaxation times and electron-phonon coupling factor are taken into account. The problem formulated has been solved by means of the interval lattice Boltzmann method using the rules of directed interval arithmetic. In the final part of the paper the results of numerical computations are shown.


Sign in / Sign up

Export Citation Format

Share Document