A general boundary integral formulation for compressible Stokes flows is theoretically described within the framework of hydrodynamic potentials. The integral equation is implemented numerically to the study of drop expansion in compressible viscous flows. Marker point positions on the drop interface are involved by using the boundary integral method for calculation of fluid velocity. Surface discretization is adaptive to the instantaneous drops shapes. The interplay between viscous and surface tension and its influence on the evolving emulsion microstructure during its expansion is fundamental to the science and technology of foam processing. In this article the method is applied for 3D simulations of emulsion densification that involves an uniform expansion of a viscous fluid containing spherical drops on a body centered cubic lattice (BCC).