The motion of spheres (bubbles) in an incompressible liquid undergoing a small‐amplitude oscillatory motion is calculated by a multipole expansion method. In the limit of small viscosity, the Stokes layer is confined to the vicinity of the surface of the bubbles, which therefore interact approximately only through the pressure field. The motion of the spheres is parametrized in terms of added mass, Basset, and drag forces; the coefficients of which are obtained from the simulation. To obtain results useful for the study of pressure wave propagation in bubbly liquids, several bubble configurations are studied for different (finite) volume fractions and the results then averaged. The effects of surface tension and bubble density are also considered. [Work supported by DOE.]