The dynamic interactions between gas bubbles, rigid particles and liquid can lead to profound nonlinearities in the aggregate behavior of a multiphase fluid. Predicting the nonlinear dynamics of the multiphase mixture hence requires understanding how the phases interact at the scale of individual interfaces, but these interactions are notoriously difficult to resolve in models. The goal of this paper is to develop and validate a computational method capable of capturing the complex flow interactions between gas bubbles and rigid particles immersed in a Newtonian liquid. We focus on multiphase systems that are dilute enough for the solid and gas components to move through and be moved by the ambient liquid. We use level sets with a topology-preserving advection scheme to track the gas interfaces. To include the motion of the rigid particles, we couple distributed Lagrange multipliers to an immersed-boundary method. The high viscosity contrast between the liquid and the gas requires both time splitting and approximate factorization to efficiently solve the governing equations consisting of the conservation of mass, momentum and energy. To resolve interactions between interfaces that vary drastically in size, we refine our mesh adaptively in the vicinity of the boundary. (C) 2017 Elsevier Ltd. All rights reserved.