Experimental study on flow over arced-plan porous weirs

Unlike conventional impermeable weirs, porous weirs without clogging the flow and passage of aquatic life with increased aeration and aerobic reactions with minimal negative effects on the environment are known as environmentally friendly structures. This study experimentally investigates the hydraulic performance of Arced-Plan Porous Weirs (APPWs) in different hydraulic and geometric conditions. For this purpose, four different porous and two solid weirs were examined. Experiments were conducted in a horizontal laboratory flume with length, width, and height of 20, 0.6, and 0.5 m, respectively, for a wide range of flow rates, particle sizes, and three arc lengths. Results showed that increasing filling material sizes increases the free discharge coefficient and reduces the submerged Discharge Reduction Factor (DRF). It was also concluded that the weirs’ effective length significantly impacts the free discharge coefficient and has no significant effect on the threshold submergence index and submerged DRF. Unlike solid weirs, the threshold submergence of porous weirs occurs at a downstream depth lower than the weir’s height. Finally, according to the dimensional analysis and Gene-Expression Programming (GEP) approach, three relations were extracted to calculate the free discharge coefficient, threshold submergence index, and submerged discharge reduction factor for APPWs.


INTRODUCTION
Hydraulic structures made of metal mesh and stone have been used in various fields of water engineering, such as stepped spillways for energy dissipation, groins for bed and riverbank erosion control, and detention dams for flood control. In the last decade, engineers have tried not to change nature and landscapes as much as possible alongside construction. Rigid concrete structures disrupt aquatic traffic, especially in rivers and natural canals, and negatively affect the environment and natural landscapes.
Weirs are usually made of impermeable materials such as concrete. These impermeable bodies prevent the longitudinal movement of aquatic organisms, physical and chemical substances in the water, and as a result, cause the accumulation of sediment behind the weir and negative effects on the water ecosystem (Mohamed 2010). Researchers have recently tried to use inexpensive and available natural materials to eliminate these destructive effects and not change the natural landscape. Among these, porous weirs with characteristics such as stability (no pressure and low boiling force), economic efficiency (use of available and natural materials), permeability, and flexibility against forces, are good alternatives (Fathi-moghadam et al. 2018;Rahmanshahi & Shafai Bejestan 2020). According to Chanson (2006), these structures are among the recommended and valuable structures for future studies. These types of structures with more aeration cause the activity of aerobic microorganisms and thus increase the self-purification capacity of the river. Increasing the self-purification capacity will purify water and eliminate bad odors in waterways to which sewage is discharged (Mohamed 2010;Fathi-moghadam et al. 2018).
So far, several studies have been conducted on the flow inside the permeable bodies. However, in most of these studies, the upstream depth is less than the height of its crest, and the interaction between the over and through flows and their effects on each other has been less considered. Kells (1993) studied the gradual variable flow over the gravel dam for two different flow conditions, including partial and complete fall from the dam. Kells (1994) showed that the ratio of through flow to overflow in the gravel dam ranges from 0.25 to 0.5. Michioku et al. (2005) studied the hydrodynamics of a rubble mound theoretically and experimentally. Performing a one-dimensional analysis on a constant uniform flow through the rubble mound, they described that the discharge is a function of parameters such as upstream and downstream flow depth, porosity, gravel particle diameter, and rubble mound length. Michioku et al. (2007) conducted laboratory studies on flow over weirs and gravel groins. Their findings showed the mass and momentum changes between the mainstream and the stream in the porous media. Leu et al. (2008) studied turbulent flow around a permeable body. They examined three models with porosities of 0, 34.9, and 47.5% and showed that the downstream rotation area is larger for porous structures due to the passage of flow through it. Mohamed (2010) investigates two models of rectangular gabion weirs and two solid weirs with the same geometry for comparison and provided two relationships for predicting the discharge coefficient in free and submerged flow conditions. Salmasi & Taghi Sattari (2017) experimentally studied the hydraulic performance of broad-crested gabion weirs in free and submerged flow conditions. Finally, using the M5 model, they extracted relationships for the discharge coefficient in both conditions. Mohammadpour et al. (2013) numerically modeled the flow on the porous broad-crested weirs. Their results showed that the k-ε model has better accuracy for modeling this type of weirs. Fathi-moghadam et al. (2018), using the Flow3D numerical model, investigated the flow pattern and the hydraulic performance of porous trapezoidal and triangular weirs in free-flow conditions. The study also discussed energy dissipation and over and through flow rates. Safarzadeh & Mohajeri (2018) numerically simulated broad-crested weirs in free-flow conditions and provided relationships to calculate the discharge coefficient. The ratio of thought flow to total flow in different hydraulic and geometric conditions was obtained in this research. Also, new relationships were extracted to calculate the discharge coefficient using a nonlinear multivariate regression model. Salmasi & Taghi Sattari (2017) and Salmasi et al. (2021) experimentally investigated porous broad-crested weirs' hydraulic performance under free and submerged flow conditions with four different porosities. Their results indicated that the discharge coefficient in submerged flow is approximately 20% lower than the free flow's discharge coefficient. Shariq et al. (2020) investigated the porous broad-crested weir discharge equation under through flow conditions. Their results reveal that the discharge through the porous broad-crested weir is more sensitive to the material size and upstream depth than the downstream depth of the porous broad-crested weir.
According to the common weir's head-discharge relationship (Q ¼ C d W ffiffiffiffiffi 2g p H 1:5 ) (Fritz & Hager 1998), to increase the hydraulic performance of the weirs, the discharge coefficient or the effective length of the weir can be increased (Crookston & Tullis 2012). So far, various works have been done in this field in impermeable weirs. Using an arched-plan and a W-shape or zigzag increases the weirs' effective length and increases its performance. The main idea of this research was that a porous weir with an arced-plan increases flow rate first because it is porous (because it has both through and overflow), and secondly, the arc plan increases weir's effective length. According to the research background, previous studies regarding permeable weirs mainly were involved with linear-plan broad-crested weirs. Thus, no comprehensive research has been conducted on Arced-Plan Porous Weir's (APPWs) hydraulic performance in free and submerged flow conditions. In addition, none of the previous studies have examined the threshold of submergence or modular limit of porous weirs. Thus, there is a lack of comprehensive research on the hydraulics of such weirs.
The main purpose of this study is to investigate the hydraulic performance of the APPW experimentally. For this purpose, experiments were conducted on different physical models of APPWs to determine the effect of various factors such as flow rate, material size, and geometrical dimensions of weirs on their free discharge coefficient, threshold submergence, and submerged Discharge Reduction Factor (DRF). Also, the performance of two porous weirs is compared with their corresponding solid weirs. In addition, two weirs with the same geometry but with different directions were tested to investigate the effect of the arced-plan weir direction on its hydraulic performance in free and submerged conditions. Finally, the Gene-Expression Programming (GEP) approach with multivariable nonlinear regression is used to develop three mathematical equations for free discharge coefficient, threshold submergence index, and submerged DRF for engineers to use in APPW design.

Experimental setup and measurements
The experiments of this research were performed in a horizontal, glass-walled laboratory flume with length, width, and height of 20, 0.6, and 0.5 m, respectively, in Khuzestan Water and Power Authority (KWPA) laboratory, Iran. The flume has a closed water flow system, and an end gate was used to control the downstream depth in submerged flow conditions. A rectangular weir at the downstream end of the system was used to measure the discharge, and a point gauge with an accuracy of +0.2 mm was used to measure water surface elevation. The weirs were made of a metal rod with a diameter of 4 mm and covered with 5 Â 5 mm fine wire mesh, smaller than the smallest desire filling stones. The model frames were filled with different gravel sizes with a coefficient of uniformity of the materials near 1. The filling materials used in porous weirs include five sizes with almost uniform diameters and sharp corners in 7.13-31.75 mm. In this research, four porous and two solid weirs were constructed. Two solid weirs were corresponded to PW-M1 and PW-M3 and were made of plexiglas. They have been constructed to compare their results with corresponding porous weirs. The physical model test matrix also included one APPW weir with inverted direction. This experiment was performed to investigate the directional effect of APPW.
It is not always practical to set the weir in only free-flow conditions or to avoid submergence. On the other hand, the control flow structures' performance depends mainly on the downstream water level. Hence, all experiments have been carried out for free and submerged flow conditions to understand better the flow status's effect on the hydraulic behaviour of flow over and through the porous weirs. In addition, the threshold submergence is determined as the water level boundary in which the free approach flow regime turns into submergence. This is defined as the water level where the upstream water depth starts to increase by þ1 mm due to the increase of downstream water depth (Hager & Schwalt 1994;Fritz & Hager 1998).
In free condition experiments, after stabilizing the flow discharge and the upstream head over the weir, the upstream water depth was measured at a distance about four times the head over the weir. In order to determine the threshold submergence boundary, the downstream water depth was slowly increased using the end gate until the upstream flow depth was increased by 1 mm. Then, using up and downstream water depths, the threshold submergence index was calculated. This procedure continued for various submergence ratios to calculate submerged DRF. All measurements were made along the centreline of the flume. Figure 1 shows definition sketches and experimental views of porous and solid weirs used in this study. Geometric details and acronyms for the constructed models are given in Table 1. In Table 1, PW refers to a Porous Weir, IPW refers to the Inverted Porous Weir, and SW refers to a Solid Weir. Also, L c is the length of the weir crest, P is weir height, Q is flow discharge, W is channel width, d 50 is gravel median size, n is gravel porosity, and L arc is weir's effective length. In this research, all models have been tested in different discharges, and a total of 175, 175, and 1,708 experiments have been performed in free, threshold submergence, and submerged flow conditions, respectively.

Dimensional analysis of APPW
Free discharge coefficient, C d and flow rate could be related as in Equation (1) (Fritz & Hager 1998).
where Q is water discharge, W is the width of the weir, g is the gravitational acceleration, H is the total upstream head , h is the head over the weir, and v is the approach velocity.
For an APPW, the discharge coefficient can be expressed by the following functional relationship: in which f is a functional symbol, s is surface tension, ρ is water density, μ is water viscosity, h is upstream water level, and n is gravel porosity.

Uncorrected Proof
Using the Π theorem of dimensional analysis (Barenblatt 1987) for the parameters in Equation (2), C d can be expressed by following dimensionless groups; where Re ¼ Reynolds number; and We ¼ Weber number.

Uncorrected Proof
At a specified boundary of the downstream depth,h TS , as the depth increases, the weir's upstream depth begins to increase. This depth's boundary is known as the threshold submergence. For porous weirs, the threshold submergence index can be explained as Y TS ¼ (h TS À P)=P. Similar to free flow conditions, by applying the dimensional analysis, the threshold submergence index can be formulated as Generally, for a specified weir's upstream head, the submerged discharge, Q s , is smaller than that in the free flow condition, Q. Hence, the discharge for the submerged flow is defined by the following equation: In addition to the governing parameters in free-flow conditions, the downstream depth, h d , affects the DRF. In this study, DRF can be described by the following functional relationship is defined as the submergence index, that for free flow Γ ,0 (which indicates that the weir is not submerged), for modular limit Γ ¼ 0, and for fully submerged flow condition Γ ¼ 1.
In all experiments, the head over the weirs was larger than 3 cm to eliminate the effect of surface tension and the Weber number from Equations (3), (4), and (6) (Horton 1907;Novak & Cabelka 1981). According to Table 2, the Reynolds number in all experiments was high enough to provide a fully turbulent flow and is not taken into account. The ranges of effective independent variables in Equations (3), (4), and (6) based on experimental data are given in Table 2. Consequently, three functional relationships are simplified as follows

RESULTS AND DISCUSSIONS
The flow over weirs can be modular or non-modular. Figure 2 shows the modular or free-flow condition on PW-M3. The water surface profile downstream of the arced-plan weirs in free-flow conditions has a convexity in the channel's centreline that increases with increasing flow rate. Figure 2(b) illustrates the non-modular or submerged flow condition over SW-M2. Like other type of weirs, the weir's downstream water surface in the submerged conditions has surface waves that decrease with increasing submergence ratio. Figure 3 depicts the discharge coefficient for PW-M1 and PW-M2 at the different flow rates in free and submerged flow conditions. Clearly, with the increase of weirs submergence, the discharge coefficient decreases, and according to Equation (5), the submerged discharge coefficient is a reduction coefficient in a free condition. The results also show that in low flow rates, the discharge coefficient decreases at a higher rate. According to these figures, the submerged discharge coefficient sometimes is 68% less than the corresponding free-flow discharge coefficient. The free conditions discharge coefficients are calculated using Equation (1). Figure 4 shows the discharge coefficient versus H/P for PW-M1 and PW-M3 for all values of d 50 =P (the ratio of the mean size of filling materials to the weir's height). The discharge coefficient of the two corresponding solid weirs (SW-M1 and SW-M2) is also presented for comparison with the porous weir. Results show that with increasing material sizes, the discharge coefficient increases. Also, experimental observation indicated that increasing particle sizes reduces the upstream water level due to increased pores size. According to the Darcy-Weisbach equation, as the size of the weir's filling material increases, the pore diameter increases, and the through-flow resistance decreases. Due to the lower flow resistance in porous weirs, the through-flow fraction and the discharge coefficient increase. Unlike porous weirs, the discharge coefficient increases with increasing H/P (or discharge) ratio in solid weirs and becomes almost constant for H/P .0.8. Results show that the difference between C d of different particle sizes is reduced with increasing H/P. This is a direct consequence of the reduced ratio of through-flow to total flow discharge. For H/P values greater than 0.9, the discharge coefficient values of solid and porous weirs are very close to each other. Hence, in low  Uncorrected Proof values of H/P ratios, smaller than 0.5, the discharge coefficient is significantly high due to the higher proportion of through flow to the total flow. The discharge coefficients for different APPWs are presented in Figure 5. It can be inferred that in weirs with an arced-plan, increasing the weir's effective length increases the overflow and increases the weir's upstream contact surface; the throughflow fraction also increases. Thus, the maximum discharge coefficient belongs to the APPW with maximum effective length (PW-M2).
Threshold submergence boundary Figure 6 shows the threshold submergence index Y TS , versus H/P for various ratios of d 50 =P. There is a relatively linear correlation of the threshold submergence index with H/P. In this study, the threshold submergence index is defined purposefully as a function of the weir's height (P). According to this definition, if the threshold submergence occurs in a downstream depth lower/higher than the weir's crest level, then the threshold submergence index would be negative/positive, respectively. For example, if Y TS is equal to À0.2 (h TS , P) or 0.2 (h TS . P), the threshold submergence occurs in a downstream depth of 0:2 P lower or higher than the weir's crest level, respectively. For Y TS ¼ 0, the threshold submergence level is equal to the weir's height, P. Threshold submergence in the conventional impermeable weirs occurs at the stage in which the tailwater is higher than the weir's crest level (Hager & Schwalt 1994;Fritz & Hager 1998;Schmocker et al. 2011;Azimi et al. 2014). On the contrary, for most of the low discharge cases in this study, the threshold submergence occurred at the stage in which the downstream depth,h d , was lower than the crest level of the porous weir. Results indicated that material sizes had no significant effect on the threshold submergence index.
The variation of threshold submergence index Y TS versus H/P for different APPWs for two material sizes is illustrated in Figure 7. Results show that, for finer filling material, with the increases of weir's effective length in the arced-plan models, the threshold submergence index increases, and the APPW is less sensitive to downstream depth. Also, Figure 7(b) depicts that the threshold submergence index in all four weirs is not significantly different for the coarse filling material.
The DRF of submerged flow Figure 8 shows the changes in the DRF versus relative submergence Γ for PW-M2 and PW-M4 for maximum (31.75 mm) and minimum (7.13 mm) filling material sizes. The relative submergence was defined as the G ¼ ((h d =h u ) À (h TS =h u ))=(1 À (h TS =h u )), where 0 , G , 1. Results show in both models, with increasing the flow rate, the DRF increases. In other words, with increasing flow rate, the sensitivity of the APPW to the downstream depth decreases. However, for a high flow rate, the DRF is not significantly different. Figure 9 compares the DRF values in each model for two different material sizes simultaneously with the corresponding solid weirs. Results show that with increasing material size, the DRF coefficient decreases, and the weir is more prone to submergence. Comparing the porous and its corresponding solid weir results indicated that the impermeable weirs are less sensitive to downstream depths than porous weirs and are submerged at higher downstream depths. Figure 10 shows the effect of different APPWs on DRF. For this purpose, the DRF values were compared in two different granulations for a constant flow rate. The results show that the weir's shape does not have a significant effect on the DRF.

Extraction of mathematical relationships
In this section, the aim is to extract the three relationships developed in the dimensional analysis section. To derive the free flow discharge coefficient, threshold submergence index, and submerged DRF, frequency relationships were first tested based on nonlinear multivariate regression using SPSS software. Considering that the user's selection of mathematical function type is performed at each step in a standard regression model, and the model determines only constant ratios, this method is very time-intensive. Further, it cannot be claimed that the best type of mathematical function is fitted due to the limitation of possible combinations of the selected mathematical functions. However, in soft computing models similar to GEP, the model intelligently selects both mathematical functions and constant ratios. Regarding the model's metaheuristic nature, the appropriate function is selected with a high speed. See Ferreira (2001aFerreira ( , 2001b and Rahmanshahi & Shafai Bejestan (2020) for more details about GEP. It is noteworthy that d 50 /P and L c //L Arc were removed from Equations (8) and (9), respectively. Because the dependence of the threshold submergence index on the particle size as well as the DRF on the APPW shape can be neglected. The optimal equations are developed as Equations (10)-(12). Also, in Figure 11, the observed and calculated   Table 3.

APPW direction effects
The following section compares the hydraulic performances of two types of APPWs with the same geometry but in the opposite direction. For this purpose, PW-M3 with the same grain size in two different directions has been tested according to   Figure 12(a) and 12(b). In Figure 12(c), the free discharge coefficient and Figure 12(d), the threshold submergence index of these two models are compared. The results show that there is no significant difference between these two weirs. Figure 12(e) and 12(f) depicts their DRF for two different flow rates. The results indicated that for a specific discharge, DRF is almost the

CONCLUSIONS
In this study, an experimental investigation was performed to determine the free discharge coefficient, threshold submergence, and submerged DRF of APPWs in different geometric and hydraulic conditions. In free-flow conditions, the water surface profile downstream of the APPWs has a convexity in the channel's centreline that increases with increasing flow rate. Like other weirs, the weir's downstream water surface in the submerged conditions has surface waves that decrease with increasing submergence ratio. The submerged discharge coefficient sometimes is 68% less than the corresponding free-flow discharge coefficient. Results show that with increasing material sizes, the discharge coefficient increases. Unlike solid weirs, the discharge coefficient decreases with increasing H/P (or discharge) ratio in porous weirs. The difference between C d of different particle sizes is reduced with increasing H/P and for H/P values greater than 0.9, the discharge coefficient values of solid and porous weirs are very close to each other. It can be inferred that in weirs with an arced-plan, increasing the weir's effective length increases the discharge coefficient. The threshold submergence of porous weirs occurred at the stage in which the downstream depth was lower than the weir's crest. Results indicated that material sizes had no significant effect on the threshold submergence index. For finer filling material, with the increases of weir's effective length, the threshold submergence index of APPWs increases and is not significantly different for the coarse filling material. The DRF increases with increasing the flow rate, and for a high flow rate, the DRF is not significantly different. The DRF coefficient decreases with increasing material sizes. The weir's shape does not have a significant effect on the DRF. Hence, it can be concluded that for similar material size and flow rate, the effects of APPWs direction do not affect weirs' hydraulic performance in both free and submerged flow conditions. Finally, three mathematical relationships were developed to estimate the free discharge coefficient, the threshold submergence index, and the submerged DRF.