Journal of Chemical Physics, Vol.109, No.1, 196-207, 1998
Direct simulation of hydrodynamic relaxation in microchannels
Simulations were performed on a fluid confined between two parallel walls. The fluid is modeled by a Lennard-Jones potential and the walls by a simple cubic lattice of harmonically bonded sites. A Lennard-Jones potential is also used to model the interactions between the wall and the fluid. The simulation consisted of over 30 000 sites arranged to form a liquid film approximately 35 Lennard-Jones diameters in thickness. This is large enough to begin approximating the range where classical hydrodynamics is expected to be applicable. Both equilibrium simulations and simulations of velocity transients were performed on the system. Two values of the wall-fluid interaction strength were examined, which appear to correspond to a wetting and nonwetting surface. Results from equilibrium simulations show that both the density and the stress tensor relax to their bulk values within a short distance of the wall. Furthermore, examination of the relaxation of spontaneous momentum fluctuations indicates that there is little change in the value of transport coefficients near the boundary compared to the bulk fluid. Nonequilibrium simulations on the decay of a parabolic velocity profile, however, suggest that the decay of the profile is faster than would be predicted from classical hydrodynamics and that the type of boundary conditions that should be used in a hydrodynamic analysis may depend on the details of the wall-fluid interaction. (C) 1993 American Institute of Physics. [S0021-9606(98)50725-8].