Boundary integral methods for sound propagation with subsonic potential mean flows