We study the dynamics of two air bubbles driven by the motion of a suspending viscous fluid in a Hele-Shaw channel with a small elevation along its centreline via physical experiment and numerical simulation of a depth-averaged model. For a single-bubble system we establish that, in general, the bubble propagation speed monotonically increases with bubble volume so that two bubbles of different sizes, in the absence of any hydrodynamic interactions, will either coalesce or separate in a finite time. However, our experiments indicate that the bubbles interact and that an unstable two-bubble state is responsible for the eventual dynamical outcome: coalescence or separation. These results motivate us to develop an edge-tracking routine and to calculate these weakly unstable two-bubble steady states from the governing equations. The steady states consist of pairs of ‘aligned’ bubbles that appear on the same side of the centreline with the larger bubble leading. We also discover, through time-dependent simulations and physical experiment, another class of two-bubble states that, surprisingly, are stable. In contrast to the ‘aligned’ steady states, these bubbles appear on either side of the centreline and are ‘offset’ from each other. We calculate the bifurcation structures of both classes of steady states as the flow rate and bubble volume ratio are varied. We find that they exhibit intriguing similarities to the single-bubble bifurcation structure, which has implications for the existence of $n$ -bubble steady states.

Unpredictable dynamics arising from a sensitivity to initial conditions is commonly associated with chaos. We demonstrate how similar unpredictability manifests in a nonlinear system that possesses a large number of long-term outcomes, namely the propagation of an air bubble within a viscous fluid-filled channel. The system under investigation supports various stable states of single-bubble propagation. In addition, bubbles can readily break up during their propagation. Upon subjecting steadily-propagating bubbles to finite-amplitude perturbations in the form of localised channel constrictions, we identify localised regions of the driving flow rate for which the resulting evolutions are unpredictable. Visibly-indistinguishable bubbles are observed to evolve towards a multitude of long-term outcomes, including each of the stable states available to the initial bubble and various states of permanently-changed bubble topology. By combining high-precision experimental results with simulations of a depth-averaged lubrication model of the system, we determine that this behaviour is driven by a sensitive dependence on initial conditions within the vicinity of an unstable periodic orbit.