Vertical mixing of oil droplets by breaking waves

Oil spilled on a sea surface can be dispersed by a variety of natural processes, of which the influence of breaking waves is dominant. Breaking waves are able to split the slick into small droplets, facilitating oil mixing in the water column. Vertical dynamics of the droplets plays a major role in the oil mass exchange between the slick and the water column. In this paper a mathematical model of oil droplet mixing by breaking waves is developed. The model uses a kinetic approach to describe the vertical exchange of the droplets at the interface between the slick and the water column. The majority of the coefficients and parameters are conveniently combined into a single "mixing factor". The model is verified using sensitivity analysis and empirical formulae of other authors. The model permits a rapid estimation of the amount of dispersed oil under the breaking waves. The ultimate goal of the research is to parameterise influence of breaking waves on vertical mixing of oil droplets to be used in a general 3-D oil spill model.     

Assessment of chemical dispersant effectiveness in a wave tank under regular non-breaking and breaking wave conditions

Current chemical dispersant effectiveness tests for product selection are commonly performed with bench-scale testing apparatus. However, for the assessment of oil dispersant effectiveness under real sea state conditions, test protocols are required to have hydrodynamic conditions closer to the natural environment, including transport and dilution effects. To achieve this goal, Fisheries and Oceans Canada and the US Environmental Protection Agency (EPA) designed and constructed a wave tank system to study chemical dispersant effectiveness under controlled mixing energy conditions (regular non-breaking, spilling breaking, and plunging breaking waves). Quantification of oil dispersant effectiveness was based on observed changes in dispersed oil concentrations and oil-droplet size distribution. The study results quantitatively demonstrated that total dispersed oil concentration and breakup kinetics of oil droplets in the water column were strongly dependent on the presence of chemical dispersants and the influence of breaking waves. These data on the effectiveness of dispersants as a function of sea state will have significant implications in the drafting of future operational guidelines for dispersant use at sea.     

Evaluating crude oil chemical dispersion efficacy in a flow-through wave tank under regular non-breaking wave and breaking wave conditions

Testing dispersant effectiveness under conditions similar to that of the open environment is required for improvements in operational procedures and the formulation of regulatory guidelines. To this end, a novel wave tank facility was fabricated to study the dispersion of crude oil under regular non-breaking and irregular breaking wave conditions. This wave tank facility was designed for operation in a flow-through mode to simulate both wave- and current-driven hydrodynamic conditions. We report here an evaluation of the effectiveness of chemical dispersants (Corexit^® EC9500A and SPC 1000^™) on two crude oils (Medium South American [MESA] and Alaska North Slope [ANS]) under two different wave conditions (regular non-breaking and plunging breaking waves) in this wave tank. The dispersant effectiveness was assessed by measuring the water column oil concentration and dispersed oil droplet size distribution. In the absence of dispersants, nearly 8-19% of the test crude oils were dispersed and diluted under regular wave and breaking wave conditions. In the presence of dispersants, about 21-36% of the crude oils were dispersed and diluted under regular waves, and 42-62% under breaking waves. Consistently, physical dispersion under regular waves produced large oil droplets (volumetric mean diameter or VMD ? 300 μm), whereas chemical dispersion under breaking waves created small droplets (VMD ? 50 μm). The data can provide useful information for developing better operational guidelines for dispersant use and improved predictive models on dispersant effectiveness in the field.     

The relationship between oil droplet size and upper ocean turbulence

Oil spilled at sea often forms oil droplets in stormy conditions. This paper examines possible mechanisms which generate the oil droplets. When droplet Reynolds numbers are large, the dynamic pressure force of turbulent flows is the cause of droplet breakup. Using dimensional analysis, Hinze (1955, A.I.Ch.E. Journal 1, 289–295) obtained a formula for the maximum size of oil droplets that can survive the pressure force. When droplet Reynolds numbers are small, however, viscous shear associated with small turbulent eddies is the cause of breakup. For the shear mechanism, we obtain estimates of droplet size as a function of energy dissipation rate, the ratio of oil-to-water viscosity and the surface tension coefficient.

The two formulae are applied to oil spills in the ocean. At dissipation rates expected in breaking waves, the pressure force is the dominant breakup mechanism and can generate oil droplets with radii of hundreds of microns. However, when chemical dispersants are used to treat an oil slick and significantly reduce the oil-water interfacial tension, viscous shear is the dominant breakup mechanism and oil droplets with radii of tens of microns can be generated. Viscous shear is also the mechanism for disintegrating water-in-oil emulsions and the size of a typical emulsion blob is estimated to be tens of millimeters.     

The Relation Between Humidity and Liquid Water Content in Fog: An Experimental Approach

Microphysical measurements of orographic fog were performed above a montane cloud forest in northeastern Taiwan (Chilan mountain site). The measured parameters include droplet size distribution (DSD), absolute humidity (AH), relative humidity (RH), air temperature, wind speed and direction, visibility, and solar short wave radiation. The scope of this work was to study the short term variations of DSD, temperature, and RH, with a temporal resolution of 3?Hz. The results show that orographic fog is randomly composed of various air volumes that are intrinsically rather homogeneous, but exhibit clear differences between each other with respect to their size, RH, LWC, and DSD. Three general types of air volumes have been identified via the recorded DSD. A statistical analysis of the characteristics of these volumes yielded large variabilities in persistence, RH, and LWC. Further, the data revealed an inverse relation between RH and LWC. In principle, this finding can be explained by the condensational growth theory for droplets containing soluble or insoluble material. Droplets with greater diameters can exist at lower ambient RH than smaller ones. However, condensational growth alone is not capable to explain the large observed differences in DSD and RH because the respective growth speeds are too slow to explain the observed phenomena. Other mechanisms play key roles as well. Possible processes leading to the large observed differences in RH and DSD include turbulence induced collision and coalescence, and heterogeneous mixing. More analyses including fog droplet chemistry and dynamic microphysical modeling are required to further study these processes. To our knowledge, this is the first experimental field observation of the anti-correlation between RH and LWC in fog.     

Effects of temperature and wave conditions on chemical dispersion efficacy of heavy fuel oil in an experimental flow-through wave tank

The effectiveness of chemical dispersants (Corexit 9500 and SPC 1000) on heavy fuel oil (IFO180 as test oil) has been evaluated under different wave conditions in a flow-through wave tank. The dispersant effectiveness was determined by measuring oil concentrations and droplet size distributions. An analysis of covariance (ANCOVA) model indicated that wave type and temperature significantly (p < 0.05) affected the dynamic dispersant effectiveness (DDE). At higher temperatures (16 °C), the test IFO180 was effectively dispersed under breaking waves with a DDE of 90% and 50% for Corexit 9500 and SPC 1000, respectively. The dispersion was ineffective under breaking waves at lower temperature (10 °C), and under regular wave conditions at all temperatures (10-17 °C), with DDE < 15%. Effective chemical dispersion was associated with formation of smaller droplets (with volumetric mean diameters or VMD ? 200 μm), whereas ineffective dispersion produced large oil droplets (with VMD ? 400 μm).     

Droplet breakup in subsurface oil releases – Part 1: Experimental study of droplet breakup and effectiveness of dispersant injection

Size distribution of oil droplets formed in deep water oil and gas blowouts have strong impact on the fate of the oil in the environment. However, very limited data on droplet distributions from subsurface releases exist. The objective of this study has been to establish a laboratory facility to study droplet size versus release conditions (rates and nozzle diameters), oil properties and injection of dispersants (injection techniques and dispersant types). This paper presents this facility (6 m high, 3 m wide, containing 40 m³ of sea water) and introductory data. Injection of dispersant lowers the interfacial tension between oil and water and cause a significant reduction in droplet size. Most of this data show a good fit to existing Weber scaling equations. Some interesting deviations due to dispersant treatment are further analyzed and used to develop modified algorithms for predicting droplet sizes in a second paper (Johansen et al., 2013).     

Droplet breakup in subsea oil releases – Part 2: Predictions of droplet size distributions with and without injection of chemical dispersants

A new method for prediction of droplet size distributions from subsea oil and gas releases is presented in this paper. The method is based on experimental data obtained from oil droplet breakup experiments conducted in a new test facility at SINTEF. The facility is described in a companion paper, while this paper deals with the theoretical basis for the model and the empirical correlations used to derive the model parameters from the available data from the test facility. A major issue dealt with in this paper is the basis for extrapolation of the data to full scale (blowout) conditions. Possible contribution from factors such as buoyancy flux and gas void fraction are discussed and evaluated based on results from the DeepSpill field experiment.     

Effects of rainfall on oil droplet size and the dispersion of spilled oil with application to Douglas Channel,British Columbia,Canada

Raindrops falling on the sea surface produce turbulence. The present study examined the influence of rain-induced turbulence on oil droplet size and dispersion of oil spills in Douglas Channel in British Columbia, Canada using hourly atmospheric data in 2011–2013. We examined three types of oils: a light oil (Cold Lake Diluent - CLD), and two heavy oils (Cold Lake Blend - CLB and Access Western Blend - AWB). We found that the turbulent energy dissipation rate produced by rainfalls is comparable to what is produced by wind-induced wave breaking in our study area. With the use of chemical dispersants, our results indicate that a heavy rainfall (rain rate > 20 mm h^? 1) can produce the maximum droplet size of 300 μm for light oil and 1000 μm for heavy oils, and it can disperse the light oil with fraction of 22–45% and the heavy oils of 8–13%, respectively. Heavy rainfalls could be a factor for the fate of oil spills in Douglas Channel, especially for a spill of light oil and the use of chemical dispersants.     

Evolution of droplets in subsea oil and gas blowouts: Development and validation of the numerical model VDROP-J

The droplet size distribution of dispersed phase (oil and/or gas) in submerged buoyant jets was addressed in this work using a numerical model, VDROP-J. A brief literature review on jets and plumes allows the development of average equations for the change of jet velocity, dilution, and mixing energy as function of distance from the orifice. The model VDROP-J was then calibrated to jets emanating from orifices ranging in diameter, D, from 0.5 mm to 0.12 m, and in cross-section average jet velocity at the orifice ranging from 1.5 m/s to 27 m/s. The d₅₀/D obtained from the model (where d₅₀ is the volume median diameter of droplets) correlated very well with data, with an R² = 0.99. Finally, the VDROP-J model was used to predict the droplet size distribution from Deepwater Horizon blowouts. The droplet size distribution from the blowout is of great importance to the fate and transport of the spilled oil in marine environment.     

Chemical dispersant effectiveness testing: influence of droplet coalescence

Thermodynamic and kinetic investigations were performed to determine the influence of coalescence of chemically dispersed crude oil droplets in saline waters. For the range of pH (4-10) and salinity (10 per thousand, 30 per thousand, 50 per thousand ) values studied, zeta-potential values ranged from -3 to -10 mV. As the interaction potential values calculated using Derjaguin-Landau-Verway-Overbeek (DLVO) theory were negative, the electrostatic barrier did not produce significant resistance to droplet coalescence. Coalescence kinetics of premixed crude oil and chemical dispersant were determined within a range of mean shear rates (Gm = 5, 10, 15, 20 s(-1)) and salinity (10 per thousand, 30 per thousand ) values. Coalescence reaction rates were modeled using Smoluchowski reaction kinetics. Measured collision efficiency values (alpha = 0.25) suggest insignificant resistance to coalescence in shear systems. Experimentally determined dispersant efficiencies (alpha = 0.35) were 10-50% lower than that predicted using a non-interacting droplet model (alpha = 0.0). Unlike other protocols in which the crude oil and dispersant are not premixed, salinity effects were not significant in this protocol. This approach allowed the effects of dispersant-oil contact efficiency eta(contact) to be separated from those of water column transport efficiency (eta(transport)) and coalescence efficiency (eta(coalescence)).     

Lab tests on the biodegradation of chemically dispersed oil should consider the rapid dilution that occurs at sea

Most crude oils spread on open water to an average thickness as low as 0.1 mm. The application of dispersants enhances the transport of oil as small droplets into the water column, and when combined with the turbulence of 1 m waves will quickly entrain oil into the top 1 m of the water column, where it rapidly dilutes to concentrations less than 100 ppm. In less than 24 h, the dispersed oil is expected to mix into the top 10 m of the water column and be diluted to concentrations well below 10 ppm, with dilution continuing as time proceeds. Over the multiple weeks that biodegradation takes place, dispersed oil concentrations are expected to be below 1 ppm. Measurements from spills and wave basin studies support these calculations. Published laboratory studies focused on the quantification of contaminant biodegradation rates have used concentrations orders of magnitude greater than this, as it was necessary to ensure the concentrations of hydrocarbons and other chemicals were higher than the detection limits of chemical analysis. However, current analytical methods can quantify individual alkanes and PAHs (and their alkyl homologues) at ppb and ppm levels. To simulate marine biodegradation of dispersed oil at dilute concentrations commonly encountered in the field, laboratory studies should be conducted at similarly low hydrocarbon concentrations.     

Development of a unified oil droplet size distribution model with application to surface breaking waves and subsea blowout releases considering dispersant effects

An oil droplet size model was developed for a variety of turbulent conditions based on non-dimensional analysis of disruptive and restorative forces, which is applicable to oil droplet formation under both surface breaking-wave and subsurface-blowout conditions, with or without dispersant application. This new model was calibrated and successfully validated with droplet size data obtained from controlled laboratory studies of dispersant-treated and non-treated oil in subsea dispersant tank tests and field surveys, including the Deep Spill experimental release and the Deepwater Horizon blowout oil spill. This model is an advancement over prior models, as it explicitly addresses the effects of the dispersed phase viscosity, resulting from dispersant application and constrains the maximum stable droplet size based on Rayleigh-Taylor instability that is invoked for a release from a large aperture.     

A note on Stokes production of turbulence kinetic energy in the oceanic mixed layer: observations in the Baltic Sea

Shear- and convection-driven turbulence coexists with wind-generated surface gravity waves in the upper ocean. The turbulent Reynolds stresses in the oceanic mixed layer can therefore interact with the shear of the wave-generated Stokes drift velocity to extract energy from the surface waves and inject it into turbulence, thus augmenting the mean shear-driven turbulence. Stokes production of turbulence kinetic energy (TKE) is difficult to measure in the field, since it requires simultaneous measurement of the turbulent stress and the Stokes drift profiles in the water column. However, it is readily inferred using second moment closure models of the oceanic mixed layer provided: (1) wave properties are available, along with the usual water mass properties, and radiative and air–sea fluxes needed to drive the mixed layer model and (2) the model skill can be assessed by comparing the model results against the observed dissipation rates of TKE. Comprehensive measurements made during the Reynolds 2002 campaign in the Baltic Sea have made the estimation of Stokes production possible, and in this paper, we report on the effort and the conclusions reached. Measurements of air–sea exchange parameters and water mass properties during the campaign allowed a mixed layer model to be run and the turbulent stress in the water column to be inferred. Simultaneous wave spectrum measurements enabled Stokes drift profile to be deduced and wave breaking to be included in the model run, and the Stokes production of TKE in the water column estimated. Direct measurements of the TKE dissipation rate from an upward traversing microstructure profiler were used to assure that the model could reproduce the turbulent dissipation rate in the water column. The model results indicate that the Stokes production of TKE in the mixed layer is of the same order of magnitude as the shear production and must therefore be included in mixed layer models.     

Turbulence structure in the upper ocean: a comparative study of observations and modeling

Observations of turbulent dissipation rates measured by two independent instruments are compared with numerical model runs to investigate the injection of turbulence generated by sea surface gravity waves. The near-surface observations are made by a moored autonomous instrument, fixed at approximately 8 m below the sea surface. The instrument is equipped with shear probes, a high-resolution pressure sensor, and an inertial motion package to measure time series of dissipation rate and nondirectional surface wave energy spectrum. A free-falling profiler is used additionally to collect vertical microstructure profiles in the upper ocean. For the model simulations, we use a one-dimensional mixed layer model based on a k–ε type second moment turbulence closure, which is modified to include the effects of wave breaking and Langmuir cells. The dissipation rates obtained using the modified k–ε model are elevated near the sea surface and in the upper water column, consistent with the measurements, mainly as a result of wave breaking at the surface, and energy drawn from wave field to the mean flow by Stokes drift. The agreement between observed and simulated turbulent quantities is fairly good, especially when the Stokes production is taken into account.     

Air–water two-phase flow modeling of turbulent surf and swash zone wave motions

Wave breaking and wave runup/rundown have a major influence on nearshore hydrodynamics, morphodynamics and beach evolution. In the case of wave breaking, there is significant mixing of air and water at the wave crest, along with relatively high kinetic energy, so prediction of the free surface is complicated. Most hydrodynamic studies of surf and swash zone are derived from single-phase flow, in which the role of air is ignored. Two-phase flow modeling, consisting of both phases of water and air, may be a good alternative numerical modeling approach for simulating nearshore hydrodynamics and, consequently, sediment transport. A two-phase flow tool can compute more realistically the shape of the free surface, while the effects of air are accounted for. This paper used models based on two-dimensional, two-phase Reynolds-averaged Navier–Stokes equations, the volume-of-fluid surface capturing technique and different turbulence closure models, i.e., k–ε, k–ω and re-normalized group (RNG). Our numerical results were compared with the available experimental data. Comparison of the employed method with a model not utilizing a two-phase flow modeling demonstrates that including the air phase leads to improvement in simulation of wave characteristics, especially in the vicinity of the breaking point. The numerical results revealed that the RNG turbulence model yielded better predictions of nearshore zone hydrodynamics, although the k–ε model also gave satisfactory predictions. The model provides new insights for the wave, turbulence and means flow structure in the surf and swash zones.     

Numerical modelling of wind effects on breaking waves in the surf zone

Wind effects on periodic breaking waves in the surf zone have been investigated in this study using a two-phase flow model. The model solves the Reynolds-averaged Navier–Stokes equations with the k ? ?? turbulence model simultaneously for the flows both in the air and water. Both spilling and plunging breakers over a 1:35 sloping beach have been studied under the influence of wind, with a focus during wave breaking. Detailed information of the distribution of wave amplitudes and mean water level, wave-height-to-water-depth ratio, the water surface profiles, velocity, vorticity, and turbulence fields have been presented and discussed. The inclusion of wind alters the air flow structure above water waves, increases the generation of vorticity, and affects the wave shoaling, breaking, overturning, and splash-up processes. Wind increases the water particle velocities and causes water waves to break earlier and seaward, which agrees with the previous experiment.     

Numerical simulation of surf–swash zone motions and turbulent flow

A two-dimensional numerical model was presented for the simulation of wave breaking, runup and turbulence in the surf and swash zones. The main components of the model are the Reynolds-Averaged Navier–Stokes equations describing the average motion of a turbulent flow, a k–ε turbulence closure model describing the transformation and dissipation processes of turbulence and a volume of fluid technique for tracking the free surface motion. Nearshore wave evolution on a sloping bed, the velocity field and other wave characteristics were investigated. First, the results of the model were compared with experimental results for different surf zone hydrodynamic conditions. Spilling and plunging breakers were simulated and the numerical model investigated for different wave parameters. The turbulence field was also considered and the spatial and time-dependent variations of turbulence parameters were discussed. In the next stage of the study, numerical results were compared with two sets of experimental data in the swash zone. Generally, there is good agreement except for turbulence predictions near the breaking point where the model does not represent well the physical processes. On the other hand, turbulence predictions were found to be excellent for the swash zone. The model provides a precise and efficient tool for the simulation of the flow field and wave transformations in the nearshore, especially in the swash zone. The numerical model can simulate the surface elevation of the vertical shoreline excursion on sloping beaches, while swash–swash interactions within the swash zone are accounted for.     

Meso-scale testing and development of test procedures to maintain mass balance

The Conrad Blucher Institute for Surveying and Science (Texas A&M University––Corpus Christi) has conducted numerous petroleum experiments at the Shoreline Environmental Research Facility (Corpus Christi, Texas, USA). The meso-scale facility has multiple wave tanks, permitting some control in experimental design of the investigations, but allowing for real-world conditions. This paper outlines the evolution of a materials balance approach in conducting petroleum experiments at the facility. The first attempt at a materials balance was during a 1998 study on the fate/effects of dispersant use on crude oil. Both water column and beach sediment samples were collected. For the materials balance, the defined environmental compartments for oil accumulation were sediments, water column, and the water surface, while the discharge from the tanks was presumed to be the primary sink. The “lessons learned” included a need to quantify oil adhesion to the tank surfaces. This was resolved by adhering strips of the polymer tank lining to the tank sides that could be later removed and extracted for oil. Also, a protocol was needed to quantify any floating oil on the water surface. A water surface (oil slick) quantification protocol was developed, involving the use of solid-phase extraction disks. This protocol was first tested during a shoreline cleaner experiment, and later refined in subsequent dispersant effectiveness studies. The effectiveness tests were designed to simulate shallow embayments which created the need for additional adjustments in the tanks. Since dispersant efficacy is largely affected by hydrodynamics, it was necessary to scale the hydrodynamic conditions of the tanks to those expected in our prototype system (Corpus Christi Bay, Texas). The use of a scaled model permits the experiment to be reproduced and/or evaluated under different conditions. To minimize wave reflection in the tank, a parabolic wave dissipater was built. In terms of materials balance, this design reduced available surface area as a sink for oil adsorption.     

Influence of sea surface wind wave turbulence upon wind-induced circulation, tide–surge interaction and bed stress

A three-dimensional finite volume unstructured mesh model of the west coast of Britain, with high resolution in the coastal regions, is used to investigate the role of wind wave turbulence and wind and tide forced currents in producing maximum bed stress in the eastern Irish Sea. The spatial distribution of the maximum bed stress, which is important in sediment transport problems, is determined, together with how it is modified by the direction of wind forced currents, tide–surge interaction and a surface source of wind wave turbulence associated with wave breaking. Initial calculations show that to first order the distribution of maximum bed stress is determined by the tide. However, since maximum sediment transport occurs at times of episodic events, such as storm surges, their effects upon maximum bed stresses are examined for the case of strong northerly, southerly and westerly wind forcing. Calculations show that due to tide–surge interaction both the tidal distribution and the surge are modified by non-linear effects. Consequently, the magnitude and spatial distribution of maximum bed stress during major wind events depends upon wind direction. In addition calculations show that a surface source of turbulence due to wind wave breaking in shallow water can influence the maximum bed stress. In turn, this influences the wind forced flow and hence the movement of suspended sediment. Calculations of the spatial variability of maximum bed stress indicate the level of measurements required for model validation.