Molecular dynamics simulations of thermal resistance at the liquid-solid interface

Heat conduction between parallel plates separated by a thin layer of liquid Argon is investigated using three-dimensional molecular dynamics (MD) simulations employing 6-12 Lennard-Jones potential interactions. Channel walls are maintained at specific temperatures using a recently developed interactive thermal wall model. Heat flux and temperature distribution in nanochannels are calculated for channel heights varying from 12.96 to 3.24 nm. Fourier law of heat conduction is verified for the smallest channel, while the thermal conductivity obtained from Fourier law is verified using the predictions of Green-Kubo theory. Temperature jumps at the liquid/solid interface, corresponding to the well known Kapitza resistance, are observed. Using systematic studies thermal resistance length at the interface is characterized as a function of the surface wettability, thermal oscillation frequency, wall temperature, thermal gradient, and channel height. An empirical model for the thermal resistance length, which could be used as the jump coefficient of a Navier boundary condition, is developed. Temperature distribution in nanochannels is predicted using analytical solution of continuum heat conduction equation subjected to the new temperature jump condition. Analytical predictions are verified using MD simulations.