Investigation of velocity distribution in open channel flows based on conditional average of turbulent structures

We report the development of a new analytical model similar to the Reynolds-averaged Navier-Stokes equations to determine the distribution of streamwise velocity by considering the bursting phenomenon. It is found that, in two-dimensional (2D) flows, the underlying mechanism of the wake law in 2D uniform flow is actually a result of up/down events. A special experiment was conducted to examine the newly derived analytical model, and good agreement is achieved between the experimental data in the inner region and the model's prediction. The obtained experimental data were also used to examine the DML-Law (dip-modifiedlog-law), MLW-Law (modified-log-wake law), and CML-Law (Cole's wake law), and the agreement is not very satisfactory in the outer region. Disciplines Engineering | Science and Technology Studies Publication Details Han, Y., Yang, S., Sivakumar, M. & Qiu, L. (2017). Investigation of velocity distribution in open channel flows based on conditional average of turbulent structures. Mathematical Problems in Engineering, 2017 1458591-1-1458591-9. This journal article is available at Research Online: http://ro.uow.edu.au/eispapers1/458 Research Article Investigation of Velocity Distribution in Open Channel Flows Based on Conditional Average of Turbulent Structures Yu Han, Shu-Qing Yang, Muttucumaru Sivakumar, and Liu-Chao Qiu 1College of Water Resources & Civil Engineering, China Agricultural University, Beijing 100083, China 2School of Civil, Mining & Environmental Engineering, University of Wollongong, NSW 2522, Australia Correspondence should be addressed to Yu Han; yh916@uowmail.edu.au and Liu-Chao Qiu; qiuliuchao@cau.edu.cn Received 4 January 2017; Accepted 24 April 2017; Published 29 May 2017 Academic Editor: Jian G. Zhou Copyright © 2017 Yu Han et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. We report the development of a new analytical model similar to the Reynolds-averaged Navier-Stokes equations to determine the distribution of streamwise velocity by considering the bursting phenomenon. It is found that, in two-dimensional (2D) flows, the underlying mechanism of the wake law in 2D uniform flow is actually a result of up/down events. A special experiment was conducted to examine the newly derived analytical model, and good agreement is achieved between the experimental data in the inner region and the model’s prediction. The obtained experimental data were also used to examine the DML-Law (dip-modifiedlog-law), MLW-Law (modified-log-wake law), and CML-Law (Cole’s wake law), and the agreement is not very satisfactory in the outer region.


Introduction
Velocity distribution in open channel flows offers a wide range of applications in the fields of hydrometry, sediment transport, river restoration, power plant design, and so forth [1][2][3].Accurate investigation of the velocity distribution in open channel flows has been conducted over the past century.The law of the wall, developed for pipe flows, has often been applied to uniform flows in open channel flows [4]; this "universal wall function" or log-law of the wall is expressed as and herein  is the von Karman constant,  * is the mean friction velocity, and  is the vertical direction.In (1), the boundary condition is  = 0 at  =  0 [5], and the parameter  0 is the constant of integration that must be determined using the boundary condition [6].It is expected that the parameter  0 in (1) should be a variable as well [7], for a smooth boundary  0 should be related to the thickness of the local viscous sublayer; that is, where  is the coefficient to be determined experimentally, ] is the kinematic viscosity, and  *  denotes the local friction velocity, which may be different from  * .The velocity vertical profile is well described by the classical log-law in the inner region, /ℎ < 0.2, but the loglaw normally deviates from the experimental data in the outer ration, /ℎ > 0.2 [8].Since this method is unable to represent the outer region of a velocity profile, Coles [9] suggested an improvement, called the wake law.Many researchers have discussed the validity of such a law [10,11] and have improved both the log and wake laws for smooth or rough flows [12,13].
Coles [9] extended the log-law by introducing a purely empirical correction function.He developed an additional term (i.e., the wake term) as a supplement to the log-law, to account for the velocity deviation caused by wake.This model is known as the wake law: where Π denotes the "profile" parameter; Coles [9] gave Π = 0.55.For open channel flows, Coleman [14] obtained an average Π value of 0.19.Nezu and Rodi [15] yielded a Π value of 0 to 0.20, while Kirkgöz [16] reported a value of 0.1.Cardoso et al. [10] obtained a Π value of −0.077 over a smooth bed.Kironoto and Graf [11] found Π value in the range of −0.08 to 0.15 for water flows over a gravel bed.However, the values of Π are established by experimentation, with large variations reported in the literatures.While significant advances have been made by using Cole's wake law, the mechanism of Cole's wake law (CWL-Law) and the associated wake strength parameter are not fully understood.These investigations often ascribe the large variations to either smooth or rough boundaries and few researchers have taken into account the velocity-dip phenomenon caused by secondary currents [17] (Guo 2015).Different from the positive deviation of measured velocity from the log-law's prediction or the CWL-Law, Yang et al. [18] discussed the negative deviation of measured velocity from the log-law's prediction or the dip phenomenon, in which the maximum longitudinal velocity occurs below the water surface.They suggested that Cole's wake law is not able to describe the entire velocity profile when the dip phenomenon exists.Yang et al. [18] proposed a dip-modified-log-law (DML-Law) based on the analysis of the Reynolds-averaged Navier-Stokes (RANS) equations.This law, involving two logarithmic distances, one from the bed (i.e., the log-law) and the other from the free surface, has the advantage that it contains only one parameter  for dip-correction.The DML-Law reverts to the classical log-law for  = 0. Yang et al. [18] modified the log-law by adding a term to express the dip phenomenon instead of Cole's wake law based on Reynolds equations: Guo and Julien [12] proposed a modified-log-wake law (MLW-Law) which fits velocity profiles with a dip phenomenon.However, this law cannot be used for predictive applications since it requires fitting the near-free-surface velocities to the parabolic law to obtain dip position and maximum velocity.As Guo and Julien indicate, the MLW-Law can be used only in flow measurements since it requires measured velocities; that is, where  = / and  is the distance from the bed to the point where the velocity is maximum.Nezu and Nakagawa [19] found that in the outer region the mean velocity data deviate systematically from the loglaw distribution; at sufficiently high Reynolds numbers, this deviation cannot be neglected in the free-surface region; that is, /ℎ > 0.6.It is thus worth investigating the underlying mechanism of velocity deviation from the log-law.Bonakdari et al. [20] also proposed a new formulation of the vertical velocity profile in open channel flows based on an analysis of the Navier-Stoke equations.On the basis of this law and of the suggestions of Absi [21], a new analytical procedure was presented by Pu et al. [22].
With the advent of nonintrusive instrumentation, visualisation techniques reveal that turbulence is dominated by forceful structures that are well organised and ordered eddies with a certain lifespan, now termed as coherent structures.Coherent structures can be placed into two categories: bursting phenomena that occur in the wall region and large-scale vertical motion in the outer region.For coherent structural analysis, the method of quadrants has been developed based on measured velocity fluctuations, that is, quadrant I,   > 0 and V  > 0; quadrant II,   < 0 and V  > 0; quadrant III,   < 0 and V  < 0; and quadrant IV,   > 0 and V  < 0 [23].Cellino and Lemmin [24] pointed out that it would be meaningful to concentrate on V  > 0 and V  < 0 in quadrants because the most significant feature of turbulence is fluctuations in the wall-normal direction.They also developed an approach to calculate conditionally averaged upward/downward velocity based on velocity fluctuations.However, the linkage of quadrant analysis and upward/downward events with velocity deviation from the log-law has not been established.
The objective of this research is to present a new analytical approach to determine the streamwise velocity profile and establish a relationship between the log wake function and up/down events.This leads to the present research aim to develop a universal model to express the velocity profile in uniform flows.Such a novel analytical approach will be used to explain that streamwise velocity profiles in the outer region depart from the universal logarithmic law of the wall.The proposed analytical approach and those previously proposed mathematical methods will be examined by using experimental data.

Theoretical Consideration
The RANS momentum and continuity equations can be written for steady state flow as where , V, and  denote the mean velocity in the streamwise (), lateral (), and vertical () directions, respectively,  is the gravitational acceleration,  is the energy slope, and −    are the Reynolds stress tensor components.For a uniform and fully developed flow, (7) becomes Alternatively, (9) becomes where  is the density of the fluid,  is the dynamic viscosity of the fluid, and   , V  , and   are the turbulent velocity fluctuation in the , , and  directions.For 2D flow, ( 10) can be simplified as Integration of ( 13) yields where  is the integration constant.Using the boundary condition,  = 0,  = V = 0,   V  = 0, and ]/ =  2 * , and thus  =  2 * .For water surface,  = ℎ, V = 0, and   V  = 0 and then the shear velocity can be rewritten as Therefore, the momentum equation has the following form by dividing (15): where  + = / * and V + = V/ * .Cellino and Lemmin [24] realised that it would be meaningful to develop an approach to calculate the distribution of streamwise velocity by considering the bursting phenomenon.However, the above time averaged method cannot be used to express the bursting phenomenon.Yang [25] proposed an approach to calculate the conditionally averaged upward/downward velocity based on velocity fluctuations.He first defines the occurrence probabilities  based on the direction of vertical velocity.For upflow (subscript ) and downflow (subscript ),  is expressed as where Δ is the time interval between two consecutive velocity acquisitions in an experiment;  is the total observed time period; that is, ∑ Δ  + ∑ Δ  =   +   = .The upflow and downflow velocity definitions are given in the following way: where   and   denote upflow and downflow events, respectively, and the symbol "∼" represents instantaneous velocities.Yang [25] obtained the Reynolds shear stress for upflow and downflow as follows, and both deviate from the standard linear distribution; that is, where    = ũ −   ,    = ũ −   ,  +  =   / *  , and  +  =   /  *  ; ℎ is water depth; and  *  and  *  are friction velocities for upflow (subscript ) and downflow (subscript ), respectively.
The expression of eddy viscosity where  = /ℎ and ( 24) can be expressed as The velocity gradient in upflow and downflow can then be expressed by Similarly, the gradient of mean velocity can be expressed by Inserting (26) into (27) yields Here, we denote Then Integration of (31) with respect to  yields In nonslip boundary conditions  = 0 and then  =  0 =  0 /ℎ and thus  1 = −(1/)ln( 0 /ℎ) Finally, we obtained the corrected log equation as The value of  used in (30) corrects the deviation of velocity from the log-law by considering upward and downward motion adequately.Thus, we can find the existence of the wake function from the newly derived log-law.Upward and downward motion can decompose turbulent flow into two distinct events, and it also indicates that if  +  =  +  , then the classic log-law can be obtained.Many researchers determine the different  and  parameters in  + = 1/ ln  * /]+ from their experiments.Equation (33) implies that the main reason for such different results is not from different experimental conditions, but it is rather from upflow and downflow effects.It is necessary to verify (33) using experimental data.

Experimental Set-Up
The experiments were undertaken in a large-scale flow loop at the University of Wollongong, Australia, in which water was supplied from a head tank.The main components of this flume were the head tanks tail tank, glass water channel, recirculation pipe system, and two pumps as shown in Figure 1.
The head tank and the tail tank were constructed using stainless steel.The head tank was aligned centrally to the flume and symmetrical about the centre line of the channel so that the flow at the entrance of the channel was as uniform as possible.A honeycomb section was also located at the entrance to the flume.Through the addition of this honeycomb, the velocity distribution within the channel cross-section became increasingly uniform.The water depth through the flume was controlled by an adjustable tailgate in the downstream tail tank.The relative height of this tailgate was adjusted by rotating the attached handle.This was installed at the entrance of the tail tank.From the data recorded at these locations, it was possible to visually analyse the development of uniform flow in the channel prior to experiments.
The experimental flume was 10.5 m long, 0.3 m wide, and 0.45 m high with a slope set to a certain value for different experimental conditions.Experimental conditions are given in Table 1.A series of velocity profiles were measured in central line with different /ℎ values.The flow depth ℎ was changed in order to examine the effect of aspect ratio in this study.Hence, a comprehensive set of high quality data, such as velocities, turbulence intensities, Reynolds stresses, and boundary shear stresses of the 3D turbulence structure, were also obtained.
Throughout the experimental process in this study the point velocities were measured using a Dantec twocomponent LDA system.It is a two-component configuration with a 60 mm optical fibre probe and a front lens with a 400 mm focal length.The system consists of a 300 mW continuous wave Argon-Ion laser and transmitting optics including a beam splitter Bragg-cell and signal processors.Green (514.5 nm) and blue (488 nm) components were used to measure the horizontal and vertical components of the velocity, respectively.The optics for 2D LDA was arranged in off-axis backscatter mode.The backscatter mode of receiving a signal from a particle refers to the reflected light being captured by the receiving optics on the same side as the incident laser beams, with reflected light going back into the prove head through the front lens.The dimension of the measurement volume was approximately 0.189 × 0.189 × 3.97 mm for both colours when measured in air.
Hence, no calibration for LDA was necessary prior to the experiments.It is important to note that, during the measurement of velocity  and V by LDA, data acquisition of the components needed to be coincidental in time for all the data pairs () and V().The signal processor is a Dantec Burst Spectrum Analyser (BSA) connected to an oscilloscope and a PC.The BSA converts the electrical signals, which are processed by the oscilloscope, into velocity data, which again are monitored online by a PC.Raw data exported from BSA can be processed into statistical values like mean velocities, root mean squared values of velocity fluctuations, and Reynolds stresses.

Results and Discussion
The velocity profiles for a smooth rectangular channel are plotted in Figures 2-7, and the velocity profiles predicted by DML-Law, MLW-Law, and CML-Law are used for comparison to experimental data.These figures show that the log-law is able to predict all experimental data in the inner region.It can be clearly found that the wake function indeed exists in uniform flow.However, DML-Law, MLW-Law, and CML-Law are not able to agree with the experimental data properly in the outer region; that is, /ℎ > 0.6.It is observed that the measured values for all runs are relatively higher than the calculated results from DML-Law in the outer region.It is impossible to improve the velocity profiles predicted by the DML-Law by adjusting parameter , so further improvement of ( 5) is recommended.The DML-Law seems to be able to predict the velocity-dip phenomenon, but with inaccurate maximum velocities.Accurate predictions of the velocity-dip phenomenon therefore require larger values of .
It can be seen from Figures 2-4 that the MLW-Law and CML-Law significantly overestimate the velocity profile in the outer region.Both of these two methods seem not to be a suitable predictor for channels with /ℎ > 3 in the outer region.In particular, the CML-Law is important for an inadequate prediction of the dip phenomenon with /ℎ > 3, and Π = 0.55 seems to be too small to improve predictions.Figures 5-7 show that the MLW-Law and CML-Law are more accurate in predicting the velocity profiles for /ℎ < 3. The effect of each term in the MLW-Law and CML-Law is not significant.In general, the average differences between the MLW-Law, CML-Law, and DML-Law theories and the current measurements give relatively large averaged  values of 7.1%, 5.9%, and 5.7%, respectively.The average differences between measurements and the current approach give a relatively small average  value of 1.0%.
Figures 2-7 present comparisons of predicted velocity profiles obtained by (33) with experimental data, and we note that the parameter  in (33) is an important value in improving the predicted velocity profiles.Thus, to plot the velocity distribution in the entire flow region using the proposed theoretical equation, it was necessary to determine the value of  in (30) accurately.The values of   and V +  in (30) can be easily calculated using (17) and (21), respectively, from a given data set of upflow velocities.We also note that  * in (30) is the mean shear stress; that is,  * = √.Furthermore, in a plot of measured upflow/downflow velocity against log  for a profile along the central line, the unknowns  *  and  *  can be evaluated from its slope, respectively.The values of   , V +  ,  * ,  *  , and  *  used in (30) are then calculated and included in Tables 2 and 3.

Mathematical Problems in Engineering
The experimental data of   and   scale with  *  and  *  , and then  +  and  +  can be calculated, respectively.In (30), the ratio of  +  and  +  is larger than 1 in the outer region, as shown in Figure 8; thus, the values of  are all negative in the outer region.Hence, the wake function (/) ln(/ℎ/(1 − /ℎ)), in (33) would also be a negative value in the outer region, implying that the velocity from (33) must be lower than that predicted by the log-law, or there is a velocity   dip near the water surface.Then, it is clearly indicated that downward motions will mainly cause the wake strength.In other words, in 2D flows there indeed exists a wake function, which is caused by variations in the wall-normal velocity.These different  values have been applied in the proposed theoretical equation and have amounted to approximately 1.5% errors from the velocity encountered in the experiments.They are much smaller than the corresponding errors for other comparisons.It will be seen that the value of  used in (30) adequately accounts for the errors in comparison.Similarly,  = ( +  −  +  ) can be determined from experimental data.The above analysis applies to (30) in which all parameters are found, and thus the relationship between  and  can then be plotted in Figure 9.The dimensionless profiles of  collapse reasonably well between the experiments, indicating that the relationship between  and  for different water depths are similar.Therefore, an empirical equation can be fitted using data  and , as shown in Figure 9.
The parameter  in (33) is an important value for improving the predicted velocity profiles.To simplify the calculation of  in practice, Figure 8 depicts the fitting of a regression equation.For a particular channel (constant aspect ratio), it is easier to check the value of  +  / +  from Figure 8 to  obtain the value of ; thus facilitating obtaining the value of  in Figure 9.

Conclusions
In this study, we theoretically investigated conditionally averaged turbulent structures for upward and downward motions in a 2D uniform flow.The connection between wake law and up/down events was also explored.The following conclusions can be drawn from the foregoing analysis: (1) In 2D flows, the wake function is caused by the additional momentum flux.By considering up/down events, it explores the underlying mechanism of Coles' wake law and establishes a relationship between the wake function and up/down events.
(2) Applying the proposed theoretical equation to experimental measurements, good agreement is achieved between the theoretical equation and the measurements.The theoretical results show that the momentum flux caused by the conditional up/down events leads to deviation of velocity from the log-law.
(3) The DML-Law, MLW-Law, and CML-Law models recorded in the literature were examined in this study.Figures 2-7 show that all of these models are able to predict all experimental data in the inner region.However, DML-Law, MLW-Law, and CML-Law are not able to agree with the experimental data properly in the outer region.It is clearly found that the wake function indeed exists in uniform flow.

Figure 1 :
Figure 1: General view of the flume.

Figure 9 :
Figure 9: Variation of  with parameter .