Yes / Self-aeration in open channel flows occurs owing to free surface air entrainment. Self-aeration development and fully cross-sectional distribution of air concentration are not thoroughly understood. In the present study, an analytical solution for the averaged cross-sectional air concentration in the gradually varying region is established using a simplified mechanism of free surface air entrainment. For a fully cross-sectional distribution of air concentration affected by the channel bottom, a model of a diffusion region without wall restraint is proposed, and two situations are classified based on averaged cross-sectional air concentration. Good agreement between measured data and calculations is obtained, and the computational accuracy of the air concentration distribution near the wall is improved. The results reveal that the channel slope determines the air entrainment quantity, while water flow discharge determines the self-aeration evolution distance. The solutions for the averaged cross-sectional air concentration and the effect of the bottom wall on air diffusion promote air–water flow applications in hydraulic engineering practices. / This work was supported by National Natural Science Foundation of China [grant numbers 51939007, 51979183]; Sichuan Province Science and Technology support program [grant number 2019JDTD0007].

The movement of granular material along a streambed has been a challenging subject for researchers for more than a century. Predicting the limiting case of nearly zero bedload transport, usually referred to as threshold of motion or critical condition, is even more challenging due to the highly fluctuating nature of turbulent flow. Numerous works have advocated that the peak turbulent forces, randomly occurring in time and space with magnitudes higher than the average, initiate the bed material motion. More recent findings have shown that not only the magnitude of the peak turbulent forces acting on individual grains but their duration as well have to be considered for determining the incipient conditions. Their product, or impulse, is better suited for specifying such conditions. The goal of this study was to investigate the mechanism responsible for initiation of sediment motion under turbulent flow conditions. The impulse concept was investigated by utilizing appropriate measurement methods in the laboratory for determining the condition of incipient motion. The experimental program included measurements of particle entrainment rates of a mobile grain and turbulence induced forces acting upon a fixed grain for a range of flow conditions. In addition, near bed flow velocities were measured synchronously with both the entrainment and pressure measurements at turbulent resolving frequencies. Results of this work covered the limitations and uncertainties associated with the experimental methods employed, and the description of the inadequacies of existing incipient motion models via the impulse framework. The extreme sensitivity of bed material activity to minute adjustments in flow conditions was explained by the associated change in the frequency of impulse events. The probability density function proposed for impulse was used together with the critical impulse to estimate the particle entrainment rate for a range of flow conditions. It was shown that the impulse events with potential to dislodge the grain were occurring mostly during sweep type of flow structures. The impulse events were also typically accompanied by positive lift forces. The force patterns showed that the positive peaks in the lift consistently occurred before and after the impulse events in the drag force. The magnitude of these lift forces were significantly higher in the wake of a cylinder compared to that of uniform flow conditions. The time average lift force in the wake of a cylinder was also observed to be positive with magnitudes reaching more than 30% of the submerged weight of the particle. The cylinder caused the downstream turbulence intensity to increase slightly but the particle entrainment rate to increase significantly. This finding provided a physically based explanation for the modification of turbulent force fluctuations and resulting changes in the particle movement rates by such unsteady flow conditions. / Ph. D.

No / This study investigates an inclusive curriculum design based on student-centred approach. This proposed design approach has been applied to Open Channel Hydraulics module (CSE6008-A) at School of Engineering, University of Bradford, United Kingdom. This paper will introduce in step-by-step manner the full curriculum design and how the student-centred approach is being adapted in each step of the design. The required criteria will be designed based on learning outcomes design, curriculum organization, assessment strategies and student achievement evaluation. Besides, a key discussion will also be allocated for the inclusive practice that allows the vastly diverse student group to benefit from this approach, and a separate section will also be utilized to fully discuss this inclusive approach in the proposed curriculum design. This paper proposes a useful student-centred curriculum design concept, which is adaptable for different engineering modules.

It is demonstrated that the method for locating a velocity profile origin, or plane of zero velocity, by fitting log profiles to streamwise velocity measurements is applicable to a larger range of roughness scales than previously expected. Five different sets of detailed, experimental velocity measurements were analyzed encompassing sediment-scale roughness elements, roughness caused by rigid vegetation, and large-scale roughness elements comprised of mobile bedforms. The method resulted in similar values of normalized zero-plane displacement for all roughness types considered. The ratios of zero-plane displacement, dh, to roughness height, ks, were 0.20 and 0.26 for the sediment- and vegetation-scale experiments, respectively. The results for the two experiments with bedform dominated roughness were 0.34 and 0.41. An estimate of dh/ks ranging from 0.2 to 0.4 is therefore recommended for a range of roughness types with the higher end of the range being more appropriate for the larger, bedform-scale roughness elements, and the lower end for the sediment-scale roughness elements. In addition, it is demonstrated that the location of the plane of zero velocity is temporally constant even when the bed height is not. The effects of roughness element packing density were also examined with the identification of a possible threshold at 4%, above which zero-plane displacement is independent of packing density. The findings can be applied to field velocity measurements under mobile bed conditions, facilitating the calculation of turbulence parameters such as shear velocity, by using point measurements and providing guidelines for the estimation of an appropriate value for zero-plane displacement. / Master of Science

Yes / A two-dimensional steady potential flow theory is applied to calculate the flow in an open channel with complex solid boundaries. The boundary integral equations for the problem under investigation are first derived in an auxiliary plane by taking the Cauchy integral principal values. To overcome the difficulties of a nonlinear curvilinear solid boundary character and free water surface not being known a priori, the boundary integral equations are transformed to the physical plane by substituting the integral variables. As such, the proposed approach has the following advantages: (1) the angle of the curvilinear solid boundary as well as the location of free water surface (initially assumed) is a known function of coordinates in physical plane; and (2) the meshes can be flexibly assigned on the solid and free water surface boundaries along which the integration is performed. This avoids the difficulty of the traditional potential flow theory, which seeks a function to conformally map the geometry in physical plane onto an auxiliary plane. Furthermore, rough bed friction-induced energy loss is estimated using the Darcy-Weisbach equation and is solved together with the boundary integral equations using the proposed iterative method. The method has no stringent requirement for initial free-water surface position, while traditional potential flow methods usually have strict requirement for the initial free-surface profiles to ensure that the numerical computation is stable and convergent. Several typical open-channel flows have been calculated with high accuracy and limited computational time, indicating that the proposed method has general suitability for open-channel flows with complex geometry.

No / This paper presents an integration procedure for the Reynolds-averaged Navier-Stokes equations for the determination of the distribution of the streamwise velocity using the vertical component. This procedure is dedicated to the outer region and central part of channels. The proposed model is applicable to both rough and smooth flow regimes, provided the velocity at the inner-outer boundary has been properly defined. To generate a simplified expansion, a number of hypotheses are proposed, focusing in particular on the analytical modeling of the vertical component by adopting a negligible viscosity. The proposed hypotheses are validated by the experimental data existing in the literature. The proposed simplified expansion is studied through a sensitivity analysis and proved consistent in regards to model experimental data. The proposed model seems capable of demonstrating different kinds of flows, including dip phenomenon flow patterns.

Yes / The Prandtl second kind of secondary current occurs in any narrow channel flow causing velocity dip in the flow velocity distribution by introducing the anisotropic turbulence into the flow. Here, a study was conducted to explain the occurrence of the secondary current in the outer region of flow velocity distribution using a universal expression. Started from the basic Navier-Stokes equation, the velocity profile derivation was accomplished in a universal way for both smooth and rough open channel flows. However, the outcome of the derived theoretical equation shows that the smooth and rough bed flows give different boundary conditions due to the different formation of log law for smooth and rough bed cases in the inner region of velocity distribution. Detailed comparison with a wide range of different measurement results from literatures (from smooth, rough and field measured data) evidences the capability of the proposed law to represent flow under all bed roughness conditions.

Yes / Aquatic vegetation can significantly affect flow structure, sediment transport, bed scour and water quality in rivers, lakes, reservoirs and open channels. In this study, the lattice Boltzmann method is applied for performing the two dimensional numerical simulation of the flow structure in a flume with rigid vegetation. A multi-relaxation time model is applied to improve the stability of the numerical scheme for flow with high Reynolds number. The vegetation induced drag force is added in lattice Boltzmann equation model with the algorithm of multi-relaxation time in order to improve the simulation accuracy,. Numerical simulations are performed for a wide range of flow and vegetation conditions and are validated by comparing with the laboratory experiments. Analysis of the simulated and experimentally measured flow field shows that the numerical simulation can satisfactorily reproduce the laboratory experiments, indicating that the proposed lattice Boltzmann model has high accuracy for simulating flow-vegetation interaction in open channel. / National Natural Science Foundation of China (grant number: 11861003 and 11761005)

Yes / Bridge pier is a common feature in hydraulic structure. Its impact to the river usually occurs in group form rather than single pier, so this challenging piers-group influence towards river hydraulics and turbulence needs to be explored. In this paper, the measurements were conducted using an Acoustic doppler velocimeter (ADV) to study velocities in three dimensions (longitudinal, transversal, and vertical). Based on the experimental data, we have observed reversed depth-averaged velocity vector after each pier in the group of three-pier. The analysis has been conducted on the contribution of each bursting event to Reynolds shear stress (RSS) generation, in order to identify the critical events and turbulence structures around the piers. In the upstream near-wake flow in the bed-wall layer, strong sweep and ejection events have been observed; while at downstream, sweeps were more dominant. The pattern of burst changed in the outer layer of flow, where ejections were more dominant. Furthermore, the contribution fractional ratio to RSS variation at hole size H = 0 indicates that sweeps and ejections were significantly generated at the near wake-flow in upstream.

Yes / This paper aims to investigate the mean velocities and turbulence characteristics, which include Reynolds shear stresses (RSS) and turbulent intensities, in order to study temporal and spatial flow patterns around a complex three square-piers group. Measurements have been performed by acoustic Doppler velocimeter (ADV) to identify turbulence flow structures alteration and vortices formation in longitudinal, transverse, and vertical directions. In total, 164 locations around the bridge piers have been measured, while each location consists of no lesser than 15 vertical depth ADV data point measurements. To investigate the flow turbulence, critical locations have been observed from high vortex-influenced regions. The impacted mean velocity profiles have been compared to that at the unaffected region to study the flow alteration caused by piers group. The result reveals that the log-law layer near wall has been enhanced by the installation of pier. Furthermore, as identified from the most critical locations of turbulent structures across different measurement rows, the accelerated flow will result in lower observed turbulent structures. Although flow after each pier has been influenced by strong reversal velocities, vortices have managed to form at near wake. Analysis has been divided into inner and outer flow regions (IFR and OFR) based on flow depth to enhance the understanding of how bridge piers contribute to the development of flow turbulence. / Shuyan Zang acknowledges the support of Chinese National Natural Science Foundation (Grant No. 22308222) to this study.