Rising bubbles draw surface patterns: a numerical study

Small bubbles rising in a chain can self-organize into regular patterns upon reaching a liquid’s free surface. This phenomenon is investigated through direct numerical simulations. By varying the bubble release period, distinct branching patterns characterized by different numbers of arms are observed. These macroscopic regular configurations arise from localized non-contact repulsion and pair collisions between bubbles as they arrive at the free-surface emergence site. A theoretical model is proposed to quantitatively relate the number of branches to the bubble release period. The model also predicts probabilities of observing specific arm counts in reality. This study provides insights into broader nonlinear pattern formation and self-organization phenomena.

💡 Research Summary

This paper presents a comprehensive numerical investigation into the self-organization of periodically released bubble chains into regular geometric patterns upon reaching a liquid free surface. Using direct numerical simulations (DNS), the study unravels the underlying physical mechanisms and establishes a predictive model linking the bubble release period to the emergent pattern.

The research employs the open-source multiphase flow solver Aphros, utilizing the Volume-of-Fluid (VOF) method with PLIC interface reconstruction and a multi-marker technique. This setup prevents bubble coalescence and bursting, allowing for the accurate capture of non-contact interactions crucial to the phenomenon. The simulation domain is a partially liquid-filled box, with bubbles of fixed equivalent diameter released periodically from a central point at the bottom. Fluid properties are chosen to match prior experimental studies for validation.

The core finding is that the bubble release period (T) is the primary control parameter governing the pattern morphology. Systematic variation of T reveals four distinct regimes: a disordered mode for large T (~100ms), a stable two-armed mode for intermediate T (~50ms), and multi-armed modes (e.g., five-armed at ~42ms, spiral three-armed at ~36ms) for smaller T. The analysis focuses on the “divergence angle” (φ), the angular difference in the emission direction of consecutive bubbles from the central emergence site. Ordered patterns correspond to a constant or near-constant φ.

The microscopic physics driving pattern formation are dissected into four key processes: (i) direct pair collisions between bubbles, (ii) non-contact repulsion mediated by the intervening fluid flow, (iii) confinement by a central surface bump created by the rising bubble chain, and (iv) advection by the radial outward flow beneath the surface. For large T, direct collisions dominate, leading to disorder. As T decreases, bubbles arrive at the surface before preceding bubbles have moved far away. This switches the dominant interaction to non-contact repulsion, where a new bubble is repelled from the position of the most recent surface bubble, not necessarily in the exact opposite direction. This allows for a constant, non-180° divergence angle and the formation of multi-armed patterns. In these modes, bubbles are observed to orbit the central bump before settling into a radial arm.

Based on these insights, the authors propose a heuristic theoretical model. The model posits that the number of arms (n) and the divergence angle (φ) are related by nφ = 2πm, where m is a small integer. They further relate φ to the bubble release period T and the average radial velocity of bubbles on the surface, providing a quantitative bridge between the control parameter (T) and the macroscopic pattern (n). This model not only predicts the number of arms but can also be extended to estimate the probability of observing specific arm counts in real experiments, accounting for natural fluctuations.

In summary, this work successfully decodes a complex self-organization phenomenon by combining high-fidelity DNS with mechanistic analysis. It demonstrates how macroscopic spatial order can emerge from localized, repetitive interactions (collisions and repulsions) governed by a single timescale. The study offers a fundamental understanding and a predictive framework that could inform research in broader areas of nonlinear pattern formation and collective dynamics.