Next Article in Journal
Assumption-Free Derivation of the Bell-Type Criteria of Contextuality/Nonlocality
Previous Article in Journal
First Integrals of Shear-Free Fluids and Complexity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Entropy-Based Shear Stress Distribution in Open Channel for All Types of Flow Using Experimental Data

Department of Civil and Environmental Engineering, Pusan National University, Busan 46241, Korea
*
Author to whom correspondence should be addressed.
Entropy 2021, 23(11), 1540; https://doi.org/10.3390/e23111540
Submission received: 28 October 2021 / Accepted: 18 November 2021 / Published: 19 November 2021
(This article belongs to the Section Information Theory, Probability and Statistics)

Abstract

:
Korean river design standards set general design standards for rivers and river-related projects in Korea, which systematize the technologies and methods involved in river-related projects. This includes measurement methods for parts necessary for river design, but does not include information on shear stress. Shear stress is one of the factors necessary for river design and operation. Shear stress is one of the most important hydraulic factors used in the fields of water, especially for artificial channel design. Shear stress is calculated from the frictional force caused by viscosity and fluctuating fluid velocity. Current methods are based on past calculations, but factors such as boundary shear stress or energy gradient are difficult to actually measure or estimate. The point velocity throughout the entire cross-section is needed to calculate the velocity gradient. In other words, the current Korean river design standards use tractive force and critical tractive force instead of shear stress because it is more difficult to calculate the shear stress in the current method. However, it is difficult to calculate the exact value due to the limitations of the formula to obtain the river factor called the tractive force. In addition, tractive force has limitations that use an empirically identified base value for use in practice. This paper focuses on the modeling of shear-stress distribution in open channel turbulent flow using entropy theory. In addition, this study suggests a shear stress distribution formula, which can easily be used in practice after calculating the river-specific factor T. The tractive force and critical tractive force in the Korean river design standards should be modified by the shear stress obtained by the proposed shear stress distribution method. The present study therefore focuses on the modeling of shear stress distribution in an open channel turbulent flow using entropy theory. The shear stress distribution model is tested using a wide range of forty-two experimental runs collected from the literature. Then, an error analysis is performed to further evaluate the accuracy of the proposed model. The results reveal a correlation coefficient of approximately 0.95–0.99, indicating that the proposed method can estimate shear-stress distribution accurately. Based on this, the results of the distribution of shear stress after calculating the river-specific factors show a correlation coefficient of about 0.86 to 0.98, which suggests that the equation can be applied in practice.

1. Introduction

Understanding fluid interaction is very important in almost all studies of open channel flows. Shear stress is used and applied in hydraulics, hydrology, fluid mechanics, and in various fields, and is one of the most important mechanical factors [1]. It is a great challenge to the river engineers and researchers working in the field to estimate the distribution of bed shear stress in open channel flows [2].
Leighly [3] proposed that the bed shear stress can be balanced by the downstream component of the weight of water contained within the bounding orthogonal. Lundgren and Jonsson [4] modified the logarithmic law to a parabolic cross-sectional open channel and suggested a method to estimate the velocity and shear stress distribution. Chiu et al. [5,6] studied the complex interaction between primary and secondary flows, shear stress distribution, channel characteristics such as roughness, slope and geometry, and other related factors in open channel flows. However, the velocity profile was required to estimate boundary shear stress.
Keulegan [7] and Johnson [8] contributed to the early development of shear stress, and Einstein’s [9] hydraulic radius separation method is still used in various studies. Following this idea, Knight and his associates [10,11,12,13,14,15,16] proposed several empirical relations that are very helpful in understanding open channel flows and sediment transport. Noutsopoulos and Hadjipanos [17], Hu [18], and Patel [19] have led to an improved understanding of the lateral distributions of wall shear stress in rectangular channels, prismatic channels, and ducts. Past literature shows that the shear stress profile in open channel flows has been studied either experimentally or theoretically using deterministic approaches.
For these reasons, this study focuses on theoretical and statistic methods, which are probability and entropy concepts. Combining these, a useful method can be developed to explain the shear stress profile in open channel flows.
A probabilistic method of shear stress distribution in open channel flows using entropy concepts has been studied. Bonakdari et al. [20] compared Shannon and Tsallis entropies for shear stress distribution prediction in open channels. Sheikh and Wan [21] used the new Tsallis-based equation to predict shear stress distribution in circular and trapezoidal channels. Then, Mirauda and Maria [22] used an entropic parameter for modeling bed shear stress distribution in rectangular channels.
For the past few decades, entropy theory has been applied in the field of river hydraulic geometry and fluvial hydraulics. The entropy concept was introduced in hydraulics by Chiu [23] with Shannon’s entropy. Chiu studied the two-dimensional velocity distribution in an open channel [24,25]. Later, Choo [26] used Chiu’s velocity equation to calculate the momentum and energy coefficients. Especially, Chiu et al. [27] modified the entropy concept to be applied in pipe flows. Here, Chiu et al. compared with the Schlichting equation and the relationship between the frictional loss coefficient and entropy coefficient M, but there was no study on the shear stress. Singh [28] studied a wide range of hydrology and water resources based on entropy theory.
Since then, various entropy-related research has been carried out. Chiu et al. [29,30] applied the maximum velocity and regularity and a one-dimensional velocity distribution in an open channel with the entropy concept. Singh and Luo [31] examined the one-dimensional velocity distribution in an open channel with entropy theory, where they used Shannon’s entropy to derive the power law and logarithmic velocity distribution. Cui and Singh [32,33,34] studied velocity distribution and sediment concentration in open channels using Tsallis entropy. Singh and Cui [35] developed sediment concentration in a debris flow by Tsallis entropy.
Shear stress has to be estimated to determine flow characteristics. Considering the importance in an open channel, the complete evaluation is highly difficult due to the complexity of the cross-section and the various hydraulic parameters. The research related to the shear stress is based on empirical outcomes. That being so, it is difficult to apply the equations generally. Therefore, the current Korean river design standards use simple and obtainable tractive force and critical tractive force, which can be obtained through empirical methods instead of shear stress.
The objective of this study, therefore, is to model shear stress distribution using entropy theory, verify the model using twenty-one experimental datasets obtained from Song’s experimental data [36], and to prove the utility of the proposed equation. Based on the proven equation, it is then suggested that the tractive force and critical tractive force in the Korean river design standards be revised to the shear stress obtained using the proposed shear stress distribution formula by presenting a method that can be easily used in the actual conditions.

2. Methodology

2.1. Entropy Theory

The variable x related to information I ( x ) is shown as Equation (1), which provides the amount of information. Here, information I ( x ) is the measurement of uncertainty related to a certain state as:
I ( x ) = l n p ( x )
Considering every state, the average value of information I ( x ) can be expressed as Equation (2). Function H ( x ) is defined using Shannon’s [37] entropy as:
H ( x ) = + p ( x ) I ( x ) d x = + p ( x ) l n p ( x ) d x
where p ( x ) and l n p ( x ) are dimensionless, but d x has a dimension; thus, H ( x ) has the same dimension as d x . In Equation (2), a probability density function (PDF) p ( x ) of a continuous variation state means maximizing the entropy of uncertainty x .

2.2. Constraint Conditions

The probability distribution of maximizing the entropy produces more information from already acquired basic knowledge. To solve the PDF, i.e., the available information of variable x in Equation (3), the constraint conditions, such as average, variance, distortion, etc., are applied as:
a b Φ i ( x , p ) d x       i = 1 ,   2 ,   3 ,   n
Therefore, the PDF p ( x ) , which maximizes the entropy, can be obtained using the method from Lagrange in Equation (4); we have:
I ( x , p ) p + i = 1 n λ i ϕ i ( x , p ) p = 0
where λ i denotes the Lagrange multipliers.

2.3. Entropy Maximization

The entropy concept can be applied to the shear stress using Equation (5) from Shannon’s entropy as:
H ( τ ) = 0 τ 0 p ( τ ) I ( τ ) d τ = 0 τ 0 p ( τ ) l n p ( τ ) d τ
where τ is point shear stress.
The available information for τ uses constraint conditions. First, the total probability must be satisfied for the PDF p ( τ ) as:
0 τ 0 p ( τ ) d τ = 1
which follows from the total probability rule.
Then, average information can be expressed as:
0 τ 0 τ · p ( τ ) d τ = τ ¯

2.4. Lagrange Method

Arranging the independent constraint conditions can be given as:
a b Φ i ( τ , p ) d τ       i = 1 ,   2
Therefore, PDF p ( τ ) , which maximizes the entropy, can be obtained using the method of Lagrange as:
I ( τ , p ) p + i = 1 2 λ i ϕ i ( τ , p ) p = 0
I ( τ , p ) = p ( τ ) l n p ( τ )
where ϕ 1 ( τ , p ) = p ( τ ) ,   ϕ 2 ( τ , p ) = τ · p ( τ ) .
ϕ 1 ( τ , p ) p = 1 ,   ϕ 2 ( τ , p ) p = τ
Substituting Equation (9) into Equations (10) and (11) can be constructed as follows:
1 l n p ( τ ) + λ 1 + λ 2 τ = 0
where λ 1 1 = a 1 ,   λ 2 = a 2 are the Lagrange multipliers and differentiating Equation (12) with respect to p ( τ ) results in the shear stress PDF as:
p ( τ ) = e x p a 1 + a 2 τ
We apply the cumulative probability function to fluid flow using Equation (14), and Equation (13) will become Equation (14) by applying the PDF as:
F ( τ ) = 0 τ p ( τ ) d τ = 0 τ e a 1 + a 2 τ d τ = 1 [ ξ ξ 0 ξ m a x ξ 0 ]
where ξ denotes the spatial coordinates ( 0 ξ 1 ), τ is the point shear stress at ξ , ξ 0 is the minimum value of ξ (occurring at the channel boundary where τ = 0 ), and ξ m a x is the maximum value of ξ (where τ is at its maximum) (see Figure 1).
The ξ - η coordinates are the isovel system, which was first developed by Chiu [25,38] to explain two-dimensional velocity distribution in the cross-section of an open channel.

2.5. Proposed Shear Stress Distribution Model

Solving Equation (14) can be written as:
τ = 1 a 2 l n [ 1 + a 2 e a 1 ( 1 [ ξ ξ 0 ξ m a x ξ 0 ] ) ]
Open channel flow also uses the constraint conditions of Equations (6) and (7). Substituting Equation (13) into Equation (6) is given as:
a 2 e a 1 = ( e a 2 τ 0 1 ) e a 1 = T ( e T 1 ) τ 0
where T = a 2 τ 0 (normally called the entropy coefficient) to represent the model in a simple-to-use form and substituting Equation (13) into Equation (7) can be expressed as:
τ ¯ = 0 τ 0 τ · p ( τ ) d τ = 0 τ 0 τ e a 1 + a 2 τ d τ τ ¯ = [ e T e T 1 1 T ] τ 0 = ϕ ( T ) τ 0
where τ 0 is the boundary shear stress, T is the entropy coefficient, and ϕ ( T ) is a function of T . Equation (17) is one of the proposed shear stress distribution equations for open channel flows. Substituting Equation (17) into Equation (16) and rearranging it in respect of the e a 1 term gives:
e a 1 = T ( e T 1 ) τ 0 = T · ϕ ( T ) ( e T 1 ) τ ¯
Substituting Equation (18) into Equation (15) and rearranging it results in another proposed mean shear stress distribution formula as:
τ = τ ¯ ( e T 1 ) ( T e T e T + 1 ) l n [ 1 + ( e T 1 ) ( 1 [ ξ ξ 0 ξ m a x ξ 0 ] ) ]
By inserting Equation (16) into Equation (15) we have:
τ = τ 0 T l n [ 1 + ( e T 1 ) ( 1 [ ξ ξ 0 ξ m a x ξ 0 ] ) ]
Equation (20) is the last proposed model for boundary shear stress distribution. Generally, ξ 0 is close to 0, ξ m a x is 1, and ξ for Equations (19) and (20) is the same, which was formulated as follows:
ξ = y D h e x p ( 1 y D h )
where D is the maximum depth, h is the depth where shear stress is 0 from the water surface (maximum velocity also occurs at this location), and y is the vertical depth from the bed for a given shear stress.

2.6. Shear Stress in Fluid

Normally when laminar and turbulent flows coincide, they can be written as:
τ = ρ ν d u d y ρ u v ¯ = ρ ( ν + ε ) d u d y
where ρ is the fluid density, ν is the fluid viscosity, d u / d y is the velocity gradient, u v ¯ is the Reynolds stress, and ε is the eddy viscosity. Reynolds stress is the shear stress caused by turbulent fluctuating velocity. The kinematic coefficient of viscosity, ν , is caused by the molecular motion of fluid, and ε y is caused by fluid particle mixing, which is much larger than molecular motion.
Boundary shear stress from Equation (22) can be expressed as follows: when there is shear stress at the bed, y = 0 , the velocity is 0, and u = 0 . This is given by:
τ 0 = ρ ( ν + ε ) [ d u d y ] y = 0 = ρ g R h S f
where g is the gravitational acceleration, R h is the shape form of the cross-section, and S f is the energy gradient.

2.7. Tractive Force Formula in Fluid

Tractive force means the running water force when the silt on the river bed is moved by water. The commonly used tractive force ( τ 0 ) formula is:
τ 0 = ω R i
where ω is the unit weight of water, i is the bed slope, and R is the hydraulic radius.

2.8. Critical Tractive Force Formula in Fluid

Critical tractive force means the tractive force at the beginning of the movement of the river bed silt due to the fact that the running water force is greater than the resistance of the river bed. The commonly used critical tractive force formula ( F s ) is:
F s = τ 0 ( ρ s ρ ) g d = u 2 * c 1 ρ ( ρ s ρ ) g d
where u * c is critical friction velocity, ρ s is the density of silt particles, ρ is the density of water, and g is gravitational acceleration.
In terms of simplicity and convenience, Equations (22) and (23) have an advantage, but their accuracy is suspect. The reason for this is that the energy gradient is actually a difficult factor to estimate. As can be seen in Equation (22), the measured point velocity of the whole cross-section is required for shear stress to reach each gradient. In other words, the velocity gradient ( d u / d y ), eddy viscosity coefficient ( ε ), and the energy gradient ( S f ) in Equations (22) and (23) are factors that are very difficult to estimate. In addition, as shown by the Equation (24), it is difficult to accurately calculate the tractive force used by the Korean river design standards due to the hard-to-find river factors such as the river bed. Critical tractive force is readily calculated, but is not certain in terms of accuracy because the formulas are empirical formulas obtained from experiments. Therefore, this paper suggests an equation that can express the shear stress distribution and boundary shear stress in open channel turbulent flow using entropy-based modeling. This study demonstrates the utility of the proposed equation by using the Song data. It also proposes to revise the tractive force and critical tractive force of the Korean river design standards to shear stress by presenting measures easily applicable in practice.

3. Experimental Data

The proposed model of shear stress distribution was validated with experimental observations available in the literature. To test the validity of the model, i.e., Equations (19) and (20) with a wide range of slope, discharge, and sediment flow conditions, experimental data from Song [36] were selected (see Table 1). Forty-six flows were used in this study: twenty-one uniform flows, twenty-one non-uniform flows, and four unsteady flows according to four slope conditions. Out of twenty-one uniform flows, six runs were experimented under sediment conditions. For non-uniform flow, twelve accelerating flows and nine decelerating flows were tested. Four unsteady flows were tested according to four slope conditions. This study considered a various range of experimental runs for verification of the shear stress distribution.

4. Parameter Estimation and Comparison with Experimental Data

4.1. Parameter Estimation

Proposed shear stress distribution was used to estimate the entropy parameter T . First, we estimate parameter T by inserting experimental point shear stress ( τ 1 ,   τ 2 ,   ,   τ n ) and point vertical depth from the bed for a given shear stress ( y 1 ,   y 2 ,   ,   y n ) into Equations (19) and (20). We estimate the best boundary shear stress τ 0 and mean shear stress τ ¯ value, which has the least error for each run. We then use τ 0 and τ ¯ from Equations (19) and (20) to estimate ϕ ( T ) from Equation (17). Lastly, we calculate the shear stress distribution for given vertical depths ( y 1 ,   y 2 ,   ,   y n ) . The parameter estimation process is shown in Figure 2.

4.2. Comparison with Experimental Data

Figure 3, Figure 4, Figure 5 and Figure 6 compare the proposed model with the experimental data of Song [36] to determine if the estimated shear stress distribution fits well with the observed shear stress distribution. Each diagram shows shear stress at the x-axis and y / D at the y-axis with the correlation coefficient of observed data and estimated shear stress.
For first (uniform), second (accelerating non-uniform), third (decelerating non-uniform) and fourth (unsteady) flow conditions (Figure 3, Figure 4, Figure 5 and Figure 6), the proposed model was applied to compare the estimated and observed values of shear stress distribution to see how well they were expressed.
In Figure 3, the proposed shear stress model showed a good agreement with experimental data despite the scattered nature of the data. As for the six sediment runs, the correlation coefficient was actually better in some cases compared to small sediment runs. For uniform flows, the entropy-based model seems to estimate shear stress accurately. In all uniform flows, the correlation coefficient showed a small range from 0.9375 to 0.9931.
From Figure 4, it is found that the proposed model simulated well even in more complex flows. For accelerating non-uniform flows, the entropy-based model seems to estimate shear stress accurately. In all accelerating non-uniform flows, the correlation coefficient showed a small range from 0.9522 to 0.9959. However, some flows have shown very scattered experimental data, especially in the bottom layer due to the difficulties in measuring.
From Figure 5, it can be seen that the entropy-based model simulated very well. For decelerating non-uniform flows, the proposed model seems to compute shear stress accurately. In all decelerating non-uniform flows, the correlation coefficient showed a small range from 0.9475 to 0.9822. The proposed model seems to express very accurate matching results.
From Figure 6, it can be seen that the entropy-based model simulated well. For unsteady flows, the proposed model appears to compute shear stress well. In an unsteady flow (S-25-931), the correlation coefficient was between 0.8262 and 0.9843. The results for the other three slopes are as follows. The correlation coefficient in S-60-933 was between 0.9034 and 0.952, the correlation coefficient in S-10-31 was between 0.8536 and 0.9899, and the correlation coefficient in S-30-932 was between 0.712 and 0.9829. In four cases, the correlation coefficient was above 0.89 on average. This demonstrates the utility of the shear-stress distribution equation using entropy.

4.3. Major Parameter Estimation Results

Table 2 shows the result of major parameters from Figure 3, Figure 4, Figure 5 and Figure 6. The entropy-based model showed a 0.9375 to 0.9959 range of correlation coefficients in steady flow conditions and a 0.712 to 0.9899 range of correlation coefficients in unsteady flow conditions. From these results, it seems that the number of measured shear stresses in one distribution does not have a large effect on simulation. Entropy parameter T seems to have a range of −1.441 to 6.405 in steady flow conditions and a range of −5.576 to 4.6124 in unsteady flow conditions. Looking at the correlation coefficients, the decelerating non-uniform flow showed the worst average, 0.9693, whereas accelerating non-uniform flow showed the best average, 0.9887 in steady flow conditions. For unsteady flow conditions, S30-932 showed the worst average, 0.89353, whereas S-10-932 showed the best average, 0.93602.

5. Proposal and Verification of Shear Stress Distribution Method

5.1. Easily Applied Shear Stress Distribution Formula for Practice

If the main parameters of Table 2 are used in expressions (17) to graph, the distribution of mean shear stress and floor shear stress, ϕ ( T ) , can be obtained based on the slope values of the graph. Using the value of the obtained ϕ ( T ) , it is possible to calculate the river-specific factor ( T ) in Equation (17). Here, T is the average value of T . The entropy parameter T can be obtained from previous equations and represents one cross-section in a river. However, there are various cross-sections in a river which have many values of T , so T is used for a representative factor for one river.
Therefore, in other words, if the river-specific factor T is calculated or known, it is much easier to calculate the boundary shear stress, which is an important hydraulic factor in river design. For example, when the river-specific factor T is used in the shear stress distribution Equations (19) and (20), less calculation and more accurate shear stress distribution can be obtained, because one of the parameters is already known.
Song data were used to prove the utility of shear stress distribution when the river-specific factor T was fixed. The river-specific factor T for each flow state was then obtained by using the floor shear stress and average shear stress, which are the parameters of Figure 7 and Figure 8. Shear stress was obtained by substituting T for Equations (19) and (20) and the actual measurement value of the Song data was compared and analyzed.

5.2. Estimation Graph of River-Specific Factors T in All Flow Conditions

Using Equation (17) based on the data, is the results were displayed graphically to calculate the river-specific factor for each flow state. Figure 7 and Figure 8 show graphs for the calculating of river-specific factors in each flow state.
Figure 7 is a graph of calculating equilibrium T after calculating equilibrium ϕ ( T ) in uniform, non-uniform accelerating, and non-uniform decelerating flow. For each flow, equilibrium ϕ ( T ) is 0.6490 (uniform), 0.5344 (non-uniform accelerating), and 0.7466 (non-uniform decelerating); and when Equation (17) is calculated using equilibrium ϕ ( T ) , equilibrium T is calculated as 1.88 (uniform), 0.41 (non-uniform accelerating), and 3.52 (non-uniform decelerating).
Unsteady flow was divided into four cases according to the gradient, and equilibrium T was calculated after selecting equilibrium ϕ ( T ) . Figure 8 is a graph that calculates equilibrium T after obtaining equilibrium ϕ ( T ) from the unsteady flow. For each case in the unsteady flow, equilibrium ϕ ( T ) is 0.3642 (S-60-933), 0.399 (S-25-931), 0.408 (S-10-932), and 0.562 (S30-932); and when Equation (17) is calculated using equilibrium ϕ ( T ) , equilibrium T is calculated as −1.70208 (S-60-933), 1.24195 (S-25-931), −1.1267 (S-10-932), and 0.74953 (S30-932).
We now put the equilibrium T values into Equations (19) and (20) and verify the utility of the equilibrium T -fixed Equations (19) and (20) through the Song data. If the effectiveness of the equation is verified, the T -fixed method easily enables the calculation of the floor shear stress by average shear stress, and the distribution of shear stress can be easily calculated.

5.3. Result and Error Analysis

The RMSE is a measure of the residual, which is the difference between the values simulated by the model and actual observed values. The RMSE enables simulative power to be integrated into a single unit of measurement. The RMSE of the model’s simulation for the estimated variable X e s t , i is defined as the square root of the mean square error Equation (26).
RMSE = i = 1 n ( X o b s , i X e s t , i ) 2 n
where X o b s , i indicates the actual observed value and X e s t , i is the simulated value obtained from the model.
The determined shear-stress distribution was compared with those measured by Song. The determined RMSE values were very low from 0.000245 to 0.001108 as shown in Table 3.
Basically, the steady flow shows a high correlation coefficient value and proves the utility of the shear-stress distribution equation. However, the unsteady flow shows a lower mean correlation value compared to the steady flow. Looking at the results of the unsteady flow in detail, note that not all areas are observed with low correlation, but with high and low correlation. This indicates that the low correlation coefficient value from the unsteady flow is the error caused by the observation of the data. In practice, using point shear stress values calculated based on accurate observations can reduce error and have higher correlation coefficient values.
Based on the advantages of being able to express and easily obtain shear stress distribution by entropy shear stress distribution using equilibrium T , it is suggested that the contents of tractive force and critical tractive force with unclear accuracy can be revised to shear stress in Korea’s river design standards.

6. Discussion

The proposed equation has proved to show reasonable results. From Figure 3, Figure 4, Figure 5 and Figure 6, in most cases, there is a good match between Song’s observed data and shear stress values from the proposed model. It was confirmed that shear stress distributions estimated are accurate showing an average 0.9780 correlation coefficient for forty-two types of steady open channel flows. In four cases where the flow was unsteady, the mean value of the correlation coefficient was found to be 0.9084. The estimated values are not from an empirical formula but from a theoretical method, which has a great meaning for open channel fields.
Using the estimated boundary shear stress and mean shear stress from all forty-two runs, entropy parameter T was analyzed as 1.629. Figure 7 and Figure 8 were plotted using the relationship between bed shear stress and mean shear stress, Equation (17), which seems to have a tendency for an equilibrium state. The equilibrium state of a velocity distribution was studied previously by Chiu et al. [38,39], which shows similar results.
The aforementioned discussion delineates that an entropy-based model on shear stress is able to describe the characteristics of shear stress from the frictional force caused by viscosity and fluctuating fluid velocity in an open channel turbulent flow. The proposed equation has proved to show reasonable results.
Based on the parameters calculated, the mean values of the correlation coefficient in the results of the shear stress distribution equation for each flow were 0.9708 for the uniform flow, 0.9867 for the non-uniform flow (accelerating), 0.9471 for the non-uniform flow(decelerating), and 0.8604–0.9212 in the unsteady flow for each gradient. Based on these results, it was confirmed that the shear stress distribution formula using river-specific factors can be available in practice.
Of course, part of the unsteady flow shows lower correlation coefficient values, but the correlation coefficient values are error generated from observations. When measuring the distribution of point shear stress for practical use, the conclusions based on accurate observations are determined to reduce the error further and obtain high accuracy shear stress distributions and results.
However, there are some minor limitations to this study. Some observed data, such as shear stress distribution and vertical depth from the bed, should be known in order to estimate the model. In other words, there needs to be some information in order to use this model. In addition, the basic shape of the model depends on the PDF of shear stress and Equation (21). This is important because in complex flows such as unsteady flows, it can be difficult to show the distribution exactly, especially if the maximum velocity occurs below the water surface. The model is based on probability and statistics and does not consider basic hydraulics in the beginning, but only depends on constraints and data. Because the accuracy of the observed data can affect the correlation value of the distribution equation, accurate point shear stress data must be obtained for use in practice.

7. Conclusions and Proposition

Although the tractive force and critical tractive force formulas used in the current Korean river design standards are simple to obtain, the tractive force is a formula using a factor that is difficult to obtain, which makes it difficult to calculate the exact value. In addition, critical tractive force is a value obtained empirically, which also has the disadvantage of difficulty in calculating the exact value. However, in order to use the existing shear stress formula, shear stress can be calculated only with the energy gradient, a factor that is very difficult to obtain, so Korea’s river design standards use tractive force and critical tractive force instead of shear stress.
Therefore, using the entropy concept, this study has proposed shear stress distribution and boundary layer shear stress, which can be applied in open channel flows. To determine how well the shear stress model fits with the observed data, Song’s data were used. From this aspect, this paper validated the model with a wide range of forty-two runs of experimental data published in the literature. The results show the utility and reliability of Equations (19) and (20), in which the mean shear stress is considered, then using the qualified shear stress distribution, the shear stress distribution over the whole flow depth of open channel turbulent flows is shown.
Furthermore, based on the equations’ proven utility and reliability, this paper has proposed a method that can easily be used in practice by obtaining river-specific factor T if point shear stress is given. We suggest that the contents of tractive force and critical tractive forces with shortcomings in terms of accuracy in the Korean river design standards are revised to the shear stress distribution method presented in this study. The results of the distribution of shear stress after the calculation of river-specific factor T were also considered to be highly correlated, and thus it was determined that the method could be used in practice.
However, there is no best entropy parameter T value for every run. This means that we need to somewhat reach an agreement with solutions. To improve this, very precise methods must be applied in order to enhance the results. One of these methods can be an optimum technique. Optimization techniques such as genetic algorithms or harmony searching can be applied for finding better solutions to these kinds of problems, where we might even be getting closer to what we are seeking.
There have been very few studies about shear stress in open channels recently. It seems that there are only few models that can estimate shear stress distribution. Some of the models require parameters that are difficult to calculate or even obtain, or even have hypotheses that eventually reach the limit of the formula. Other models have many parameters to estimate, but the features of the proposed model have only two parameters, the entropy parameter (T) and the unknown boundary shear stress or mean shear stress.
In addition, if the river-specific factors T are calculated for practical use based on the proposed model, the boundary shear stress, which is an important river factor for river design, can be calculated immediately, and shear stress distribution can easily be calculated. The proposed model can be used regardless of the shape or flow of the river (except for unsteady flow). Nevertheless, in future research, our model has to be compared with some shear stress models.
The analysis has limitations but the results appear to be useful. If the point shear stress in the open channel and vertical depth from the bed are given, the shear stress distribution can be estimated simply from the model, which will show high availability when designing or managing the open channel. In addition, the boundary shear stress can be estimated easily without the energy gradient when calculating the boundary shear stress in an open channel.

Author Contributions

Y.-M.C. and J.-C.S. conceptualized the research. Y.-M.C., H.-S.J. and J.-C.S. carried out the investigation. Y.-M.C. wrote the manuscript. J.-C.S. conceived the methodology of the proposed method. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Institute of Industrial Technology, Pusan National University.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This research was supported by the Institute of Industrial Technology, Pusan National University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Day, M.A. The no-slip condition of fluid dynamics. Erkenntnis 1990, 33, 285–296. [Google Scholar] [CrossRef]
  2. Patra, K.; Sahoo, N.; Khatua, K.K. Distribution of boundary shear in compound channel with rough floodplains. River Basin Manag. VII 2013, 172, 99–110. [Google Scholar] [CrossRef] [Green Version]
  3. Leighly, J.B. Toward a Theory of the Morphologic Significance of Turbulence in the Flow of Water in Streams; University of California Press: Berkeley, CA, USA, 1932; Volume 6, pp. 1–9. [Google Scholar]
  4. Lundgren, H.; Jonsson, I.G. Shear and Velocity Distribution in Shallow Channels. J. Hydraul. Div. 1964, 90, 1–21. [Google Scholar] [CrossRef]
  5. Chiu, C.; Lin, G. Computation of 3-D Flow and Shear in Open Channels. J. Hydraul. Eng. 1983, 109, 1424–1440. [Google Scholar] [CrossRef]
  6. Chiu, C.; Chiou, J. Structure of 3-D Flow in Rectangular Open Channels. J. Hydraul. Eng. 1986, 112, 1050–1067. [Google Scholar] [CrossRef]
  7. Keulegan, G.H. Laws of turbulent flow in open channels. J. Res. Natl. Inst. Stand. Technol. 1938, 21, 707. [Google Scholar] [CrossRef]
  8. Johnson, J.W. The Importance of Side Wall Friction in Bed-load Investigation. Proc. Civ. Eng. USA 1942, 12, 331–339. [Google Scholar]
  9. Einstein, H.A. Formulas for the Transportation of Bed Load. Trans. Am. Soc. Civ. Eng. 1942, 107, 561–577. [Google Scholar] [CrossRef]
  10. Ghosh, S.; Knight, D. Discussion. Boundary Shear Distribution in Channels with Varying Wall Roughness. Proc. Inst. Civ. Eng. 1973, 55, 503–511. [Google Scholar] [CrossRef] [Green Version]
  11. Knight, D.W.; Macdonald, J.A. Hydraulic Resistance of Artificial Strip Roughness. J. Hydraul. Div. 1979, 105, 675–690. [Google Scholar] [CrossRef]
  12. Knight, D.W.; Macdonald, J.A. Open Channel Flow with Varying Bed Roughness. J. Hydraul. Div. 1979, 105, 1167–1183. [Google Scholar] [CrossRef]
  13. Knight, D.W. Boundary Shear in Smooth and Rough Channels. J. Hydraul. Div. 1981, 107, 839–851. [Google Scholar] [CrossRef]
  14. Knight, D.; Patel, H.; Demetriou, J.; Hamed, M. Boundary shear stress distributions in open channel and closed conduit flows. In Mechanics of Sediment Transport; CRC Press: London, UK, 2020; Volume 156, pp. 33–40. [Google Scholar] [CrossRef]
  15. Knight, D.W.; Demetriou, J.D. Flood Plain and Main Channel Flow Interaction. J. Hydraul. Eng. 1983, 109, 1073–1092. [Google Scholar] [CrossRef]
  16. Knight, D.W.; Patel, H.S. Boundary Shear in Smooth Rectangular Ducts. J. Hydraul. Eng. 1985, 111, 29–47. [Google Scholar] [CrossRef]
  17. Noutsopoulos, G.C.; Hadjipanos, P.A. Discussion of “Boundary Shear in Smooth and Rough Channels”. J. Hydraul. Div. 1982, 108, 809–812. [Google Scholar] [CrossRef]
  18. Hu, C.H. The Effect of the Width/Depth Ratio and Side-Wall Roughness on Velocity Profile and Resistance in Rectangular Open-Channels. Ph.D. Thesis, Tsinghua University, Beijing, China, 1985. [Google Scholar]
  19. Patel, H.S. Boundary Shear in Rectangular and Compound Ducts. Ph.D. Thesis, Birmingham University, Birmingham, UK, 1984. Available online: https://ethos.bl.uk/OrderDetails.do?uin=uk.bl.ethos.566061 (accessed on 1 June 2021).
  20. Bonakdari, H.; Sheikh, Z.; Tooshmalani, M. Comparison between Shannon and Tsallis entropies for prediction of shear stress distribution in open channels. Stoch. Environ. Res. Risk Assess. 2014, 29, 1–11. [Google Scholar] [CrossRef]
  21. Sheikh Khozani, Z.; Wan Mohtar, W.H.M. Investigation of New Tsallis-Based Equation to Predict Shear Stress Distribution in Circular and Trapezoidal Channels. Entropy 2019, 21, 1046. [Google Scholar] [CrossRef] [Green Version]
  22. Mirauda, D.; Russo, M.G. Modeling Bed Shear Stress Distribution in Rectangular Channels Using the Entropic Parameter. Entropy 2020, 22, 87. [Google Scholar] [CrossRef] [Green Version]
  23. Chiu, C.-L. Entropy and Probability Concepts in Hydraulics. J. Hydraul. Eng. 1987, 113, 583–599. [Google Scholar] [CrossRef]
  24. Chiu, C. Entropy and 2-D Velocity Distribution in Open Channels. J. Hydraul. Eng. 1988, 114, 738–756. [Google Scholar] [CrossRef]
  25. Chiu, C.-L. Application of Entropy Concept in Open-Channel Flow Study. J. Hydraul. Eng. 1991, 117, 615–628. [Google Scholar] [CrossRef]
  26. Choo, T.H. Estimation of Energy Momentum Coefficients in Open Channel Flow by Chiu’s Velocity Distribution Equation. Master’s Thesis, Pittsburgh University, Pittsburgh, PA, USA, 1990. [Google Scholar]
  27. Chiu, C.; Lin, G.; Lu, J. Application of Probability and Entropy Concepts in Pipe? Flow Study. J. Hydraul. Eng. 1993, 119, 742–756. [Google Scholar] [CrossRef]
  28. Singh, V.P. The use of entropy in hydrology and water resources. Hydrol. Process. 1997, 11, 587–626. [Google Scholar] [CrossRef]
  29. Chiu, C.-L.; Tung, N.-C. Maximum Velocity and Regularities in Open-Channel Flow. J. Hydraul. Eng. 2002, 128, 390–398. [Google Scholar] [CrossRef]
  30. Chiu, C.-L.; Hsu, S.-M. Probabilistic approach to modeling of velocity distributions in fluid flows. J. Hydrol. 2006, 316, 28–42. [Google Scholar] [CrossRef]
  31. Singh, V.P.; Luo, H. Entropy Theory for Distribution of One-Dimensional Velocity in Open Channels. J. Hydrol. Eng. 2011, 16, 725–735. [Google Scholar] [CrossRef]
  32. Cui, H.; Singh, V.P. Two-Dimensional Velocity Distribution in Open Channels Using the Tsallis Entropy. J. Hydrol. Eng. 2013, 18, 331–339. [Google Scholar] [CrossRef]
  33. Cui, H.; Singh, V.P. One-Dimensional Velocity Distribution in Open Channels Using Tsallis Entropy. J. Hydrol. Eng. 2014, 19, 290–298. [Google Scholar] [CrossRef]
  34. Cui, H.; Singh, V.P. Suspended Sediment Concentration in Open Channels Using Tsallis Entropy. J. Hydrol. Eng. 2014, 19, 966–977. [Google Scholar] [CrossRef]
  35. Singh, V.P.; Cui, H. Modeling sediment concentration in debris flow by Tsallis entropy. Phys. A Stat. Mech. Its Appl. 2015, 420, 49–58. [Google Scholar] [CrossRef]
  36. Song, T. Velocity and turbulence distribution in non-uniform and unsteady open-channel flow. EPFL 1995, 4, 78. [Google Scholar] [CrossRef]
  37. Shannon, C.E. The mathematical theory of communications, I and II. Bell Syst. Tech. J. 1948, 27, 379–423. [Google Scholar] [CrossRef] [Green Version]
  38. Chiu, C.-L. Velocity Distribution in Open Channel Flow. J. Hydraul. Eng. 1989, 115, 576–594. [Google Scholar] [CrossRef]
  39. Chiu, C.-L.; Said, C.A.A. Maximum and Mean Velocities and Entropy in Open-Channel Flow. J. Hydraul. Eng. 1995, 121, 26–35. [Google Scholar] [CrossRef]
Figure 1. ξ - η coordinates in an open channel flow (Chiu [25,38]).
Figure 1. ξ - η coordinates in an open channel flow (Chiu [25,38]).
Entropy 23 01540 g001aEntropy 23 01540 g001b
Figure 2. Parameter estimation flowchart.
Figure 2. Parameter estimation flowchart.
Entropy 23 01540 g002
Figure 3. Verification of the proposed shear stress distribution model with six uniform flows.
Figure 3. Verification of the proposed shear stress distribution model with six uniform flows.
Entropy 23 01540 g003
Figure 4. Verification of the proposed shear-stress distribution model with six accelerating non-uniform flows.
Figure 4. Verification of the proposed shear-stress distribution model with six accelerating non-uniform flows.
Entropy 23 01540 g004
Figure 5. Verification of the proposed shear-stress distribution model with six decelerating non-uniform flows.
Figure 5. Verification of the proposed shear-stress distribution model with six decelerating non-uniform flows.
Entropy 23 01540 g005aEntropy 23 01540 g005b
Figure 6. Verification of the proposed shear-stress distribution model with unsteady flow of one slope (S-25-931).
Figure 6. Verification of the proposed shear-stress distribution model with unsteady flow of one slope (S-25-931).
Entropy 23 01540 g006
Figure 7. Equilibrium ( T ) in steady flow.
Figure 7. Equilibrium ( T ) in steady flow.
Entropy 23 01540 g007
Figure 8. Equilibrium ( T ) in unsteady flow.
Figure 8. Equilibrium ( T ) in unsteady flow.
Entropy 23 01540 g008
Table 1. Experimental data.
Table 1. Experimental data.
Flow TypeSlope
(%)
Discharge
(cm/s)
Sediment
(cm)
UniformS25-Q400.254012.2
S50-Q500.55011.8
S75-Q1200.7512018.5
S90-Q700.907012.2
S100-Q551.005510.5
S125-Q501.25509.1
Non-uniform
(accelerating)
AS00-Q8008014.5
AS-25-Q1000.2510016.9
AS-25-Q800.258015.1
AS-50-Q1100.5011017.9
AS-75-Q800.758016.5
AS-93-Q800.938017.2
Non-uniform
(deaccelerating)
DS25-Q700.257016.0
DS25-Q900.259020.0
DS50-Q700.507016.5
DS50-Q900.509018.5
DS75-Q800.758020.5
DS90-Q700.907018.0
UnsteadyS-25-9310.631–0.6974.056.26–7.04
S-60-9330.564–0.6127.015.11–6.14
S-10-9320.642–0.7035.745.94–7.63
S30-9320.821–0.9445.706.34–7.66
Table 2. Summary of major parameter results for steady (right) and unsteady (left) flow conditions.
Table 2. Summary of major parameter results for steady (right) and unsteady (left) flow conditions.
Data Set τ 0 τ ¯ T R 2 Data Set τ 0 τ ¯ T R 2
S-60-933(t1)69.577312.2146−5.57600.9491S25-Q3120.53114.4582.7420.9375
S-60-933(t3)58.262311.4943−4.87930.9340S25-Q4021.51916.1343.5880.9834
S-60-933(t57)47.918013.5993−2.94910.9512S150-Q40102.85058.9270.8870.9809
S-60-933(t59)51.224812.9439−3.53830.9520S150-Q50123.53578.0791.6560.9677
S-10-932(t2)39.577523.74301.22880.9255AS00-Q14592.14452.4840.8450.9957
S-10-932(t6)39.212025.30551.84040.9176AS00-Q10071.48141.4640.9760.9934
S-10-932(t54)41.615516.2094−1.36670.9867AS-93-Q10041.03623.8460.9890.9829
S-10-932(t78)37.957420.70720.54920.9267AS-93-Q8038.19419.6740.1820.9902
S-25-931(t2)28.011121.04773.62270.8393DS25-Q9027.40819.0212.5730.9769
S-25-931(t58)37.655613.5880−1.750.9800DS25-Q7023.98317.5533.2360.9822
S-25-931(t134)32.233722.90282.85010.8262DS90-Q8048.08740.6596.4050.9734
S-25-931(t138)30.191622.26383.34150.8382DS90-Q7033.71127.7515.5320.9683
S-30-932(t1)39.660523.97801.28930.9098
S-30-932(t17)44.397833.93483.90190.7628
S-30-932(t57)44.830520.6737−0.46790.9607
S-30-932(t59)48.567617.8320−1.66650.9461
Table 3. Simulation results R2 and RMSE.
Table 3. Simulation results R2 and RMSE.
Flow TypeR2RMSE
UniformS25-Q400.98340.000491
S50-Q500.98480.000431
S75-Q1200.99080.000394
S90-Q700.98120.000424
S100-Q550.97210.000564
S125-Q500.98710.000417
Non-uniform
(accelerating)
AS00-Q800.99590.000284
AS-25-Q1000.99330.000304
AS-25-Q800.99140.000346
AS-50-Q1100.99570.000245
AS-75-Q800.99310.000364
AS-93-Q800.99020.000384
Non-uniform
(deaccelerating)
DS25-Q700.98220.000436
DS25-Q900.97690.000513
DS50-Q700.96700.000685
DS50-Q900.96900.000613
DS75-Q800.98140.000531
DS90-Q700.96830.000678
UnsteadyS-25-931 2 s0.83930.001108
S-25-931 26 s0.86800.001085
S-25-931 42 s0.93820.000712
S-25-931 50 s0.96330.000675
S-25-931 66 s0.96430.000645
S-25-931 98 s0.86660.000907
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Choo, Y.-M.; Jeon, H.-S.; Seo, J.-C. Entropy-Based Shear Stress Distribution in Open Channel for All Types of Flow Using Experimental Data. Entropy 2021, 23, 1540. https://doi.org/10.3390/e23111540

AMA Style

Choo Y-M, Jeon H-S, Seo J-C. Entropy-Based Shear Stress Distribution in Open Channel for All Types of Flow Using Experimental Data. Entropy. 2021; 23(11):1540. https://doi.org/10.3390/e23111540

Chicago/Turabian Style

Choo, Yeon-Moon, Hae-Seong Jeon, and Jong-Cheol Seo. 2021. "Entropy-Based Shear Stress Distribution in Open Channel for All Types of Flow Using Experimental Data" Entropy 23, no. 11: 1540. https://doi.org/10.3390/e23111540

APA Style

Choo, Y. -M., Jeon, H. -S., & Seo, J. -C. (2021). Entropy-Based Shear Stress Distribution in Open Channel for All Types of Flow Using Experimental Data. Entropy, 23(11), 1540. https://doi.org/10.3390/e23111540

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