Experimental Study of the Criteria of Flow Laminarization in Two-Dimensional Dense Gas Plumes

Laminarization of flow in a two-dimensional dense gas plume was experimentally investigated in this study. The plume was created by releasing CO₂ through a ground-level line source into a simulated turbulent boundary layer over an aerodynamically rough surface in a meteorological wind tunnel. The bulk Richardson number (Ri_*), based on negative plume buoyancy, plume thickness, and friction velocity, was varied over a wide range so that the effects of stable stratification on plume laminarization could be observed. A variety of ambient wind speeds as well as three different sizes of roughness arrays were used so that possible effects of roughness Reynolds number (Re_*) on plume laminarization could also be identified. Both flow visualization methods and quantitative measurements of velocity and intermittency of turbulence were used to provide quantitative assessments of plume laminarization.Flow visualization provided an overall picture of how the plume was affected by the negative buoyancy. With increasing Ri_*, both the plume depth and the vertical mixing were significantly suppressed, while upstream propagation of the plume from the source was enhanced. The most important feature of the flow revealed by visualization was the laminarization of flow in the lower part of the plume, which appeared to be closely related to both Ri_* and Re_*.Measurements within the simulated dense gas plumes revealed the influence of the stable stratification on mean velocity and turbulence intensity profiles. Both the mean velocity and turbulence intensity were significantly reduced near the surface; and these reductions systematically depended on Ri_*. The roughness Reynolds number also had considerable influence on the mean flow and turbulence structure of the dense gas plumes.An intermittency analysis technique was developed and applied to the digitized instantaneous velocity signals. It not only confirmed the general flow picture within the dense plume indicated by the flow visualization, but also clearly demonstrated the changes of flow regime with variations in Ri_* and Re_*. Most importantly, based on this intermittency analysis, simple criteria for characterizing different flow regimes are formulated; these may be useful in predicting when plume laminarization might occur.     

Oil spreading in instantaneous and continuous spills on rotating earth

The effect of the Coriolis force on the oil spill spreading in the gravity-viscous regime is examined. A new shallow water model for the transport and spreading of oil slick of arbitrary shape is described in which the Coriolis force is included in the momentum equations and the oil–water friction is parameterized in a frame of the boundary layer theory including the Ekman friction. The numerical Lagrangian method based on smoothed particle dynamics is described. New similarity solutions of the model equations are obtained for unidirectional and axisymmetric spreading in gravity-viscous, gravity-turbulent and gravity-viscous-rotational regimes for instantaneous as well as continuous releases. The numerical simulation extends these results for the case of continuous release in the presence of currents. It was shown that Coriolis term in the momentum equation can be omitted if slick thickness is much less of the laminar Ekman layer thickness. However, the Ekman friction should be retained for slicks of any thickness for larger times. The Ekman friction results in the essential slowdown of the spreading as well as in the deflection of the oil spreading velocity at 45° from the direction of velocity in the non-rotation case. Numerical simulations of large-scale spills showed that after the 2?days the slick area with the Coriolis effect was approximately less than half of that without rotation. Therefore, the earth rotation can be also important in the oil weathering.     

Dredged Material Disposal at the Edge of the Florida Current

A field data collection project was undertaken to investigate the short-term fate of dredged material discharged in the designated Miami Ocean Dredged Material Disposal Site (ODMDS) before dredging of the Miami River and the Miami Harbor Turning Basin begins. the designated ODMDS is located in relatively deep water for discharge sites with a typical bottom depth of 150 metres and is also located in the western boundary region of the Gulf Stream current off Miami. Acoustical backscattering, current, particulate, temperature and salinity data were gathered over a three day period from April 24, 1990 through April 26, 1990. the major generic features of shallow-water discharge plumes were observed to be present: (a) the presence of a rapid convective descending plume portion; (b) impact of that plume portion with the ocean bottom and concomitant generation of a bottom surge; (c) rapid horizontal width growth of the descending plume through entrainment; and (d) retention of a residual plume portion within the water column. A well-mixed upper water column layer extending to a depth of 40 to 60 metres below the surface of the ocean permitted measurements of the plume entrainment coefficient free from bottom boundary, water column density gradient, and vertical current shear effects which are usually present in relatively shallow, e.g. less than 40 metres bottom depth, coastal ocean discharge studies. Entrainment coefficient estimates obtained in this study were between 0.5 to 0.7. the residual water plume material was tracked over one-half hour during each of eight discharge events and was transported in a north-northeast direction.     

Dredged Material Disposal at the Edge of the Florida Current

A field data collection project was undertaken to investigate the short-term fate of dredged material discharged in the designated Miami Ocean Dredged Material Disposal Site (ODMDS) before dredging of the Miami River and the Miami Harbor Turning Basin begins. the designated ODMDS is located in relatively deep water for discharge sites with a typical bottom depth of 150 metres and is also located in the western boundary region of the Gulf Stream current off Miami. Acoustical backscattering, current, particulate, temperature and salinity data were gathered over a three day period from April 24, 1990 through April 26, 1990. the major generic features of shallow-water discharge plumes were observed to be present: (a) the presence of a rapid convective descending plume portion; (b) impact of that plume portion with the ocean bottom and concomitant generation of a bottom surge; (c) rapid horizontal width growth of the descending plume through entrainment; and (d) retention of a residual plume portion within the water column. A well-mixed upper water column layer extending to a depth of 40 to 60 metres below the surface of the ocean permitted measurements of the plume entrainment coefficient free from bottom boundary, water column density gradient, and vertical current shear effects which are usually present in relatively shallow, e.g. less than 40 metres bottom depth, coastal ocean discharge studies. Entrainment coefficient estimates obtained in this study were between 0.5 to 0.7. the residual water plume material was tracked over one-half hour during each of eight discharge events and was transported in a north-northeast direction.     

Determination of the entrainment coefficient of a pure plume using the salt-bath technique

Environmental Fluid Mechanics - The entrainment coefficient of the pure plume, the ratio of the radial velocity of the entraining fluid at the edge of the plume to the axial velocity, is...     

Blockage of saline intrusions in restricted,two-layer exchange flows across a submerged sill obstruction

Results are presented from a series of large-scale experiments investigating the internal and near-bed dynamics of bi-directional stratified flows with a net-barotropic component across a submerged, trapezoidal, sill obstruction. High-resolution velocity and density profiles are obtained in the vicinity of the obstruction to observe internal-flow dynamics under a range of parametric forcing conditions (i.e. variable saline and fresh water volume fluxes; density differences; sill obstruction submergence depths). Detailed synoptic velocity fields are measured across the sill crest using 2D particle image velocimetry, while the density structure of the two-layer exchange flows is measured using micro-conductivity probes at several sill locations. These measurements are designed to aid qualitative and quantitative interpretation of the internal-flow processes associated with the lower saline intrusion layer blockage conditions, and indicate that the primary mechanism for this blockage is mass exchange from the saline intrusion layer due to significant interfacial mixing and entrainment under dominant, net-barotropic, flow conditions in the upper freshwater layer. This interfacial mixing is quantified by considering both the isopycnal separation of vertically-sorted density profiles across the sill, as well as calculation of corresponding Thorpe overturning length scales. Analysis of the synoptic velocity fields and density profiles also indicates that the net exchange flow conditions remain subcritical (G < 1) across the sill for all parametric conditions tested. An analytical two-layer exchange flow model is then developed to include frictional and entrainment effects, both of which are needed to account for turbulent stresses and saline entrainment into the upper freshwater layer. The experimental results are used to validate two key model parameters: (1) the internal-flow head loss associated with boundary friction and interfacial shear; and (2) the mass exchange from the lower saline layer into the upper fresh layer due to entrainment.     

A model of vertical distribution of suspended matter in an open channel flow

A quasi-stationary model of vertical distribution of concentration of suspended particular matter in the bottom layer of 1D open channel with a sloped bottom and varying free surface slope is discussed. The model proceeds from the balance between the turbulent diffusion and settling with the buoyancy flux effects on the medium turbulence neglected. The model outcome is formulated in the form of an analytic formula for the vertical distribution of concentration. It is shown that the derived formula embraces two basic types of vertical distribution of concentration, one with a monotonic decrease of concentration gradient and the other with a gradient maximum (lutocline) located at some distance from the bottom. The first distribution type realizes for a relatively large settling velocity or low intensity of turbulence and the second type for a small settling velocity or high intensity of turbulence. The skill of the model to mimic realistic situations is demonstrated on data measured in the Jiaojiang Estuary (China).     

A simplified computational analysis of turbulent plumes and jets

The application of computational fluid dynamics (CFD), particularly Large Eddy Simulation, for the modelling of buoyant turbulent plumes, has been demonstrated to be very accurate, but computationally expensive. Here a more basic, and therefore more generally practicable, approach is presented. Commercial CFD software is used to model such plumes using Reynolds-Averaged Navier-Stokes (RANS) turbulence models. A careful comparison is made between the numerical predictions and well-established results regarding the bulk properties of plumes. During this process, we are able to observe the well-known approximate Gaussian nature of the plume and achieve quantitative agreement with empirical plume spread coefficients. The use of numerical modelling allows for the investigation of the flow field and turbulence in those regions of the plume of most interest—the plume edge and near source regions. A comprehensive sensitivity study is conducted to identify the limits of applicability of this modelling approach. It is shown that the standard modelling approach of Morton, Taylor and Turner, which introduced the well-known entrainment assumption, pertains in a region well above the source region. At the plume edge, the levels of turbulence are contrasted with the value of the entrainment parameter. Finally, the effects of forcing the plumes with additional momentum at the source are considered, including the case of a pure jet. We show how these forced plumes eventually lose their momentum excess and tend to the behaviour of a pure, buoyant plume.     

The Information Content of a Scalar Plume – A Plume Tracing Perspective

The ability of many animals and insects to track a plume to its source is a particularly impressive feat when the fluid dynamics is considered. Inspired by this observation this research seeks to identify the information in a passive scalar plume suitable for developing robust and efficient plume tracing algorithms. The subject of this study is a scalar plume emanating from a point source in a turbulent boundary layer which has been modeled in a laboratory facility built specifically for this purpose. A coupled PIV-LIF technique is used to measure the velocity and scalar field in a time resolved fashion. This data set is analyzed and the convergence rates of five single-point statistics, suitable as kernels of plume tracing algorithms, are investigated. The experimental data shows that the scalar fluctuations over long downstream distances from the source are characterized by filamentary structures that lead to relatively slow convergence rates for any statistic that is based on mean concentrations. The scalar intermittency, however, converges rapidly toward its true value, in fact converging to a testable hypothesis for source location direction faster than the time scale of the larger scale plume meander.     

Enrichments of Trace Metals in Particulate Matter From the Southern East China Sea,North of Taiwan

Abstract

This study investigates the distributions and enrichments of trace metals in suspended and sinking particulate matter from southern East China Sea (ECS) north of Taiwan during the period April 1992 to April 1993. According to these results, concentration of suspended particulate matter in the inner shelf of southern China Sea, the upwellinginfluenced shelf break, and in Kuroshio water are 1.30 (surface)–4.2 (bottom) mg1^?1, ca. 0.4 mg1^?1 and 0.1–0.2 mg1^?1, respectively, reflecting various influences of terrestrial inputs. A benthic nepheloid layer (BNL), apparently owing to resuspension of local and/or remote bottom sediments, formed over the shelf region. Temporal variations in trace metal contents and enrichments in suspended matter from the shelf region reflect the variation of metal inputs from Chinese rivers, particularly from the Changjiang runoff. the enriched metals are more likely to be derived from anthropogenic input, rather than from biological accumulation. in addition, a decrease in metal contents and an increase in salinity confirm the transport of suspended particulate metals from the East China Sea shelf to the open ocean. the feature of metal plume in the intermediate layer (550–800m) of Kuroshio water also verifies this occurrence. Moreover, the sinking particles collected from a sediment trap on the upper slope are relatively enriched in lithogenic matter and trace metals, suggesting the deposit of anthropogenic metals in the slope area.     

Enrichments of Trace Metals in Particulate Matter From the Southern East China Sea, North of Taiwan

This study investigates the distributions and enrichments of trace metals in suspended and sinking particulate matter from southern East China Sea (ECS) north of Taiwan during the period April 1992 to April 1993. According to these results, concentration of suspended particulate matter in the inner shelf of southern China Sea, the upwellinginfluenced shelf break, and in Kuroshio water are 1.30 (surface)-4.2 (bottom) mg1^-1, ca. 0.4 mg1^-1 and 0.1-0.2 mg1^-1, respectively, reflecting various influences of terrestrial inputs. A benthic nepheloid layer (BNL), apparently owing to resuspension of local and/or remote bottom sediments, formed over the shelf region. Temporal variations in trace metal contents and enrichments in suspended matter from the shelf region reflect the variation of metal inputs from Chinese rivers, particularly from the Changjiang runoff. the enriched metals are more likely to be derived from anthropogenic input, rather than from biological accumulation. in addition, a decrease in metal contents and an increase in salinity confirm the transport of suspended particulate metals from the East China Sea shelf to the open ocean. the feature of metal plume in the intermediate layer (550-800m) of Kuroshio water also verifies this occurrence. Moreover, the sinking particles collected from a sediment trap on the upper slope are relatively enriched in lithogenic matter and trace metals, suggesting the deposit of anthropogenic metals in the slope area.     

Spreading of oil on water in the surface-tension regime

The third stage of oil spreading on water, in which surface-tension force promotes spreading against the resisting viscous effect, is investigated using a similarity solution in combination with an integral boundary-layer technique to solve the unidirectional oil-spreading dynamics problem in the last stage of spreading. The thin layer is assumed to be supplied by oil from a bulk boundary. Using a constitutive equation for oil-film surface tension versus oil-film thickness, analytical solutions near the bulk boundary and near the edge are developed. Using the asymptotic solutions to initiate integration, the differential equations for the oil thickness, oil velocity, and boundary-layer profiles are integrated starting from the leading edge and bulk boundary, which after matching provide a complete solution. The results for the spreading-law prefactors are found to differ by about 10 % from published theoretical results using the same constitutive equation. Using an empirical constitutive equation for oil-film surface tension versus distance from the bulk boundary leads to a spreading-law prefactor that is in excellent agreement with the published experimental result and published theoretical work providing and using the same empirical constitutive equation.     

Integral models for bubble,droplet, and multiphase plume dynamics in stratification and crossflow

We present the development and validation of a numerical modeling suite for bubble and droplet dynamics of multiphase plumes in the environment. This modeling suite includes real-fluid equations of state, Lagrangian particle tracking, and two different integral plume models: an Eulerian model for a double-plume integral model in quiescent stratification and a Lagrangian integral model for multiphase plumes in stratified crossflows. Here, we report a particle tracking algorithm for dispersed-phase particles within the Lagrangian integral plume model and a comprehensive validation of the Lagrangian plume model for single- and multiphase buoyant jets. The model utilizes literature values for all entrainment and spreading coefficients and has one remaining calibration parameter \( \kappa \), which reduces the buoyant force of dispersed phase particles as they approach the edge of a Lagrangian plume element, eventually separating from the plume as it bends over in a crossflow. We report the calibrated form \( \kappa = [(b - r) / b]^4 \), where b is the plume half-width, and r is the distance of a particle from the plume centerline. We apply the validated modeling suite to simulate two test cases of a subsea oil well blowout in a stratification-dominated crossflow. These tests confirm that errors from overlapping plume elements in the Lagrangian integral model during intrusion formation for a weak crossflow are negligible for predicting intrusion depth and the fate of oil droplets in the plume. The Lagrangian integral model has the added advantages of being able to account for entrainment from an arbitrary crossflow, predict the intrusion of small gas bubbles and oil droplets when appropriate, and track the pathways of individual bubbles and droplets after they separate from the main plume or intrusion layer.     

Gravity currents in rotating,wedge-shaped,adverse channels

Results are presented from a series of parametric experimental and analytical studies of the behaviour of dense gravity currents along rotating, up-sloping, wedge-shaped channels. High resolution density profile measurements at fixed cross- and along-channel locations reveal the outflowing bottom gravity currents to adjust to quasi-steady, geostrophically-balanced conditions along the channels, with the outflow layer thickness and cross-channel interface slope shown to scale with the inlet Burger number for all experimental conditions tested. A general analytical solution to the classic rotating hydraulics problem has been developed under the assumption of inviscid, zero-potential-vorticity conditions to model dense water flow through a triangular constriction and thus simulate the vee-channel configurations under consideration. Predictions from this zero-PV model are shown to provide good overall quantitative agreement with experimental measurements obtained both under hydraulically-controlled conditions at the channel exit and for subcritical conditions generated along the channel length. Quantitative discrepancies between measurements and analytical predictions are attributed primarily to assumptions and limitations associated with the zero-PV modelling approach adopted, as well as the to the rapid adjustment in outflow characteristics as the channel exit is approached, as characterised by the along-channel variation in densimetric Froude number for the outflows.     

Streamwise velocity profiles in coastal currents

The ability to model marine currents can be a powerful device for many planning activities, for which the knowledge of the velocity field is of pivotal importance, such as the evaluation of current-induced loading on maritime structures or the diffusion and dispersion of polluted flow discharges. Observations of time-averaged velocity profiles, taken with a vessel mounted acoustic Doppler current profiler during a monitoring survey program in the seas of Southern Italy, are analysed in this paper. The measurements were taken under non-breaking conditions, offshore the surf zone, with the aim of reproducing the vertical trends of the streamwise velocity by means of standard theoretical laws. To evaluate also the possible influence of stratification on the current velocity profile shape, together with velocity measurements water temperature and salinity were also measured at the same time and locations, by means of a CTD recorder. The examined surveys referred to different time periods and sites, to guarantee a general validity of deductions. On the basis of the experiments, we verified the actual existence of a log-layer and concluded that the upper limit of the region in which the log law is applicable extends well beyond the inner region. Moreover, the deviations of the measured velocity from the logarithmic profiles above the height of the log layer is consistent with the effects of stratification. The parameters of the log law were estimated, depending on both flow dynamics and stratification in the target area. As a second step, in the most superficial and stratified layer, the velocity profiles were modelled by means of a power law, which fitted the measured data well. According to previous studies, the power law parameters result Reynolds number dependent by means of a new proposed formulation. Finally, the bottom stress and the bottom drag coefficient were investigated.     

Effect of initial excess density and discharge on constant flux gravity currents propagating on a slope

The effect of the upstream conditions on propagation of gravity current over a slope is investigated using three-dimensional numerical simulations. The current produced by constant buoyancy flux, is simulated using a large eddy simulation solver. The dense saline solution used at the inlet is the driving force of the flow. Higher replenishment of the current is possible either by a high inflow discharge or high initial fractional density excess. In the simulations, it is observed that these two parameters affect the flow in different ways. Results show that the front speed of the descending current is proportional to the cube root of buoyancy flux, $(g_o^{\prime } Q)^{1/3}$ , which agrees with the previous experimental and numerical observations. The height of the tail of the current grows linearly in the streamwise direction. Formation of a strong shear layer at the boundary of mixed upper layer and dense lower layer is observed within the body and the tail of the current. Over the tail of the current far enough from the inlet, the vertical velocity and density profiles are compared to the ones from an experimental study. Distance from the bed to the point of maximum velocity increases with an increase in inflow discharge, while it remains practically unchanged with increasing initial fractional excess density in the simulations. Even though the velocity profiles are in good agreement, some discrepancies are observed in fractional excess density profiles among experimental and numerical results. Possible reasons for these discrepancies are discussed. Generally, gravity current type of flows could be expressed in layer-integrated formulation of governing equations. However, layer integration introduces several constants, commonly known as shape factors, to the equations of motion. The values of these shape factors are calculated based on simulation results and compared to the values from experiments and to the favorably used ‘top hat’ assumption.     

Hydrodynamics and turbulence in emergent and sparsely vegetated open channel flow

This present study reports the results of an experimental study characterizing thorough variation of turbulent hydrodynamics and flow distribution in emergent and sparsely vegetated open channel flow. An emergent and rigid sparse vegetation patch with regular spacing between stems along the flow and transverse directions was fixed in the central region of the cross-section of open channel. Experiments were conducted in subcritical flow conditions and velocity measurements were obtained with an acoustic Doppler Velocimetry system. Large variations of the turbulence intensities, Reynolds shear stress, turbulent kinetic energy and vortical motions are found in and around the vegetation patch. At any cross-section through the interior of the vegetation patch, streamwise velocity decreases with increase in streamwise length and the velocity profiles converge from the log-law to a linear profile with increasing slope. Time-averaged lateral and vertical velocities inside the vegetation patch increase with increasing streamwise distance and converge from negative values to positive values. Turbulence intensities interior of the sparse vegetation patch are more than those of without the vegetation patch. Similar to the trend of streamwise velocity profiles inside the vegetation, turbulence intensities and longitudinal-normal Reynolds shear stress profile decreases with streamwise direction. In the interior of the vegetation patch and downstream of the trailing edge, turbulent kinetic energy profiles are exhibiting irregular fluctuations and the maximum values are occurring in the outer layer. Analysis of flow distribution confirms sparse vegetation patch is inducing a serpentine flow pattern in its vicinity. At the leading edge, flow is rushing towards the right hand sidewall, and at the trailing edge, flow is turning to the left hand sidewall. In between the leading and trailing edges, the streamlines are following a zig-zag fashion at varied degree along the streamwise and lateral directions. Immediate upstream of the leading edge and in the interior of the vegetation patch, vortex motion is clearly visible and the vortices are stretched along the width of the channel with streamwise direction.     

Diel vertical migrations and nocturnal feeding of a dense coastal krill scattering layer (Thysanoessa raschi and Meganyctiphanes norvegica) in stratified surface waters

To study the nocturnal feeding of euphausiids during vertical migrations and its impact on the phytoplankton, a phytoplankton-rich water mass (drogue marked) drifting over a dense krill scattering layer (acoustic 104kHz) in the lower St. Lawrence estuary was monitored over 46 h in July 1982. Phytoplankton >20 m was abundant and mostly concentrated at the bottom of the photic layer above the pycnocline. Less than 42% of the particulate carbon was due to phytoplankton. The krill scattering layer was about 2 to 3 km in width, elongated along the 100-m bathymetric contour, and absent when the bottom was shallower than 50 m. Its upper day depth was 50 m. At deeper depths, its vertical distribution frequently changed from unimodal to polymodal shapes and viceversa, often with large concentrations of zooplankton just above the bottom. Typical vertical migrations were observed on both days. At night the scattering layer had a lower scattering strength. Most of it was below the thermocline but net catches showed that large concentrations of euphausiids (up to 57 individuals m^-3) crossed it. Stomach pigment content of Thysanoessa raschi was generally low, but mean stomach fullness was always high. They were more opportunistic than herbivorous. From stomach fullness and the presence of a food bolus in mouth parts, feeding in surface waters appeared to be intensive, but gut content indicated that food was not processed there. It is therefore suggested that individuals underwent vertical interchanges across the thermocline while feeding during the night. Meganyctiphanes norvegica had significant herbivorous activity during the night. The grazing pressure impact of the scattering layer on phytoplankton was negligible.Contribution to the program of GIROQ (Groupe interuniversitaire de recherches océanographiques du Québec)     

Model of turbulence penetration into a suspension layer on a sediment bed and effect on vertical solute transfer

The vertical diffusional mass (solute) transfer through a suspended sediment layer, e.g. at the bottom of a lake, reservoir or estuary, by the propagation of velocity fluctuations from above was investigated. The attenuation of the velocity fluctuations in the suspension layer and the associated effect on solute transfer through the suspension layer was simulated. To represent large eddies traveling downstream in water over a high-concentration suspended sediment layer, a streamwise velocity fluctuation moving in downstream direction was imposed along the upper boundary of the suspension layer. Velocity fluctuations and downstream velocity were normalized by the shearvelocity (U_*) at the top of the suspension layer. Streamwise and vertical velocity components inside the suspension layer, were obtained from the 2-D continuity and the Navier–Stokes equations. The persistence of turbulence with depth—as it penetrates from the overlying water into the suspension layer—was found to depend on its amplitude, its period, and on the apparent viscosity of the suspension. The turbulence was found to propagate efficiently into the suspension layer when its frequency is low, and the apparent viscosity of the suspension is high. Effects on vertical mass transfer were parameterized by penetration depth and effective diffusion coefficient, and related to apparent viscosity of the suspension, Schmidt number and shear velocity on top of the suspension layer. The enhancement of turbulence penetration by viscosity is similar to the flow near an oscillating flat plate (Stokes’ second problem), but is opposite to turbulence penetration into a stationary porous and permeable sediment bed. The information is applicable to water quality modeling mear the sediment/water interface of lakes, river impoundments and estuaries.