Shock wave lithotripsy and high intensity focused ultrasound produce a variety of cavitation bubble interactions in confined spaces. Densely packed bubble clusters are generated that exhibit strong collective interactions. In lithotripsy, the bubble clusters themselves interact with the surface of a kidney stone. Comminution of the stone generates an accumulation of small particles that are interspersed with, and dynamically coupled to, the cavitation bubbles. Within tissues surrounding the stone, cavitation activity can occur in urinary ducts and blood vessels, constrained by channel walls. We are developing a systematic approach to modeling coupled bubble dynamics in clusters, amidst particles, near surfaces, and in tubes. The model is based on Lagrangian and Hamiltonian formulations of the dynamical equations for interacting bubbles. The initial formulations model the acoustically driven pulsation and translation of clustered spherical bubbles in a free field [Ilinskii et al., J. Acoust. Soc. Am. 121, 786–795 (2007)]. This model has been augmented to include interaction with spherical elastic particles. Interaction with rigid surfaces, including confinement by tubes with triangular, square, and hexagonal cross sections, has been modeled via the method of images. The relevant model equations will be presented and discussed, and representative simulations will be shown.