Proposed Equation to Estimate the Length of a Hydraulic Jump in a Rectangular Open Channel Hydraulic with Variable Slope, and Its Comparison with Seven Theoretical Equations of Specialized Literature ()
1. Introduction
The hydraulic jump is a localized phenomenon that occurs in open channels hydraulics (OCH). Giorgio Bidone, an Italian investigator, studied it for the first time in a laboratory in Parella, Turin, during the '40s in the XIX century. He named it intumescence of limited length [1]. The high loss of energy and the turbulence generated in a hydraulic jump can provoke erosion and scouring [2] in an OCH and great damage to hydraulic structures, so to avoid this issue, the hydraulic jump must be contained in a stilling basin [3] whose length (L) must be calculated correctly to avoid any damage. Whoever, because of the phenomenon’s high energy loss and turbulence there is no mathematical equation that can explain or prove what exactly occurs during this intumescence or calculate length with precision, therefore length in a stilling basin is currently calculated with an empiric equation. Then, according to [4] state that: [5] - [13] have been developed various analytical and empirical models from their experimental study for hydraulic jump characteristics, some methods as [7] design of stilling basins or calculate de L.
Then, this work presents the results of the generated hydraulic jumps and the measure of its length with the methods: Smetana, Safranez, Einwatcher, Woycicki and Chertusov [14], Pavlovsky and Sieñchi methods.
Contemporary science explains Bidone’s intumescence formation with the critical regime theory: a hydraulic jump forms when water flows change from a supercritical to a subcritical regime [15]. The critical regime theory comes from the Bernoulli equation for a constant and incomprehensible flow [16] :
(1)
where C is a constant, P is pressure, ρ is fluid density, g is earth gravity, z is a horizontal comparison plane and v velocity of the water flow.
When Equation (1) is modified and written this way, it is called specific energy [17] [18] :
(2)
The y variable of Equation (2) is the hydraulic depth, Q is the flow rate, g is earth’s gravity and A is the hydraulic area. Giorgio Bidone reported the intumescence phenomenon to the Science Academy of Turin in 1824, that is to say, 86 years after Daniel Bernoulli published his book about hydrodynamics where he mentions for the first time what we now know as kinetic energy. He named it “vis viva” in Latin (which he took from an original terminology from Leibnitz). During that time there was not a solid theory about total energy or specific energy of water in open channels hydraulics, it developed over time with contributions from Euler, D’Alembert, and Lagrange [1].
On the other hand, the critical regime flow theory only explains the conditions that produce a hydraulic jump, while the hydrostatic laws and momentum law allows us to estimate the conjugate depth, but none of them determines with full mathematical precision the length of a hydraulic jump. Then, many specialized publications or textbooks, for example [14], propose several empirical equations to calculate length in a hydraulic jump, but each one of them is only valid under certain conditions. Table 1 presents seven empirical equations to calculate L and compares them against the estimations made with H-Canales and the results from direct measurements on an experimental prototype.
2. Materials and Method
We carried out the experiments in a variable slope open channel hydraulic (VSOCH) which we designed in the Civil Engineering Faculty of the Autonomous University of Chiapas, Mexico. The materials that were used to build the VSOCH (Figure 1) were the following [19] [20] : tubular metal with a rectangular profile (PTR by its acronym in Spanish), 10 cm wide and 5 m long, to build the base; hydraulic jack to change the slopes of the VSOCH; acrylic sheets to build the walls and the bottom of the channel; two tripods constructed with three pieces of steels of 3 inches diameter in order to support the weight of the channel; volumetric tank that serves not only as “reference model” for validate Q but also as a deposit “water intake” for the circulation of the flow rate (Q) into the VSOCH; 4 HP gasoline pump with 2 inches diameter of output to drive water from the deposit “water intake” to VSOCH and to circulate the flow in a closed system.
Figure 1. Transportable hydraulic channel (VSOCH) of Engineering Faculty of Chiapas State University, Mexico.
Once we installed the VSOCH, we input four different flows rates (Q) through it. Each Q had an arbitrary increase and decrease of Q (Figure 2) in order to avoid statistical tendencies in the L calculation, while the gate opening “a” remained constant.
We estimated the flow rates (Q) with volumetric method, using the circulation tank of the VSOCH which was previously calibrated. The Manning “n” used is 0.009 (acrylic), and the slope used is S = 0.0035. Each flow rate generated a hydraulic jump and for each jump; we measured the following variables directly in the experimental model: Length (L), major hydraulic depth (y1) and minor hydraulic (y2), which provided the data to calculate Froude’s numbers:
and
.
Later, we calculated L length with the equations shown in Table 1, and finally we estimated L with the proposed method and H-Canales program which uses the Sieñchi equation. The H-Canales software for channel design is the most used in Mexico and Latin America; in its 3.1.1 version it has a menu with 8 tabs [21]. One of them is named hydraulic jump, which has five options for calculating open channels hydraulics with different geometry: rectangular section, trapezoidal section, rectangular section, parabolic section, and circular section.
3. Results
The hydraulic jump that formed from each Q was directly measured in the VSOCH (Figure 3).
The results of L are shown in Table 2.
The best method in order to calculate L is the proposed method, because it has the lowest MSE and bias and the upper ME (see Table 3). Table 3 shows the statistical comparisons from the eight methods.
As seen in the results, the proposed method, the Sieñchi method (H-Canales) and the Pavlosky method are the best practices for estimating hydraulic jump
Table 3. Statistical comparation between eight methods.
Figure 2. Waterflow (arbitraryincrease and decrease).
Figure 3. Hidraulic jump under a spillgate.
length under a constant gate opening a = 2 cm, an arbitrary increase and decrease in volume and a slope of S = 0.0035 (see Figure 4).
Figure 4. Comparation lenght of hidraulic jump.
The proposed equation in order to calculate L in a rectangular open channel hydraulic shown in Figure 1 is:
(3)
where L is hydraulic jump length in cm, y2 and y1 are major and minor hydraulic depth in cm. Whoever, all the empirical equations are formulated under certain limits of operation. The limits of operation of the proposed equation are:
.
4. Conclusions
We executed a series of hydraulic experiments in a variable slope open channel hydraulic (VSOCH)from the Engineering Faculty in the Autonomous University of Chiapas, Mexico, with the objective to propose an equation to determine the hydraulic jump length of different water flows with an arbitrary increase and decrease, a constant gate opening “a” and a slope of S = 0.0035.
We then compared the proposed equation with seven equations presented in specialized literature, including the Sieñchi equation, which is used in the software H-Canales. The results show that the equation proposed in this investigation is exact and precise, and compared with the results of measurements in the experimental model it has an MSE = 0.1337, bias = −0.0049, ME = 0.9991 and R2 = 0.9993. Meanwhile, the calculation of L with the Sieñchi equation, compared with the experimental model, shows and MSE = 0.1741, bias = −0.0437, ME = 0.9984 and R2 = 0.9997. The results demonstrate that both equations are highly recommended to estimate L in rectangular open channel hydraulic as shown in Figure 1. The Pavlosky equation is also quite precise, under the experimental conditions presented in this investigation.
Acknowledgements
Thank Architect Laura Montserrat Mundo-Navarro for reviewing this paper.