A simplified approach to address turbulence modulation in turbidity currents as a response to slope breaks and loss of lateral confinement

Turbidity currents traversing canyon-fan systems flow over bed slopes that decrease in the downstream direction. This slope decrease eventually causes turbidity currents to decelerate and enter a net-depositional mode. When the slope decrease is relatively rapid in the downstream direction, the turbidity current undergoes a concomitantly rapid and substantial transition. Similar conditions are found when turbidity currents debouch to fan systems with loss of lateral confinement. In this work a simplified approach to perform direct numerical simulation of continuous turbidity currents undergoing slope breaks and loss of lateral confinement is presented and applied to study turbulence modulation in the flow. The presence of settling sediment particles breaks the top–bottom symmetry of the flow, with a tendency to self-stratify. This self-stratification damps turbulence, particularly near the bottom wall, affecting substantially the flow’s ability to transport sediment in suspension. This work reports results on two different situations: turbidity currents driven by fine and coarser sediment flowing through a decreasing slope. In the case of fine sediment, after the reduction in the slope of the channel, the flow remains turbulent with only a modest influence on turbulence statistics. In the case of coarse sediments, after the change in slope, turbulence is totally suppressed.     

Froude scaling limitations in modeling of turbidity currents

The scaling problem associated with the modeling of turbidity currents has been recognized but is yet to be explored systematically. This paper presents an analysis of the dimensionless governing equations of turbidity currents to investigate the scale effect. Three types of flow conditions are considered: (i) conservative density current; (ii) purely depositional turbidity current; and (iii) mixed erosional/depositional turbidity current. Two controlling dimensionless numbers, the Froude number and the Reynolds number, appear in the non-dimensional governing equations. When densimetric Froude similarity is satisfied, the analysis shows that the results would be scale-invariant for conservative density current under the rough turbulent condition. In the case of purely depositional flows, truly scale-invariant results cannot be obtained, as the Reynolds-mediated scale effects appear in the bottom boundary conditions of the flow velocity and sediment fall velocity. However, the scale effect would be relatively modest. The Reynolds effect becomes more significant for erosional or mixed erosional/depositional turbidity currents as Reynolds-mediated scale effects also appear in the sediment entrainment relation. Numerical simulations have been conducted at three different scales by considering densimetric Froude scaling alone as well as combined densimetric Froude and Reynolds similarity. Simulation results confirm that although the scaling of densimetric Froude number alone can produce scale-invariable results for conservative density currents, variations occur in the case of turbidity currents. The results become scale invariant when densimetric Froude and Reynolds similarities are satisfied simultaneously.     

Flow structure of unconfined turbidity currents interacting with an obstacle

Driven by a growing importance to engineered structures, investigating the flow characteristics of turbidity currents interacting with a basal obstruction has become popular over the last three decades. However, research has focused on confined studies or numerical simulations, whereas in situ turbidity currents are typically unconfined. The present study investigates experimentally the velocity and turbulence structure of an unconfined turbidity current, in the immediate regions surrounding a rectangular obstacle. Initial density of the current, and substrate condition is varied. Through a novel technique of installing ultrasonic probes within the obstacle, the presence of a velocity recirculation region immediately upstream and downstream of the obstacle is revealed and confirmed with high-resolution imagery. This was found to be comparable to previous confined studies, suggesting that stream-wise velocity profile structure is somewhat independent of confinement. The obstacle was found to reduce velocity and turbulence intensity maxima downstream of the obstacle when compared with unobstructed tests.     

Numerical simulation of particle-driven gravity currents

Particle-driven gravity currents frequently occur in nature, for instance as turbidity currents in reservoirs. They are produced by the buoyant forces between fluids of different density and can introduce sediments and pollutants into water bodies. In this study, the propagation dynamics of gravity currents is investigated using the FLOW-3D computational fluid dynamics code. The performance of the numerical model using two different turbulence closure schemes namely the renormalization group (RNG) ${k-\epsilon}$ scheme in a Reynold-averaged Navier-Stokes framework (RANS) and the large-eddy simulation (LES) technique using the Smagorinsky scheme, were compared with laboratory experiments. The numerical simulations focus on two different types of density flows from laboratory experiments namely: Intrusive Gravity Currents (IGC) and Particle-Driven Gravity Currents (PDGC). The simulated evolution profiles and propagation speeds are compared with laboratory experiments and analytical solutions. The numerical model shows good quantitative agreement for predicting the temporal and spatial evolution of intrusive gravity currents. In particular, the simulated propagation speeds are in excellent agreement with experimental results. The simulation results do not show any considerable discrepancies between RNG ${k-\epsilon}$ and LES closure schemes. The FLOW-3D model coupled with a particle dynamics algorithm successfully captured the decreasing propagation speeds of PDGC due to settling of sediment particles. The simulation results show that the ratio of transported to initial concentration C _o /C _i by the gravity current varies as a function of the particle diameter d _s . We classify the transport pattern by PDGC into three regimes: (1) a suspended regime (d _s is less than about 16 μm) where the effect of particle deposition rate on the propagation dynamics of gravity currents is negligible i.e. such flows behave like homogeneous fluids (IGC); (2) a mixed regime (16 μm < d _s <40 μm) where deposition rates significantly change the flow dynamics; and (3) a deposition regime (d _s ?> 40 μm) where the PDGC rapidly loses its forward momentum due to fast deposition. The present work highlights the potential of the RANS simulation technique using the RNG ${k-\epsilon}$ turbulence closure scheme for field scale investigation of particle-driven gravity currents.     

Feedbacks between bivalve density, flow, and suspended sediment concentration on patch stable states

We explore the role of biophysical feedbacks occurring at the patch scale (spatial scale of tens of meters) that influence bivalve physiological condition and affect patch stability by developing a numerical model for the pinnid bivalve, Atrina zelandica, in cohesive sediments. Simulated feedbacks involve bivalve density, flow conditions (assumed to be primarily influenced by local water depth and peak current speed), suspended sediment concentration (evaluated through a balance between background concentration, deposition, and erosion), and changes in the physiology of Atrina derived from empirical study. The model demonstrates that high bivalve density can lead to skimming flow and to a concomitant decrease in resuspension that will affect suspended sediment concentration over the patch directly feeding back on bivalve physiology. Consequently, for a given flow and background suspended sediment load, the stability of a patch directly depends on the size and density of bivalves in the patch. Although under a range of conditions patch stability is ensured independently of bivalve density, simulations clearly indicate that sudden changes in bivalve density or suspended sediment concentration can substantially affect patch structure and lead to different stable states. The model highlights the role of interactions between organisms, flow, and broader scale environmental conditions in providing a mechanistic explanation for the patchy occurrence of benthic suspension feeders.     

Evaluation of spatiotemporal differences in suspended sediment concentration derived from remote sensing and numerical simulation for coastal waters

Remote sensing and numerical models are often used to monitor the suspended sediment concentration (SSC) in coastal waters; however, the derived SSC varies between the two methods in both space and time. In this study, a method was proposed to assess the spatiotemporal differences in SSC derived from Moderate Resolution Imaging Spectroradiometer (MODIS) images and numerical simulation for coastal waters, using the Bohai Sea in China as an example. An empirical model for SSC retrieval from remote sensed images was initially established. A comparison of the temporal synchronicity over a single day period was performed between the observed data and the numerically simulated results. The range in the SSC at different observation sites was significantly different. Both the SSC values and their daily variation ranges were larger near the estuary of the Yellow River compared with the open area due to the concurrence of tidal flow and the introduction of fresh river water with high turbidity near the estuary. The areas that exhibited spatial differences were defined according to their differences in remotely sensed and numerically simulated SSC distribution patterns. Finally, the reasons for these spatiotemporal differences were discussed. The results provided understanding into the spatiotemporal differences that were introduced when multi-source data were used, thus improving the accuracy of the results when monitoring coastal environments for the management of coastal conservation.     

Laboratory-numerical studies of stratified spin-up flows

The stability of stratified rotating flows is investigated by means of laboratory experiments and numerical simulations in axisymmetric cylindrical and annular containers with both horizontal and sloping bottoms. The baroclinic current is initiated via incremental spin-up/down of a linearly stratified fluid by an abrupt change in the rotation rate of the system. Particular attention is given to the non-linear flow regime (finite Rossby numbers). It is found that axisymmetric spin-up current loses its azimuthal symmetry when the Burger number drops below unity, and breaks into a system of large-scale cyclonic and anticyclonic vortices with predominantly vertical axis of rotation. Eddies always develop at the density fronts formed by the corner regions adjacent to the sidewalls of the container. It is shown that the stability of the spin-up flow is largely affected by the bottom slope and the structure of the bottom boundary layer.     

Simulation-based optimization of in-stream structures design: rock vanes

We employ a three-dimensional coupled hydro-morphodynamic model, the Virtual Flow Simulator (VFS-Geophysics) in its Unsteady Reynolds Averaged Navier–Stokes mode closed with \(k-\omega\) model, to simulate the turbulent flow and sediment transport in large-scale sand and gravel bed waterways under prototype and live-bed conditions. The simulation results are used to carry out systematic numerical experiments to develop design guidelines for rock vane structures. The numerical model is based on the Curvilinear Immersed Boundary approach to simulate flow and sediment transport processes in arbitrarily complex rivers with embedded rock structures. Three validation test cases are conducted to examine the capability of the model in capturing turbulent flow and sediment transport in channels with mobile-bed. Transport of sediment materials is handled using the Exner equation coupled with a transport equation for suspended load. Two representative meandering rivers, with gravel and sand beds, respectively, are selected to serve as the virtual test-bed for developing design guidelines for rock vane structures. The characteristics of these rivers are selected based on available field data. Initially guided by existing design guidelines, we consider numerous arrangements of rock vane structures computationally to identify optimal structure design and placement characteristics for a given river system.     

Mass-based depth and velocity scales for gravity currents and related flows

Gravity driven flows on inclines can be caused by cold, saline or turbid inflows into water bodies. Another example are cold downslope winds, which are caused by cooling of the atmosphere at the lower boundary. In a well-known contribution, Ellison and Turner (ET) investigated such flows by making use of earlier work on free shear flows by Morton, Taylor and Turner (MTT). Their entrainment relation is compared here with a spread relation based on a diffusion model for jets by Prandtl. This diffusion approach is suitable for forced plumes on an incline, but only when the channel topography is uniform, and the flow remains supercritical. A second aspect considered here is that the structure of ET’s entrainment relation, and their shallow water equations, agrees with the one for open channel flows, but their depth and velocity scales are those for free shear flows, and derived from the velocity field. Conversely, the depth of an open channel flow is the vertical extent of the excess mass of the liquid phase, and the average velocity is the (known) discharge divided by the depth. As an alternative to ET’s parameterization, two sets of flow scales similar to those of open channel flows are outlined for gravity currents in unstratified environments. The common feature of the two sets is that the velocity scale is derived by dividing the buoyancy flux by the excess pressure at the bottom. The difference between them is the way the volume flux is accounted for, which—unlike in open channel flows—generally increases in the streamwise direction. The relations between the three sets of scales are established here for gravity currents by allowing for a constant co-flow in the upper layer. The actual ratios of the three width, velocity, and buoyancy scales are evaluated from available experimental data on gravity currents, and from field data on katabatic winds. A corresponding study for free shear flows is referred to. Finally, a comparison of mass-based scales with a number of other flow scales is carried out for available data on a two-layer flow over an obstacle. Mass-based flow scales can also be used for other types of flows, such as self-aerated flows on spillways, water jets in air, or bubble plumes.     

Turbulence suppression by suspended sediment within a geophysical flow

Experiments are performed in a mixing box to evaluate the effect of suspended sediment on turbulence generated by an oscillating grid. Quartz-density sand of varying sizes and concentrations is used, and particle image velocimetry is employed to quantify only the fluid phase. Results show that (1) while a relatively large secondary flow field is present in the box, turbulence is a maximum near the grid and it decreases systematically toward the water surface; (2) relatively high concentrations of fine sediment can markedly alter this secondary flow field and significantly decrease both the time-mean and turbulent kinetic energy within the flow, yet these same sediment concentrations have little effect on the integral time and length scales derived for each velocity component; and (3) the overall turbulence suppression observed can be related to the transfer of energy from the fluid to the sediment and the maintenance of a suspended sediment load rather than commonly employed turbulence modulation criteria. These experimental data demonstrate unequivocally that the presence of a suspended sediment load can significantly reduce overall turbulent kinetic energy, and these results should be applicable to a range of sediment-laden geophysical flows.     

Hierarchical modeling of the dilute transport of suspended sediment in open channels

We propose, discuss and validate a theoretical and numerical framework for sediment-laden, open-channel flows which is based on the two-fluid-model (TFM) equations of motion. The framework models involve mass and momentum equations for both phases (sediment and water) including the interactive forces of drag, lift, virtual mass and turbulent dispersion. The developed framework is composed by the complete two-fluid model (CTFM), a partial two-fluid model (PTFM), and a standard sediment-transport model (SSTM). Within the umbrella of the Reynolds-Averaged Navier-Stokes (RANS) equations, we apply K–ε type closures (standard and extended) to account for the turbulence in the carrier phase (water). We present the results of numerical computations undertaken by integrating the differential equations over control volumes. We address several issues of the theoretical models, especially those related to coupling between the two phases, interaction forces, turbulence closure and turbulent diffusivities. We compare simulation results with various recent experimental datasets for mean flow variables of the carrier as well as, for the first time, mean flow of the disperse phase and turbulence statistics. We show that most models analyzed in this paper predict the velocity of the carrier phase and that of the disperse phase within 10% of error. We also show that the PTFM provides better predictions of the distribution of sediment in the wall-normal direction as opposed to the standard Rousean profile, and that the CTFM is by no means superior to the PTFM for dilute mixtures. We additionally report and discuss the values of the Schmidt number found to improve the agreement between predictions of the distribution of suspended sediment and the experimental data.     

Modeling effects of a tidal barrage on water quality indicator distribution in the Severn Estuary

In this study, emphasis has focused on assessing the potential hydro-environmental impacts of a barrage across the Severn Estuary, with a numerical model being developed and applied to the estuary to assess the impacts of proposed Severn Barrage on the hydrodynamic, sediment transport and faecal indicator organism distribution within the estuary. The results show that the Severn Barrage has the potential to reduce the tidal currents in a highly dynamic estuary. This leads to the reduction of suspended sediment concentrations, which in turn affects the bacterial transport processes which is highly related to the sediment transport processes.     

A model examining hierarchical wetland networks for watershed stormwater management

There is increasing awareness that solutions to degraded quality and excessive quantity of stormwater and resulting impacts on downstream water bodies may require a watershed approach to management rather that the incremental approach that is now common. Examination of low-relief watersheds characteristic of the southeastern coastal plain reveals common hierarchical patterns of surface water convergence that may be emulated in developed watersheds to enhance the efficacy of peak-flow attenuation and pollutant removal. A dynamic systems model was developed to compare stormwater management using a hierarchical network of treatment wetlands with the standard incremental approach wherein treatment systems are designed considering only site-level effluent criteria. The model simulates watershed hydrology, suspended sediment transport and phosphorus removal and transformation. Results indicate that watershed planning of stormwater collection and treatment systems using hierarchical networks can greatly enhance overall effectiveness (annual retention improvements of 31% for flow, 36% for sediment and 27% for phosphorus) with respect to an equal area of uniformly sized wetlands. Further, network proportions can be adjusted to specific runoff characteristics. Distinct roles were observed for each wetland size class: small headwater wetlands effectively removed sediment, medium-sized mid-reach wetlands retained phosphorus, while large wetlands primarily stored and attenuated long-period hydrologic flows.     

Effects of simulated sublethal predation on the growth and regeneration rates of a spionid polychaete in laboratory flumes

Most spionid polychaetes switch from deposit feeding to suspension feeding as current speed and the flux of suspended food increase. Growth rates of juvenile Polydora cornuta are strongly affected by flow and can be as rapid as 60% day⁻¹ in moderate currents. Feeding palps that extend above the sediment–water interface during suspension feeding are especially vulnerable to sublethal predation, but individuals with damaged posteriors are also common. We performed a series of laboratory flume experiments to test the effects of sublethal tissue damage on the growth and regeneration rates of P. cornuta juveniles. Replicated experiments were conducted at three flow speeds in counter-rotating annular flumes containing field-collected sediment and a nonliving algal slurry as deposited and suspended food. In the first set of experiments, we removed 2, 1, or 0 of worms’ two feeding palps and measured the relative growth rates of worm bodies and palps after 3 days in the flumes. Worms that lost both palps grew significantly slower than the other two groups, but the growth rate of worms that had one undamaged palp was not significantly different from worms that had two undamaged palps. Faster flow speeds significantly increased rates of body growth, and there was a significant interaction between flow and the effect of palp loss. During the 3-day experiments, damaged palps fully regenerated and often grew larger than they were prior to being removed. Damaged palps also grew significantly faster than undamaged palps. The second set of experiments tested the effects of removing a worm’s posterior region (~18% of body volume). The growth rates of these damaged and undamaged worms did not differ significantly. By the end of a 3-day flume experiment, damaged worms had grown 6× larger than they were prior to the posterior damage. The rapid regeneration of damaged palps and posterior tissue in moderate flows that allow suspension feeding suggests that sublethal predation on spionids might be more frequent than previously estimated and will have little impact on the growth of juvenile recruits.     

Structure of a dense release produced by varying initial conditions

Buoyancy driven flows such as gravity currents, present in nature and human made applications, are conveyors of particles or dissolved substances for long distances with clear implications for the environment. This transport depends on the triggering conditions of the current. Gravity currents are experimentally simulated under varying initial conditions by combining three different initial buoyancies and five volumes of dense fluid released. The horizontal and vertical structures of the gravity currents are analysed and it is shown that the variation on the initial configuration is conditioning for these. Vertical transport through the gravity current is influenced at the bottom by the solid wall over which the current flows, and at the upper interface by the contact with the ambient water. The relative contribution of shear stress at the bottom and at the upper interfaces are estimated and analysed in terms of the initial triggering conditions of the current; these two compete with the buoyancy, the driver of the current, determining mixing and entrainment. By using a proper parametrization, which accounts for both initial volume of release and location of the observation position relative to the lock, a relation between the resistance at the bottom and at the upper interfaces with the initial conditions of release (i.e. the lock-length) has been found; this is found to be independent of the initial density in the lock. The present study shows that the variation of the initial conditions have consequences on (1) the configuration of the currents and on (2) the hydrodynamics of the currents, including mass and momentum exchanges, which are in addition mutually dependent.     

Heavy Metals in the Bed and Suspended Sediments of Anyang River, Korea: Implications for Water Quality

The objective of this study is to compare Anyang River bed sediments with water chemical composition and to assess the anthropogenic chemical inputs into the river system. Eight sampling locations were chosen along the river channel. Bed and suspended river sediments and water samples were collected, and analyzed for their chemical and physical composition. Data revealed that trace element concentrations in the river water were generally below world average, except for As, Mn, Ni and Cr. Among the three phases: water, bed and suspended sediment, more than 99% of the trace elements was associated with the bed sediment. Concentrations of trace elements in the sediment were a function a particle size distribution and organic content. The calculated degrees of enrichment based on the least influenced sample (ASD 1) indicated the river sediments were enriched with respect to background. The enrichment factors for Pb, Zn and As were relatively lower than for Cr, Co, Ni and Zn. The difference in the enrichment seems to reflect the human activities influence in the basin, and specially for Cd. Speciation of the elements in the five different chemical forms in the sediment by sequential extraction indicated that the reducible fraction was predominant for Fe, Zinc and Cu showed an irregular variation among the different fractions; whereas, Cd and Pb were more regular. Zinc and Cu highly existed mostly in exchangeable forms. Acid soluble and reducible forms were also important for most metals. The speciation implies that the metals associated with the sediment are subject to release into water bodies as goechemical variables (pH and Eh) change. Currently, the introduced metals are deposited near the source area and are mostly associated with the sediment, implying that the river bed sediment acts mainly as a sink, rather than a pool. The accumulated and enriched toxic trace elements can pose a potential pollution of river water.     

A hydromechanical principle for particle retention in Mytilus edulis and other ciliary suspension feeders

It is generally thought that the laterofrontal cirri of the bivalve gill act as filters that retain suspended particles in the through current and transfer the particles onto the frontal surface of the gill filaments. In Mytilus edulis calculations indicated that if water passed between the branching cilia of the cirri that are assumed to constitute the filter the pressure drop needed would amount to about 10 times the actual pressure drop across the whole gill. Thus, instead of acting as filters the laterofrontal cirri seem to move water. Presumably, the cirri together with the frontal cilia produce the water currents along the frontal surface of the gill filaments. Particle retention in the bivalve gill implies the transfer of suspended particles from the current of water about to enter an interfilamentar space into a neighbouring frontal surface current. The complex three-dimensional pattern of flow that arises where the 2 systems of current meet is characterized by steep velocity gradients. Particles that enter such steep, steady velocity gradients become exposed to transverse forces that cause the particles to migrate perpendicularly to the direction of flow. Whether particles enter the surface current, i.e. are retained, or they stay within the through current andescape, depends primarily upon particle size, and upon the steepness and height of the gradients within the boundary zone between the surface current and through current. Further studies are needed to evaluate the capacities and relative importance of this hydromechanical particle-trapping mechanism in suspension feeding bivalves. It is suggested that in downstream particle-retaining systems, e.g. on the tentacles of polychaetes and entoprocts, velocity gradients between through currents and surface currents also act as the particle-collecting mechanism.     

On the fine structure of the thermal bar front

The thermal bar—a hydrodynamic phenomenon, arising in natural basins due to successive changes of the water temperature across the temperature of maximum density (T_m, which is close to 4°C)—has been studied in laboratory experiments and by numerical simulations. The experiments were performed in a rectangular tank with an inclined bottom, filled with water with initial temperature T₀ < T_m and then heated at the surface. During the heating a basin-wide circulation develops, consisting of down-slope cascades in regions where T < T_m, a subsurface off-shore jet in the region where T > T_m, and a compensating flow at intermediate depths towards the shallow part of the tank, supplying both off-shore flows with waters from deeper regions. Analysis of the water temperature and density fields as well as the currents has revealed that the location of the convergence zone of the surface current (when formed) does not coincide with that of the Tm-isotherm. The thermal bar front is typically understood as a convergence zone near the 4°C-isotherm, formed due to the effect of cabbeling. Our experiments demonstrate, however, that the front is associated with the leading edge of the subsurface current. The increasing distance between the 4°C-isotherm and the subsurface jet has been recorded in the laboratory experiments. Numerical simulation results corroborate the laboratory experiments. A scaling analysis predicts the speed of propagation of this frontal zone to be U ~ [g × Δρ/ρ × H]^1/2, where H is the depth (increasing with time) of the upper thermo-active layer, ρ₀ a reference density, and Δρ is the characteristic horizontal density difference across the front. A combined analysis of laboratory, field and numerical data has corroborated this law.     

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.