E ﬀ ect of Height and Geometry of Stepped Spillway on Inception Point Location

: Air entrainment in a stepped spillway is very important to protect the spillway from cavitation damage. The inception point is the location where air starts entering the non-aerated ﬂow zone. The inception point location depends on di ﬀ erent parameters, such as the discharge, step height, and step shape. In this paper, various stepped spillways, including ﬂat steps, pooled steps, and round steps with di ﬀ erent step heights were numerically simulated using the volume of ﬂuid and realizable k- ε models. The results indicate that the inception point location moves downwards with the increase of the discharge of the stepped spillways. The length of the non-aerated ﬂow zone increases with the discharge. The inception point location moves downwards as the step height decreases and the step number increases at the same discharge. The inception point location of the round stepped spillway model is much closer to the spillway crest than that of the ﬂat stepped spillway with the same number of steps. The inception point location of the pooled stepped spillway is closer to the spillway crest than that of the ﬂat stepped spillway, but more downstream than that of the round stepped spillway.


Introduction
A stepped spillway is an important hydraulic structure widely used in hydraulic projects owing to its high energy dissipation function. Due to the increase of the energy dissipation rate, the stilling basin at the toe of the dam requires less size than that of a smooth spillway [1,2]. Aeration is a process in which air enters the surface of a spillway. Without aeration, the spillway may be subjected to cavitation damage. Thus, aeration is necessary to prevent cavitation damage in stepped spillways. Air entrainment starts when the turbulent boundary layer coincides with the water depth. The point at which air entrainment starts is called the inception point, i.e., the location where a non-aerated flow zone is converted into an aerated flow zone [3,4].
The performance of a stepped spillway depends on the air entrainment. Since cavitation often occurs in the non-aerated flow zone, measurement of the exact location of air entrainment is necessary to determine the area affected by cavitation. The inception point location indicates the length of the non-aerated flow zone. Thus, by measuring the exact location of the inception point, the aerated and non-aerated flow zones can be predicted [5][6][7]. The location of the inception point is also very important for measuring the length of the stilling basin, because energy dissipation in a stepped spillway depends on the location of the inception point. The energy dissipation downstream of the inception point is much larger than that upstream. In a smooth spillway, the location of the inception point is further downstream than in a stepped spillway. The location of the inception point in a stepped spillway depends on the discharge, the step height and step geometry [8,9].
In a stepped spillway, three kinds of flows occur: nappe flow regime, transition flow regime, and skimming flow regime. The nappe flow regime occurs at very low discharge, where free falling In a stepped spillway, three kinds of flows occur: nappe flow regime, transition flow regime, and skimming flow regime. The nappe flow regime occurs at very low discharge, where free falling jets fall in the form of nappes from one step to another. Skimming flow occurs at high discharge. In the skimming flow regime, the water moves above the pseudo-bottom as a coherent stream, where the pseudo-bottom is formed by connecting the edges of the steps. Beneath the pseudo-bottom, recirculating vortices are formed causing energy dissipation in the stepped spillway [10,11].
Skimming flow in a stepped spillway consists of two regions: non-aerated flow zone and aerated flow zone. The line connected to the edges of the step is called the pseudo-bottom. Beneath the pseudo-bottom, recirculating vortices are formed [12,13]. In Figure 1, Li is the distance from the spillway crest to the inception point. The less Li is, the greater is the amount of air that enters the water and vice versa. At the inception point, the pseudo-bottom air concentration is 0.01 [14] and the air phase suddenly changes from entrapped to entrained. Entrained air is transported within the flow, while entrapped air is transported above the water surface and between the crest. The sum of the entrained and entrapped air in a stepped spillway is called the total air concentration. A 5%-8% pseudo-bottom air concentration is required to prevent cavitation damage [15]. The point at which the pseudo-bottom air concentration reaches 5%-8% is called the critical point [16]. Beyond the inception point, the air concentration increases rapidly at the pseudo-bottom and quickly reaches the critical point, and then increases to the maximum value. Air entrainment can prevent cavitation damage caused by the occurrence of low pressure [17]. On increasing the discharge, the length of the non-aerated flow zone and the risk of cavitation both increase. The zone near the inception point is closely related to cavitation damage. The inception point in a stepped spillway is observed when the boundary layer thickness reaches 80% of the flow depth. The flow upstream of the inception point is smooth, while it is very turbulent downstream [18,19]. A mathematical formula was developed to measure the length of the inception point from the spillway crest (Li) for determining the start of air entrainment with slope of 22° [20]. Alghazali and Jasmin [21] experimentally measured the inception point location of different stepped spillway configurations. In their study, 12 different empirical equations were derived based on the experimental results for different configurations of stepped spillways. Hua et al. [22] experimentally measured the inception point of a stepped spillway model with a slope of 39.29°. Different models of different step heights were used and an empirical formula was derived based on their experiments.
Energy dissipation in a round-shaped stepped spillway increases by 3% compared with that in a flat stepped spillway. A round stepped spillway dissipates more energy due to very quick air entrainment of the skimming flow. Hence, the length of the training wall height in a round stepped spillway decreases by approximately 20%. The risk of cavitation in a round stepped spillway is also less than that of a flat stepped spillway [23].
Felder and Chanson [24] conducted experiments on pooled and flat stepped spillways with slopes of 8.9° and 26.6°. The results showed that energy dissipation in the pooled stepped spillway is larger than in the flat stepped spillway with a smaller slope. In the pooled stepped spillway, the mean air concentration is larger than that in the flat stepped spillway. Munta and Otun [25]  On increasing the discharge, the length of the non-aerated flow zone and the risk of cavitation both increase. The zone near the inception point is closely related to cavitation damage. The inception point in a stepped spillway is observed when the boundary layer thickness reaches 80% of the flow depth. The flow upstream of the inception point is smooth, while it is very turbulent downstream [18,19]. A mathematical formula was developed to measure the length of the inception point from the spillway crest (Li) for determining the start of air entrainment with slope of 22 • [20]. Alghazali and Jasmin [21] experimentally measured the inception point location of different stepped spillway configurations. In their study, 12 different empirical equations were derived based on the experimental results for different configurations of stepped spillways. Hua et al. [22] experimentally measured the inception point of a stepped spillway model with a slope of 39.29 • . Different models of different step heights were used and an empirical formula was derived based on their experiments.
Energy dissipation in a round-shaped stepped spillway increases by 3% compared with that in a flat stepped spillway. A round stepped spillway dissipates more energy due to very quick air entrainment of the skimming flow. Hence, the length of the training wall height in a round stepped spillway decreases by approximately 20%. The risk of cavitation in a round stepped spillway is also less than that of a flat stepped spillway [23].
Felder and Chanson [24] conducted experiments on pooled and flat stepped spillways with slopes of 8.9 • and 26.6 • . The results showed that energy dissipation in the pooled stepped spillway is larger than in the flat stepped spillway with a smaller slope. In the pooled stepped spillway, the mean air concentration is larger than that in the flat stepped spillway. Munta and Otun [25] performed 40 experiments on three different stepped spillway models to study the relationship among the inception length, the discharge, and the chute angle. It was found that the inception length increases with the discharge but decreases with the chute angle.
Physical modelling of a stepped spillway is very difficult and does not always yield accurate results to develop physical models for experiments. Fortunately, with the development of high-performance computers, the computational fluid dynamics (CFD) technique has been developed. Chen et al. [26] found that CFD is a very authentic and reliable source to simulate the flow over a stepped spillway. The k-ε turbulence model was used to simulate turbulent flow in a stepped spillway. Imam and Mehdi [27] studied the energy dissipation in a stepped spillway by numerical modelling with different parameters and investigated the effects of these parameters on energy dissipation. The parameters involved in energy dissipation include the step height, the number of steps, the spillway slope and the unit discharge. Afshin and Mitra [28] used the well-known commercial software ANSYS Fluent with different numerical models. They used the volume of fluid (VOF) and the mixture model to simulate the flow in a stepped spillway and evaluated the performance of each model to determine which one accurately simulated the skimming flow over a stepped spillway. Parsaie et al. [29] performed numerical simulation with the help of turbulent modelling and concluded that turbulent modelling is a very good and efficient method to determine the complexity of stepped spillway flow. In their study, Flow 3D software was used for numerical simulation. Qian et al. [30] simulated a spillway model with the help of four different turbulence models: realizable k-ε, k-ω, shear stress transport k-ω, and large eddy simulation (LES). Spillway models were simulated using these four turbulence models and the results compared. The results obtained with the realizable k-ε model were more accurate than the other three turbulence models. Benmamar et al. [31] conducted a numerical simulation and developed a numerical model for determining the boundary layer in a stepped spillway with a steep slope. The numerical model was based on the implicit finite difference scheme. Cheng et al. [32] measured the air-water volume fraction of two-phase flow using k-ε models. Dong and Lee [33] used the VOF multiphase model to simulate the stepped spillway model with slope of 10 • and to determine the characteristics of the skimming flow. Their results included velocity distribution, air concentration, and pressure distribution. It was concluded that the VOF model accurately simulates the flow pattern over a stepped spillway. Bombardelli et al. [34] used the VOF and renormalized group (RNG) turbulence model to simulate the stepped spillway with a slope of 53 • . They obtained the velocity distribution, and turbulent kinetic energy. Tabbara et al. [35] numerically simulated the flow over a stepped spillway using different step configurations, using the ADINA software with the k-ε flow model. In the prediction of the water surface profile and energy dissipation values, the results showed good agreement between the numerical and experimental values. Cheng et al. [36] used ANSYS Fluent software to simulate the flow over a stepped spillway. A mixture flow model was used as the numerical model and the RNG k-ε model was used as the turbulence model. They studied the skimming flow region and determined the interaction between cavity recirculation and air bubbles present in the skimming flow zone. The velocity distribution on the steps was also determined in their investigation. Mohammad et al. [37] used the VOF model and the turbulence models RNG and LES to simulate the interaction between air and water using stepped spillway models with a steep slope. Bai and Zhang [38] simulated a different type of stepped spillway called the V-type stepped spillway. They analyzed the pressure distribution by five different turbulence models and compared it with the physical values. By comparing the values with the physical model, it was suggested that the realizable k-ε model gave better results for simulating the pressure distribution than the others turbulence models. Li and Zhang [39] simulated pooled stepped spillways with four types of pool weirs (full pool, full pool and two-sided pool, full pool and central pool, two-sided and central pool) using the VOF and RNG turbulence models. It was suggested that the energy dissipation rate in the different types of pooled weirs followed the order, two-sided and central pool, full pool and central pool, full pool and two-sided pools, and full pool. Cheng et al. [40] measured the logarithmic velocity profile in a stepped spillway using numerical simulation and compared it with experiments. According to their researches, the air inception location over stepped spillways exits until the boundary layer thickness is between the range of 0.72 and 0.79. Ljubicic et al. [41] studied the adverse slope stilling basins with stepped chutes. They studied the momentum, the length of the hydraulic jump, and the energy dissipation by the numerical model, and then compared the results with the experimental data. Wan et al. [42] predicted the cavitation damage in a high-speed smooth spillway using numerical simulation. VOF and the standard k-turbulence model were used to simulate the high-speed flow and the potential cavitation region in the high-speed smooth spillway. The numerical results and the experimental measurement showed a good agreement.
The main aim of this research was to study the effect of different parameters, i.e., the discharge, step height, and step geometry, on the inception point location by numerical simulations. Different numerical experiments were conducted on different step spillway models using multiphase VOF combined with the realizable k-ε turbulence model. The different step spillway models differ in the various step heights and step geometry configurations, including flat, pooled and round. The data shows how the inception point location varies with the discharge, step height, and step geometry. The results improve the understanding of stepped spillways with respect to air entrainment, and provide a reference for optimal design of proper air entrainment for the stepped spillway. The results also identify the step geometry that provides more aeration and earlier air entrainment into the water to reduce cavitation damage in stepped spillways.

Materials and Methods
Numerical modelling is a widely applied technique performed with the aid of a computer to solve the Navier-Stokes equation, which is based upon the conservation of mass, energy, and momentum. CFD software is used to simulate the stepped spillway models. Fluent software uses the finite volume method (FVM) to discretize the Navier-Stokes equation. The FVM method is the most commonly used method to discretize the dynamic equations in CFD. The control volume technique converts the governing equation into algebraic equations rather than solving the equations numerically. Integration equations for the control volume are solved in an implicit form.
Six kinds of stepped spillway models with different step sizes and configurations used for numerical simulation are listed in Table 1. The models varying according to step height with flat, pooled, and round steps, are shown in Figure 2. All the six models were simulated respectively, using VOF as the multiphase model and the realizable k-ε model as the turbulence model. The slope of all the models is held constant.
The inlet boundary condition of the stepped spillway model is represented by the discharge. In the initial state, there is no water in the spillway and the spillway surface is directly connected to the atmosphere. The outlet boundary condition of the spillway is set by pressure. All the walls are stationary with no slip velocity. The three-dimensional grid models and boundary conditions used for computation are shown in Figure 3. The results of the three-dimensional model are better than those of the two-dimensional model because the three-dimensional model consider the air phase more accurately than the two-dimensional model.
The calculation domain is discretized into an unstructured grid with 194,447 mesh elements and the number of nodes used for computation is 215,760. The inflation layer is applied to the steps to yield better results. The area around the steps is particularly considered to achieve better results. Therefore, 10 layers of inflation are applied on the steps to better study the flow of the spillway. Table 1. Design parameters of different stepped spillways for numerical simulation.

Model
Step Geometry Step Height(m) Step

VOF (Volume of Fluid)
VOF is the multiphase model developed by Hirt and Nichols [43]. It is usually used when two or more phases are involved. In this model, it is assumed that the different phases will not intermingle with each other. In free surface flow, as dealt with in the present study, the purpose is to track the interface between air and water. Therefore, it is appropriate to use VOF especially when the position of the interface between the different phases is the point of interest. The VOF method tracks the interface as a mixture cell, while other multiphase models only focus on bubbles tracking. Previous literature indicates [44] that the results from other multiphase models are unsatisfactory when the interface is considered. Thus, in the present study, the VOF model is used to accurately predict the interface between water and air. When the VOF model involves two phases, water and air, then the volume fraction of fluid can be regarded as 0 if the cell does not contain any water, while 1 represents it is fully filled with water. If the value is between 0 and 1, it indicates that it is a mixture cell, seen as an interface between water and air, and a free surface between water and air. Here, the factor α is introduced, which indicates the value of the phases in each cell. The factor α a is the volume fraction of air and α w is the volume fraction of water.
Equation (1) indicates that if the value of α is between 0 and 1, then there is an interface between water and air, and it is a free surface that can be tracked by the VOF method.
The properties of any cell are either representative of one phase or the mixture. All the properties of a cell such as velocity, pressure, and temperature are shared due to the concept of volume fraction. The total volume fraction of all the phases in a cell is equal to unity. The density ρ of the cell is calculated from Equation (2). The meaning of symbols in the equations is shown at the end of the text.
The VOF method solves the continuity and momentum equations and measures the volume fraction of each phase by tracking Equations (3)- (5).

Realizable k-ε Model
The steps in a spillway have a significant effect on increasing the roughness. Owing to the effect of roughness, the flow over a stepped spillway is highly turbulent. Two-phase flow over a stepped spillway can generate an even higher turbulence rate. Thus, a turbulent model is used to predict the turbulent behavior. The realizable k-ε model was introduced by Shih et al. [45] for high Reynolds number turbulent flow, which is based upon the realizable constraints. The realizable k-ε model gives better results than other turbulence models because a new equation for the dissipation rate (ε) is proposed, which is based upon the large turbulent Reynolds number. The performance of the realizable k-ε model is improved and provides good results for the recirculating flow. Other Reynolds-averaged Navier-Stokes (RANS) models such as the standard k-ε model, RNG k-ε, standard k-ω, and shear stress transport k-ω models perform well when used for simulating other types of flows like pipe or laminar flows. Because of the unsatisfactory results from other RANS models for recirculating flow, considering that the flow in the stepped spillway is recirculating, realizable k-ε is used in the present study. The transport equations for turbulent kinetic energy (k) and turbulent dissipation rate (ε) are expressed as Equations (6) and (7), respectively. The meaning of all the symbols in Equations (6) and (7) is shown at the end of the text.

Results and Discussions
Various stepped spillway models with different step heights and configuration were simulated to predict the location of the inception point in the stepped spillways. All the stepped spillway models were simulated at different discharge rates (0.8, 1, 1.2, and 1.4).

Effect of Step Height on Inception Length Using Flat Stepped Spillways
Four stepped spillway models with different step heights (0.10, 0.09, 0.05, and 0.04 m) were simulated to predict the changes in inception length with the decrease in step height. The inception point is closer to the spillway crest with a larger step height under the same discharge. The results also reveal that the inception point length increases with the discharge in all stepped spillway models. The comparison between the VOF results of the inception point lengths with different step heights when the discharge is 0.8 m 2 /s is shown in Figure 4. It is seen clearly that with the decrease in step height, the inception point moves downward under the same discharge.
Reynolds-averaged Navier-Stokes (RANS) models such as the standard k-ε model, RNG k-ε, standard k-ω, and shear stress transport k-ω models perform well when used for simulating other types of flows like pipe or laminar flows. Because of the unsatisfactory results from other RANS models for recirculating flow, considering that the flow in the stepped spillway is recirculating, realizable k-ε is used in the present study. The transport equations for turbulent kinetic energy (k) and turbulent dissipation rate (ε) are expressed as Equations (6) and (7), respectively. The meaning of all the symbols in Equations (6) and (7) is shown at the end of the text.

Results and Discussions
Various stepped spillway models with different step heights and configuration were simulated to predict the location of the inception point in the stepped spillways. All the stepped spillway models were simulated at different discharge rates (0.8, 1, 1.2, and 1.4).

Effect of Step Height on Inception Length Using Flat Stepped Spillways
Four stepped spillway models with different step heights (0.10, 0.09, 0.05, and 0.04 m) were simulated to predict the changes in inception length with the decrease in step height. The inception point is closer to the spillway crest with a larger step height under the same discharge. The results also reveal that the inception point length increases with the discharge in all stepped spillway models. The comparison between the VOF results of the inception point lengths with different step heights when the discharge is 0.8 m 2 /s is shown in Figure 4. It is seen clearly that with the decrease in step height, the inception point moves downward under the same discharge.  The inception lengths (Li) with respect to discharge, step height, and Froude surface roughness (F) are listed in Table 2. With increase of discharge, the inception length increases under the conditions of the same number of steps. With the increase in the number of steps, the inception The inception lengths (Li) with respect to discharge, step height, and Froude surface roughness (F) are listed in Table 2. With increase of discharge, the inception length increases under the conditions of the same number of steps. With the increase in the number of steps, the inception length also increases under the same discharge. The Froude surface roughness, defined with respect to surface roughness (ks = hcosθ), which depends on unit discharge (q) and step height, is illustrated as Equation (8). The symbols in the equation are listed at the end of the text.  The variation of the inception length (Li) with respect to the discharge at different step heights is shown in Figure 5. The trend of the figure shows that at the highest discharge, the inception length is the maximum under all the conditions of different step heights. With the decrease in step height, the inception length increases under the same discharge. Table 2. Summary of unit discharge (q), step height, Froude surface roughness (F), inception point length (Li), surface roughness (ks) and normalized (Li The variation of the inception length (Li) with respect to the discharge at different step heights is shown in Figure 5. The trend of the figure shows that at the highest discharge, the inception length is the maximum under all the conditions of different step heights. With the decrease in step height, the inception length increases under the same discharge. Step Height 0.1m Step Height 0.09m Step Height 0.05m Step Height 0.04m  The relationship between inception length (Li) and Froude surface roughness (F) at different step heights is shown in Figure 6. With the increase in Froude surface roughness, the inception length increases. According to Equation (8), Froude surface roughness (F) increases with the decrease in step height and the increase in discharge. It is also apparent from Figure 6 that the Froude surface roughness (F) values increase with the decrease of the step height. The relationship between inception length (Li) and Froude surface roughness (F) at different step heights is shown in Figure 6. With the increase in Froude surface roughness, the inception length increases. According to Equation (8), Froude surface roughness (F) increases with the decrease in step height and the increase in discharge. It is also apparent from Figure 6 that the Froude surface roughness (F) values increase with the decrease of the step height.

Li F
Step Height 0.1m Step Height 0.09m Step Height 0.05m Step Height 0.04m Figure 6. Variation of inception length (Li) with Froude surface roughness (F) at different stepped heights of flat stepped spillway models.
The relationship between normalized Li and Froude surface roughness is shown in Figure 7. Normalized Li is the ratio of the inception length (Li) and Froude surface roughness (ks). Figure 7 shows that the normalized Li increases with the Froude surface roughness.

Li/ks F
Step Height 0.1m Step Height 0.09m Step Height 0.05m Step Height 0.04m

Effect of Step Geometry on Inception Length Using Different Step Shapes
In order to analyze the effect of step geometry on inception length, three different step geometry models (flat, pooled, and round) with the same number of steps were simulated to obtain the inception length in these three different stepped spillways. Conditions of the water volume fraction contours of the flat, pooled, and round stepped spillway models with the same number of The relationship between normalized Li and Froude surface roughness is shown in Figure 7. Normalized Li is the ratio of the inception length (Li) and Froude surface roughness (ks). Figure 7 shows that the normalized Li increases with the Froude surface roughness.
The relationship between inception length (Li) and Froude surface roughness (F) at different step heights is shown in Figure 6. With the increase in Froude surface roughness, the inception length increases. According to Equation (8), Froude surface roughness (F) increases with the decrease in step height and the increase in discharge. It is also apparent from Figure 6 that the Froude surface roughness (F) values increase with the decrease of the step height.

Li F
Step Height 0.1m Step Height 0.09m Step Height 0.05m Step Height 0.04m Figure 6. Variation of inception length (Li) with Froude surface roughness (F) at different stepped heights of flat stepped spillway models.
The relationship between normalized Li and Froude surface roughness is shown in Figure 7. Normalized Li is the ratio of the inception length (Li) and Froude surface roughness (ks). Figure 7 shows that the normalized Li increases with the Froude surface roughness. Step Height 0.09m Step Height 0.05m Step Height 0.04m

Effect of Step Geometry on Inception Length Using Different Step Shapes
In order to analyze the effect of step geometry on inception length, three different step geometry models (flat, pooled, and round) with the same number of steps were simulated to obtain the inception length in these three different stepped spillways. Conditions of the water volume fraction contours of the flat, pooled, and round stepped spillway models with the same number of

Effect of Step Geometry on Inception Length Using Different Step Shapes
In order to analyze the effect of step geometry on inception length, three different step geometry models (flat, pooled, and round) with the same number of steps were simulated to obtain the inception length in these three different stepped spillways. Conditions of the water volume fraction contours of the flat, pooled, and round stepped spillway models with the same number of steps and unit discharge rate of 0.8 m 2 /s and 1.0 m 2 /s were simulated. The inception point locations are shown in Figure 8. The location of the inception point in the round stepped spillway model is near the spillway crest. In the pooled stepped spillway model, the location of the inception point is further away from the spillway crest compared with the round stepped, but nearer than in the flat stepped spillway model. The location of the inception point in the flat stepped spillway model is far away from the spillway crest compared with the round or the pooled stepped spillway models under the same discharge. Thus, the location of the inception point moves toward the spillway crest with changes in step geometry from the flat to the pooled, to the round stepped spillway. The distances of the inception point from the crest Li of the stepped spillway models with different step geometry under the same number of steps and discharge are listed in Table 3    Li increases when the stepped geometry changes from round to pooled, and further increases when the geometry changes from pooled to flat, as shown in Figure 9. The round stepped spillway has the lowest Li compared with both the pooled stepped and the flat stepped spillways. The Li of the flat stepped spillway is greater than that of the pooled stepped spillway under the same unit discharge. For all the stepped spillway geometries, Li increases with unit discharge.
Li increases when the stepped geometry changes from round to pooled, and further increases when the geometry changes from pooled to flat, as shown in Figure 9. The round stepped spillway has the lowest Li compared with both the pooled stepped and the flat stepped spillways. The Li of the flat stepped spillway is greater than that of the pooled stepped spillway under the same unit discharge. For all the stepped spillway geometries, Li increases with unit discharge.

Conclusions
The inception point location of the stepped spillways depends on the discharge, step height, and step geometry. The distance of the inception point from the crest (Li) has a positive correlation with the discharge. Under the same discharge, the step height has a negative correlation with Li, while the number of steps shows a positive correlation with it. The inception point location is closer to the spillway crest in the stepped spillways with higher step height. The value of Li changes with the Froude surface roughness (F). Li increases with F, while F increases with the unit discharge but decreases with the step height.
Various simulations and comparison among the three kinds of spillway forms (round, pooled, and flat) were conducted and analyzed. It was shown that under conditions of the same discharge and number of steps, the location of the inception point in a round stepped spillway is closest to the spillway crest, while it is farthest in a flat stepped spillway. As for the non-aerated flow zone, it decreases with the stepped height, as does the probability of cavitation. Under the same boundary conditions, the round stepped spillway shows the least non-aerated flow zone, while the flat form shows the most. Therefore, the minimum probability of cavitation damage may occur in the round stepped spillway, and thus it is better to use the round stepped spillway to prevent cavitation.

Conflicts of Interest:
The authors declare no conflict of interest.

Conclusions
The inception point location of the stepped spillways depends on the discharge, step height, and step geometry. The distance of the inception point from the crest (Li) has a positive correlation with the discharge. Under the same discharge, the step height has a negative correlation with Li, while the number of steps shows a positive correlation with it. The inception point location is closer to the spillway crest in the stepped spillways with higher step height. The value of Li changes with the Froude surface roughness (F). Li increases with F, while F increases with the unit discharge but decreases with the step height.
Various simulations and comparison among the three kinds of spillway forms (round, pooled, and flat) were conducted and analyzed. It was shown that under conditions of the same discharge and number of steps, the location of the inception point in a round stepped spillway is closest to the spillway crest, while it is farthest in a flat stepped spillway. As for the non-aerated flow zone, it decreases with the stepped height, as does the probability of cavitation. Under the same boundary conditions, the round stepped spillway shows the least non-aerated flow zone, while the flat form shows the most. Therefore, the minimum probability of cavitation damage may occur in the round stepped spillway, and thus it is better to use the round stepped spillway to prevent cavitation.