Prediction of fluid slip in cylindrical nanopores using equilibrium molecular simulations
Alan Sam (Indian Institute of Technology Madras)
Remco Hartkamp (TU Delft - Process and Energy)
Sridhar Kumar Kannam (Swinburne University of Technology, Royal Melbourne Institute of Technology University)
Sarith P. Sathian (Indian Institute of Technology Madras)
More Info
expand_more
Other than for strictly personal use, it is not permitted to download, forward or distribute the text or part of it, without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license such as Creative Commons.
Abstract
We introduce an analytical method to predict the slip length (L s) in cylindrical nanopores using equilibrium molecular dynamics (EMD) simulations, following the approach proposed by Sokhan and Quirke for planar channels [39]. Using this approach, we determined the slip length of water in carbon nanotubes (CNTs) of various diameters. The slip length predicted from our method shows excellent agreement with the results obtained from nonequilibrium molecular dynamics (NEMD) simulations. The data show a monotonically decreasing slip length with an increasing nanotube diameter. The proposed EMD method can be used to precisely estimate slip length in high slip cylindrical systems, whereas, L s calculated from NEMD is highly sensitive to the velocity profile and may cause large statistical errors due to large velocity slip at the channel surface. We also demonstrated the validity of the EMD method in a BNNT-water system, where the slip length is very small compared to that in a CNT pore of similar diameter. The developed method enables us to calculate the interfacial friction coefficient directly from EMD simulations, while friction can be estimated using NEMD by performing simulations at various external driving forces, thereby increasing the overall computational time. The EMD analysis revealed a curvature dependence in the friction coefficient, which induces the slip length dependency on the tube diameter. Conversely, in flat graphene nanopores, both L s and friction coefficient show no strong dependency on the channel width.