The turbulent bubble break-up cascade. Part 2. Numerical simulations of breaking waves

Breaking waves generate a distribution of bubble sizes that evolves over time. Knowledge of how this distribution evolves is of practical importance for maritime and climate studies. The analytical framework developed in Part 1 (Chan, Johnson & Moin, J. Fluid Mech., vol. 912, 2021, A42) examined how this evolution is governed by the bubble-mass flux from large- to small-bubble sizes which depends on the rate of break-up events and the distribution of child bubble sizes. These statistics are measured in Part 2 as ensemble-averaged functions of time by simulating ensembles of breaking waves, and identifying and tracking individual bubbles and their break-up events. The large-scale break-up dynamics is seen to be statistically unsteady, and two intervals with distinct characteristics were identified. In the first interval, the dissipation rate and bubble-mass flux are quasi-steady, and the theoretical analysis of Part 1 is supported by all observed statistics, including the expected as small bubbles are also depleted more quickly. This suggests the emergence of different physical mechanisms during different phases of the breaking-wave evolution, although size-local break-up remains a dominant theme. Parts 1 and 2 present an analytical toolkit for population balance analysis in two-phase flows.