---

In the present paper, Molecular Dynamics Simulation (MDS) is performed for Poiseuille flow of liquid Argon in a nanochannel by embedding the fluid particles in an external force with different potential functions. Three types of Lennard-Jones (LJ) potentials are used as interatomistic or molecular models for evaluations of interactions and density, velocity profiles across the channel are investigated. The interatomic potentials are LJ 12-6 potential, LJ 9-6 potential and LJ-Smooth potential. Density and velocity profiles across the channel are investigated. Obtained results show that hydrodynamic characteristics and behavior of flow depends on the type of interaction potential. It is shown that the LJ 9-6 predictions for velocity and temperature are larger than those of LJ12-6 and LJ-Smooth potentials. Also, applying LJ 9-6 results in further calculations time. The results show the effect of interaction force model on the understanding and analyzing of nanoscale flows.

---

Fluid flow through channels and ducts in nano scales is an important issue which needs numerical simulations for better analysis of fluid behavior because of the limitations of experimental methods. Hence, in the present study Molecular Dynamics simulation is used as a precise method for molecular scale problems to investigate fluid behavior. This method which is based on Newton’s second law, is applied to investigate liquid Argon flow in steady Couette flows through smooth and rough nanochannels. Using LAMMPS software, were performed simulation. In the present study, the fluid velocity and fluid slip in steady Couette flows were obtained to investigate various effects including: wall velocity, channel height, wall density, fluid-wall interaction, and surface roughness with different shapes such as rectangular and triangular in different dimensions. Based on the results, an increase in wall velocity increases the fluid slip velocity. For velocity constant values, an increase of channel height will decrease the fluid slip velocity. In steady Couette flow, decrease of wall density will result in decrease of fluid slip velocity. Reducing the energy parameter between fluid and wall will increase the fluid slip velocity and on the other hand, decreasing the fluid-wall length parameter will decrease the fluid slip velocity. The rectangular and triangular roughness at the bottom wall reduces the fluid slip velocity, and an increase of roughness height will further decrease the fluid slip velocity.