Drag on a nanotube in uniform liquid argon flow

In this work, nonequilibrium molecular dynamics (MD) simulations were performed to investigate uniform liquid argon flow past a carbon nanotube. In the simulation, nanotubes were modeled as rigid cylinders of carbon atoms. Both argon-argon and argon-carbon interactions were calculated based on Lennard-Jones potential. Simulated drag coefficients were compared with (i) published empirical equation which was based on experiments conducted with macroscale cylinders and (ii) finite element (FE) analyses based on Navier-Stokes equation for flow past a circular cylinder using the same dimensionless parameters used in MD simulations. Results show that classical continuum mechanics cannot be used to calculate drag on a nanotube. In slow flows, the drag coefficients on a single-walled nanotube calculated from MD simulations were larger than those from the empirical equation or FE analysis. The difference increased as the flow velocity decreased. For higher velocity flows, slippage on the surface of the nanotube was identified which resulted in lower drag coefficient from MD simulation. This explains why the drag coefficient from MD dropped faster than those from the empirical equation or FE simulation as the flow velocity increased. It was also found that the drag forces are almost equal for single- and double-walled nanotubes with the same outer diameter, implying that inner tubes do not interact with fluid molecules.