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.     

Modelling hydraulic jump using the bubbly two-phase flow method

Hydraulic jumps have complex flow structures, characterised by strong turbulence and large air contents. It is difficult to numerically predict the flows. It is necessary to bolster the existing computer models to emphasise the gas phase in hydraulic jumps, and avoid the pitfall of treating the phenomenon as a single-phase water flow. This paper aims to improve predictions of hydraulic jumps as bubbly two-phase flow. We allow for airflow above the free surface and air mass entrained across it. We use the Reynolds-averaged Navier–Stokes equations to describe fluid motion, the volume of fluid method to track the interface, and the k–ε model for turbulence closure. A shear layer is shown to form between the bottom jet flow and the upper recirculation flow. The key to success in predicting the jet flow lies in formulating appropriate bottom boundary conditions. The majority of entrained air bubbles are advected downstream through the shear layer. Predictions of the recirculation region’s length and air volume fraction within the layer are validated by available measurements. The predictions show a linear growth of the shear layer. There is strong turbulence at the impingement, and the bulk of the turbulence kinetic energy is advected to the recirculation region via the shear layer. The predicted bottom-shear-stress distribution, with a peak value upstream of the toe of the jump and a decaying trend downstream, is realistic. This paper reveals a significant transient bottom shear stress associated with temporal fluctuations of mainly flow velocity in the jump. The prediction method discussed is useful for modelling hydraulic jumps and advancing the understanding of the complex flow phenomenon.     

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.     

Visualization of the internal flow properties and the material exchange interface in an entraining viscous Newtonian gravity current

In order to simulate a simple entraining geophysical flow, a viscous Newtonian gravity current is released from a reservoir by a dam-break and flows along a rigid horizontal bed until it meets a layer of entrainable material of finite depth, identical to the current. The goal is to examine the entrainment mechanisms by observing the interaction between the incoming flow and the loose bed. The sole parameter varied is the initial volume of the gravity current, thus altering its height and velocity. The gravity current plunges or spills into the entrainable bed and the velocity of the flow front becomes linear with time. The bed material is directly affected: motion is generated in the fluid far downstream of, and in that lying beneath the encroaching front. Shear bands are identified, separating horizontal flow downstream from flow with a strong vertical component close to the step. Downstream of the step the flow is horizontal and stratified, with no slip on the bottom boundary and very low shear near the surface. Between these two regions may lie transitional zones with linear velocity profiles, separated by horizontal bands of high shear; the number of transitional zones in the cross-section varies with the initial volume of the dam-break.     

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.     

Dressler’s theory for curved topography flows: iterative derivation,transcritical flow solutions and higher-order wave-type equations

The Dressler equations are a system of two non-linear partial differential equations for shallow fluid flows over curved topography. The theory originated from an asymptotic stretching method formulating the equations of motion in terrain-fitted curvilinear coordinates. Apparently, these equations failed to produce a transcritical flow profile changing from sub- to supercritical flow conditions. Further, wave-like motions over a flat bottom are excluded because the bed-normal velocity component is not accounted for. However, the theory was found relevant for several environmental flow problems including density currents over mountains and valleys, seepage flow in hillslope hydrology, the development of antidunes, the formation of geological deposits from hyper-concentrated flows, and shallow-water flow modeling in hydraulics. In this work, Dressler’s theory is developed in an alternative way by a systematic iteration of the stream and potential functions in terrain-fitted coordinates. The first iteration was found to be the former Dressler’s theory, whereas a second iteration of the governing equations results in velocity components generalizing Dressler’s theory to wave-like motion. Dressler’s first-order theory produces a transcritical flow solution over topography only if the total head is fixed by a minimum value of the specific energy at the transition point. However, the theory deviates from measurements under subcritical flow conditions, given that the bed-normal velocity component is significant. A second iteration to the velocity field was used to produce a second-order differential equation that resembles the cnoidal-wave theory. It accurately produces flow over an obstacle including the critical point and the minimum specific energy as part of the numerical solution. The new cnoidal-wave model compares well with the theory of a Cosserat surface for directed fluid sheets, whereas the Saint-Venant theory appears to be poor.     

Two-phase modeling of sediment clouds

A sediment cloud release in stagnant ambient fluid occurs in many engineering applications. Examples include land reclamation and disposal of dredged materials. The detailed modeling of the distinct characteristics of both the solid and fluid phases of the sediment cloud is hitherto unavailable in the literature despite their importance in practice. In this paper, the two-phase mixing characteristics of the sediment cloud are investigated both experimentally and theoretically. Experiments were carried out to measure the transient depth penetration and the lateral spread of the sediment cloud and its entrained fluid using the laser induced fluorescence technique, with a range of particle sizes frequently encountered in the field (modeled at laboratory scale). A two-phase model of the sediment cloud that provides detailed predictions of the mixing characteristics of the individual phases is also proposed. The entrained fluid characteristics are solved by an integral model accounting for the buoyancy loss (due to particle separation) in each time step. The flow field induced by the sediment cloud is approximated by a Hill’s spherical vortex centered at the centroid and with the size of the entrained fluid. The particle equation of motion under the effect of the induced flow governs each computational particle. A random walk model using the hydrodynamic diffusion coefficient is used to account for the random fluctuation of particles in the dispersive regime. Overall, the model predictions of the two-phase mixing characteristics are in good agreement with the experimental data for a wide range of release conditions.     

Characterization of bedload intermittency near the threshold of motion using a Lagrangian sediment transport model

At the smallest scales of sediment transport in rivers, the coherent structures of the turbulent boundary layer constitute the fundamental mechanisms of bedload transport, locally increasing the instantaneous hydrodynamic forces acting on sediment particles, and mobilizing them downstream. Near the critical threshold for initiating sediment motion, the interactions of the particles with these unsteady coherent structures and with other sediment grains, produce localized transport events with brief episodes of collective motion occurring due to the near-bed velocity fluctuations. Simulations of these flows pose a significant challenge for numerical models aimed at capturing the physical processes and complex non-linear interactions that generate highly intermittent and self-similar bedload transport fluxes. In this investigation we carry out direct numerical simulations of the flow in a rectangular flat-bed channel, at a Reynolds number equal to Re = 3632, coupled with the discrete element method to simulate the dynamics of spherical particles near the bed. We perform two-way coupled Lagrangian simulations of 48,510 sediment particles, with 4851 fixed particles to account for bed roughness. Our simulations consider a total of eight different values of the non-dimensional Shields parameter to study the evolution of transport statistics. From the trajectory and velocity of each sediment particle, we compute the changes in the probability distribution functions of velocities, bed activity, and jump lengths as the Shields number increases. For the lower shear stresses, the intermittency of the global bedload transport flux is described by computing the singularity or multifr actal spectrum of transport, which also characterizes the widespread range of transport event magnitudes. These findings can help to identify the mechanisms of sediment transport at the particle scale. The statistical analysis can also be used as an ingredient to develop larger, upscaled models for predicting mean transport rates, considering the variability of entrainment and deposition that characterizes the transport near the threshold of motion.     

Large-scale turbulent structure of uniform shallow free-surface flows

The hydrodynamics of super- and sub-critical shallow uniform free-surface flows are assessed using laboratory experiments aimed at identifying and quantifying flow structure at scales larger than the flow depth. In particular, we provide information on probability distributions of horizontal velocity components, their correlation functions, velocity spectra, and structure functions for the near-water-surface flow region. The data suggest that for the high Froude number flows the structure of the near-surface layer resembles that of two-dimensional turbulence with an inverse energy cascade. In contrast, although large-scale velocity fluctuations were also present in low Froude number flow its behaviour was different, with a direct energy cascade. Based on our results and some published data we suggest a physical explanation for the observed behaviours. The experiments support Jirka’s [Jirka GH (2001) J Hydraul Res 39(6):567–573] hypothesis that secondary instabilities of the base flow may generate large-scale two-dimensional eddies, even in the absence of transverse gradients in the time-averaged flow properties.     

Numerical modelling of horizontal sediment-laden jets

Sediment-laden turbulent flows are commonly encountered in natural and engineered environments. It is well known that turbulence generates fluctuations to the particle motion, resulting in modulation of the particle settling velocity. A novel stochastic particle tracking model is developed to predict the particle settling out and deposition from a sediment-laden jet. Particle velocity fluctuations in the jet flow are modelled from a Lagrangian velocity autocorrelation function that incorporates the physical mechanism leading to a reduction of settling velocity. The model is first applied to study the settling velocity modulation in a homogeneous turbulence field. Consistent with basic experiments using grid-generated turbulence and computational fluid dynamics (CFD) calculations, the model predicts that the apparent settling velocity can be reduced by as much as 30 % of the stillwater settling velocity. Using analytical solution for the jet mean flow and semi-empirical RMS turbulent velocity fluctuation and dissipation rate profiles derived from CFD predictions, model predictions of the sediment deposition and cross-sectional concentration profiles of horizontal sediment-laden jets are in excellent agreement with data. Unlike CFD calculations of sediment fall out and deposition from a jet flow, the present method does not require any a priori adjustment of particle settling velocity.     

Vertical dense jet in flowing current

The discharge of brackish water, as a dense jet in a natural water body, by the osmotic power plants, undergoes complex mixing processes and has significant environmental impacts. This paper focuses on the mixing processes that develop when a dense round jet outfall perpendicularly enters a shallow flowing current. Extensive experimental measurements of both the salinity and the velocity flow fields were conducted to investigate the hydrodynamic jet behavior within the ambient current. Experiments were carried out in a closed circuit flume at the Coastal Engineering Laboratory (LIC) of the Technical University of Bari (Italy). The salinity concentration and velocity fields were analyzed, providing a more thorough knowledge about the main features of the jet behavior within the ambient flow, such as the jet penetration, spreading, dilution, terminal rise height and its impact point with the flume lower boundary. In this study, special attention is given to understand and confirm the conjecture, not yet experimentally demonstrated, of the development and orientation of the jet vortex structures. Results show that the dense jet is almost characterized by two distinct phases: a rapid ascent phase and a gradually descent phase. The measured flow velocity fields definitely confirm the formation of the counter-rotating vortices pair, within the jet cross-section, during both the ascent and descent phases. Nevertheless, the experimental results show that the counter-rotating vortices pair of both phases (ascent and descent) are of opposite rotational direction.     

The Structure of the Shear Layer in Flows over Rigid and Flexible Canopies

Flume experiments were conducted with rigid and flexible model vegetation to study the structure of coherent vortices (a manifestation of the Kelvin–Helmholtz instability) and vertical transport in shallow vegetated shear flows. The vortex street in a vegetated shear layer creates a pronounced oscillation in the velocity profile, with the velocity near the top of a model canopy varying by a factor of three during vortex passage. In turn, this velocity oscillation drives the coherent waving of flexible canopies. Relative to flows over rigid vegetation, the oscillation in canopy geometry has the effect of decreasing the amount of turbulent vertical momentum transport in the shear layer. Using a waving plant to determine phase in the vortex cycle, each vortex is shown to consist of a strong sweep at its front (during which the canopy is most deflected), followed by a weak ejection at its rear (when the canopy height is at a maximum). Whereas in unobstructed mixing layers the vortices span the entire layer, they encompass only 70% of the flexibly obstructed shear layer studied here.     

Flows across high aspect ratio street canyons: Reynolds number independence revisited

The Reynolds number for flow in a street canyon, Re?=?U_refH/ν (where U_ref is a reference velocity, H the street canyon height, and ν the kinematic viscosity), cannot be matched between reduced-scale experiments and full-scale field measurements. This mismatch is often circumvented by satisfying the Re independence criterion, which states that above a critical Re (Re_c), the flow field remains invariant with Re. Re_c?=?11,000 is often adopted in reduced-scale experiments. In deep street canyons with height-to-width aspect ratio ≥?1.5, reduced-scale experiments have shown two recirculation vortices induced by the mean flows, but full-scale field measurements have observed only one vortex. We investigated this discrepancy by conducting water channel experiments with Re between 10⁴ and 10⁵ at three aspect ratios. The canyons with aspect ratio 1.0 have Re_c?=?11,000, the canyons with aspect ratio 1.5 have Re_c between 31,000 and 58,000, while the canyons with aspect ratio 2.0 have Re_c between 57,000 and 87,000. Therefore, the widely adopted Re_c?=?11,000 is not applicable for canyons with aspect ratio greater than 1.5. Our results also confirm that there is only one vortex in deep canyons at high Re. This single-vortex flow regime could change our fundamental understanding of deep canyons, which are often assumed to exhibit multiple-vortex flow regimes. Applications such as numerical model validation based on the multiple-vortex regime should be revisited. Our experimental data with Re up to 10⁵ could be used to validate numerical models at high Re.     

How do larvae attach to a solid in a laminar flow?

A hydrodynamic model explaining the mechanism of contact of marine larvae in vertical flows is presented. Two hydrodynamic factors—flow vorticity and larval self-propulsion—are the key components in the mathematical model. It is shown that flow vorticity causes a larva to rotate and change the direction of self-thrust, thus leading to its migration across the mean flow. The latter motion is of an oscillatory nature. Contact will be enabled only for sufficiently large amplitudes of oscillations. Simple expressions for the probability of initial contact are obtained for two-dimensional Couette and Poiseuille flows. The three-dimensional motion of a larva in a tube is studied using the Monte-Carlo simulations. It is shown that contact probability depends mainly on the ratio of the characteristic flow velocity and the larva’s swimming speed. The theoretical results compare favorably with available experimental data. Possible applications of the method and results presented here to the classical problem of larval attachment to bodies of general geometry are briefly discussed in the concluding section.     

Mobility of free surface in different liquids and its influence on water striders locomotion

Dynamics of the surface layer in different liquids is examined by means of infrared thermography of the surface and simultaneous velocity fields measurements using surface and infrared Particle Image Velocimetry. This technique allows measurements and comparison of two velocity fields—at the surface and at small depth about 50–200 μm. In distilled water the velocity fields at the surface and at small depth exhibit significant dissimilarity. The flow field below the surface is essentially 3D, whereas the surface flow is characterized by vanishing 2D divergence of velocity, indicating predominantly planar motion. In contrast, in ethanol–butanol mixture two velocity fields are well correlated, both corresponding to 3D flow with continuous surface renewal. Thermal patterns, observed at the surface, and the flow field structure in different liquids are associated with different boundary conditions for velocity at the surface. Water surface is seldom renewed, which inhibits heat and mass exchange between the liquid and atmosphere. However, absence of vertical advection also enables organisms to live within the surface layer, to stand and walk on the free surface. This is illustrated by the difficulties a water strider faces on the surface of ultrapure water, which exhibits Marangoni convection.     

Two-phase flow insights into open-channel flows with suspended particles of different densities

The effect of particle density on the turbulent open-channel flow carrying dilute particle suspensions is investigated using two specific gravities and three concentrations of solid particles. The particles, identical in size and similar in shape, were natural sand and a neutrally buoyant plastic. The particles were fully suspended, and formed no particle streaks on the channel’s bed. Accordingly, the changes in the flow are attributed to the interactions between suspended particles and flow turbulence structures. Measurements were obtained by means of image velocimetry enabling simultaneous, but distinct, measurement of liquid and particle velocities. The experimental results show that, irrespective of particle specific gravity, particle suspension influences bulk velocity of flow and the Kármán coefficient, while friction velocity essentially remains constant. The results also show that particles in suspension modify local water turbulence over the flow depth, but in ways not accurately predicted using the customary parameters for characterizing turbulence modification.     

Fourier analysis of the roll-up and merging of coherent structures in shallow mixing layers

The current study investigates the role of nonlinearity in the development of two-dimensional coherent structures (2DCS) in shallow mixing layers. A nonlinear numerical model based on the depth-averaged shallow water equations is used to investigate temporal shallow mixing layers, where the mapping from temporal to spatial results is made using the velocity at the center of the mixing layer. The flow is periodic in the stream-wise direction and the transmissive boundary conditions are used in the cross-stream boundaries to prevent reflections. The numerical results are examined with the aid of Fourier decomposition. Results show that the previous success in applying local linear theory to shallow mixing layers does not imply that the flow is truly linear. Linear stability theory is confirmed to be only valid within a short distance from the inflow boundary. Downstream of this linear region, nonlinearity becomes important for the roll-up and merging of 2DCS. While the energy required for the merging of 2DCS is still largely provided by the velocity shear, the merging mechanism is one where nonlinear mode interaction changes the velocity field of the subharmonic mode and the gradient of the along-stream velocity profile which, in turn, changes the magnitude of the energy production of the subharmonic mode by the velocity shear implicitly. The nonlinear mode interaction is associated with energy up-scaling and is consistent with the inverse energy cascade which is expected to occur in shallow shear flows. Current results also show that such implicit nonlinear interaction is sensitive to the phase angle difference between the most unstable mode and its subharmonic. The bed friction effect on the 2DCS is relatively small initially and grows in tandem with the size of the 2DCS. The bed friction also causes a decrease in the velocity gradient as the flow develops downstream. The transition from unstable to stable flow occurs when the bed friction balances the energy production. Beyond this point, the bed friction is more dominant and the 2DCS are progressively damped and eventually get annihilated. The energy production by the velocity shear plays an important role from the upstream end all the way to the point of transition to stable flow. The fact that linear stability theory is valid only for a short distance from the inflow boundary suggests that some elements of nonlinearity is incorporated in the mean velocity profile in experiments by the averaging process. The implicit nature of nonlinear interaction in shallow mixing layers and the sensitivity of the nonlinear interaction to phase angle difference between the most unstable mode and its subharmonic allows local linear theory to be successful in reproducing features of the instability such as the dominant mode of the 2DCS and its amplitude.     

Experimental investigation of meandering jets in shallow reservoirs

Meandering flows in rectangular shallow reservoirs were experimentally investigated. The characteristic frequency, the longitudinal wave length and the mean lateral extension of the meandering jet were extracted from the first paired modes, obtained by a proper orthogonal decomposition of the surface velocity field measured by large scale PIV. The depth-normalised characteristic lengths and the Strouhal number were then compared to the main dimensionless numbers characterizing the experiments: Froude number, friction number and reservoir shape factor. The normalised wave length and mean lateral extension of the meandering jet are neither correlated with the Froude number nor with the reservoir shape factor; but a clear relationship is found with the friction number. Similarly, the Strouhal number is found proportional to a negative power of the friction number. In contrast, the Froude number and the reservoir shape factor enable to predict the occurrence of a meandering flow pattern: meandering jets occur for Froude number greater than 0.21 and for a shape factor smaller than 6.2.     

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.