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.

In this paper we present an experimental and theoretical study of weak bubble plumes in unstratified and stationary water. We define a weak bubble plume as one that spreads slower than the linear rate of a classic plume. This work focuses on the characteristics of the mean flow in the plume, including centreline velocity, plume spreading and entrainment of ambient water. A new theory based on diffusive spreading instead of an entrainment hypothesis is used to describe the lateral spreading of the bubbles and the associated plume. The new theory is supported by the experimental data. With the measured data of liquid volume fluxes and the new theory, we conclude that the weak bubble plume is a decreasing entrainment process, with the entrainment coefficient $\unicode[STIX]{x1D6FC}$ in the weak bubble plume decreasing with height $z$, following $\unicode[STIX]{x1D6FC}\sim z^{-1/2}$, and taking on values much smaller than those in a classic bubble plume. An additional non-dimensional diffusion coefficient, $\hat{E_{t}}\sim E_{t}U_{s}^{2}/B_{0}$, is also needed to describe the evolution of the volume and kinematic momentum fluxes for the mean flow in the weak bubble plume. Here, $E_{t}$ is the effective turbulent diffusion coefficient, $U_{s}$ is the terminal rise velocity of the bubbles, and $B_{0}$ is the kinematic buoyancy flux of the source. Finally, we provide a unified framework for the mean flow characteristics, including volume flux, momentum flux and plume spreading for the classic and weak bubble plumes, which also provides insight on the transition from classic to weak bubble plume behaviour.

We present the turbulent kinetic energy (t.k.e.) budget of a dilute bubble plume in its asymptotic state. The budget is derived from an experimental dataset of bubble plumes formed inside an unstratified water tank. The experiments cover both the adjustment phase and asymptotic state of the plume. The diameters $d$ of air bubbles are in the range 1–4 mm and the air void fraction $\unicode[STIX]{x1D6FC}_{g}$ is between 0.7 % and 1.8 %. We measured the three components of the instantaneous liquid velocity vector with a profiling acoustic Doppler velocimeter. From the experiments, we found the following inside the heterogeneous bubble core of the plume: (i) the probability density functions of the standardized liquid fluctuations are very similar to those of homogeneous bubble swarms rising with and without background liquid turbulence; (ii) the characteristic temporal frequency $f_{cwi}$ at which bubbles inject t.k.e. into the liquid agrees with the prediction $f_{cwi}=0.14u_{s}/d$ observed and theoretically derived for homogeneous bubble swarms ($u_{s}$ is the bubble slip velocity); (iii) the liquid turbulence is anisotropic with the ratio of turbulence intensities between the vertical and horizontal components in the range 1.9–2.1; (iv) the t.k.e. production by air bubbles is much larger than that by liquid mean shear; and (v) an increasing fraction of the available work done by bubbles is deposited into liquid turbulence as one moves away from the plume centreline. Together with the existing knowledge of homogeneous bubble swarms, our results of the heterogeneous bubble plume support the view that millimetre-sized bubbles create specific patterns of liquid fluctuations that are insensitive to flow conditions and can therefore be possibly modelled by a universal form.

Characteristics of laboratory-scale bubble-driven buoyant plumes in a stably stratified quiescent fluid are studied using large-eddy simulation (LES). As a bubble plume entrains stratified ambient water, its net buoyancy decreases due to the increasing density difference between the entrained and ambient fluids. A large fraction of the entrained fluid eventually detrains and falls along an annular outer plume from a height of maximum rise (peel height) to a neutral buoyancy level (trap height), during which less buoyant scalars (e.g. small droplets) are trapped and dispersed horizontally, forming quasi-horizontal intrusion layers. The inner/outer double-plume structure and the peel/intrusion process are found to be more distinct for cases with small bubble rise velocity, while weak and unstable when the slip velocity is large. LES results are averaged to generate distributions of mean velocity and turbulent fluxes. These distributions provide data for assessing the performance of previously developed closures used in one-dimensional integral plume models. In particular, the various LES cases considered in this study yield consistent behaviour for the entrainment coefficients for various plume cases. Furthermore, a new continuous peeling model is derived based on the insights obtained from LES results. Comparing to previous peeling models, the new model behaves in a more self-consistent manner, and it is expected to provide more reliable performance when applied in integral plume models.