The interaction between stable colloidal particles arriving at a pore entrance was studied using a numerical method for the case where the particle size is smaller than but of the same order as the pore size. The numerical method was adapted from a front-tracking technique developed for studying incompressible, multifluid flow by S. O. Unverdi and G. Tryggvason (J. Comp. Phys. 100, 25, 1992). The method is based on the finite difference solution of Navier-Stokes equation on a stationary, structured, Cartesian grid and the explicit representation of the particle-liquid interface using an unstructured grid that moves through the stationary grid. The simulations are in two dimensions, considering both deformable and nondeformable particles, and include interparticle colloidal interactions. The interparticle and particle-pore hydrodynamic interactions, which are very difficult to determine using existing analytical and semi-numerical, semi-analytical techniques in microhydrodynamics, are naturally accounted for in our numerical method and need not be explicity determined. Two- and three-particle motion toward a pore has been considered in our simulations. The simulations demonstrate how the competition between hydrodynamic forces and colloidal forces acting on particles dictate their flow behavior near the pore entrance. The predicted dependence of the particle flow behavior on the flow velocity and the ratio of pore size to particle size are qualitatively consistent with the experimental observations of V. Ramachandran and H. S. Fogler (J. Fluid Mech. 385, 129, 1999).