Thermal conductivity predictions of herringbone graphite nanofibers using molecular dynamics simulations

Non-equilibrium molecular dynamics (NEMD) simulations are used to investigate the thermal conductivity of herringbone graphite nanofibers (GNFs) at room temperature by breaking down the axial and transverse conductivity values into intralayer and interlayer components. The optimized Tersoff potential is used to account for intralayer carbon-carbon interactions while the Lennard-Jones potential is used to model the interlayer carbon-carbon interactions. The intralayer thermal conductivity of the graphene layers near room temperature is calculated for different crease angles and number of layers using NEMD with a constant applied heat flux. The edge effect on a layer's thermal conductivity is investigated by computing the thermal conductivity values in both zigzag and armchair directions of the heat flow. The interlayer thermal conductivity is also predicted by imposing hot and cold Nosé-Hoover thermostats on two layers. The limiting case of a 90° crease angle is used to compare the results with those of single-layer graphene and few-layer graphene. The axial and transverse thermal conductivities are then calculated using standard trigonometric conversions of the calculated intralayer and interlayer thermal conductivities, along with calculations of few-layer graphene without a crease. The results show a large influence of the crease angle on the intralayer thermal conductivity, and the saturation of thermal conductivity occurs when number of layers is more than three. The axial thermal conductivity, transverse thermal conductivity in the crease direction, and transverse thermal conductivity normal to the crease for the case of a five-layer herringbone GNF with a 45° crease angle are calculated to be 27 W∕m K, 263 W∕m K, and 1500 W∕m K, respectively, where the axial thermal conductivity is in good agreement with experimental measurements.