A unified Monte Carlo treatment of the transport of electromagnetic energy, electrons, and phonons in absorbing and scattering media

The scalar Boltzmann transport equation (BTE) is often applicable to radiative energy transfer, electron-beam propagation, as well as thermal conduction by electrons and phonons provided that the characteristic length of the system is much larger than the wavelength of energy carriers and that certain interference phenomena and the polarization nature of carriers are ignored. It is generally difficult to solve the BTE analytically unless a series of assumptions are introduced for the particle distribution function and scattering terms. Yet, the BTE can be solved using statistical approaches such as Monte Carlo (MC) methods without simplifying the underlying physics significantly. Derivations of the MC methods are relatively straightforward and their implementation can be achieved with little effort; they are also quite powerful in accounting for complicated physical situations and geometries. MC simulations in radiative transfer, electron-beam propagation, and thermal conduction by electrons and phonons have similar simulation procedures; however, there are important differences in implementing the algorithms and scattering properties between these simulations. The objective of this review article is to present these simulation procedures in detail and to show that it is possible to adapt an existing MC computer code, for instance, in radiative transfer, to account for physics in electron-beam transport or phonon (or electronic thermal) conduction by sorting out the differences and implementing the correct corresponding steps. Several simulation results are presented and some of the difficulties associated with different applications are explained.