Influence of the secondary motions on pollutant mixing in a meandering open channel flow

This paper presents large eddy simulation of turbulent flow in a meandering open channel with smooth wall and rectangular cross-section. The Reynolds number based on the channel height is 40,000 and the aspect ratio of the cross-section is 4.48. The depth-averaged mean stream-wise velocity agree well to experimental measurements. In this specific case, two interacting cells are formed that swap from one bend to the other. Transport and mixing of a pollutant is analysed using three different positions of release, e.g. on the inner bank, on the outer bank and on the centre of the cross section. The obtained depth-average mean concentration profiles are reasonably consistent with available experimental data. The role of the secondary motions in the mixing processes is the main focus of the discussion. It is found that the mixing when the scalar is released on the centre of the cross-section is stronger and faster than the mixing of the scalar released on the sides. When the position of release is close to a bank side, the mixing is weaker and a clear concentration of scalar close to the corresponding side-wall can be observed in both cases.     

A numerical model of wave- and current-driven nutrient uptake by coral reef communities

We developed a numerical model capable of simulating the spatial zonation of nutrient uptake in coral reef systems driven by hydrodynamic forcing (both from waves and currents). Relationships between nutrient uptake and bed stress derived from flume and field studies were added to a four-component biogeochemical model embedded within a three-dimensional (3-D) hydrodynamic ocean model coupled to a numerical wave model. The performance of the resulting coupled physical-biogeochemical model was first evaluated in an idealized one-dimensional (1-D) channel for both a pure current and a combined wave-current flow. Waves in the channel were represented by an oscillatory flow with constant amplitude and frequency. The simulated nutrient concentrations were in good agreement with the analytical solution for nutrient depletion along a uniform channel, as well as with existing observations of phosphate uptake across a real reef flat. We then applied this integrated model to investigate more complex two-dimensional (2-D) nutrient dynamics, firstly to an idealized coral reef-lagoon morphology, and secondly to a realistic section of Ningaloo Reef in Western Australia, where nutrients were advected into the domain via alongshore coastal currents. Both the idealized reef and Ningaloo Reef simulations showed similar patterns of maximum uptake rates on the shallow forereef and reef crest, and with nutrient concentration decreasing as water flowed over the reef flat. As a result of the cumulative outflow of nutrient-depleted water exiting the reef channels and then being advected down the coast by alongshore currents, both reef simulations exhibited substantial alongshore variation in nutrient concentrations. The coupled models successfully reproduced the observed spatial-variability in nitrate concentration across the Ningaloo Reef system.     

Characteristics of flow structure of free-surface flow in a partly obstructed open channel with vegetation patch

Free-surface flows over patchy vegetation are common in aquatic environments. In this study, the hydrodynamics of free-surface flow in a rectangular channel with a bed of rigid vegetation-like cylinders occupying half of the channel bed was investigated and interpreted by means of laboratory experiments and numerical simulations. The channel configurations have low width-to-depth aspect ratio (1.235 and 2.153). Experimental results show that the adjustment length for the flow to be fully developed through the vegetation patch in the present study is shorter than observed for large-aspect-ratio channels in other studies. Outside the lateral edge of the vegetation patch, negative velocity gradient (\(\partial \overline{u}/\partial z < 0\)) and a local velocity maximum are observed in the vertical profile of the longitudinal velocity in the near-bed region, corresponding to the negative Reynolds stress (\(- \overline{{u^{\prime}w^{\prime}}} < 0\)) at the same location. Assuming coherent vortices to be the dominant factor influencing the mean flow field, an improved Spalart–Allmaras turbulence model is developed. The model improvement is based on an enhanced turbulence length scale accounting for coherent vortices due to the effect of the porous vegetation canopy and channel bed. This particular flow characteristic is more profound in the case of high vegetation density due to the stronger momentum exchange of horizontal coherent vortices. Numerical simulations confirmed the local maximum velocity and negative gradient in the velocity profile due to the presence of vegetation and bed friction. This in turn supports the physical interpretation of the flow processes in the partly obstructed channel with vegetation patch. In addition, the vertical profile of the longitudinal velocity can also be explained by the vertical behavior of the horizontal coherent vortices based on a theoretical argument.     

Experimental study of a positive surge. Part 1: basic flow patterns and wave attenuation

A positive surge results from a sudden change in flow that increases the depth. It is the unsteady flow analogy of the stationary hydraulic jump and a geophysical application is the tidal bore. Positive surges are commonly studied using the method of characteristics and the Saint-Venant equations. The article presents the results from new experimental investigations conducted in a large rectangular channel. Detailed unsteady velocity measurements were performed with a high temporal resolution using acoustic Doppler velocimetry and non-intrusive free-surface measurement devices. Several experiments were conducted with the same initial discharge (Q = 0.060 m³/s) and six different gate openings after closure resulting in both non-breaking undular and breaking bores. The analysis of undular surges revealed wave amplitude attenuation with increasing distance of surge propagation were in agreement with Ippen and Kulin theory. Also, undular wave period and wave length data were relatively close to the values predicted by the wave dispersion theory for gravity waves in intermediate water depths.     

Turbulent velocity profile in fully-developed open channel flows

The determination of velocity profile in turbulent narrow open channels is a difficult task due to the significant effects of the anisotropic turbulence that involve the Prandtl’s second type of secondary flow occurring in the cross section. With these currents the maximum velocity appears below the free surface that is called dip phenomenon. The well-known logarithmic law describes the velocity distribution in the inner region of the turbulent boundary layer but it is not adapted to define the velocity profile in the outer region of narrow channels. This paper relies on an analysis of the Navier–Stokes equations and yields a new formulation of the vertical velocity profile in the center region of steady, fully developed turbulent flows in open channels. This formulation is able to predict time averaged primary velocity in the outer region of the turbulent boundary layer for both narrow and wide open channels. The proposed law is based on the knowledge of the aspect ratio and involves a parameter C_Ar depending on the position of the maximum velocity (ξ_dip). ξ_dip may be derived, either from measurements or from an empirical equation given in this paper. A wide range of longitudinal velocity profile data for narrow open channels has been used for validating the model. The agreement between the measured and the computed velocities is rather good, despite the simplification used.     

Morphology of submarine channel-levee systems in the eastern Bay of Bengal near Andaman region

Multibeam echo-sounder (MBES), sub-bottom profiler, singlebeam tracks and public-domain bathymetry geophysical datasets were analyzed in a Geographic Information System (GIS) for accuracy assessment and further used to update the previous map (Curray et al. Mar Pet Geol 19: 1191–1223, 2003) of submarine channel-levee systems in the eastern Bengal Fan near the Andaman region. Comparative analysis of singlebeam with MBES indicated 2.8 % of root-mean-square difference in bathymetry, and could clearly bring out distinct submarine channel profiles. The channel map revealed that the previously named E6 [updated as E7(1a)] and E7 [updated as E7(1b)] channels emanate from the middle and lower Bengal Fan, respectively and merge together at 91.28°E, 10.64°N where they form a second order channel [E7(2)] that flows southward along the east of Ninety-east Ridge. Similarly, two first order channels [E6(1a) and E6(1b)] emanate from the lower Bengal Fan which merge together at 88.78°E, 06.92°N to form a second order channel [E6(2)] that flows along the west of Ninety-east Ridge towards the equator. This channel was previously thought to be linked to the middle Bengal Fan. High-resolution MBES tracks over the channel shows extensive meandering patterns with the channel gently deepening and widening towards south. Variations in the morphology and stratigraphy indicate that the E7 channel is characterized by a curvilinear course, with higher observed sinuosity wherever gradient decreases.     

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.     

Influence of planform geometry and momentum ratio on thermal mixing at a stream confluence with a concordant bed

The effects of planform geometry and momentum flux ratio on thermal mixing at a stream confluence with concordant bed morphology are investigated based on numerical simulations that can capture the dynamics of large-scale turbulence. In two simulations, the bathymetry and asymmetrical planform geometry are obtained from field experiments and the momentum flux ratio is set at values of one and four. These two conditions provide the basis for studying differences in thermal mixing processes at this confluence when the wake mode and the Kelvin–Helmholtz mode dominate the development of coherent structures within the mixing interface (MI). The effects of channel curvature and angle between the two incoming streams on thermal mixing processes are investigated based on simulations conducted with modified planform geometries. Two additional simulations are conducted for the case where the upstream channels are parallel but not aligned with the downstream channel and for the zero-curvature case where the upstream channels are parallel and aligned with the downstream channel. The simulations highlight the influence of large-scale coherent structures within the MI and of streamwise-oriented vortical (SOV) cells on thermal mixing processes within the confluence hydrodynamics zone. Simulation results demonstrate the critical role played by the SOV cells in promoting large-scale thermal mixing for cases when such cells form in the immediate vicinity of the MI and in modifying the shape of the thermal MI within cross sections of the downstream channel—predictions consistent with empirical measurements of thermal mixing at the confluence. The set of numerical simulations reveal that the degree of thermal mixing occurring within the confluence hydrodynamic zone varies dramatically with planform geometry and incoming flow conditions. In some cases thermal mixing at the downstream end of the confluence hydrodynamic zone is limited to the MI and its immediate vicinity, whereas in others substantial thermal mixing has occurred over most of the cross-sectional area of the flow. Overall, the simulations highlight the flow conditions and the controls of these conditions that influence mixing within the immediate vicinity of a confluence.     

Alternating skimming flow over a stepped spillway

The study of stepped spillways in laboratory scales has been essentially focused on two separated sub-regimes within skimming flow. In this paper we investigate the appearance of an unclassified alternating skimming flow regime in a 0.5 m wide stepped spillway which does not fit on these earlier definitions, and which does not occur in a 0.3 m wide spillway. Our aim is to explain the genesis of this unclassified flow which is visualised in the physical stepped spillway, by using 3D numerical modelling. Flow depths and velocities are measured using an ultrasonic sensor and Bubble Image Velocimetry in the wider flume (0.5 m). The numerical model is validated with the experimental data from the 0.5 m wide spillway. After validation, the channel width of the same numerical model is reduced to 0.3 m wide spillway in order to characterise (compare) the case without (with) alternating skimming flow. Both cases are solved using Reynolds-Averaged Navier–Stokes equations together with the Volume-of-Fluid technique and SST k-\(\omega\) turbulence model. The experimental results reveal that the alternating skimming flow regime is characterised by an evident seesaw pattern of flow properties over consecutive steps. In turn, the numerical modelling clarified that this seesaw pattern is due to the presence of a complex system of cross waves along the spillway. These cross waves are also responsible for a mass and momentum exchange in the transversal direction and for the formation of the alternating skimming flow in the spillway.     

Trajectory of a jet in crossflow in a channel bend

Mixing in rivers is an important issue with many applications in water quality and water resource management. Mixing of effluents with ambient river water is especially important, particularly in river bends, where secondary circulation complicates the mixing process. By comparing measured trajectories from dye tests to velocimetry data measured with an acoustic Doppler velocimeter, this paper models the trajectory of a jet in an open channel bend using a modified formula for a jet trajectory in a straight crossflow. The original formula is shown to be insufficient for modeling the trajectory in the bend. Modifications are proposed using the position of the centre of the main secondary circulation cell to account for the bend effects. In the absence of secondary circulation, the modified formula reduces to the original formula. Once the secondary circulation has developed, the proposed formula is shown to have better residuals, lower root mean squared error, and higher \(R^2\) than the original formula.     

Formation of tidal starting-jet vortices through idealized barotropic inlets with finite length

This paper presents a surface particle image velocimetry study to investigate the dynamics of shallow starting-jet dipoles formed by tidal flow through inlets and their interaction with vorticity formed at the inlet channel lateral boundaries. Vortical structure in the flow field is identified using a local swirl strength criterion evaluated from the two-dimensional flow field. The starting jet dipole vortices and vortices formed as the lateral boundary layers are expelled during flow reversal are characterized by their trajectory, size, and circulation. Using these quantities, a model is developed to predict the size and strength of the expelled lateral boundary layer vortices based on the inlet velocity, channel length, and width of the lateral boundary layer. The expelled boundary layer vortices are found to disrupt the formation of the primary tidal jet dipole through two mechanisms. First, because the boundary layer vortices themselves form a dipole with each half of the starting-jet dipole, the starting-jet vortices are pulled apart and advected away from the inlet mouth early in the tidal cycle, resulting in a reduction in the spin-up time and the amount of vorticity input during starting-jet vortex formation. Second, the advection of each dipole away from the inlet disconnects each starting-jet vortex from the starting jet; hence, the vortices are not fed by fluid in the jet or energized by shear in the jet boundary layers. These influences of the lateral boundary layer on the starting-jet vortices’ formation and propagation are found to be a function of the channel length L, maximum velocity U, and tidal period T, resulting in a predictive value to characterize their trajectory, strength, and evolution.     

流动注射-化学发光流通池芯片及其在环境污染物测定中的应用

1　IntroductionMiniaturizationofanalyticalsystemhasbeenpaidmuchattentiontointherecentyearsforitsobviousadvantages:fastanalysis,portability,reducedconsumptionofreagentsandimprovingper formance .Suchdeviceisso calledaslaboratory on a chip[1 ],whichistheinte…     

Compound channel flow with a longitudinal transition in hydraulic roughness over the floodplains

Flows in a compound open-channel (two-stage geometry with a main channel and adjacent floodplains) with a longitudinal transition in roughness over the floodplains are experimentally investigated in an 18 m long and 3 m wide flume. Transitions from submerged dense vegetation (meadow) to emergent rigid vegetation (wood) and vice versa are modelled using plastic grass and vertical wooden cylinders. For a given roughness transition, the upstream discharge distribution between main channel and floodplain (called subsections) is also varied, keeping the total flow rate constant. The flows with a roughness transition are compared to flows with a uniformly distributed roughness over the whole length of the flume. Besides the influence of the downstream boundary condition, the longitudinal profiles of water depth are controlled by the upstream discharge distribution. The latter also strongly influences the magnitude of the lateral net mass exchanges between subsections, especially upstream from the roughness transition. Irrespective of flow conditions, the inflection point in the mean velocity profile across the mixing layer is always observed at the interface between subsections. The longitudinal velocity at the main channel/floodplain interface, denoted \(U_{int}\), appeared to be a key parameter for characterising the flows. First, the mean velocity profiles across the mixing layer, normalised using \(U_{int}\), are superimposed irrespective of downstream position, flow depth, floodplain roughness type and lateral mass transfers. However, the profiles of turbulence quantities do not coincide, indicating that the flows are not fully self-similar and that the eddy viscosity assumption is not valid in this case. Second, the depth-averaged turbulent intensities and Reynolds stresses, when scaled by the depth-averaged velocity \(U_{d,int}\) exhibit two plateau values, each related to a roughness type, meadow or wood. Lastly, the same results hold when scaling by \(U_{d,int}\) the depth-averaged lateral flux of momentum due to secondary currents. Turbulence production and magnitude of secondary currents are increased by the presence of emergent rigid elements over the floodplains. The autocorrelation functions show that the length of the coherent structures scales with the mixing layer width for all flow cases. It is suggested that coherent structures tend to a state where the magnitude of velocity fluctuations (of both horizontal vortices and secondary currents) and the spatial extension of the structures are in equilibrium.     

An ecological risk assessment model for a pulsed contaminant emission into a wetland channel flow

As a continuation of the modelling on ecological degradation and hydraulic dispersion of pollutant emission into an idealized two-dimensional free-surface wetland flow (Zeng, L., Chen, G.Q., 2009b. Ecological degradation and hydraulic dispersion of contaminant in wetland. Ecol. Model., doi:10.1016/j.ecolmodel.2009.10.024), an ecological risk assessment model for the typical case of a pulsed contaminant emission into a realistic three-dimensional wetland channel flow is presented in this paper for the fate of cross-sectional mean concentration under environmental dispersion. An environmental dispersion model for the mean concentration is devised as an extension of Taylor’s classical analysis on dispersion in fluid flows. The velocity distribution and the environmental dispersivity in the fully developed steady flow through the wetland is found and illustrated with limiting cases covering various known solutions for the porous media flow between parallel plates, flow in a shallow wetland, sweeping flow in a densely vegetated wetland, and single phase flow in a channel. Obtained by Aris’s method of moments, the environmental dispersivity is shown characterized with multi-scale asymptotic time variations with stem dominated stage, transitional stage, and width-depth-stem dominated stage. Based on the solution for the evolution of contaminant cloud in the wetland channel flow, critical length and duration of the contaminant cloud with concentration beyond given environmental standard level are concretely illustrated for typical pollutant constituents in wastewater emission. Under the same emission intensity and environmental standard, the duration of contaminant cloud in the wetland channel is revealed shorter than that in a free surface wetland, due to the lateral effect.     

Modeling the buoyant flow of heated water discharged from surface and submerged side outfalls in shallow and deep water with a cross flow

In this study, a three-dimensional model was used to numerically study the buoyant flow, along with its mixing characteristics, of heated water discharged from the surface and submerged side outfalls in shallow and deep water with a cross flow. Hydraulic experimental data were used to evaluate the applicability of the model. The simulation results agree well with the experimental results, particularly, the jet trajectories, the dimensions of the recirculating zone, and the distribution of the dimensionless excess temperature. The level of accuracy of the simulation results of the present study is nearly identical to that of the results conducted by McGuirk and Rodi (1978). If the heated water is discharged into shallow water where the momentum flux ratio and the discharge densimetric Froude number are high, the submerged discharge method is better than the surface discharge method in terms of the scale of the recirculating zone and the minimum dilution. In deep water, where the momentum flux ratio and discharge densimetric Froude number are low, however, the submerged discharge method had few advantages. In shallow water, the discharge jet is deflected by the ambient cross flow, while forcing the ambient flow to bend towards the far bank for the full depth. For a submerged discharge in shallow water, the recirculating zone is the largest in the lowest layer but becomes smaller in the upper layer. As the water depth increases, the ambient flow goes over the jet and diminishes the blocking effect, thereby decreasing the bending of the jet.     

Experimental study of recirculating flows generated by lateral shock waves in very large channels

Due to the lack of data on hydraulic-jump dynamics in very large channels, the present paper describes the main characteristics of the velocity field and turbulence in a large rectangular channel with a width of 4 m. Although a hydraulic jump is always treated as a wave that is transversal to the channel wall, in the case of this study it has a trapezoidal front shape, first starting from a point at the sidewalls and then developing downstream in an oblique manner, finally giving rise to a trapezoidal shape. The oblique wave front may be regarded as a lateral shockwave that arises from a perturbation at a certain point of the lateral wall and travels obliquely toward the centreline of the channel. The experimental work was carried out at the Coastal Engineering Laboratory of the Water Engineering and Chemistry Department of the Technical University of Bari (Italy). In addition to the hydraulic jump formation, a large recirculating flow zone starts to develop from the separating point of the lateral shock wave and a separate boundary layer occurs. Intensive measurements of the streamwise and spanwise flow velocity components along one-half width of the channel were taken using a bidimensional Acoustic Doppler Velocimeter (ADV). The water surface elevation was obtained by means of an ultrasonic profiler. Velocity vectors, transversal velocity profiles, turbulence intensities and Reynolds shear stresses were all investigated. The experimental results of the separated boundary layer were compared with numerical predictions and related work presented in literature and showed good agreement. The transversal velocity profiles indicated the presence of adverse pressure gradient zones and the law of the wall appears to govern the region around the separated boundary layer.     

Using discriminant analysis to determine the breaking criterion for an ISW propagating over a ridge

This study aims to develop methods that are capable of deciding the breaking criterion for an internal solitary wave (ISW) propagating over a submarine ridge. Laboratory experiments were conducted in a wave tank to measure ISWs propagating over a submarine ridge. The results suggest that the ISW-ridge interaction can be grouped according to three degrees of magnitude based on the blockage parameter ζ and the degree of blocking B. For classification reasons, we first present an alternative decision model for evaluating the interaction of ISWs with an underwater ridge in a two-layer system. This approach is based on a multivariate statistical method and discriminant analysis. Information obtained from the eigenvalues is used to combine different ratio measures which are defined according to every single input and output. The discriminant model effectively classifies units into distinct predefined groups. An experimental simulation is conducted to demonstrate the practical implementation of the ISW-ridge interaction. The results of the method applied in this example are statistically significant, demonstrating the effectiveness of the ISW-ridge interaction classification method.     

Advent of sheet flow in suction affected alluvial channels

The presence of suction (flow of water from channel to ground water) affects the channel hydrodynamics and increases the bed shear stress. At high bed shear stress in alluvial channels made of the non-cohesive material, sediment transport occurs as sheet flow layer of high sediment concentration. The sediment transport in the form of sheet flow has been observed in the present study when suction was applied to the non-transporting channels designed on incipient motion condition. The erosion of the channel banks contributed to the sheet flow because of the increased channel bed shear stress. An empirical relation for the thickness of sheet flow layer has been developed which includes suction as independent parameter along with others.     

Discussion of “Tang,H., Zhang,H., &amp; Yuan,S. (2018). Hydrodynamics and contaminant transport on a degraded bed at a 90-degree channel confluence. Environmental Fluid Mechanics, 18(2), 443–463”

The authors measured and computed the hydrodynamics and passive scalar dispersion in 90-degree open channel confluences over flat and degraded beds with a dominant upstream or tributary inflow. The present discussion essentially deals with the direction of rotation of the secondary currents, reported for the flat bed configuration with dominant tributary inflow. This rotation direction is indeed surprisingly opposite to the ones reported in the literature, both from calculations and measurements, even if present geometry slightly differs from literature geometries.     

Tsunami generation by submarine landslides: comparison of physical and numerical models

An idealised two-dimensional laboratory model of tsunamis generated by submarine landslides is described. The experimental configuration corresponds to the benchmark configuration suggested by other researchers in the international tsunami community. It comprises a semi-elliptical rigid landslide with a height to length ratio of 0.052 sliding down a 15° slope. The initial landslide submergence and specific gravity are varied, the second of which primarily determines the initial landslide acceleration. In these experiments the landslide motion is generally well approximated as consisting of two periods of constant acceleration. The first phase of positive acceleration finishes as the landslide reaches the base of the slope, while the second period of a slower deceleration continues until the landslide comes to rest along the horizontal base of the tank. A novel experimental technique, which utilises laser-induced fluorescence (LIF), is employed to measure the free surface displacement over the entire space and time domains. This enables the wave potential energy field to be computed directly and provides a vivid picture of the wave generation and development process. Particle tracking velocimetry provides detailed information on the landslide motion and also some data on the sub-surface velocity field. Experimental runs require multiple repeats (typically 35–50) of the same setup in order to capture the entire wave field with the desired resolution. Thus high level experimental repeatability is required, and this is demonstrated. A range of parameters relevant to hazard management are presented and discussed. Maximum crest and trough amplitudes of the offshore propagating waves are shown to be approximately proportional to the initial landslide acceleration and somewhat less strongly dependent on the initial landslide submergence. The maximum wave run-up experienced at the shoreline is shown to depend almost linearly on the magnitude of a high deceleration that occurs for a short period when the landslide nears the toe of the slope. The initial submergence and initial acceleration do not directly determine the maximum wave run-up, although for these experiments they impact indirectly on the magnitude of the deceleration. The efficiency of the energy transfer from the landslide potential energy to the wave field potential energy reaches values of up to 6% and is found to be strongly dependent on the initial submergence. However because of the link between the landslide mass and its acceleration, this efficiency is almost completely independent of the initial acceleration. The results from a numerical model based on linear, inviscid and irrotational wave theory, and solved with the boundary element method, are compared with the data from the experimental program. The numerical model accurately produces the generated sequence of wave crests and troughs, but slightly overpredicts their phase speed by between 2 and 4%. For all other parameters the numerical model predictions are within 25% of the experimental values, although this includes both under- and overprediction for the range of independent parameters covered.