Search for: All records

Abstract Microplastics are globally ubiquitous in marine environments, and their concentration is expected to continue rising at significant rates as a result of human activity. They present a major ecological problem with well-documented environmental harm. Sea spray from bubble bursting can transport salt and biological material from the ocean into the atmosphere, and there is a need to quantify the amount of microplastic that can be emitted from the ocean by this mechanism. We present a mechanistic study of bursting bubbles transporting microplastics. We demonstrate and quantify that jet drops are efficient at emitting microplastics up to 280μm in diameter and are thus expected to dominate the emitted mass of microplastic. The results are integrated to provide a global microplastic emission model which depends on bubble scavenging and bursting physics; local wind and sea state; and oceanic microplastic concentration. We test multiple possible microplastic concentration maps to find annual emissions ranging from 0.02 to 7.4—with a best guess of 0.1—mega metric tons per year and demonstrate that while we significantly reduce the uncertainty associated with the bursting physics, the limited knowledge and measurements on the mass concentration and size distribution of microplastic at the ocean surface leaves large uncertainties on the amount of microplastic ejected. 

We describe the rising trajectory of bubbles in isotropic turbulence and quantify the slowdown of the mean rise velocity of bubbles with sizes within the inertial subrange. We perform direct numerical simulations of bubbles, for a wide range of turbulence intensity, bubble inertia and deformability, with systematic comparison with the corresponding quiescent case, with Reynolds number at the Taylor microscale from 38 to 77. Turbulent fluctuations randomise the rising trajectory and cause a reduction of the mean rise velocity$$\tilde {w}_b$$compared with the rise velocity in quiescent flow$$w_b$$. The decrease in mean rise velocity of bubbles$$\tilde {w}_b/w_b$$is shown to be primarily a function of the ratio of the turbulence intensity and the buoyancy forces, described by the Froude number$$Fr=u'/\sqrt {gd}$$, where$$u'$$is the root-mean-square velocity fluctuations,$$g$$is gravity and$$d$$is the bubble diameter. The bubble inertia, characterised by the ratio of inertial to viscous forces (Galileo number), and the bubble deformability, characterised by the ratio of buoyancy forces to surface tension (Bond number), modulate the rise trajectory and velocity in quiescent fluid. The slowdown of these bubbles in the inertial subrange is not due to preferential sampling, as is the case with sub-Kolmogorov bubbles. Instead, it is caused by the nonlinear drag–velocity relationship, where velocity fluctuations lead to an increased average drag. For$$Fr > 0.5$$, we confirm the scaling$$\tilde {w}_b / w_b \propto 1 / Fr$$, as proposed previously by Ruthet al.(J. Fluid Mech., vol. 924, 2021, p. A2), over a wide range of bubble inertia and deformability. 

Thin-film equations are utilised in many different areas of fluid dynamics when there exists a direction in which the aspect ratio can be considered small. We consider thin free films with Marangoni effects in the extensional flow regime, where velocity gradients occur predominantly along the film. In practice, because of the local deposition of surfactants or input of energy, asymmetric distributions of surfactants or surface tension more generally, are possible. Such examples include the surface of bubbles and the rupture of thin films. In this study, we consider the asymmetric thin-film equations for extensional flow with Marangoni effects. Concentrating on the case of small Reynolds number$$ Re $$, we study the deposition of insoluble surfactants on one side of a liquid sheet otherwise at rest and the resulting thinning and rupture of the sheet. The analogous problem with a uniformly thinning liquid sheet is also considered. In addition, the centreline deformation is discussed. In particular, we show analytically that if the surface tension isotherm$$\sigma = \sigma (\varGamma )$$is nonlinear (surface tension$$\sigma$$varies with surfactant concentration$$\varGamma$$), then accounting for top–bottom asymmetry leads to slower (faster) thinning and pinching if$$\sigma = \sigma (\varGamma )$$is convex (concave). The analytical progress reported in this paper allows us to discuss the production of satellite drops from rupture via Marangoni effects, which, if relevant to surface bubbles, would be an aerosol production mechanism that is distinct from jet drops and film drops. 

We present experiments on large air cavities spanning a wide range of sizes relative to the Hinze scale $$d_{H}$$ , the scale at which turbulent stresses are balanced by surface tension, disintegrating in turbulence. For cavities with initial sizes $$d_0$$ much larger than $$d_{H}$$ (probing up to $$d_0/d_{H} = 8.3$$ ), the size distribution of bubbles smaller than $$d_{H}$$ follows $$N(d) \propto d^{-3/2}$$ , with $$d$$ the bubble diameter. The capillary instability of ligaments involved in the deformation of the large bubbles is shown visually to be responsible for the creation of the small bubbles. Turning to dynamical, three-dimensional measurements of individual break-up events, we describe the break-up child size distribution and the number of child bubbles formed as a function of $$d_0/d_{H}$$ . Then, to model the evolution of a population of bubbles produced by turbulent bubble break-up, we propose a population balance framework in which break-up involves two physical processes: an inertial deformation to the parent bubble that sets the size of large child bubbles, and a capillary instability that sets the size of small child bubbles. A Monte Carlo approach is used to construct the child size distribution, with simulated stochastic break-ups constrained by our experimental measurements and the understanding of the role of capillarity in small bubble production. This approach reproduces the experimental time evolution of the bubble size distribution during the disintegration of large air cavities in turbulence. 

Breaking waves modulate the transfer of energy, momentum, and mass between the ocean and atmosphere, controlling processes critical to the climate system, from gas exchange of carbon dioxide and oxygen to the generation of sea spray aerosols that can be transported in the atmosphere and serve as cloud condensation nuclei. The smallest components, i.e., drops and bubbles generated by breaking waves, play an outsize role. This fascinating problem is characterized by a wide range of length scales, from wind forcing the wave field at scales of [Formula: see text](1 km–0.1 m) to the dynamics of wave breaking at [Formula: see text](10–0.1 m); air bubble entrainment, dynamics, and dissolution in the water column at [Formula: see text](1 m–10 μm); and bubbles bursting at [Formula: see text](10 mm–1 μm), generating sea spray droplets at [Formula: see text](0.5 mm–0.5 μm) that are ejected into atmospheric turbulent boundary layers. I discuss recent progress to bridge these length scales, identifying the controlling processes and proposing a path toward mechanistic parameterizations of air–sea mass exchange that naturally accounts for sea state effects.