Yusuke Masao*, Masaya Okano, Mitsuhiro Matsumoto
The Journal of Mechanical Science and Technology, vol. 25, no. 1, pp.21-26, 2011

Abstract : In order to solve the Boltzmann transport equation (BTE) of phonons for investigating heat conduction in non-metallic solids, we propose employing a DSMC (direct simulation Monte Carlo) scheme to simulate the dynamics of phonons analogous to rarefied gas. In contrast to treating the BTE with conventional linear approximation, this scheme requires no relaxation times as input parameters. We can directly investigate couplings among phonons with different modes, although we have to assume an appropriate scattering model for phonon-phonon interactions. In this paper, we describe the DSMC scheme for phonon dynamics and present some results with our prototype codes for a simple solid model. In the first case of single-branch four-phonon processes, we carried out simulations of nonequilibrium thin films with a temperature gradient. We found that the temperature jump at the boundaries can be successfully achieved. For the second case of three phonon processes, we developed a simulation code that takes into consideration the different acoustic branches to evaluate the mode-dependent relaxation time and mean free path. This type of DSMC scheme for phonons enables us to include other relevant factors, such as optical branches and phonon-electron interactions.

Keyword : Boltzmann transport equation; DSMC; Heat conduction; Phonon dynamics; Solid thin film