We present a new framework for solving volume integral equation (VIE) formulations of Poisson, Stokes and Helmholtz equations in complex geometries. Our method is based on a high-order discretization of the data on adaptive octrees and uses a kernel independent fast multipole method along with precomputed quadratures for evaluating volume integrals efficiently. We will present new performance optimizations and demonstrate scalability of our method on distributed memory platforms with thousands of compute nodes.