We analytically studied the flow field in multiple fluid sphere systems in low Knudsen number regime. The expanded zero vorticity cell model, based on Kuwabara's theory (1959), and the zero shear stress model of Happel (1958) were applied to bubbles and liquid spheres; solid spheres with the effects of gas slippage at the collecting surface were also considered. The analytical solution obtained for the finite Knudsen number region in a multiple fluid sphere was shown to converge with that of the existing solutions for the flow field around fluid sphere systems under proper boundary conditions. The flow field inside the collector shows increases in the velocity of flow with an increasing Knudsen number for cases with a viscosity ratio of less than 1, while decreases in the velocity are observed with a viscosity ratio greater than 1. Based on the resolved flow field, the drag force and terminal velocity around the collecting fluid spheres were obtained and compared with existing theory. Subsequently, this study evaluated the most general solution for the flow field around multiple fluid sphere systems, including the low Knudsen number region.