A central part of integral equation methods are the quadrature methods used to evaluate boundary integrals that are singular or have nearby singularities. A relatively new such method is quadrature by expansion (QBX), which evaluates the solution using local expansions formed at points inside the domain. I will in this talk discuss how QBX can be used for Stokes flow simulations in both two and three dimensions. In particular, I will discuss techniques for accelerating the method such that it can be combined with existing fast methods, and also how parameter selection can be simplified through accurate estimation of quadrature errors.