ABSTRACT
In this research, the effect of sudden non-continuous contraction on the energy dissipation of supercritical flows was numerically investigated. This study focuses on energy dissipation in sudden contractions in supercritical flows. The numerical models were studied using FLOW-3D software and the random number generator (RNG) turbulence model. The laboratory tests were performed on sudden contractions of 10–15 cm and used shapes such as geometric trapezoidal and semicircular. Numerical simulations were carried out at a distance of 1.5 m from the supercritical flow generator gate, at a fixed opening of 2.6 cm, with the Froude number in the range of 2.5–7 and the relative contraction in the range of 8.9–9.73. This laboratory study found that as the Froude number of the upstream flow increases, the energy dissipation increases in the upstream (ΔE/E0) and downstream (ΔE/E1). For the 15 cm contraction, the results indicated that the energy loss compared with section A is 48.25% and compared with section B, it is 69.5% more than a free hydraulic jump in this channel. According to the conclusion, this value in trapezoidal constriction is 45.73 and 63.6% higher than in free hydraulic jump, respectively.
HIGHLIGHTS
The Froude number in the range of 2.5–7 and the relative contraction in the range of 8.9–9.73, trapezoidal and semicircular shapes were tested.
It was observed that as the Froude number increased, energy consumption increased both upstream and downstream.
According to the hydraulic jump, the energy loss in section A is 48.25% and in section B it is 69.5%.
The energy loss is 45.73% higher upstream and 63.6% higher downstream than in the free hydraulic jump.
INTRODUCTION
The design of channels downstream of hydraulic structures is very important to control hydraulic jumps (Ead & Rajaratnam 2002; Harrison 2015; Abbaszadeh et al. 2023; Norouzi et al. 2023). Chutes are one of the useful structures to reduce the energy downstream (Chinnarasri & Wongwises 2006). For the first time, the discharge coefficient in these structures was investigated by Hager & Dupraz (1985). Discharge equations for contraction flumes were developed by Samani et al. (2021). Results indicated that this function is useful for the design of hydraulic structures. Shock waves are created by the vertical side walls of the chute contract and cause the supercritical flow to deflect. In recent years, many studies have been completed in this field (Molinas 2000; Kaushik et al. 2012; Das et al. 2014). Hydraulic jumps were investigated in different roughened beds by Hughes & Flack (1984). Trapezoidal wavy beds were investigated by Izadjoo et al. (2005). Channels that are often too narrow to reduce the destructive energy of the flow are used and the energy of supercritical flow is higher below the gate. Therefore, solutions must be found to reduce the high flow rate, because high speed causes erosion, foaming, and damage to structures. Parsaie et al. (2018) numerically investigated the flow pattern in spillway approach channels. The results show that the current geometry of the left wall causes instability in the flow pattern and makes the secondary and vortex flows at the beginning approach channel. Daneshfaraz et al. (2020) investigated numerically and experimentally the effect of sudden narrowing on hydraulic parameters. Yamini et al. (2022) investigated the hydraulic performance of the seawater inlet system using computational fluid dynamics (CFD) modeling. The results show that the designer should make every effort to avoid small entrance and filtration areas from the basin to the intake forebay bottom, which could result in jet outlet and/or supercritical flow. Kyriakopoulos et al. (2022) investigated the hydraulic performance of Howell-Bunger and butterfly gates used for downstream outlets in large dams under flood risk. Results showed that based on the numerical study, the parameters of flow patterns, velocity profiles and vectors, turbulence kinetic energy, and formation of flow vorticity were investigated as criteria to determine the appropriate longitudinal distance.
Afaridegan et al. (2023) numerically investigated the hydraulic parameters of modified semi-cylindrical weirs. The analysis of streamline patterns from upstream to downstream of the weirs showed that the curve of the crest provided an opportunity for the flow to harmonize with the surface of the crest.
Al-bedyry et al. (2023) investigated the moving hydraulic jump in a trapezoidal channel in a laboratory and showed the connections between discharge and shifting hydraulic jump variables. Parsaie et al. (2023) investigated the flow rate coefficient in Samani, Magallanez, Baiamonte, Ferro (SMBF) flumes under free and submerged flow conditions. The results showed that the submerged condition does not reduce the Cd, but it reduces the discharge capacity.
The converging structures that are used downstream of the channels cause the formation of hydraulic jumps. Primary and secondary jump depth changes are among the most important topics in hydraulic engineering. Energy dissipation and the parameters affecting it have been investigated in various types of research.
In recent years, with the development of three-dimensional models and the advancement of computers, it has become possible to extract relationships with the desired accuracy from these models, but three-dimensional models require considerable time to be implemented. Three-dimensional models of flow simulation, including FLOW-3D software, have been widely used in the simulation of flow in open-channel structures. Comparison with experimental results shows that the models in this software can predict flow characteristics with high accuracy. In this model, there are five parental mixing length turbulence models, a kinetic energy equation of turbulence, a two-equation k-ε model, a random number generator (RNG) normalized group model, and a large eddy simulation (LES) model. Therefore, in this research, the hydraulic parameters were investigated using FLOW-3D software. Theoretical equations for the sequential-area and sequential-depth ratios for hydraulic jumps in contraction were developed considering the effects of a contracting width and a sloping bottom. Previous studies show that research conducted in the field of contracting channels has provided effective and useful results. However, further research is needed on the influence of the geometry of the flow path narrowing elements and their arrangement in the channel. Here, the influence of the geometry of the limiting elements of the transverse path of the flow with two semicircular and trapezoidal geometries has been investigated.
METHODS
Numerical modeling
In the above relationships, ΔE is the energy loss between sections A and B in meters. EA is the flow energy in section A in meters, n: is the number of datapoints, and Equations (5)–(7) are used to calculate absolute mistake, relative mistake, and root mean square mistake, respectively. The results are provided in Table 1.
Row . | Total cell . | ΔE (exp) . | ΔE (num) . | RE . | AE . | RMSE . |
---|---|---|---|---|---|---|
1 | 454,253 | 0.54 | 0.586 | 14.176 | 0.073 | 0.0788 |
2 | 577,845 | 0.54 | 0.531 | 5.212 | 0.034 | 0.0458 |
3 | 858,450 | 0.54 | 0.505 | 0.734 | 0.008 | 0.0119 |
4 | 920,955 | 0.54 | 0.502 | 0.682 | 0.008 | 0.0108 |
Row . | Total cell . | ΔE (exp) . | ΔE (num) . | RE . | AE . | RMSE . |
---|---|---|---|---|---|---|
1 | 454,253 | 0.54 | 0.586 | 14.176 | 0.073 | 0.0788 |
2 | 577,845 | 0.54 | 0.531 | 5.212 | 0.034 | 0.0458 |
3 | 858,450 | 0.54 | 0.505 | 0.734 | 0.008 | 0.0119 |
4 | 920,955 | 0.54 | 0.502 | 0.682 | 0.008 | 0.0108 |
The results showed that the reduction in the size of the meshes has reduced error. Because the error rates of rows 3 and 4 are almost the same, so row 3 was used to increase the speed of the numerical models. Mesh sizes are between 0.05 and 0.032 m and usually the number of meshes is taken as 3. In this study, two mesh blocks were used. Figure 1(b) shows the boundary conditions and networks in this study. In mesh block 1, the pressure boundary condition for Xmin, outflow for Xmax, wall for Ymin, Ymax, and Zmin and symmetry for Zmax was used, and in mesh block 2, the symmetry boundary condition for Xmin and Xmax, wall for Ymin, Ymax, and Zmin and symmetry for Zmax was used.
1. LES model,
2. k-ω model,
3. k-ε model, and
4. RNG normalized group model.
Validation
RESULTS AND DISCUSSION
Longitudinal profile and velocity
Energy dissipation
For a contraction of 15 cm, results show that the energy loss is 48.25% higher than the energy in section A and 69.5% higher than the free hydraulic jump compared with section B. Figures 5 and 6 have been prepared to show the changes in fluid energy compared with the energy of the fluid at sections A and B (Figure 5). One of the most important results obtained from the current research is that the narrowness elements have caused an increase in energy dissipation. With the result obtained, it can be suggested that instead of creating stilling basins with large dimensions and high construction costs, the dimensions of the stilling basins should be reduced or removed altogether. This will save a lot of project costs.
With a general and more detailed examination of the obtained results, it was observed that with the narrowing of the channel width, the energy dissipation also increased. This result can somehow indicate that the change in the geometry of the main channel can also play an effective role in the management and design of hydraulic structures.
Downstream Froude number
The results show that as the flow velocity at the inlet increases, the downstream Froude number increases strongly with the constriction elements. The reason for this is that when the streamlines collide with the narrow section, the hydraulic jump causes an increase in the energy loss, and then, the Froude number in the contraction section decreases. This is even though the Froude number drop compared with the simple case is 68 and 66% for semicircular and trapezoidal constrictions, respectively (Figure 7).
Downstream depth
The downstream depth parameter is one of the most important parameters that must be considered when using additional elements. The use of supplementary elements has the aim of increasing the depth downstream and reducing the energy of that section.
CONCLUSION
In this research, the effect of sudden non-continuous contraction on the energy dissipation of supercritical flow was numerically investigated. The numerical simulation model was performed with FLOW-3D for two contractions of 10 and 15 cm. Semicircular and trapezoidal geometries were used to create narrowing. The most important results from this research are given below:
1. Investigating the disturbance models indicated that the RNG disturbance model provides superior results compared with other models (Daneshfaraz et al. 2018).
2. When the models are compared with each other, the relative energy dissipation increases with the narrowing of the channel width.
3. Some of the energy loss is related to flow turbulence and swirling currents between the gaps in the narrowing elements. The results showed that the energy loss in the non-continuous constriction of 10 cm is 39.39% higher than the energy at location A and 62.5% higher than the free hydraulic jump compared with section B.
4. For a 15 cm narrowing, the results showed that the energy loss compared with section A is 48.25%, and compared with section B, it is 69.5% more than a free hydraulic jump. This value in trapezoidal constriction is 45.73 and 63.6% higher than the free hydraulic jump, respectively.
5. In general, the amount of energy loss in a 15 cm contraction exceeded that of a 10 cm contraction.
6. This result can somehow indicate that the change in the geometry of the main channel can also play an effective role in the management and design of hydraulic structures.
DATA AVAILABILITY STATEMENT
All relevant data are included in the paper or its Supplementary Information.
CONFLICT OF INTEREST
The authors declare there is no conflict.