Jean-Philippe M. Péraud

An alternative approach to efficient simulation of micro/nanoscale phonon transport

Starting from the recently proposed energy-based deviational formulation for solving the Boltzmann equation [J.-P. Peraud and N. G. Hadjiconstantinou, Phys. Rev. B 84, 205331 (2011)], which provides significant computational speedup compared to standard Monte Carlo methods for small deviations from equilibrium, we show that additional computational benefits are possible in the limit that the governing equation can be linearized. The proposed method exploits the observation that under linearized conditions (small temperature differences) the trajectories of individual deviational particles can be decoupled and thus simulated independently; this leads to a particularly simple and efficient algorithm for simulating steady and transient problems in arbitrary three-dimensional geometries, without introducing any additional approximation.

On efficient simulations of multiscale kinetic transport.

We discuss a new class of approaches for simulating multiscale kinetic problems, with particular emphasis on applications related to small-scale transport. These approaches are based on a decomposition of the kinetic description into an equilibrium part, which is described deterministically (analytically or numerically), and the remainder, which is described using a particle simulation method. We show that it is possible to derive evolution equations for the two parts from the governing kinetic equation, leading to a decomposition that is dynamically and automatically adaptive, and a multiscale method that seamlessly bridges the two descriptions without introducing any approximation. Our discussion pays particular attention to stochastic particle simulation methods that are typically used to simulate kinetic phenomena; in this context, these decomposition approaches can be thought of as control-variate variance-reduction formulations, with the nearby equilibrium serving as the control. Such formulations can provide substantial computational benefits in a broad spectrum of applications because a number of transport processes and phenomena of practical interest correspond to perturbations from nearby equilibrium distributions. In many cases, the computational cost reduction is sufficiently large to enable otherwise intractable simulations.

Monte Carlo study of non-diffusive relaxation of a transient thermal grating in thin membranes

The impact of boundary scattering on non-diffusive thermal relaxation of a transient grating in thin membranes is rigorously analyzed using the multidimensional phonon Boltzmann equation. The gray Boltzmann simulation results indicate that approximating models derived from previously reported one-dimensional relaxation model and Fuchs-Sondheimer model fail to describe the thermal relaxation of membranes with thickness comparable with phonon mean free path. Effective thermal conductivities from spectral Boltzmann simulations free of any fitting parameters are shown to agree reasonably well with experimental results. These findings are important for improving our fundamental understanding of non-diffusive thermal transport in membranes and other nanostructures.

Thermal transport in suspended silicon membranes measured by laser-induced transient gratings

Studying thermal transport at the nanoscale poses formidable experimental challenges due both to the physics of the measurement process and to the issues of accuracy and reproducibility. The laser-induced transient thermal grating (TTG) technique permits non-contact measurements on nanostructured samples without a need for metal heaters or any other extraneous structures, offering the advantage of inherently high absolute accuracy. We present a review of recent studies of thermal transport in nanoscale silicon membranes using the TTG technique. An overview of the methodology, including an analysis of measurements errors, is followed by a discussion of new findings obtained from measurements on both “solid” and nanopatterned membranes. The most important results have been a direct observation of non-diffusive phonon-mediated transport at room temperature and measurements of thickness-dependent thermal conductivity of suspended membranes across a wide thickness range, showing good agreement with first-princ...

Extending the range of validity of Fourier's law into the kinetic transport regime via asymptotic solution of the phonon Boltzmann transport equation

We derive the continuum equations and boundary conditions governing phonon-mediated heat transfer in the limit of small but finite mean free path from asymptotic solution of the linearized Boltzmann equation in the relaxation time approximation. Our approach uses the ratio of the mean free path to the characteristic system lengthscale, also known as the Knudsen number, as the expansion parameter to study the effects of boundaries on the breakdown of the Fourier descrition. We show that, in the bulk, the traditional heat conduction equation using Fouriers law as a constitutive relation is valid at least up to second order in the Knudsen number for steady problems and first order for time-dependent problems. However, this description does not hold within distances on the order of a few mean free paths from the boundary; this breakdown is a result of kinetic effects that are always present in the boundary vicinity and require solution of a Boltzmann boundary-layer problem to be determined. Matching the inner, boundary layer, solution to the outer, bulk, solution yields boundary conditions for the Fourier description as well as additive corrections in the form of universal kinetic boundary layers; both are found to be proportional to the bulk-solution gradients at the boundary and parametrized by the material model and the phonon-boundary interaction model (Boltzmann boundary condition). Our derivation shows that the traditional no-jump boundary condition for prescribed temperature boundaries and no-flux boundary condition for diffusely reflecting boundaries are appropriate only to zeroth order in the Knudsen number; at higher order, boundary conditions are of the jump type. We illustrate the utility of the asymptotic solution procedure by demonstrating that it can be used to predict the Kapitza resistance (and temperature jump) associated with an interface between two materials.

Geometry and wetting of capillary folding

Capillary forces are involved in a variety of natural phenomena, ranging from droplet breakup to the physics of clouds. The forces from surface tension can also be exploited in industrial applications provided the length scales involved are small enough. Recent experimental investigations showed how to take advantage of capillarity to fold planar structures into three-dimensional configurations by selectively melting polymeric hinges joining otherwise rigid shapes. In this paper we use theoretical calculations to quantify the role of geometry and fluid wetting on the final folded state. Considering folding in two and three dimensions, studying both hydrophilic and hydrophobic situations with possible contact-angle hysteresis, and addressing the shapes to be folded to be successively infinite, finite, curved, kinked, and elastic, we are able to derive an overview of the geometrical parameter space available for capillary folding.

Reconstruction of the phonon relaxation times using solutions of the Boltzmann transport equation

We present a method for reconstructing the phonon relaxation time distribution

On the Equations and Boundary Conditions Governing Phonon-Mediated Heat Transfer in the Small Mean Free Path Limit: An Asymptotic Solution of the Boltzmann Equation

{\ensuremath{\tau}}_{\ensuremath{\omega}}=\ensuremath{\tau}(\ensuremath{\omega})

Deviational Phonons and Thermal Transport at the Nanoscale

(including polarization) in a material from thermal spectroscopy data. The distinguishing feature of this approach is that it does not make use of the effective thermal conductivity concept and associated approximations. The reconstruction is posed as an optimization problem in which the relaxation times

MONTE CARLO METHODS FOR SOLVING THE BOLTZMANN TRANSPORT EQUATION

{\ensuremath{\tau}}_{\ensuremath{\omega}}=\ensuremath{\tau}(\ensuremath{\omega})