Prediction of Kapitza resistance at fluid-solid interfaces

Understanding the interfacial heat transfer and thermal resistance at an interface between two dissimilar materials is of great importance in the development of nanoscale systems. This paper introduces a new and reliable linear response method for calculating the interfacial thermal resistance or Kapitza resistance in fluid-solid interfaces with the use of equilibrium molecular dynamics (EMD) simulations. The theoretical predictions are validated against classical molecular dynamics (MD) simulations. MD simulations are carried out in a Lennard-Jones (L-J) system with fluid confined between two solid slabs. Different types of interfaces are tested by varying the fluid-solid interactions (wetting coefficient) at the interface. It is observed that the Kapitza length decreases monotonically with an increasing wetting coefficient as expected. The theory is further validated by simulating under different conditions such as channel width, density, and temperature. Our method allows us to directly determine the Kapitza length from EMD simulations by considering the temperature fluctuation and heat flux fluctuations at the interface. The predicted Kapitza length shows an excellent agreement with the results obtained from both EMD and non-equilibrium MD simulations.