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

TL;DR

This paper derives extended continuum heat transfer equations from the phonon Boltzmann equation, revealing the limits of Fourier's law near boundaries and providing boundary conditions that include kinetic effects for small but finite mean free paths.

Contribution

It introduces an asymptotic solution method that extends Fourier's law validity into the kinetic regime, incorporating boundary layer effects and higher-order boundary conditions.

Findings

Fourier's law holds up to second order in Knudsen number in the bulk.

Boundary effects require kinetic boundary layer solutions for accurate modeling.

The method predicts interface thermal resistance and temperature jumps.

Abstract

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 Fourier's 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.