Showering in Flint, MI: Is there a DBP problem?

Lead contamination in the City of Flint, MI has been well documented over the past two years, with lead levels above the EPA Action Level until summer 2016. This resulted from an ill-fated decision to switch from Detroit water (Lake Huron) with corrosion control, to Flint River water without corrosion control. Although lead levels are now closer to normal, reports of skin rashes have sparked questions surrounding tap water in some Flint homes. This study investigated the presence of contaminants, including disinfection by-products (DBPs), in the hot tap water used for showering in the homes of residents in Flint. Extensive quantitative analysis of 61 regulated and priority unregulated DBPs was conducted in Flint hot and cold tap water, along with the analysis of 50 volatile organic compounds and a nontarget comprehensive, broadscreen analysis, to identify a possible source for the reported skin rashes. For comparison, chlorinated hot and cold waters from three other cities were also sampled, including Detroit, which also uses Lake Huron as its source water. Results showed that hot water samples generally contained elevated levels of regulated and priority unregulated DBPs compared to cold water samples, but trihalomethanes were still within regulatory limits. Overall, hot shower water from Flint was similar to waters sampled from the three other cities and did not have unusually high levels of DBPs or other organic chemicals that could be responsible for the skin rashes observed by residents. It is possible that an inorganic chemical or microbial contaminant may be responsible.     

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.     

Buoyancy-Driven Two-Layer Exchange Flows Across a Slowly Submerging Barrier

Results are presented from a combined analytical and laboratory study of unsteady, two-layer, density-driven, sub-maximal exchange over a slowly-descending estuarine barrier located very close to the junction of the river mouth and the near-shore coastal zone. As in the precursor study (Cuthbertson et al. 2004, Environ. Fluid Mech. 4, 127–155) that the present investigation extends, the rate of descent of the barrier is assumed to be sufficiently slow for the unsteady exchange flow to adjust continuously to the appropriate quasi-steady conditions at every stage of the descent. The results demonstrate that the thickness of each layer at the barrier crest can be predicted satisfactorily by a hydraulic analysis that (i) assumes the existence of a single control point at the barrier crest and (ii) incorporates the hydraulic losses arising from the sudden expansion and contraction of the upper and lower layers, respectively, at the channel exit. Predictions of the normalised elevations of the interface at the barrier and exit for the “inviscid” maximal exchange case are shown to coincide with the maximal exchange predictions of Zhu and Lawrence (2000, J. Hydraul. Eng. ASCE 126(12), 921–928).     

A Modelling Study of Transient,Buoyancy-Driven Exchange Flow over a Descending Barrier

Results are presented from a series of model studies of the transient exchange flow resulting from the steady descent of an impermeable barrier separating initially-quiescent fresh and saline water bodies having density ₀ and ₀+()₀, respectively. A set of parametric laboratory experiments has been carried out (i) to determine the characteristic features of the time-dependent exchange flow over the barrier crest and (ii) to quantify the temporal increase in the thickness and spatial extent of the brackish water reservoir formed behind the barrier by the outflowing, partly-mixed saline water. The results of the laboratory experiments have been compared with the predictions of a theoretical model adapted from the steady, so-called maximal exchange flow case and good qualitative agreement between theory and experiment has been demonstrated. The comparisons indicate that head losses of between 7% and 3% are applicable to the flow over the ridge crest in the early and late stages, respectively, of the barrier descent phase, with these losses being attributed to mixing processes associated with the counterflowing layers of fresh and saline water in the vicinity of the ridge crest. The experimental data show (and the theoretical model predictions confirm) that (i) the dimensionless time of detection t _det(g/H _b)^1/2 of the brackish water pool fed by the dense outflow increases (at a given distance from the barrier) with increasing values of the descent rate parameter g'H _b/(dh _b/dt)² and (ii) the normalised thickness (x,t)/H _b of the pool at a given reference station increases monotonically with increasing values of the modified time (t–t _det)/(H _b/g)^1/2, with the rate of thickening decreasing with increasing values of the descent rate parameter g'H _b(dh _b/dt)². Here, g=(g/₀)()₀ is the modified gravitational acceleration, H _b is the mean depth of the water and dh _b/dt denotes the rate of descent of the barrier height h _b with elapsed time t after the two water bodies are first brought into contact.     

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.     

Rotational effects on exchange flows across a submerged sill

Environmental Fluid Mechanics - This paper presents new laboratory-scale numerical simulations of density-driven exchange flows generated across an idealised, submerged sill obstruction under both...