Evaluating performance of various methods in predicting triangular sharp-crested side weir discharge

This paper investigated the efficiency of the traditional weir equation (TWE), Domínguez, adjusted Domínguez, and Schmidt approaches, as an alternative to the De Marchi procedure, for computing discharge of a sharp-crested triangular side weir. Comprehensive experimental data were used for the analysis, including 342 data from the present study and 140 data from other sources. The effects of approach Froude number Fr1, the ratio of weir height to upstream flow depth p/y1, and weir apex angle θ on the discharge coefficients obtained from different methods were studied. Sensitivity analysis using the partial swarm optimization-support vector regression method indicated that Fr1, p/y1, and θ affect the discharge coefficients. It was found that Fr1 with sensitivity indices equal to 1.89, 3.74, and 4.04 has the most substantial effect on the De Marchi coefficient, TWE coefficient, and adjusted Domínguez coefficient; meanwhile, p/y1 has the most significant impact on Domínguez coefficient and Schmidt coefficient with sensitivity index equal to 1.57. In addition, it was found that θ had the lowest sensitivity indices in estimating discharge coefficients. New equations for forecasting sharp-crested triangular side weir discharge coefficient were presented based on dimensional analysis. The new De Marchi coefficient executed better for calculating triangular side weir discharge than earlier De Marchi coefficients. Moreover, TWE, Domínguez, adjusted Domínguez, and Schmidt methods performed better than the De Marchi procedure (with MSE = 4.581) in calculating sharp-crested triangular side weir discharge. However, considering the simplicity of the TWE approach compared to other methods, this approach with R2 = 0.975, NSE = 0.975, MSE = 3.610, MRE = 0.097, and CP10% = 71.36 was introduced as the superior procedure.


List of symbols B
Main channel width C m De Marchi coefficient C d Discharge coefficient obtained from traditional weir equation

Introduction
A side weir is an overflow inserted into the main channel laterally to divert part of the flow from the main channel into a side channel. Side weirs are generally used in irrigation, land drainage, urban sewage systems, and sanitary engineering and are also widely used for storm relief and head regulators of distributaries. Like normal weirs, side weirs have diverse styles (labyrinth, sharp, and broad-crested) and forms such as rectangular, trapezoidal, and triangular. Depicted in Fig. 1 is the subcritical flow along the length of a triangular side weir. Where Q 1 and Q 2 are upstream and downstream discharges, y 1 and y 2 are the depth of water at upstream and downstream sections, p is weir crest height, B is main channel width, E is the specific energy ((E = y + Q/(2gB 2 y 2 )), g is the gravity acceleration, and θ is side weir apex angle. The hydraulic behavior of side weirs in channels has been studied since the turn of the twentieth century. However, many investigations have been based on empirical and experimental work. Other studies implemented theoretical approaches such as specific energy and momentum principle. Almost all experimental work and theoretical investigation are limited to prismatic rectangular channels with a horizontal overflow weir crest.
The flow along a side weir is a typical spatially varied flow with decreasing discharge. The energy equation generally emanates the governing equation for flow over side weirs. The general differential equation of spatially varied flow along a side weir ( Fig. 1) with decreasing discharge is expressed as (Henderson 1966): here α is the kinetic energy correction coefficient, S 0 is the main channel slope, S f is the friction slope, x is the distance along the side weir from the upstream end, dQ/dx is the discharge per unit length of the side weir, A is the crosssectional area of the flow, and T is top width of the channel. Equation (1) reveals that the longitudinal water surface profile along a side weir under subcritical flow is an ascending curve (Fig. 1). For a horizontal prismatic rectangular main channel, considering the kinetic energy correction coefficient α as unity and ignoring friction losses, Eq. (1) is rewritten as follows: De Marchi (1934), by considering constant specific energy E along the length of the side weir, solved the above equation for a rectangular side weir.
The discharge over a triangular side weir Q s is provided by (Kumar 1985): Differentiating Eq. (3) and considering x (= 2(y − p) tan(θ/2)) leads to (Kumar and Pathak 1987): where C m is the De Marchi coefficient. For a triangular side weir located in a rectangular channel, the following equation was introduced by Kumar and Pathak (1987) for C m : in which ϕ is the De Marchi function calculated using: where L is the effective length of the triangular side weir calculated by (Balahang & Ghodsian 2021): Side weir Side channel here, h 1 (= y 1 − p) and h 2 (= y 2 − p) are the depths of water on the weir crest at upstream and downstream sections, respectively. Presuming that E is constant along the side weir, Eqs. (5, 6, 7) are combined to form:

Main channel
By the above equation, y 2 is obtained by trial-and-error method, provided the values of C m , B, y 1 , p, E, and θ are known. Then, the downstream discharge Q 2 is calculated by using the energy equation. Then after, the discharge of triangular side weir Q s is calculated by: Various researchers have investigated the effect of different variables on C m of rectangular side weirs. Subramanya andAwasthy (1972), Nandesamoorthy andThomson (1972), Yu-Tek (1972), Prasad (1976), andRanga Raju et al. (1979) correlated C m to upstream Froude number for rectangular side weirs. While Singh et al. (1994) and Jalili and Borghei (1996) correlated C m with the upstream Froude number and the relative flow head for rectangular side weirs. Additional studies by Borghei et al. (1999) indicated that the ratio of side weir length to main channel width indirectly impacts the discharge of rectangular side weirs.
Unlike rectangular side weirs, only a few studies have been conducted on the discharge coefficient of triangular side weirs. Kumar and Pathak (1987) related the De Marchi coefficient to the upstream Froude number for (9) Q s = Q 1 − Q 2 triangular side weirs, while Ghodsian (2004) linked the De Marchi coefficient to the upstream Froude number and relative head for triangular side weirs. Various equations for the De Marchi coefficient obtained for triangular side weirs with various values of apex angles by Kumar and Pathak (1987) and Ghodsian (2004) are given in Table 1.
As mentioned earlier, De Marchi's approach is based on the trial-and-error method and hence, complicated to obtain triangular side weir discharge. So, it is necessary to examine the capability of other alternatives as introduced in the following: The traditional triangular weir equation (TWE) is expressed as follows (French, 1985): where C d is the discharge coefficient of a sharp-crested triangular weir. Unlike normal weirs, the flow over a side weir is affected by the velocity head (V 1 2 /2g) in the main channel. Therefore, TWE for triangular side weirs is defined as follows: here H 1 (= h 1 + V 1 2 /2g) is the total upstream head (i.e., in section 1 in Fig. 1). The variation of the flow surface profile along the side weir is not considered in TWE. Domínguez (1935), based on the following assumptions, introduced a simple method for estimating lateral discharge and water surface profile variations along a side weir: (1) Specific energy is constant along the side weir.
(2) The water surface varies linearly along the side weir (h(x) = h 1 + (h 2 -h 1 (x/L))). Based on Domínguez, the side weir discharge per unit length of a triangular side weir dQ s /dx is expressed as follows: here C dD is Domínguez coefficient. By integrating Eq. (23) with respect to x (x = 0 to x = L), the following equation is obtained for the discharge of a triangular side weir Q s : To consider the velocity head, Eq. (24) is rewritten by substituting h 1 2.5 with H 1 2.5 as follows: here C dD * is the adjusted Domínguez coefficient. Bagheri et al. (2014a) and Bagheri et al. (2014b) showed that TWE, Domínguez, and adjusted Domínguez approaches perform better than the De Marchi method for rectangular side weir discharge. They also reported that TWE and adjusted Domínguez coefficients strongly correlate with the approach Froude number for a rectangular side weir, while the Domínguez coefficient weakly correlates with this parameter. Schmidt (1954), by assuming a linear variation of water surface along a side weir, used (h 1 + h 2 + h 3 )/3 instead of h 1 to calculate side weir discharge. Based on Schmidt's approach, Eq. (26) is presented to calculate the discharge of triangular side weirs: here C sc is the Schmidt discharge coefficient, and h 3 is the average flow depths on the side weir crest at upstream and downstream sections. Emiroglu and Ikinciogullari (2016) reported the dependability of the Schmidt procedure for calculating the discharge coefficient of a rectangular side weir for Froude number in the range of 0.75 to 1. Balahang and Ghodsian (2021) proposed the following straightforward equation for calculating discharge Q s of a triangular side weir: where L is the effective length of the side weir, which is calculated using Eq. (7). It is clear from Eq. (27) that h 2 is essential for calculating Q s . Researchers have recently focused on using machine learning algorithms to solve engineering problems. Due to its high accuracy, the support vector machine (SVM) has been one of the most prevalent machine learning methods in solving problems of side weirs. Azamathulla et al. (2016) showed that the SVM method is more precise in estimating discharge coefficients than artificial neural network (ANN) and adaptive neuro-fuzzy inference system (ANFIS) techniques. Roushangar et al. (2016) reported that the combined support vector machine with the genetic algorithm (GA-SVR) method has better performance than the geneexpression programming (GEP) method in calculating the discharge coefficient of rectangular and trapezoidal side weirs. Zaji and Bonakdari (2017) illustrated that the SVR technique produces more precise results than the nonlinear regression (NLR) method for estimating a rectangular side weir discharge coefficient. Li et al. (2021) stated that the SVM algorithm produces minor errors in forecasting the rectangular side weir discharge coefficient compared to ANN and extreme learning machine (ELM) methods. Balahang and Ghodsian (2021) reported that the SVR method calculates the discharge of triangular side weirs better than ANFIS, ANN, and gradient-boosted regression trees (GBRT) techniques.
The review of prior studies shows that most investigations have concentrated on forecasting C m for sharp-crested rectangular side weirs. At the same time, less concentration has been paid to investigate the flow through triangular side weirs. The capabilities of the TWE, Domínguez, adjusted Domínguez, and Schmidt procedures, as an alternative to the De Marchi method, are analyzed to compute the discharge of sharp-crested triangular side weirs, by using a more comprehensive range of values of influencing parameters, compared to earlier researches. Thus, the goals and novelties of the present research are: (1) For the first time, the capabilities of the TWE, Domínguez, adjusted Domínguez, Schmidt, and De Marchi methods are evaluated for calculating triangular sharp-crested side weir discharge.
(2) The sensitivity analysis for influencing parameters on the discharge coefficients of triangular sharp-crested side weir (C m , C d , C dD , C dD * , and C sc ) using a hybrid machine learning approach (PSO-SVR) was implemented for the first time.
(3) By using a more comprehensive range of data sets, more accurate equations are introduced for calculating the discharge of triangular sharp-crested side weirs.

Experiments
The experiments were conducted in a prismatic horizontal channel by the second author. The length of the main channel was 9.0 m, its width was 0.5, and its depth was 0.5 m. At the end of the channel, a sluice gate was installed to regulate the water depth. The side channel was perpendicular to the main ones. The sharp-crest side weirs are made of the mild steel plate and installed at the upstream end of the side channel. A supply pipe provided the main channel discharge from an overhead tank with a constant head. A calibrated sharpcrested weir measured the discharges in the main and side channels. Point gauge with ± 0.1 mm accuracy measured the flow depths y 1 and y 2 ( Fig. 1) at the center line of the main channel. Figure 2 shows the laboratory setup used in the present study. Experiments were carried out for various discharges, flow depths, weir heights, and apex angles. All the experiments were performed under subcritical flow conditions. In addition, data obtained by Kumar (1985) and Mohan (1987) are also used for the analysis. Table 2 summarizes the data utilized.

Dimensional analysis
The discharge coefficient of a sharp-crested triangular side weir is a function of the subsequent geometric and hydraulic variables: , V 1 is the mean flow velocity in the main channel at the upstream section, μ is the dynamic viscosity of water, σ is the surface tension coefficient, ρ is the water density, and S 0 is the main channel slope. Using the Buckingham Π-theorem Eq. (28) is written as follows: where Fr 1 (= V 1 /((gy 1 ) 0.5 )), Re (= V 1 y 1 /ν), and W e (= py 1 V 1 2 /σ) are the Froude number, Reynolds number, and Weber number, respectively, and υ is the kinematic viscosity of water. When the flow is turbulent, the viscosity effect can be neglected compared to the inertial force. For (y 1 − p) > 30 mm, surface tension influence on the flow over a weir is insignificant (Novák & Čabelka 1981). El-Khashab and Smith (1976) and Borghei et al. (1999) stated that the effect of the main channel slope S 0 on the discharge coefficient is nominal. Thus, by ignoring the insignificant variables, Eq. (29) is written as follows:

Statistical indices
The following statistical indices were utilized to compare the performance of different equations in forecasting the discharge coefficient and the discharge of sharp-crested triangular side weir:   (1) Pearson correlation coefficient (R 2 ) measures the linear correlation between two random variables and is computed by: where o i and e i are the ith observed and estimated values of discharge coefficient or side weir discharge, o i and e i are the average of observed and estimated values of the discharge coefficient or side weir discharge, respectively, and N stands for the number of data used.
(2) Mean squared error (MSE) for determining the error value and difference between the measured and estimated values expressed as follows: (3) Mean relative error (MRE) is calculated using the following equation: (4) Ratio of data in the scope of less than ± 10% error (CP 10 %) is obtained from the following equation: here N 10 expresses the number of data with error

Particle swarm optimization: support vector regression
A support vector regression (SVR) is a kind of support vector machine for solving regression issues. The main purpose of the SVR model is to discover a function that provides a connection between dependent variable f(x) and independent variables {[x 1 , …, x n ], which is expressed as the following equation (Raschka 2015): where x i is n input vector, w i is the weight vector, b is bias, l is the number of samples, and K is the kernel function, which maps x i to linear space if the relationship between f(x) and x i be nonlinear. In the present study, radial basis function is used as kernel, which is expressed as follows: here γ determines the radius of the impact of support vectors. Two important SVR parameters besides γ that must be optimized for data training are: • C: This parameter adjusts the ratio between the complexity of the model and the required accuracy of the training data and is always greater than zero. • ε: This parameter determines the allowable error of the model, which can be a decimal number.
The PSO method is used in this study to get optimal SVR parameters. The details behind of hybrid PSO-SVR algorithm, as shown in Fig. 3, are as follows (Hu et al. 2015): Input: train dataset, number of particles n, inertia w, the cognitive element that models the direction of particles to replace to formerly discovered most satisfactory position φ 1 , a social element that quantifies execution of a particle close to optimal global particle φ 2 , and the highest iteration number T.
Output: an optimal hybrid of C, ε and γ. Initialization: Let t = 0, followed by: (1) Initialize (4) Update the position and velocity of the n particles, the current historical optimal position of the particles. For this objective, randomly generate ⇀ U 1 and ⇀ U 2 (n × n diagonal matrices in which entries on the major diagonal are arbitrary numbers uniformly dispersed in the gap [0,1]), then, i . Finally, when one particle passes out of the searching space, it can be controlled by stopping updating the particles fitness value (the MSE of the SVR's train consequence utilizing this particle) and going to the next step; 5) if MSE ( The present study implements the modeling process using Sklearn, Pandas, and Numpy library in Python 3.8.1. The radial basis function is considered the SVR kernel, and ε is assumed as 0.01 in the modeling process. PSO algorithm is used by setting w = 1.2 and n = 3 to obtain the optimal values of C and γ.

Sensitivity analysis
The sensitivity analysis was used to study the effect of independent parameters Fr 1 , p/y 1 , and θ on C m , C d , C dD , C dD * , and C sc . The sensitivity analysis was accomplished using the PSO-SVR technique. Eighty percent of the data is employed to train the models, and the rest is utilized to test the models. Based on Eq. (30), four models are represented for sensitivity analysis by skipping each variable simultaneously, as indicated in Table 3. The sensitivity index is computed for each model to execute the sensitivity analysis. The sensitivity index is the proportion of the model MSE with a skipped parameter to the model MSE in the existence of all the variables. A sensitivity index > 1 means the parameter's significance in the model. The results of the sensitivity analysis are given in Table 3. The number of iterations, values of C and γ, and statistical indices, including R 2 , MSE, MRE, and CP 10 % for each model, are also presented in Table 3.
The sensitivity indices in Table 3 indicate that Fr 1 is the most influencing variable in determining C m , C d , and C dD * . By ignoring Fr 1, the MSE index due to C m , C d , and C dD * increased by 89.28%, 273.91%, and 304.16%, respectively.
Unlike C m , C d , and C dD * , C dD and C sc are more influenced by p/y 1 than Fr 1 . In calculating C dD and C sc , the sensitivity index of p/y 1 is 1.29 times the sensitivity index of Fr 1 . By omitting the impact of p/y 1 on C dD and C sc , the MSE index increased by 56.76%. According to Table 3, in all the approaches, the side weir angle θ has the least effect on determining the discharge coefficient. By ignoring the effect of θ on C m , C d , C dD , C dD * , and C sc , the MSE index increases by 21.43%, 13.04%, 13.51%, 4.69%, and 13.51%, respectively.

Analysis of results
This part aims to study the effect of dimensionless independent variables of Eq. (30) on discharge coefficients obtained from the De Marchi, TWE, Domínguez, adjusted Domínguez, and Schmidt approaches using experimental data introduced in Table 2. The variations of C m , C d , C dD , C dD * , and C sc with approach Froude number are shown in Fig. 4.
The results displayed in Fig. 4a-e confirm that the triangular side weir discharge coefficient decrease by increasing Fr 1 . Increasing the value of Fr 1 is due to increased longitudinal velocity or decreased flow depth. The flow divagation angle and the outflow velocity of a side weir decrease  (Hager 1987). As a result, the discharge coefficient decreases. The linear water surface profile in terms of the approach Froude number is shown in Fig. 5. Due to the direct relationship between Fr 1 and Δh/L, it can be stated that the approach Froude number indirectly affects C dD and C sc . The discharge coefficients C dD and C sc have a similar correlation with Fr 1 due to considering the impact of the water surface profile and disregarding the influence of V 1 in Domínguez and Schmidt's approaches. The scattering of data in Fig. 4 proves that in addition to Fr 1, other dimensionless parameters such as p/y 1 and θ may affect the discharge coefficients C m , C d , C dD , C dD * , and C sc . The lowest percentage of data scatter is observed in Fig. 4b, d, which indicates that C d and C dD * have the most correlation with Fr 1 . However, the correlation of C dD * with Fr 1 is higher than C d . It is due to considering the water surface profile variation and V 1 in the adjusted Domínguez approach. The data points in Fig. 4c, e show more scatter, which means C dD , C sc , and Fr 1 depend on other dimensionless parameters. As observed in the sensitivity analysis results, the dependency of C dD and C sc on p/y 1 is greater than the dependency of these two coefficients on Fr 1 . According to Fig. 4b, d, the variations of C d and C dD * versus Fr 1 is almost linear. The absolute rate of these discharge coefficients versus Fr 1 equals 0.511 and 0.554, respectively. In contrast, the variations of C m , C dD , and C sc versus Fr 1 are nonlinear. According to Fig. 4a, c, and e, as Fr 1 increase, the absolute rate of C m , C dD , and C sc increases. Figure 6 shows the variations of discharge coefficients C m , C d , C dD , C dD * , and C sc versus Fr 1 for various side weir apex angles. Figure 6a shows that for Fr 1 ≥ 0.2, the discharge coefficients C m , C d , C dD , C dD * , and C sc increase slightly with the side weir apex angle from 60° to 120°, which confirms the findings of Kumar and Pathak (1987). It is also observed that with increasing Fr 1 , the variation of C m , C d , C dD , C dD * , and C sc increases slightly with the side weir angle. Comparison of Figs. 6a with 6b-e shows that the dependency of C m with θ at the higher Froude numbers is slightly more than other discharge coefficients. Based on trend lines in Fig. 6, when Fr 1 = 0.8, with increasing the side weir apex angle from 60° to 120°, the De Marchi coefficient increases by 0.09, while C d , C dD , C dD * , and C sc increase by only 0.06, 0.05, 0.03, and 0.05, respectively. According to Eq. (7), the effective length increases with increasing the side weir apex angle. Increasing the side weir's effective length intensifies the flow divagation angle and the outflow velocity. As a result, the discharge coefficient increases. Therefore, it is inferred that with increasing Fr 1 , the incremental impact of the weir length on discharge coefficients is more prevalent compared to decreasing effect of Fr 1 .
Another parameter influencing the discharge coefficient is p/y 1 . Depicted in Fig. 7a-e are the variations of C m , C d , C dD , C dD * , and C sc versus p/y 1 , respectively. It is deduced that the variation of p/y 1 with discharge coefficients can be considered linear. According to Fig. 7, C dD and C sc are slightly more correlated with p/y 1 than other discharge coefficients, which is logical based on the sensitivity analysis results (Table 3). In contrast, the discharge coefficients C m , C d , and C dD * are more influenced by Fr 1 than p/y 1 . Based on trend lines demonstrated in Fig. 7, with increasing p/y 1 , the discharge coefficients C m , C d , C dD , C dD * , and C sc increase with the rates of 0.197, 0.113,  192, 0.113, and 0.192, respectively. This trend of variations of discharge coefficient with p/y 1 is also reported by Singh et al. (1994) and Ghodsian (2004) for rectangular and triangular side weirs, respectively. Other investigators, such as Borghei et al. (1999) and Jalili and Borghei (1996), reported the decreasing effect of p/y 1 on the discharge coefficient of rectangular side weir. Therefore, more data are needed to explore this contradiction. The above analysis indicates a linear variation of p/y 1 with the discharge coefficients, while Fr 1 has a linear effect on C d and C dD * and a nonlinear effect on C m , C dD and C sc . So with these interpretations and employing the dimensionless parameters, Eqs. (37-41) are obtained to predict the discharge coefficients C m , C d , C dD , C dD * , and C sc for sharp-crested triangular side weirs. Equations (37-41) are valid for 0.03 ≤ Fr 1 ≤ 0.97, 0 < p/y 1 < 0.79, and 30 ≤ θ ≤ 120. Table 4 shows the statistical indices R 2 , MSE, and CP 10 % due to the above equations and those introduced by previous studies for predicting the discharge coefficients. The MSE indices due to Eq. (37) for different apex angles indicate the improved performance of this equation compared to previous equations. For the apex angle of 120°, owing to Eq. (37) for calculating C m , the MSE index is 0.0042, while Eqs. (12) and (10) provide MSE indices equal to 0.0058 and 0.0178, respectively. The values of statistical indices R 2 , MSE, MRE, and CP 10 % due to Eq. (37) are 0.410, 0.0045, 0.103, and 66.6 in estimating the De Marchi coefficient and, despite its comprehensiveness, this equation has a good performance compared to earlier equations. The results presented in Table 4 shows that the general Eqs. (38-41) produce suitable statistical indices in estimating the discharge coefficients C d , C dD , C dD * , and C sc . The MSE indices due to Eqs.  According to Table 4, Eq. (37) for apex angle 120° performs better in estimating the De Marchi coefficient than other equations. Table 5 shows improved MSE and CP 10 % indices due to Eq. (12) in calculating side weir discharge Q s . Therefore, a more accurate estimation of the De Marchi coefficient does not necessarily lead to a more precise calculation of Q s . Equation (37) leads to a lower MSE index in calculating Q s for other values of the side weir apex angle. For example, the MSE index due to Eq. (37), compared to that of Eqs. (13) and (15), in calculating the discharge of a side weir with 60° apex angle, is reduced by 67.10% and 3.88%, respectively. While equations presented by Ghodsian (2004) give a better CP 10 % in calculating the side weir discharge with apex angles of 60° and 30°. Equation (37) provides the MSE index of 4.581 in calculating the discharge of a triangular side weir. Table 5 shows that the statistical indices due to the De Marchi coefficients  and Eq. (37)] do not differ much. Thus, due to its comprehensiveness, Eq. (37) is preferred for computing the discharge of a triangular side weir. Table 5 reveals that the MSE index due to Eqs. (27,38,39,40 and 41) reduced by 5.26,21.20%,18.40%,21.22%,and 21.20%,respectively,compared with MSE due to Eq. (37). Also CP 10 % index due to the above equations increased by 1. 04, 8.50%, 8.3%, 10.17%, and 8.3%, respectively. To compare the accuracy of Eqs. (37-41) in computing the side weir discharge, the Nash-Sutcliffe Efficiency criteria (NSE) were also computed using Eq. (42) and compared with R 2 , MSE, MRE, and CP 10 %. The obtained NSE due to  and (27) in predicting Qs is 0. 968, 0.975, 0.974, 0.975, 0.974, and 0.970, respectively. Therefore, due to the statistical indices R 2 , NSE, MSE, MRE, CP 10 %, and NSE, it is clear that the TWE, Dominguez, adjusted Dominguez, Schmidt, and straightforward approaches have better performance than the De Marchi approach. However, it is necessary to mention that in the Dominguez, adjusted Dominguez, Schmidt, and straightforward approaches, the values of flow depths at Sects. "Introduction" and "Material and Methods" and apex angle is required for calculating side weir discharge. While in the TWE approach, only the flow condition at Sect. "Introduction" (Fig. 1) and apex angle are sufficient for calculating

Conclusion
This research investigated the effects of dimensionless parameters Fr 1 , p/y 1 , and θ on the triangular sharp-crested side weir discharge coefficients obtained from De Marchi (C m ), TWE (C d ), Domínguez (C dD ), adjusted Domínguez (C dD * ), and Schmidt (C sc ) approaches. Sensitivity analysis performed by the PSO-SVR method showed that Fr 1 with sensitivity indices equal to 1.89, 3.74, and 4.04 is the most significant parameter for estimating C m , C d , and C dD * , respectively. While p/y 1 with sensitivity index equal to 1.22 is the most important parameter for predicting C dD and C sc .
The results revealed that C m , C d , C dD , C dD * , and C sc decrease with increasing Fr 1 and increase with increasing p/y 1 and θ. The dependency of C m on θ was slightly more than the other discharge coefficients (C d , C dD , C dD * , and C sc ). At a constant value of p/y 1 , Δh/L has a significant correlation with Fr 1 .
New equations were proposed to estimate C m , C d , C dD , C dD * , and C sc for a sharp-crested triangular side weir. As an alternative to the De Marchi approach, all the methods used in this study (TWE, Domínguez, adjusted Domínguez, Schmidt, and straightforward approaches) produced better statistical indices R 2 , NSE, MSE, MRE (except Balahang and Ghodsian 2021), and CP 10 %. However, due to the lack of downstream flow depth y 2 in the TWE approach (Eq. 38), this method is more practical and introduced as the superior model by producing statistical indices R 2 = 0.975, NSE = 0.975, MSE = 3.610, MRE = 0.097, and CP 10 % = 71.36.
In the present study, it was assumed that the water surface varies linearly along the side weir, while the water surface in the subcritical condition varies as an ascending curve. It is suggested that the performance of the Schmidt approach in calculating triangular side weir outflow be evaluated by considering the nonlinearity of the water surface profile. Although the developed equations estimate the triangular side weir outflow with high precision, the obtained equations must be confirmed at the prototype scale.
Funding The authors received no specific funding for this work.
Data availability Some or all data, models, or code that support the findings of this study are available from the first author upon reasonable request.

Conflict of interest
The authors declare that they have no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.