To model the imbibition of liquids into porous solids, use is often made of the Lucas-Washburn equation, which relates the distance of penetration of a liquid at a given time to the pore radius, the viscosity and surface tension of the liquid, and the effective contact angle between the liquid and the solid. In this paper, we extend previous large-scale molecular dynamics simulations to show how this tool can be used to study the details of liquid imbibition, including the impact of the contact angle on the dynamics of penetration and the evolution of the internal flow field. In particular, we show that the asymptotic behavior of the contact angle versus time for a completely wetting liquid is given by similar tot(-1/4). (C) 2003 Elsevier Inc. All rights reserved.