Next Article in Journal
Influence of the Evolution of Marine Industry Structure on the Green Total Factor Productivity of Marine Economy
Next Article in Special Issue
Velocity Field and Turbulence Structure around Spur Dikes with Different Angles of Orientation under Ice Covered Flow Conditions
Previous Article in Journal
‘Hidden Hot Springs’ as a Source of Groundwater Fluoride and Severe Dental Fluorosis in Malawi
Previous Article in Special Issue
Assessment of Critical Shear Stress and Threshold Velocity in Shallow Flow with Sand Particles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analytical Models of Velocity, Reynolds Stress and Turbulence Intensity in Ice-Covered Channels

1
State Key Laboratory of Eco-Hydraulics in Northwest Arid Region of China, Xi’an University of Technology, Xi’an 710048, China
2
Yellow River Engineering Consulting Co., Ltd., Zhengzhou 450003, China
3
Southwest Branch of China Construction Third Engineering Bureau Group Co., Ltd., Chengdu 610041, China
4
PowerChina Northwest Engineering Corporation Limited, Xi’an 710065, China
*
Author to whom correspondence should be addressed.
Water 2021, 13(8), 1107; https://doi.org/10.3390/w13081107
Submission received: 28 February 2021 / Revised: 11 April 2021 / Accepted: 14 April 2021 / Published: 17 April 2021
(This article belongs to the Special Issue Fluvial Hydraulics Affected by River Ice and Hydraulic Structures)

Abstract

:
Ice cover in an open channel can influence the flow structure, such as the flow velocity, Reynolds stress and turbulence intensity. This study analyzes the vertical distributions of velocity, Reynolds stress and turbulence intensity in fully and partially ice-covered channels by theoretical methods and laboratory experiments. According to the experimental data, the vertical profile of longitudinal velocities follows an approximately symmetry form. Different from the open channel flow, the maximum value of longitudinal velocity occurs near the middle of the water depth, which is close to the channel bed with a smoother boundary roughness compared to the ice cover. The measured Reynolds stress has a linear distribution along the vertical axis, and the vertical distribution of measured turbulence intensity follows an exponential law. Theoretically, a two-power-law function is presented to obtain the analytical formula of the longitudinal velocity. In addition, the vertical profile of Reynolds stress is obtained by the simplified momentum equation and the vertical profile of turbulence intensity is investigated by an improved exponential model. The predicted data from the analytical models agree well with the experimental ones, thereby confirming that the analytical models are feasible to predict the vertical distribution of velocity, Reynolds stress and turbulence intensity in ice-covered channels. The proposed models can offer an important theoretical reference for future study about the sediment transport and contaminant dispersion in ice-covered channels.

1. Introduction

Most rivers at high altitude in cold northern regions always freeze in the winter and form ice sheets. The ice sheet in some Canadian rivers is at least 0.6 m thick and lasts for at least 4 months [1]. The wetted perimeter of the cross section and flow resistance in ice-covered flows increases with the presence of the ice sheet, which significantly affects the hydraulic characteristics of river and topographical features and greatly change the flow velocity distribution, flow transport capacity and sediment transport rate [2,3,4,5,6,7,8,9,10,11]. Therefore, it is necessary to study the ice-covered flows.
Unlike the open channel flow, flows in the ice-covered channel have asymmetric forms. The presence of the ice sheet makes the maximum streamwise velocity appear at the inner center of the flow, and the location of the maximum velocity is generally considered the division point in the asymmetrically distributed flow [12,13]. The asymmetric distribution mainly depends on the roughness of the ice sheet and the riverbed. The division point of the velocity tends to be away from the rougher surface [14]. Previous studies have shown that the main effects of the ice sheets on alluvial channels can be summarized as: they increase the water level (compared to open channels at the same flow rate), reduce the average flow velocity, increase the channel drag force and reduce the bed sediment transport rate [15].
Previous investigations on the flow characteristics of ice-covered flows are mainly obtained through experiments [16,17,18]. Many laboratory experimental data and field observations show that the vertical distribution of streamwise velocity forms a double-layer, which is characterized by the plane of maximum velocity. Shen and Harden (1978) [19] and Lau and Krishnappan (1981) [20] applied this double-layer theory in the ice-covered flows. They divided the ice-covered flows into two separate layers: upper ice layer and lower channel bed. Parthasarathy and Muste (1994) [13] found that the zero shear stress plane in the ice-covered flows was inconsistent with the maximum velocity plane. Chen et al. (2015) [21] proposed that the horizontal plane of zero shear stress should be determined as the dividing plane of sublayers and modified the double-layer assumption.
Some researchers have studied the vertical profiles of the longitudinal velocity by numerical simulation methods and proposed two-dimensional and three-dimensional models [22,23], which may require sensitive hydraulic parameters. However, these hydraulic parameters do not have explicit expressions to be calculated, and they have uncertainty that cannot be ignored. Except the numerical methods, previous researchers also attempted to obtain the vertical distribution of streamwise velocity by theoretical methods [24,25]. Uzuner (1975) [26] separated the ice-covered flow into two layers based on the position of the maximum velocity and assumed that the velocity distributions in the upper ice layer and lower bed layer were consistent with the logarithmic distribution law, and Manning formula could be independently applied to each layer. In addition, a two-power-law function was adopted to calculate the vertical distribution of the longitudinal velocity. Teal et al. (1994) [27] reported a reasonable fit of streamwise velocity data to the two-power-law function. Compared to the two-power-law function, the logarithmic law appears to overestimate velocities near the location of maximum velocity. The two-power-law function is a reasonable extension of the power-law expression that is used to describe velocity profiles of open channel flow. The advantage for the two-power-law function is that it describes the entire flow with a single continuous curve. Hence, the two-power-law function deserves further study.
The objective of this study is to obtain the vertical profiles of longitudinal velocity, Reynolds stress and turbulence intensity in ice-covered channels. Hence, this study focuses on (1) adopting a two-power-law function to calculate the vertical profile of the longitudinal velocity, (2) simplifying the time-averaged momentum equation to analyze the vertical distributions of Reynolds stress, (3) improving the exponential model to calculate the turbulence intensity and (4) employing experimental data to validate the theoretical models and giving detail discussion on the coefficients in the theoretical models.

2. Material and Methods

Since the free water surface is covered with ice, the open channel flow changes to a closed conduit flow but retains the flow characteristics of the open channel flow. The vertical distribution of longitudinal velocity significantly changes in the open channel flow with ice-cover, which reflects that the plane of maximum longitudinal velocity shifts from the free water surface to the inner water (see Figure 1). Here, we used a Cartesian coordinate system with the x axis in the main flow direction and z axis in the water depth direction. Corresponding to the x and z axes, u and w are the longitudinal and vertical velocities, respectively.

2.1. Vertical Distribution of Longitudinal Velocity

The presence of the ice cover causes the increase of the wetted perimeter, which increases the composite flow resistance. Owing to the different roughness of ice cover and channel bed, the velocity profile is vertically asymmetric (see Figure 1). In Figure 1, the vertical flow structure can be divided into two independent layers at the plane of maximum velocity, i.e., the upper ice layer and lower bed layer. The flow in the upper ice layer is mainly affected by the ice cover, and the lower flow is primarily affected by the channel bed. The vertical location of the maximum velocity is determined by the roughness of the ice cover and channel bed, i.e., the maximum velocity will not occur at the free water surface but near the middle of the water depth.
Tsai and Ettema (1994) [28] adopted a power-law function to predict the velocity profile in the open channel flow. For the ice-covered flow, a two-power-law function was used to obtain the vertical profile of the longitudinal velocity [27]. The first advantage of the two-power-law function is that it describes the flow just using a single continuous curve. The second is that the analytical solution of velocity is adjustable in accordance with the roughness change of channel bed and ice cover. The two-power-law is shown as
u = K 0 ( z H ) 1 / m b ( 1 z H ) 1 / m i
where z is the vertical axis with z = 0 at the channel bed; H is the water depth; K 0 is the flow parameter based on a given flow discharge per unit flow width and m b and m i are the parameters corresponding to the boundary roughness of the channel bed and ice cover. When parameter m i approaches infinity, the velocity profile becomes equivalent to a single power law expression, which indicates that the ice cover disappears.
The velocity gradient obtained from Equation (1) is
d u d z = K 0 H [ 1 m b ( z H ) 1 m b 1 ( 1 z H ) 1 / m i 1 m i ( z H ) 1 m b ( 1 z H ) 1 m i 1 ]
By setting the velocity gradient u z to zero, we deduced the position of the maximum velocity as
( z H ) u m a x = h m H = m i m i + m b
where u m a x is the maximum longitudinal velocity and h m is the height of the maximum velocity from the channel bed.

2.2. Vertical Distribution of Reynolds Stress

To predict the vertical distribution of Reynolds stress in ice-covered channels for a steady uniform flow, the time-averaged momentum equation in the longitudinal direction can be simplified to
u w ¯ + ν d u d z + u i 2 z H = u b 2 ( 1 z H )
where u and w are turbulent fluctuations of the longitudinal and vertical velocities, u w ¯ is Reynolds stress; ν is the flow kinematic viscosity; u i is the shear velocity at the ice cover and u b is the shear velocity at the channel bed.
Since the viscosity shear stress in the flow is much smaller than the Reynolds stress, the viscosity shear stress can be neglected in Equation (4) [29,30]. Then, there is a linear relationship between the shear stress ( τ ) and the vertical axis ( z ) . Equation (4) has another form and is shown as
τ ( z ) = τ b ( τ b τ i ) z H
where τ ( z ) is the shear stress at distance z ; τ b and τ i are the Reynolds stresses at the channel bed and ice cover; z is the vertical distance from the channel bed and H is the water depth.
Following the method of Rowinski and Kubrak (2002) [31], who combine the eddy viscosity and flow conditions by a mixing length concept, the shear stress can be given as
τ ( z ) = ρ l 2 u z | u z |
where ρ is the flow density and l is the proposed mixing length. Chen et al. (2015) [21] demonstrate that the mixing length theory can be applied at a certain distance from the fixed boundaries. Two mixing lengths, l b and l i , are proposed corresponding to the two fixed boundaries. The first one is the mixing length of channel bed and the second one is the mixing length of the ice cover. Hence, the linear relationships of two mixing lengths are approximately written as
{ l b = κ z ,   0 z δ b l i = κ ( H z ) ,   H δ i z H  
where δ b and δ i are the distances from the channel bed and ice cover, where the mixing length theory is valid. κ = 0.41 is the von Karman constant.
By substituting Equation (7) into Equation (6), we updated the expression of shear stress as
{ τ ( z ) = ρ κ 2 z 2 u z | u z | ,   0 z δ b τ ( z ) = ρ κ 2 ( H z ) 2 u z | u z | ,   H δ i z H  
The gradient of shear stress from Equation (8) is
{ τ z = 2 ρ κ 2 z 2 u z ( 1 z u z + 2 u z 2 ) ,   0 z δ b τ z = 2 ρ κ 2 ( H z ) 2 u z ( 1 ( H z ) u z 2 u z 2 ) ,   H δ i z H  
By substituting the velocity gradient (Equation (2)) into the Reynolds stress gradient (Equation (9)) and letting τ z = 0 , we obtained the two maximum Reynolds stresses, τ b m a x and τ i m a x , near the channel bed and ice cover, respectively, which are presented as
{ τ b m a x = ρ κ 2 z b m a x 2 u z b m a x | u z b m a x | ,   0 z δ b τ i m a x = ρ κ 2 ( H z i m a x ) 2 u z i m a x | u z i m a x | ,   H δ i z H
where z b m a x and z i m a x denote the positions where the corresponding maximum shear stresses τ b m a x and τ i m a x occurred, and they can be calculated as
{ z bmax = A b A b 2 4 B b C b 2 B b H z imax = A i A i 2 4 B i C i 2 B i H
where A b , B b , C b , A i , B i , and C i are the constants relevant to the roughness of the channel bed and ice cover. Their expressions are
{ { A b = m b 2 m i m b 2 B b = ( m b + m i ) ( 1 m b + 1 m i ) C b = m i m b { A i = m i + 2 m i m b + 2 B i = ( m b + m i ) ( 1 m b + 1 m i ) C i = m i m i m b
The shear stresses at the fixed boundaries, i.e., the channel bed and ice cover, can be considered the maximum shear stresses in the lower bed layer and upper ice layer. Hence, the linear relationship (Equation (5)) between shear stress and vertical position can be rewritten as
τ ( z ) = τ b m a x ( τ b m a x τ i m a x ) z H

2.3. Vertical Distribution of Turbulence Intensity

The flow turbulence characteristics can be presented by the turbulence intensity, which can be calculated by the root mean square of the fluctuating longitudinal velocity, i.e., u r m s = ( u ) 2 . In the open channel flow, the term u r m s reaches its maximum value near the channel bed and it has a linear relationship with the vertical distance. Nezu and Rodi (1986) [32] established an exponential model to describe the vertical distribution of the turbulence intensity. However, in the ice-covered channel flow, it first decreases from the channel bed and reaches its minimum value near the middle flow depth, after which it gradually increases towards to the ice cover. For the profile of the turbulence intensity in the ice-covered channel, we adopted the division scheme to investigate the profile of the turbulence intensity. Based on the location of zero shear stress, its profile is separated into two regions. Then, the exponential model can be applied in each region as follows
{ u r m s u b = D b e E b z / h τ , 0 < z < h τ u r m s u i = D i e E i z / h τ , h τ < z < H
where D b , E b , D i and E i are constant parameters.

3. Experimental Verification

Experiments were performed in a rectangular glass flume with a total length of 20 m, a width of 1 m and a depth of 0.5 m in the State Key Laboratory of Water Resources and Hydropower Engineering Science in Wuhan University to investigate the longitudinal velocity profile in the ice-covered flow (Figure 2). The flow is circulated by a pump system and adjusted to be steady by a tailgate at the end of the flume. A plastic foam board was used to simulate the ice with the length of 15 m and a width of 1 m [33].
The rectangular coordinates are as follows: x is the main flow direction with x = 0 at the beginning of the ice cover; y is the lateral direction with y = 0 at the sidewall and z is the vertical direction counting from the channel bed. The velocity was measured at the cross section with x = 9 m, which is sufficiently far from the ice cover entry to form a fully developed flow. In the measured cross section, 14 measured lines are arranged, where the measured points are also arranged with a uniform vertical distance 1 cm. The detail layout is shown in Figure 3.
A 3D acoustic Doppler velocimeter (ADV) with a precision of ±0.25 cm/s was adopted to measure the instantaneous velocity. The maximum sampling frequency was 50 Hz, and the sampling time for every measured point was 120 s, which resulted in 6000 velocity data for each point. Five experimental conditions were considered here, and the variables among these five cases were the cover conditions and water depth H, which changes from 15 to 20 cm. The details of the experimental parameters are listed in Table 1.

4. Model Parameters

To predict the profiles of velocity, Reynolds stress and turbulence intensity in the ice-covered flow, the model parameters, i.e., m b and m i , n b and n i , K 0 , etc., should be given first.

4.1. m b and m i

In the single-power law for the free surface flow, exponent m is affected by Darcy–Weisbach resistance coefficient f , i.e., m = κ ( 8 f ) 0.5 [30]. For the two-power expression, exponents m b and m i can also be related to Darcy–Weisbach resistance coefficients of channel bed f b and ice cover f i . Specifically, the value of one exponent m can be related to the resistance coefficient of one fixed boundary. Hence, exponents m b and m i can be expressed as
{ m b = κ ( 8 f b ) 0.5 m i = κ ( 8 f i ) 0.5
The resistance coefficients can be calculated as [34]
{ f b = 8 n b 2 g R b 1 / 3 f i = 8 n i 2 g R i 1 / 3
where n denotes Manning’s roughness coefficients; R indicates the hydraulic radiuses and the subscript b and i in them represent the channel bed and ice cover. g is the gravity acceleration. Here, the ratio of channel width to water depth was approximately greater than 5, which indicates that the flow could be considered a shallow flow. Hence, the hydraulic radius of each sublayer could be simplified as the corresponding sublayer depth, i.e., R b = h b and R i = h i . The flow depths of the lower bed layer and upper ice layer have been defined from the location of the maximum velocity and are represented by h m and H h m . Hence, we easily obtained
{ h m = h b H h m = h i
Substituting Equation (16) into Equation (15), we obtained
{ m b = κ ( R b 1 / 3 n b 2 g ) 0.5 m i = κ ( R i 1 / 3 n i 2 g ) 0.5
Then, we substituted Equations (17) and (18) into Equation (3) ( h m H = m i m i + m b ) and obtained
h b = H n b ( H h b ) 1 / 6 n b ( H h b ) 1 / 6 + n i h b 1 / 6
Once Manning’s coefficients n b and n i are known, the depth of the lower bed layer ( h b ) is iteratively solved. Then, exponents m b and m i were calculated by Equation (18).

4.2. n b and n i

The basic calculation model of Darcy–Weisbach friction coefficient f is
f = 8 u 2 U 2
where U ( = Q / B H ) is the averaged velocity in a cross section. Substituting Equation (16) into Equation (20), we obtained
n = u U g R 1 / 6
For the ice-covered flow, the flow is divided as two layers, and two depth-averaged velocities, U b and U i , are averaged from the channel bed and ice cover to the position of the maximum velocity. Hence, Manning’s coefficient in each sublayer can be expressed as
{ n b = u b U b g h b 1 / 6 n i = u i U i g h i 1 / 6
These two Manning’s parameters can be determined by the measured velocity profile.
For the fully developed asymmetric flow, the method of logarithmic law can be adopted to predict the vertical distribution of the longitudinal velocity. However, because the velocity gradient of the maximum velocity is discontinuous, the method of logarithmic law may not be completely feasible for the entire water depth. For the ice-covered flow, the method of logarithmic law can be applied in each sublayer, i.e., the lower bed layer and upper ice layer. Bonakdari et al. (2008) [35] further demonstrate that the turbulent boundary layer consists of the inner region near the sidewall and outer region far from the sidewall and propose that the flow velocity inside the inner region can better satisfy the distribution form of the logarithmic law. In this study, we define that the inner region inside the lower bed layer starts from the channel bed ( z = 0) to 0.2 H , and the inner region inside the upper ice layer is from 0.8 H to the ice cover ( H ). Therefore, we used the following logarithmic law to describe the velocity profile in the two inner regions [36]
{ U b u b = 1 κ l n ( 30 z k s b ) U i u i = 1 κ l n ( 30 z k s i )
where k s b and k s i are the roughness heights of the turbulence boundary in the lower bed layer and upper ice layer. The logarithmic law can be simplified to
{ U b = a b l n ( z ) + F b U i = a i l n ( z ) + F i
where a b = u b / κ , a i = u i / κ , F b = a b ( l n 30 k s b ) and F i = a i ( l n 30 k s i ) . These parameters can be obtained by regression analysis based on the measured longitudinal velocity data; then, shear velocities u b and u i and roughness heights k s b and k s i can be obtained.

4.3. K 0

According to the previous work [28,37], K 0 is the flow model parameter for a known flow rate per unit flow width and is
K 0 = U / K 1
where K 1 is the normalized depth-averaged velocities for total flow and is calculated by
K 1 = 0 1 ( z H ) 1 / m b ( 1 z H ) 1 / m i d ( z / H )

5. Results and Discussion

5.1. Model Verification

With the above model parameters, the theoretical model can be applied to the experimental data. Figure 4 compares the predicted velocity from the model with the measured ones. Both profiles followed the assumed one, where the maximum velocity occurs near the middle of the water depth, and the minimum velocity appears near the fixed boundaries. The analytical results were consistent with the experimental velocities, so the proposed model and corresponding parameters are reasonable and reliable to be used to predict the velocity profile in the straight open channel flow covered by the ice cover.
Figure 5 presents the profiles of measured Reynolds stress and analytical Reynolds stress for all experimental cases. The Reynolds stress in the longitudinal direction demonstrates a linear distribution along the water depth. Zero shear stress occurred at the location of h τ near the middle of the water depth, after which the absolute value of Reynolds stress gradually increased and reached each peak near the fixed boundaries. The measured and predicted data in Figure 5 were basically consistent and had identical trends, which indicates that our proposed model for Reynolds stresses was feasible.
Figure 6 shows the results of the turbulence intensity from the experiments and analytical model. The turbulence intensity reached its minimum near the middle water depth, from which it had an increasing trend towards the channel bed and ice cover. The profiles of turbulence intensity are consistent with the results from Papanicolaou et al. (2007) [38], who demonstrate that the turbulence production near the central region is small, and its diffusion effect is significant, but the turbulence production near the fixed boundaries reaches its maximum value. Overall, the analytical model can catch the general trend of the measured turbulence intensity, although the measured one has some fluctuations, which confirms that the model can be applied to calculate the turbulence intensity.
To find the difference between analytical and measured data, the error analysis was conducted. An absolute error is defined as the difference between analytical and measured time-averaged velocities. Hence, the average absolute error is calculated as
Δ a ¯ = 1 N 1 N | ( I ) a n a l y t i c a l ( I ) m e a s u r e d |
where N is the number of measured points in each measured line for each case, ( I ) a n a l y t i c a l and ( I ) m e a s u r e d are the analytical and measured values and I represents the variables, i.e., longitudinal velocity, Reynolds stress and turbulence intensity.
The average relative error is defined as
Δ r ¯ = 1 N 1 N | ( I ) a n a l y t i c a l ( I ) m e a s u r e d ( I ) m e a s u r e d | × 100 %
As shown in Table 2, the time-averaged longitudinal velocities obtained from the proposed model were reliable within an accuracy of 0.019 m/s in terms of the average absolute error Δ a ¯ . The average relative error Δ r ¯ was 2.86–10.97%. The average absolute error and average relative error of the Reynolds stress were within 0.093 m/s and 13.63%, respectively. The average absolute error and relative error of the turbulence intensity were within 0.0019 m/s and 12.54%, respectively. All values in Table 2 were below 20%, which further confirmed that the proposed analytical model was reliable and feasible to predict the velocity and turbulence structure in the ice-covered flow.

5.2. Discussion

5.2.1. Manning’s Coefficients n b and n i

Table 3 lists the calculated Manning’s roughness coefficients of the channel bed and ice cover for each case, i.e., n b and n i . Manning’s roughness coefficient of the channel bed was 0.012–0.015, and its mean value was 0.0138 with the standard deviation of 0.0012, which verified that n b hardly changed among all studied cases. Manning’s roughness coefficients of the ice cover in all cases were 0.017–0.02, which had small fluctuations and were lightly larger than those of the channel bed. The mean value of n i was 0.0182 with a standard deviation of 0.0012. For the experimental channel, Manning’s roughness coefficients of the channel bed and ice cover should be considered two specific constants. We took the mean Manning’s roughness coefficient in each layer as the final Manning’s roughness coefficient. Hence, n b was set to be 0.0138, and n i was equal to 0.0182.

5.2.2. Flow Parameters m b and m i

The values of exponents m b and m i in this study are shown in Table 3. All values of m b were 5.4 m b 7.1 , and all m i were 4.3 m i 5.3 . Teal et al. (1994) [27] estimated m b and m i by nonlinear regression for more than measured 2300 vertical velocity profiles, and they found that these two parameters were 1.5–8.5, which includes the theoretical range (see Table 3). For an open channel flow, only exponent m b is considered and is approximately 6–7, which is unlike the values used here. Hence, both m b and m i in the covered flow are affected by the roughness characteristics of both the channel bed and ice cover.
The shape of the vertical profile of the longitudinal velocity is determined by exponents m b and m i . The ratio of them can be given from Equation (17) as
m i m b = ( h i h b ) 1 / 6 n b n v
According to Equation (29), the ratio of m b and m i is mainly influenced by the roughness coefficients of the channel bed and ice cover. This result is also confirmed by the measured data of Li et al. (2020) [29]. They demonstrate that the vertical distribution of the velocity remains constant with changing water depth and flow rate under the same m b and m i or the same ratio of m b and m i . In the asymmetric flow, the maximum velocity tended to be closer to the smooth boundary with smaller roughness coefficient. In Table 3, the channel bed had a smaller roughness coefficient than the ice cover, which corresponded to the close maximum velocity to the channel bed, which was verified by h m < H h m . In general, both open channel flow and symmetry flow can be considered as the special cases of asymmetric flow. For the open channel flow, the exponent m i tends to infinity. For the symmetry flow, m b = m i .

5.2.3. Comparison of h m and h τ

Table 3 presents the locations ( h m and h τ ) where the maximum velocity and zero shear stress occur for all cases. These two locations were not consistent and had distinct difference, i.e., calculated h m / H > h τ / H . The location of zero shear stress is closer to the channel bed than that of maximum velocity [13,29,38]. By contrast, the locations of the maximum velocity and zero shear stress for the symmetry and open channel flows were the same. Specifically, the location of maximum velocity and zero shear stress for the symmetry flows was at the middle water depth, and the location for the open channel flows was at the free water surface.
Considering the locations of the maximum velocity and zero shear stress for cases 1, 3 and 5, when the ratio of n i to n b increased, h m / H decreased, so the location of the maximum velocity approached the channel bed, and the vertical inhomogeneity of the velocity profile was strengthened. Meanwhile, τ i / τ b increased with the increase in n i / n b ; then, h τ / H decreased, which indicates the location of the zero shear stress gets closer to the channel bed.

5.2.4. Empirical Constants D b , E b , D i and E i

The empirical constants of turbulence intensity for all cases are listed in Table 4. No remarkable changes of D b and D i were observed in any case. The mean D b was 2.22 with a standard deviation of 0.05, and the mean D i was 2.16 with a standard deviation of 0.06. Hence, it is reasonable to consider that D b = D i . The difference between E b and E i was not negligible, similar to the results of Li et al. (2020) [29].

6. Conclusions

The existence of ice cover dramatically changed the flow velocity and turbulence structure. We here proposed theoretical models to describe the vertical distribution of longitudinal velocity, shear stress and turbulence intensity. By dividing the ice-covered flow into an ice-affected layer and a channel bed-affected layer, a two-power-law function was adopted to predict the vertical profile of velocity. The calculated velocity distribution presents that the maximum velocity occurred near the middle of the water depth close to the channel bed with smooth boundary. Theoretical analysis shows that the shear stress had a linear distribution form in the vertical direction, with the positive values in the lower bed layer and negative values in the upper ice layer. Moreover, the Manning’s roughness coefficient of the ice cover was larger than that of the channel bed. The two exponents m b and m i were influenced by the roughness coefficients of the channel bed and ice cover. The location of zero shear stress was not the same as that of maximum velocity and was closer to the smooth fixed boundary than the plane of maximum velocity, namely, h m > h τ . A comparison of the analytical and experimental velocities, the Reynolds stress and turbulence intensity displays that the theoretical models can provide satisfied predictions of the vertical distribution of these flow characteristics. This study expands our understanding of the effects of ice cover on the hydraulic characteristics in the open channels. However, we still need to do more research to explore the application of the proposed models in other conditions, like compound channels or confluence channels and we will involve comprehensive experiments to reveal detailed flow characteristics, such as vortex structure.

Author Contributions

Conceptualization, J.Z. and W.W.; methodology, J.Z.; software, Z.X., Y.Z.; resources, Q.L., Y.Z.; data curation, H.Q.; writing—original draft preparation, J.Z.; writing—review and editing, J.Z.; visualization, H.Q.; supervision, W.W.; project administration, Z.L.; funding acquisition, Z.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key R&D Program of China (2017YFC0504704), the National Natural Science Foundation of China (Grant No. 51609198), the Science and Technology Project Funded by Shaanxi Provincial Department of Water Resources (2020slkj-10) and the Technology Project Funded by Clean Energy and Ecological Water Conservancy Engineering Research Center (QNZX-2019-03).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Davar, K.S.; Elhadi, N.A. Management of ice-covered rivers: Problems and perspectives. J. Hydrol. 1981, 51, 245–253. [Google Scholar] [CrossRef]
  2. Chen, Y.; Wang, Z.; Zhu, D.; Liu, Z. Longitudinal dispersion coefficient in ice-covered rivers. J. Hydraul. Res. 2016, 54, 558–566. [Google Scholar] [CrossRef]
  3. Chen, G.; Gu, S.; Li, B.; Zhou, M.; Huai, W. Physically based coefficient for streamflow estimation in ice-covered channels. J. Hydrol. 2018, 563, 470–479. [Google Scholar] [CrossRef]
  4. Knack, I.; Shen, H.-T. Sediment transport in ice-covered channels. Int. J. Sediment Res. 2015, 30, 63–67. [Google Scholar] [CrossRef]
  5. Lee, M.; Moser, R.D. Direct numerical simulation of turbulent channel flow up to Reτ ≈ 5200. J. Fluid Mech. 2015, 774, 395–415. [Google Scholar] [CrossRef]
  6. Lotsari, E.; Tarsa, T.; Mri, M.K.; Alho, P.; Kasvi, E. Spatial variation of flow characteristics in a subarctic meandering river in ice-covered and open-channel conditions: A 2d hydrodynamic modelling approach. Earth Surf. Process. Landf. 2019, 44, 1509–1529. [Google Scholar] [CrossRef] [Green Version]
  7. Turcotte, B.; Morse, B.; Bergeron, N.E.; Roy, A.G. Sediment transport in ice-affected rivers. J. Hydrol. 2011, 409, 561–577. [Google Scholar] [CrossRef]
  8. Wang, F.; Huai, W.; Liu, M.; Fu, X. Modeling depth-averaged streamwise velocity in straight trapezoidal compound channels with ice cover. J. Hydrol. 2020, 585, 124336. [Google Scholar] [CrossRef]
  9. Smith, B.T.; Ettema, R. Flow Resistance in Ice-Covered Alluvial Channels. J. Hydraul. Eng. 1997, 123, 592–599. [Google Scholar] [CrossRef]
  10. Wang, J.; Wu, Y.; Sui, J.; Karney, B. Formation and movement of ice accumulation waves under ice cover—An experimental study. J. Hydrol. Hydromech. 2019, 67, 171–178. [Google Scholar] [CrossRef] [Green Version]
  11. Namaee, M.R.; Sui, J. Velocity profiles and turbulence intensities around side-by-side bridge piers under ice-covered flow condition. J. Hydrol. Hydromech. 2020, 68, 70–82. [Google Scholar] [CrossRef] [Green Version]
  12. Hanjalić, K.; Launder, B.E. Fully developed asymmetric flow in a plane channel. J. Fluid Mech. 1972, 51, 301–335. [Google Scholar] [CrossRef]
  13. Parthasarathy, R.N.; Muste, M. Velocity Measurements in Asymmetric Turbulent Channel Flows. J. Hydraul. Eng. 1994, 120, 1000–1020. [Google Scholar] [CrossRef]
  14. Tatinclaux, J.; Gogus, M. Asymmetric Plane Flow with Application to Ice Jams. J. Hydraul. Eng. 1983, 109, 1540–1554. [Google Scholar] [CrossRef]
  15. Lau, Y.L.; Krishnappan, B.G. Sediment Transport Under Ice Cover. J. Hydraul. Eng. 1985, 111, 934–950. [Google Scholar] [CrossRef]
  16. Muste, M.; Braileanu, F.; Ettema, R. Flow and sediment transport measurements in a simulated ice-covered channel. Water Resour. Res. 2000, 36, 2711–2720. [Google Scholar] [CrossRef]
  17. Robert, A.; Tran, T. Mean and turbulent flow fields in a simulated ice-covered channel with a gravel bed: Some laboratory observations. Earth Surf. Process. Landforms 2012, 37, 951–956. [Google Scholar] [CrossRef]
  18. Tao, L. Experimental Study on vertical velocity distribution of water flow under ice sheet. Eng. Constr. 2015, 29, 370–371. [Google Scholar]
  19. Shen, H.T.; Harden, T.O. The effect of ice cover on vertical transfer in stream channels. J. Am. Water Resour. Assoc. 1978, 14, 1429–1439. [Google Scholar] [CrossRef]
  20. Lau, Y.L.; Krishnappan, B.G. Ice Cover Effects on Stream Flows and Mixing. J. Hydraul. Div. 1981, 107, 1225–1242. [Google Scholar] [CrossRef]
  21. Chen, G.; Gu, S.; Huai, W.; Zhang, Y. Boundary Shear Stress in Rectangular Ice-Covered Channels. J. Hydraul. Eng. 2015, 141, 06015005. [Google Scholar] [CrossRef]
  22. Attar, S.; Li, S. Data-fitted velocity profiles for ice-covered rivers. Can. J. Civ. Eng. 2012, 39, 334–338. [Google Scholar] [CrossRef]
  23. Sui, J.; Wang, J.; Yun, H.E.; Krol, F. Velocity profiles and incipient motion of frazil particles under ice cover. Int. J. Sediment Res. 2010, 25, 39–51. [Google Scholar] [CrossRef]
  24. Larsen, P.A. Head losses caused by an ice cover on open channels. J. Boston Soc. Civ. Eng. 1969, 56, 45–67. [Google Scholar]
  25. Sayre, W.W.; Song, G.B. Effects of Ice Covers on Alluvial Channel Flow and Sediment Transport Processes. IIHR Report No. 218, University of Lowa. Available online: https://apps.dtic.mil/dtic/tr/fulltext/u2/a066991.pdf (accessed on 16 April 2021).
  26. Uzuner, M.S. The composite roughness of ice covered streams. J. Hydraul. Res. 1975, 13, 79–102. [Google Scholar] [CrossRef]
  27. Teal, M.J.; Ettema, R.; Walker, J.F. Estimation of Mean Flow Velocity in Ice-Covered Channels. J. Hydraul. Eng. 1994, 120, 1385–1400. [Google Scholar] [CrossRef]
  28. Tsai, W.; Ettema, R. Modified Eddy Viscosity Model in Fully Developed Asymmetric Channel Flows. J. Eng. Mech. 1994, 120, 720–732. [Google Scholar] [CrossRef]
  29. Li, Q.; Zeng, Y.-H.; Bai, Y. Mean flow and turbulence structure of open channel flow with suspended vegetation. J. Hydrodyn. 2020, 32, 314–325. [Google Scholar] [CrossRef]
  30. Zhang, J.; Lei, J.; Huai, W.; Nepf, H. Turbulence and Particle Deposition Under Steady Flow Along a Submerged Seagrass Meadow. J. Geophys. Res. Oceans 2020, 125, 2019–015985. [Google Scholar] [CrossRef]
  31. Rowiński, P.M.; Kubrak, J. A mixing-length model for predicting vertical velocity distribution in flows through emergent vegetation. Hydrol. Sci. J. 2002, 47, 893–904. [Google Scholar] [CrossRef] [Green Version]
  32. Nezu, I.; Rodi, W. Open-channel Flow Measurements with a Laser Doppler Anemometer. J. Hydraul. Eng. 1986, 112, 335–355. [Google Scholar] [CrossRef]
  33. Zhong, Y.; Huai, W.; Chen, G. Analytical Model for Lateral Depth-Averaged Velocity Distributions in Rectangular Ice-Covered Channels. J. Hydraul. Eng. 2019, 145, 04018080. [Google Scholar] [CrossRef]
  34. Huai, W.-X.; Xu, Z.-G.; Yang, Z.-H.; Zeng, Y.-H. Two dimensional analytical solution for a partially vegetated compound channel flow. Appl. Math. Mech. 2008, 29, 1077–1084. [Google Scholar] [CrossRef]
  35. Bonakdari, H.; Larrarte, F.; Lassabatere, L.; Joannis, C. Turbulent velocity profile in fully-developed open channel flows. Environ. Fluid Mech. 2008, 8, 1–17. [Google Scholar] [CrossRef]
  36. Zare, S.G.A.; Moore, S.A.; Rennie, C.D.; Seidou, O.; Ahmari, H.; Malenchak, J. Estimation of composite hydraulic resistance in ice-covered alluvial streams. Water Resour. Res. 2016, 52, 1306–1327. [Google Scholar] [CrossRef] [Green Version]
  37. Han, L.; Zeng, Y.; Chen, L.; Li, M. Modeling streamwise velocity and boundary shear stress of vegetation-covered flow. Ecol. Indic. 2018, 92, 379–387. [Google Scholar] [CrossRef]
  38. Papanicolaou, A.N.; Elhakeem, M.; Hilldale, R. Secondary current effects on cohesive river bank erosion. Water Resour. Res. 2007, 43, 497–507. [Google Scholar] [CrossRef]
Figure 1. Schematic profile of the longitudinal velocity in an iced-covered open channel.
Figure 1. Schematic profile of the longitudinal velocity in an iced-covered open channel.
Water 13 01107 g001
Figure 2. (a) Schematic of the experimental flume, the dark blue denoting the flow in the flume and (b) images of the experimental site.
Figure 2. (a) Schematic of the experimental flume, the dark blue denoting the flow in the flume and (b) images of the experimental site.
Water 13 01107 g002
Figure 3. Layout of the measured cross section, lines and points.
Figure 3. Layout of the measured cross section, lines and points.
Water 13 01107 g003
Figure 4. Comparison of the measured and analytical velocities. Black squares denote the measured data and red lines denote the analytical ones.
Figure 4. Comparison of the measured and analytical velocities. Black squares denote the measured data and red lines denote the analytical ones.
Water 13 01107 g004
Figure 5. Comparison of the measured and analytical Reynolds stresses. Black squares denote the measured data and red lines denote the analytical ones.
Figure 5. Comparison of the measured and analytical Reynolds stresses. Black squares denote the measured data and red lines denote the analytical ones.
Water 13 01107 g005aWater 13 01107 g005b
Figure 6. Comparison of the measured and analytical turbulence intensity. Black squares denote the measured data and red lines denote the analytical ones.
Figure 6. Comparison of the measured and analytical turbulence intensity. Black squares denote the measured data and red lines denote the analytical ones.
Water 13 01107 g006aWater 13 01107 g006b
Table 1. Basic details of the characteristic parameters. B is the width of the flume, H is the water depth, S is the bed slope and R e is the Reynolds number.
Table 1. Basic details of the characteristic parameters. B is the width of the flume, H is the water depth, S is the bed slope and R e is the Reynolds number.
CasesCover Condition B (m) H (m) S R e
1Full ice cover10.150.00125590
2Full ice cover10.1850.00133725
3Symmetrical shore cover10.160.00132480
4Symmetrical shore cover10.200.00149160
5Asymmetrical shore cover10.160.00131888
Table 2. Error statistics for the longitudinal velocity, Reynolds stress and turbulence intensity. Δ a ¯ is the average absolute error and Δ r ¯ is the average relative error.
Table 2. Error statistics for the longitudinal velocity, Reynolds stress and turbulence intensity. Δ a ¯ is the average absolute error and Δ r ¯ is the average relative error.
CasesVelocityReynolds StressTurbulence Intensity
Δ a ¯ ( m / s ) Δ r ¯ ( % ) Δ a ¯ ( m / s )   Δ r ¯ ( % ) Δ a ¯ ( m / s )   Δ r ¯ ( % )
10.00653.990.0367.010.001310.52
20.00573.180.03810.000.00118.09
30.00592.860.0509.110.00128.05
40.0115.250.0828.420.00126.99
50.01910.970.09313.630.001912.54
Table 3. Values of the model parameters to predict the velocity and Reynolds stress.
Table 3. Values of the model parameters to predict the velocity and Reynolds stress.
Cases n b n i m b m i h m / H h τ / H
CalculatedMeasuredCalculatedMeasured
10.0130.0186.354.840.430.560.370.35
20.0120.0177.135.310.430.580.430.44
30.0150.0175.414.960.480.420.430.43
40.0140.026.574.300.400.440.390.44
50.0150.0195.634.590.450.540.390.39
Table 4. Parameters to predict the turbulence intensity in each case.
Table 4. Parameters to predict the turbulence intensity in each case.
Cases D b E b D i E i
12.221.442.061.82
22.211.462.131.87
32.221.732.241.50
42.151.002.211.83
52.311.812.161.51
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, J.; Wang, W.; Li, Z.; Li, Q.; Zhong, Y.; Xia, Z.; Qiu, H. Analytical Models of Velocity, Reynolds Stress and Turbulence Intensity in Ice-Covered Channels. Water 2021, 13, 1107. https://doi.org/10.3390/w13081107

AMA Style

Zhang J, Wang W, Li Z, Li Q, Zhong Y, Xia Z, Qiu H. Analytical Models of Velocity, Reynolds Stress and Turbulence Intensity in Ice-Covered Channels. Water. 2021; 13(8):1107. https://doi.org/10.3390/w13081107

Chicago/Turabian Style

Zhang, Jiao, Wen Wang, Zhanbin Li, Qian Li, Ya Zhong, Zhaohui Xia, and Hunan Qiu. 2021. "Analytical Models of Velocity, Reynolds Stress and Turbulence Intensity in Ice-Covered Channels" Water 13, no. 8: 1107. https://doi.org/10.3390/w13081107

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop