Experimental surface flow visualization and numerical investigation of flow structure around an oblong stockpile

Emission factors are largely used to quantify particle emissions from industrial open storage piles. These factors are based on the knowledge of velocity distribution and flow patterns over the stockpile surface which still requires further research. The aim of the present work is to investigate the airflow characteristics over a single typical oblong pile and in its near-ground surroundings for various wind flow directions. Wind tunnel experiments using an oil-film surface coating technique were carried out for near-wall flow visualization. Numerical simulation results, favorably compared to PIV measurements, were used to allow comparison analysis of flow features. For the stockpile oriented 90° to the wind main direction, typical topology of flow around wall-mounted obstacles were observed, notably a wake zone downstream the pile including two main counter-rotating vortices. Further analysis of numerical wall shear stress distribution and streamlines indicates that two complex three-dimensional vortical flow structures develop downstream the pile. For other incoming wind flow directions (30 and 60°), the flow characteristics over the storage pile greatly differ as a single helical main vortex develops from the pile’s crest. Corresponding high values of wall shear stress are noticed downstream the storage pile. For each configuration studied, downwash and upwash zones are induced by the vortical structures developed. This near-wall flow topology combined with areas of high friction levels may be linked to potential dust emission from the ground surface surrounding industrial stockpiles.     

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.     

Analysis of the velocity field in a large rectangular channel with lateral shockwave

In this work the authors describe the main characteristics of the velocity field of hydraulic jumps in a very large channel where lateral shockwaves occur. Experiments were carried out at the Coastal Engineering Laboratory of the Water Engineering and Chemistry Department of the Technical University of Bari (Italy). Extensive flow velocity measurements were investigated in order to have a clearer understanding of both hydraulic jump development and lateral shockwave formation in a very large channel. Eight experiments were performed in a 4m wide rectangular channel; the experiments differed in the inlet Froude number F ₀ and the jump type. Seven tests were carried out with undular jumps and one with a roller jump. The flow velocity and the flow free surface measurements were taken using a two-dimensional Acoustic Doppler Velocimeter (ADV) and an ultrasonic profiler, respectively. The experimental results can be summarized as follow: (i) the formation of well developed lateral shockwaves similar to those of oblique jumps were observed; (ii) the comparison of the experimental and theoretical data shows that the classic shockwave theory is sufficiently confirmed in the analyzed range of Reynolds number, taking into account the experimental errors and the difference between the theoretical and experimental assumptions; (iii) the transversal flow velocity profiles in the recirculating zone show a good agreement with the numerical simulations presented in literature in the case of a separated turbulent boundary layer over a flat plate. This conclusion enables us to confirm the hypothesis that the lateral shockwaves in the channel are the result of a boundary layer which, as observed, forms on the channel sidewalls.     

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.     

Hydrodynamics and contaminant transport on a degraded bed at a 90-degree channel confluence

Channel confluences at which two channels merge have an important effect on momentum exchange and contaminant diffusion in both natural rivers and artificial canals. In this study, a three-dimensional numerical model, which is based on the Reynolds Averaged Navier–Stokes equations and Reynolds Stress Turbulence model, is applied to simulate and compare flow patterns and contaminant transport processes for different bed morphologies. The results clearly show that the distribution of contaminant concentrations is mainly controlled by the shear layer and two counter-rotating helical cells, which in turn are affected by the discharge ratio and the bed morphology. As the discharge ratio increases, the shear flow moves to the outer bank and the counter-clockwise tributary helical cell caused by flow deflection is enlarged, leading the mixing happens near the outer bank and the mixing layer distorted. The bed morphology can induce shrinkage of the separation zone and increase of the clockwise main channel helical cell, which is initiated by the interaction between the tributary helical cell and the main channel flow and strengthened by the deep scour hole. The bed morphology can also affect the distortion direction of the mixing layer. Both a large discharge ratio and the bed morphology could lead to an increase in mixing intensity.     

Experimental study of a positive surge. Part 2: comparison with literature theories and unsteady flow field analysis

A positive surge is a unsteady open channel flow resulting from the rapid rise of the free-surface. The phenomenon may be observed in water supply canals and channels as well as in some estuaries during spring tidal conditions. The formation and development of positive surges can be predicted using the method of characteristics and shallow water equations. The paper is the second part of a study presenting 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. A comparison between main features of the undular surges with literature theories demonstrated that the experimental data were mostly in agreement with Andersen’s theory. The analysis of unsteady flow field including Reynolds stresses confirmed and extended previous findings about positive surge hydrodynamics.     

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.     

Front velocity and structure of bottom gravity currents with a low volume of release propagating in a porous medium

The paper reports results of large eddy simulations of lock exchange compositional gravity currents with a low volume of release advancing in a horizontal, long channel. The channel contains an array of spanwise-oriented square cylinders. The cylinders are uniformly distributed within the whole channel. The flow past the individual cylinders is resolved by the numerical simulation. The paper discusses how the structure and evolution of the current change with the main geometrical parameters of the flow (e.g., solid volume fraction, ratio between the initial height of the region containing lock fluid and the channel depth, ratio between the initial length and height of the region containing lock fluid) and the Reynolds number. Though in all cases with a sufficiently large solid volume fraction the current transitions to a drag-dominated regime, the value of the power law coefficient, α, describing the front position’s variation with time (x _f  ~ t ^α , where t is the time measured from the removal of the lock gate) is different between full depth cases and partial depth cases. The paper also discusses how large eddy simulation (LES) results compare with findings based on shallow-water equations. In particular, LES results show that the values of α are not always equal to values predicted by shallow water theory for the limiting cases where the current height is comparable, or much smaller, than the channel depth.     

LE of shallow mixing interfaces: A review

Eddy-resolving techniques have become a powerful tool to investigate shallow flows at both laboratory and field scale. In this paper several examples are given where high-resolution 3D numerical simulation are used to investigate the spatial development of mixing interfaces (MIs) forming in shallow environments like open channels with idealized and natural bathymetry where the bed friction plays a major role in the spatial development of the MI and associated large-scale turbulence. The focus is on the coherent structures forming within the MI and in its vicinity that control the momentum and mass exchange and heat transfer between the two sides of the MI. Examples include: (1) a MI developing in a flat-bed open channel downstream of a splitter wall separating two parallel fully-turbulent streams of different velocities, (2) a MI developing in a flat-bed open channel downstream of a 60 \(^{\circ }\) wedge separating two non-parallel fully turbulent streams of different velocities, (3) a MI developing downstream of a river confluence for cases with a large and, respectively, a small difference between the mean velocities of the two streams. Stratification effects due to unequal densities of the two incoming streams are also discussed, (4) a MI developing between a main rectangular straight channel and a series of shallow embayments present at one of the channel banks. Besides using available experimental data to demonstrate that eddy resolving techniques can accurately predict the structure of the MI and its development, the paper discusses new insights into the physics of these flows obtained based on the simulations. The paper also provides an overview of the main numerical approaches that can be used to simulate the unsteady dynamics of the large scale turbulence in flows containing shallow MIs.     

Turbulence and aeration in hydraulic jumps: free-surface fluctuation and integral turbulent scale measurements

In an open channel, a change from a supercritical to subcritical flow is a strong dissipative process called a hydraulic jump. Herein some new measurements of free-surface fluctuations of the impingement perimeter and integral turbulent time and length scales in the roller are presented with a focus on turbulence in hydraulic jumps with a marked roller. The observations highlighted the fluctuating nature of the impingement perimeter in terms of both longitudinal and transverse locations. The results showed further the close link between the production and detachment of large eddies in jump shear layer, and the longitudinal fluctuations of the jump toe. They highlighted the importance of the impingement perimeter as the origin of the developing shear layer and a source of vorticity. The air–water flow measurements emphasised the intense flow aeration. The turbulent velocity distributions presented a shape similar to a wall jet solution with a marked shear layer downstream of the impingement point. The integral turbulent length scale distributions exhibited a monotonic increase with increasing vertical elevation within 0.2 < L_z/d₁ < 0.8 in the shear layer, where L_z is the integral turbulent length scale and d₁ the inflow depth, while the integral turbulent time scales were about two orders of magnitude smaller than the period of impingement position longitudinal oscillations.     

One-dimensional modelling of sediment deposition in reservoirs with sinking streams

It is common in karst regions that rivers are occasionally cut by mountains, resulting in the alternate appearances of open channel flow and pressurized flow. With more and more reservoirs being built in this region, the complicated sediment transport processes of such mixture flow are urgently needed to study. In this paper a one-dimensional numerical model with free surface and pressurized flows coupled together is presented. The simulated results are analyzed to explore deposition process in reservoirs with sinking streams; impacts of different hydraulic conditions on the sedimentation are also studied. To verify the computed results, a flume experiment is also conducted. The results show that deposition of sediments mainly occurs in open channel reaches where the longitudinal profile is similar to that of conventional reservoirs, i.e. typical delta has formed, demonstrating characteristics of deltaic deposition morphology in that the crest of delta moves to the downstream direction over time. The model provided by this paper is not only proved to simulate the characteristics of deposition in karst reservoirs successfully, but also reveal the impacts of hydraulic conditions in such circumstances.     

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.     

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.     

Experimental investigation of bubbly flow and turbulence in hydraulic jumps

Many environmental problems are linked to multiphase flows encompassing ecological issues, chemical processes and mixing or diffusion, with applications in different engineering fields. The transition from a supercritical flow to a subcritical motion constitutes a hydraulic jump. This flow regime is characterised by strong interactions between turbulence, free surface and air–water mixing. Although a hydraulic jump contributes to some dissipation of the flow kinetic energy, it is also associated with increases of turbulent shear stresses and the development of turbulent eddies with implications in terms of scour, erosion and sediment transport. Despite a number of experimental, theoretical and numerical studies, there is a lack of knowledge concerning the physical mechanisms involved in the diffusion and air–water mixing processes within hydraulic jumps, as well as on the interaction between the free-surface and turbulence. New experimental investigations were undertaken in hydraulic jumps with Froude numbers up to Fr = 8.3. Two-phase flow measurements were performed with phase-detection conductivity probes. Basic results related to the distributions of void fraction, bubble frequency and mean bubble chord length are presented. New developments are discussed for the interfacial bubble velocities and their fluctuations, characterizing the turbulence level and integral time scales of turbulence representing a “lifetime” of the longitudinal bubbly flow structures. The analyses show good agreement with previous studies in terms of the vertical profiles of void fraction, bubble frequency and mean bubble chord length. The dimensionless distributions of interfacial velocities compared favourably with wall-jet equations. Measurements showed high turbulence levels. Turbulence time scales were found to be dependent on the distance downstream of the toe as well as on the distance to the bottom showing the importance of the lower (channel bed) and upper (free surface) boundary conditions on the turbulence structure.     

Hydrodynamics and turbulence in emergent and sparsely vegetated open channel flow

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

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.     

Structural design of a natural windbreak using computational and experimental modeling

Computational and experimental approaches are employed for the structural design of a natural windbreak. It is intended to find the optimum tree shelterbelt to obviate the uneven wind speed distribution across the width dimension of a high-level competition rowing channel. The experimental results, obtained in a wind tunnel, and consisting of erosion-technique images and local wind-speed measurements, are used to benchmark the computational model. A good agreement between the two sets of results is obtained. Several windbreak configurations, considering one or two rows of different cross-sectional shape and porosity, are computationally modeled. For the shortest row a rectangular shape, with porosity of 35%, is considered; for the tallest row, which aims the modeling of a poplar tree, a porosity of 87% is assumed at the trunk level, and 60% at the crown. The optimum shelterbelt consists of two rows, composed by bamboo and poplar trees, which allows the attainment of a low and nearly uniform wind flow across the width of the channel.     

Evaluation of numerical simulations of CO2 transport in a city block with field measurements

Studying urban air-transport phenomena is highly complex, because of the heterogenous flow patterns that can arise. The main reason for these is the variable topology of urban areas, however, there is a large number of influencing variables such as meteorological conditions (e.g., wind situation, temperature) and anthropogenic factors such as traffic emissions. During a one-year CO₂ measurement campaign in the city of Basel, Switzerland, steep CO₂ gradients were measured around a large building. The concentration differences showed a strong dependency on the local flow regimes. Analysis of the field data alone did not provide a complete explanation for the mechanisms underlying the observed phenomena. The key numerical parameters were defined and the influence of turbulent kinetic energy dependency on the time interval for the Reynolds decomposition was studied. A Reynolds-Average Navier-Stokes Computational Fluid Dynamics (CFD) approach was applied in the study area and the CO₂ concentrations were simulated for six significant meteorological situations and compared to the measured data. Two flow regimes dependent on the wind situation, which either enhanced or suppressed the concentration of CO₂ in the street canyon, were identified. The enhancement of CO₂ in the street canyon led to a large difference in CO₂ concentration between the backyard- and street-sides of a building forming the one wall of the canyon. The specific characteristics of the flow patterns led to the identification of the processes determining the observed differences in CO₂ concentrations. The combined analysis of measurement and modeling showed the importance of reliable field measurements and CFD simulations with a high spatial resolution to assess transport mechanisms in urban areas.