Search

Understanding the size of oil droplets released from a jet in crossflow is crucial for estimating the trajectory of hydrocarbons and the rates of oil biodegradation/dissolution in the water column. We present experimental results of an oil jet with a jet-to-crossflow velocity ratio of 9.3. The oil was released from a vertical pipe 25 mm in diameter with a Reynolds number of 25 000. We measured the size of oil droplets near the top and bottom boundaries of the plume using shadowgraph cameras and we also filmed the whole plume. In parallel, we developed a multifluid large eddy simulation model to simulate the plume and coupled it with our VDROP population balance model to compute the local droplet size. We accounted for the slip velocity of oil droplets in the momentum equation and in the volume fraction equation of oil through the local, mass-weighted average droplet rise velocity. The top and bottom boundaries of the plume were captured well in the simulation. Larger droplets shaped the upper boundary of the plume, and the mean droplet size increased with elevation across the plume, most likely due to the individual rise velocity of droplets. At the same elevation across the plume, the droplet size was smaller at the centre axis as compared with the side boundaries of the plume due to the formation of the counter-rotating vortex pair, which induced upward velocity at the centre axis and downward velocity near the sides of the plume.

We investigate and quantify stirring due to chaotic advection within a steady, three-dimensional, Ekman-driven, rotating cylinder flow. The flow field has vertical overturning and horizontal swirling motion, and is an idealization of motion observed in some ocean eddies. The flow is characterized by strong background rotation, and we explore variations in Ekman and Rossby numbers, $E$ and ${R}_{o} $, over ranges appropriate for the ocean mesoscale and submesoscale. A high-resolution spectral element model is used in conjunction with linear analytical theory, weakly nonlinear resonance analysis and a kinematic model in order to map out the barriers, manifolds, resonance layers and other objects that provide a template for chaotic stirring. As expected, chaos arises when a radially symmetric background state is perturbed by a symmetry-breaking disturbance. In the background state, each trajectory lives on a torus and some of the latter survive the perturbation and act as barriers to chaotic transport, a result consistent with an extension of the KAM theorem for three-dimensional, volume-preserving flow. For shallow eddies, where $E$ is $O(1)$, the flow is dominated by thin resonant layers sandwiched between KAM-type barriers, and the stirring rate is weak. On the other hand, eddies with moderately small $E$ experience thicker resonant layers, wider-spread chaos and much more rapid stirring. This trend reverses for sufficiently small $E$, corresponding to deep eddies, where the vertical rigidity imposed by strong rotation limits the stirring. The bulk stirring rate, estimated from a passive tracer release, confirms the non-monotonic variation in stirring rate with $E$. This result is shown to be consistent with linear Ekman layer theory in conjunction with a resonant width calculation and the Taylor–Proudman theorem. The theory is able to roughly predict the value of $E$ at which stirring is maximum. For large disturbances, the stirring rate becomes monotonic over the range of Ekman numbers explored. We also explore variation in the eddy aspect ratio.