A simplified mathematical model and numerical simulations of the governing Navier-Stokes equations are used to predict the shape evolution, rupture distance, and liquid distribution of stretching pendular liquid bridges between two equal-sized spherical solid particles. In the simplified model, the bridge shape is approximated with a parabola, and it is assumed that the surface tension effects dominate the viscous, inertial, and gravitational effects. For the numerical simulations, a commercial Computational Fluid Dynamics (CFD) software package-FLUENT-isused. The rupture distance predictions obtained with both models are compared with experimental data and a reasonable agreement is found. The results of the numerical investigations show that for simulations with negligible viscous, inertial, and gravitational effects, the rupture distance approaches an asymptotic value, which is close to the value predicted by the simplified model. The bridge profiles predicted using the simplified model and the numerical simulation are compared. It is found that a second-order polynomial appropriately represents the stable bridge shape for particles with identical contact angles; however, for liquid bridges between particles with different contact angles, the numerical simulations of the governing Navier-Stokes equations should be used. (C) 2010 Elsevier Ltd. All rights reserved.