AbstractA new method for the simulation of the translational and rotational motions of a system containing a sedimenting particle interacting with a neutrally buoyant particle has been developed. The method is based on coupling the quasi‐static Stokes equations for the fluid with the rigid body equations of motion for the particles. The Stokes equations are solved at each time step with the boundary element method. The stresses are then integrated over the surface of each particle to determine the resultant forces and moments. These forces and moments are inserted into the rigid body equations of motion to determine the translational and rotational motions of the particles. Unlike many other simulation techniques, no restrictions are placed on the shape of the particles. Superparametric boundary elements are employed to achieve accurate geometric representations of the particles. The simulation method is able to predict the local fluid velocity, resolve the forces and moments exerted on the particles, and track the particle trajectories and orientation