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.

Heat transfer across fluid–solid interfaces in nanoconfinement has received significant attention due to its relevance in nanoscale systems. In this study, we investigate the Kapitza resistance at the water–graphene interface with the help of classical molecular dynamics simulation techniques in conjunction with our recently proposed equilibrium molecular dynamics (EMD) method [S. Alosious et al., J. Chem. Phys. 151, 194502 (2019)]. The size effect of the Kapitza resistance on different factors such as the number of graphene layers, the cross-sectional area, and the width of the water block was studied. The Kapitza resistance decreases slightly with an increase in the number of layers, while the influence of the cross-sectional area and the width of the water block is negligible. The variation in the Kapitza resistance as a function of the number of graphene layers is attributed to the large phonon mean free path along the graphene cross-plane. An optimum water–graphene system, which is independent of size effects, was selected, and the same was used to determine the Kapitza resistance using the predicted EMD method. The values obtained from both the EMD and the non-equilibrium molecular dynamics (NEMD) methods were compared for different potentials and water models, and the results are shown to be in good agreement. Our method allows us to compute the Kapitza resistance using EMD simulations, which obviates the need to create a large temperature gradient required for the NEMD method.

The Kapitza resistance (R k ) at the water−carbon nanotube (CNT) interface, with water on the inside of the nanotube, was investigated using molecular dynamics simulations. We propose a new equilibrium molecular dynamics (EMD) method, also valid in the weak flow regime, to determine the Kapitza resistance in a cylindrical nanoconfinement system where nonequilibrium molecular dynamics (NEMD) methods are not suitable. The proposed method is independent of the correlation time compared to Green−Kubo-based methods, which only work in short correlation time intervals. R k between the CNT and the confined water strongly depends on the diameter of the nanotube and is found to decrease with an increase in the CNT diameter, the opposite to what is reported in the literature when water is on the outside of the nanotube. R k is furthermore found to converge to the planar graphene surface value as the number of water molecules per unit surface area approaches the value in the graphene surface and a higher overlap of the vibrational spectrum. A slight increase in R k with the addition of the number of CNT walls was observed, whereas the chirality and flow do not have any impact.