A Review on Existing Methods to Assess Hydraulic Erodibility Downstream of Dam Spillways

: In recent years, rock scouring or erosion downstream of dams has become an increasing dam safety concern. Several theoretical, semi-theoretical, semi-analytical and numerical methods can be used to assess the rock erosion in hydraulic structures. Semi-theoretical approaches determine the correlation between the erosive intensity of ﬂuid ﬂow and the resistive capacity of rock. Such approaches establish the scour thresholds as a function of erosive intensity of water and several rock mass indices by using in situ data and a curve-ﬁtting approach. In some studies, the excavatability index, initially developed for rock mass stability analysis, was used to analyse the rock mass resistance in hydraulic erodibility analysis. The effectivity and weight of the geomechanical parameters used are yet to be determined on the basis of the erodibility phenomenon. The semi-analytical methods are developed on the basis of the mechanical and hydraulic interaction of rock mass and water. Four methods developed by Bollaert et al. are important in determining the erodibility in the plunge pool, but they are not applicable in the case of spillways. They used the comprehensive fracture mechanics for closed-end joints, quasi-steady impulsion, and dynamic impulsion ( DI ) for blocky rock erosion. The application of these methods to each site is necessary to identify constants that are difﬁcult to determine. Few numerical methods are available to assess the rock mass erosion in hydraulic structures. In the case of numerical methods, the erosive agent is indistinct, and the hydraulic hazard parameter on the spillway surface is almost challenging to apply. This study comprehensively reviews the mechanism of erosion and the methods for assessing the risk of potential rock mass erosion downstream of dams and hydraulic structures. The advantages and disadvantages of all methods are discussed and the potential future research directions in this domain are proposed.


•
Existing methodologies for investigating the hydraulic erodibility of rock in dam spillways and plunge pools are reviewed comprehensively.

•
Hydraulic erodibility or scouring downstream of dams can be assessed by theoretical, semi-theoretical, semi-analytical, and numerical methods.

•
The advantages and disadvantages of existing methods for evaluating erosion and scouring are summarized.

•
Future directions for assessing the hydraulic erodibility of rock in hydraulic structures (dams) are discussed.

Introduction
In recent years, the risk of rock scouring or erosion in hydraulic structures increased the concern about dam safety.The terms 'erodibility', 'scour', and 'hydraulic erosion' are considered synonymous technical words to explain critical centralised erosion that occurs when the erosive intensity of fluid surpasses the resistive capacity of the rock mass.The

Introduction
In recent years, the risk of rock scouring or erosion in hydraulic structures increased the concern about dam safety.The terms 'erodibility', 'scour', and 'hydraulic erosion' are considered synonymous technical words to explain critical centralised erosion that occurs when the erosive intensity of fluid surpasses the resistive capacity of the rock mass.The design of dam spillways that can discharge a wide range of floods with the minor scouring of material is one of the important challenges in designing hydraulic structures [1,2].
Rock scouring is a complicated mechanism.Rock scour basically occurs by four principal mechanical processes as revealed by its occurrence near engineering structures: 1. Rock block abrasion (long term); 2. Intact rock fracturing (brittle failure and fatigue failure); 3. Single block deletion (uplift pressures in fractures and shear forces); 4. Rock block peeling off.These mechanisms of rock erosion are presented in Figure 1.The instantaneous processes consist of rock block removal and brittle failure.In the case of block removal mechanism, single rock blocks can be lifted by horizontal and vertical pressures.Along these lines, the stream turbulence, block size, block dimensions, bulge of the blocks, surrounding rock mass, joint roughness, and cohesion of filling as a shear destabilised force are relevant and effective parameters that straightforwardly characterise the lifting forces.
The breaking of rock in fractured media is scientifically explained by the hypothesis of linear elastic fracture mechanics and may occur by abrupt or progressive hydraulic-driven fracturing.The stress intensity at the tip of shut-end cracks and fracture The instantaneous processes consist of rock block removal and brittle failure.In the case of block removal mechanism, single rock blocks can be lifted by horizontal and vertical pressures.Along these lines, the stream turbulence, block size, block dimensions, bulge of the blocks, surrounding rock mass, joint roughness, and cohesion of filling as a shear destabilised force are relevant and effective parameters that straightforwardly characterise the lifting forces.
The breaking of rock in fractured media is scientifically explained by the hypothesis of linear elastic fracture mechanics and may occur by abrupt or progressive hydraulic-driven fracturing.The stress intensity at the tip of shut-end cracks and fracture toughness of the rock play important roles in identifying the types of fracturing; if the stress intensity surpasses the fracture toughness, then brittle failure will occur; otherwise, dynamic fracturing of rock or fatigue failure will occur [4].Fatigue failure spreads a current fracture on the basis of the intensity and the number of pressure pulses and stops when fracture formation is finished.
Water 2021, 13, 3205 3 of 42 In the case of the block peeling-off mechanism, quasi-steady pressure forces play an important role.This type of erosion ordinarily occurs in thin rock layers with protruded open-end rock blocks.The block protrusion causes the stream deviation.The deviation of the stream creates the lifting and drag forces on the block because of the quasi-steady stream velocity, which causes block ejection [3].
Figure 1 demonstrates the other time-dependent rock mass scouring mechanism that occurs on the rock surface by abrasion phenomenon.Abrasion, occasionally referred to as ball-milling or bedrock wear, alludes to gradual grinding resulting from the repeated impacts due to different particles, such as sands or cobbles that are carried by flowing fluid.The occurrence of scouring because of abrasion commonly requires a large amount of time.Thus, dam overtopping or spillway erodibility evaluations are rarely conducted.
These mechanisms are used for the development of the methods for the prevision of rock mass erosion.
This review paper attempts to provide an extensive survey of the existing methods to assess the rock mass erosion in several types of hydraulic structures.We extract the advantages and disadvantages of the existing semi-theoretical, semi-analytical, and numerical approaches.A comparison between the semi-theoretical methods is conducted on the basis of the committed error of each method by considering real datasets and the prediction of each methodology.To efficiently explain the existing methods, we present the algorithm of each methodology in various flowcharts.Some potential future research directions are presented on the basis of this review.Finally, the limitations of the application of some methods in erosion prediction are defined, and the new necessary research in this domain is established.
In this work, the semi-theoretical methods are first explained.Thereafter, the semianalytical methodologies are explained by mainly focusing on the comprehensive scour model (CSM).The scouring at the multiblock system with numerical methods is explained.The summary of the paper and conclusion of the authors are presented in the last two sections of this work.

Scour Prediction Approaches
Many authors have studied scour hole development below plunging spillway jets [5][6][7][8].The disintegration capability of soil along inclines or unlined spillway channels is also considered in several studies [9][10][11][12].Amongst the approaches used for assessing rock mass erosion, the most common are semi-theoretical methods according to the correlation between the resistive capacity of rock and the erosive intensity of fluid, such as Annandale's and Pells' approach [1].Some classical equations can be used for computing the final scour depth in plunge pools.The flowchart in Figure 2 shows the classifications of the empirical, theoretical, and semi-theoretical methods.

Comparative Methods Based on Resistance of Materials for Erosion and Rate of Energy Dissipation
A certain form is considered in the design of a dam spillway to evaluate the potential of hydraulic erosion of dams.The scour threshold is an erodibility index and a function of P that is widely used in the industry [13][14][15].This function is calculated for different comparative methods on the basis of the interpreted erosion obtained from various case studies.A threshold line with a specific scour condition that distinguishes case studies can be identified in the plotted data.The correlation between the resistance of material for erosion and the energy dissipation rate P is calculated as in Equation (1): where P is the total water erosive intensity that is released on or over the material.Material resistance to erosion depends on its strength and can be explained by existing methods that determine the overall component commitment.In P < f (K h ), the limit of erodibility is not surpassed, while P > f (K h ) indicates that the values above the erodibility threshold are exceeded; accordingly, erosion will occur, and the material will scour [16][17][18].
Water 2021, 13, x FOR PEER REVIEW 4 of 44

Comparative Methods based on Resistance of Materials for Erosion and Rate of Energy Dissipation
A certain form is considered in the design of a dam spillway to evaluate the potential of hydraulic erosion of dams.The scour threshold is an erodibility index and a function of P that is widely used in the industry [13][14][15].This function is calculated for different comparative methods on the basis of the interpreted erosion obtained from various case studies.A threshold line with a specific scour condition that distinguishes case studies can be identified in the plotted data.The correlation between the resistance of material for erosion and the energy dissipation rate P is calculated as in Equation ( 1): where P is the total water erosive intensity that is released on or over the material.Material resistance to erosion depends on its strength and can be explained by existing methods that determine the overall component commitment.In P < f(Kh), the limit of erodibility is not surpassed, while P > f(Kh) indicates that the values above the erodibility threshold are exceeded; accordingly, erosion will occur, and the material will scour [16][17][18].Equation ( 1) is a two-part equation representing the erosive force and the material resistance threshold to scouring, which is determined by the erodibility index, rock mass erosion index (RMEI), or geological strength index (GSI) (eGSI).In this work, an improved strategy is needed to find an agent that is easy to be determined while explaining the overall strength of the fluctuating turbulence.Hydraulic water energy (kW/m 2 ) denotes the erosive force of stream water and is represented as hydraulic stream power (Pa) Equation ( 1) is a two-part equation representing the erosive force and the material resistance threshold to scouring, which is determined by the erodibility index, rock mass erosion index (RMEI), or geological strength index (GSI) (eGSI).In this work, an improved strategy is needed to find an agent that is easy to be determined while explaining the overall strength of the fluctuating turbulence.Hydraulic water energy (kW/m 2 ) denotes the erosive force of stream water and is represented as hydraulic stream power (P a ) or unit stream power dissipation Π UD .The favoured agent for this purpose is the rate of energy dissipation.Hydraulic turbulence is the main cause of energy loss and pressure fluctuations.Increased turbulence intensity leads to a great energy dissipation rate and high levels of water pressure fluctuations.Erosive force and pressure fluctuations of flowing water are evaluated and calculated, as in the following equations, in energy dissipation rate analysis: where Π UD is the unit stream power dissipation and the change in unit stream power between two locations along a flow path (kW/m 2 ); γ = unit weight of water (kN/m 3 ); q = unit discharge rate (m 2 /s); ∆E = energy loss (J = kg•m 2 •s −2 ); h l = head loss; β f = width of a flow section measured at the water surface (m); and S f = dh l dx and presents the friction slope or the gradient of the total hydraulic energy line.Figure 3a shows the variation of properties over a differential flow section in an open channel.Figure 3b shows the cross-sectional view of the flow channel.

𝛽 𝑑𝑥
= ( ) where  is the unit stream power dissipation and the change in unit stream power between two locations along a flow path (kW/m 2 );  = unit weight of water (kN/m 3 ); q = unit discharge rate (m 2 /s); ΔE = energy loss (J = kg•m 2 •s −2 ); hl = head loss; x = horizontal direction of flow; y = depth of water; Zb = the vertical distance of the bottom surface from the reference datum;  = water density (kg/m 3 ); V = velocity of flowing water; g = gravitational acceleration (m/s 2 ); Q = total discharge amount (m 3 /s);  = width of a flow section measured at the water surface (m); and Sf = and presents the friction slope or the gradient of the total hydraulic energy line.Figure 3a shows the variation of properties over a differential flow section in an open channel.Figure 3b shows the cross-sectional view of the flow channel.This equation was used for evolving equations to assess the energy dissipation rate for a characteristic manner of erosion [1,20].Annandale provided alternate forms of Equation (3) to decide the erosive limit concerning an assortment of flow cases, including open channels, plunge pools, hydraulic jumps, knick-points, and head cuts [21].
Part two of Equation (1) shows the rock resistance capacity, which is calculated by using the rock mass parameters and forms an index when inserted into the equation.In this regard, different categories of engineering rock mass that are used in evaluating the rock mechanical excavatability apply a large number of agents that modify the rock hydraulic erodibility [22,23].Van Schalkwyk used some rock mass indices for this purpose, including the rock mass rating (RMR) classification of Bieniawski; the Q-classification of This equation was used for evolving equations to assess the energy dissipation rate for a characteristic manner of erosion [1,20].Annandale provided alternate forms of Equation (3) to decide the erosive limit concerning an assortment of flow cases, including open channels, plunge pools, hydraulic jumps, knick-points, and head cuts [21].
Part two of Equation (1) shows the rock resistance capacity, which is calculated by using the rock mass parameters and forms an index when inserted into the equation.In this regard, different categories of engineering rock mass that are used in evaluating the rock mechanical excavatability apply a large number of agents that modify the rock hydraulic erodibility [22,23].Van Schalkwyk used some rock mass indices for this purpose, including the rock mass rating (RMR) classification of Bieniawski; the Q-classification of Barton, which was primarily suggested for underground space designs; and the index suggested by Kirsten for evaluating the excavatability of earth materials [24][25][26][27][28].The authors reported similar results for all these indices.However, the highest accuracy was found for the Kirsten index (N), which is used to evaluate the erosion of rock mass [29][30][31].
Equation ( 5) is determined by using certain rock mass parameters, such as uniaxial compressive strength (UCS), mass strength (M s ), joint shear strength (K d ), relative block structure (J s ), and rock block size (K b ).The proposed system applies across the full range of natural materials from the weakest soils to the hardest rocks.
In 1982, Kirsten used the classifications and charts of Jennings to propose a descriptive chart with M s ranging from 0.87 to 280 [27,32].In 1970, Cecil combined the RQD index with the joint set (J n ) and introduced the K b factor ranging from 1 to 100 [33].The corresponding changes for RQD and J n were 5-100% and 1-5, which were proposed by Barton and Kirsten, respectively [26][27][28].The K b factor was later adopted in the N index.The J r /J a ratio represents the joint shear strength for the K d factor, where J r is the joint roughness, ranging from 0.5 to 4, and J a is the joint surface change ranging from 0.75 to 18 [27].The K d changes Water 2021, 13, 3205 6 of 42 from 0.03 to 5.33.J s is the effect of rock block shape and the path by which related joints are arranged to flow.In 1982, Kirsten proposed the rating for orthogonal fractured systems, which ranges from 0.37 to 1.5.The author also used the orthogonal case as the measuring standard because it was simple and effective, maintaining that the jointing angle was correct in most cases, even though slight errors exist.However, this notion was rejected by Boumaiza et al. in 2019.The finding of Boumaiza et al. in 2019 encouraged a proposed new rating from 0.09 to 1.38 for nonorthogonal fractured systems.In 2018, Kamali stressed that the number and connectivity of blocks increase and have considerable effects on erodibility under orthogonal joint conditions [34,35].
The following are the two other indices recently proposed by Pells to evaluate the limit of rock resistance to erosive power of water: geological strength index (GSI) (eGSI) and RMEI.The eGSI was developed on the basis of the GSI of Hoek and previously suggested to define and classify the rock mass [1,36].
Bieniawski concluded that the joint orientation parameter is not considered in the RMR system, when the GSI system is identified utilising the RMR system [37].In this regard, Pells suggested the eGSI index for illustrating the importance of the joint orientation relative to the flow direction and a rock block shape by considering the recently proposed discontinuity orientation adjustment parameter (E doa ) (Equation ( 6)) [1].
The RMEI index is the other index suggested by Pells [1].This index may be noted on the basis of the likelihood factor (LF) and relative significance factor (RF) (Equation ( 7)).Prefixes P 1 to P 5 in Equation ( 7) are various sets of parameters introduced in the classification system.P 1 to P 5 indicate the kinematically viable system for isolation, the feature of the potentially eroding surface, the feature of joints, the joint spacing, and the rock block shape, respectively [1].
Pells argued that the erosion mechanism in this field is poorly represented by existing indices, including Kirsten's index.The RMEI systems represent the rock mass parameters that gain significance from studies on eroded rocky spillways to control this mechanism.A high RF is used to measure the relative importance of the rock mass parameters.The RF measures the isolation and eroding surface by using three values: joints (RF = 2), joint spacing (RF = 1), and rock block shape (RF = 1).
The RMEI system designed by Pells for hydraulic erosion is structually similar to the Q-system of Barton that was initially designed for field investigaions.The size of the rock blocks determined by joint spacing and dip is indirectly included in RMEI.However, the joint shear strength is not considered in the RMEI system.The joint nature that maintains the natural conditions of joints can be used as an alternative for determining the rock block size.Joint spacing is a less significant factor (RF = 1) in the RMEI system compared with joint factor (RF = 2).The K b rating factor as an indicator of the rock block size is given more importance in the Q-system compared with the K d factor representing the joint shear strength [26].This phenomenon shows that the field evaluation to determine the importance of the rock mass parameters is highly influenced by the individual judgement of the analyst.
Pells assumed that the existing indices, including Kirsten's index, do not indicate the erosion mechanism in the field.Accordingly, the RMEI system attempts to represent the rock mass parameters controlling the erosion mechanism.The relative importance of these parameters is assumed to be due to the field observations of the eroded rocky spillways.A high RF weights the important rock mass parameters.The kinematically viable system for isolation and the feature of the potentially eroding surface are weighted with a high RF value of three compared with the nature of joints (RF = 2), joint spacing (RF = 1), and rock block shape (RF = 1).The structure of the RMEI system for hydraulic erosion is considered similar to that of the Q-system, which was also developed on the basis of the field investigation.In the RMEI system of Pells, the rock block size is indirectly included.However, the joint spacing and dip could provide an idea about the rock block size given that the wide spacing of joints produces a more important rock block volume than the close spacing of joints.The joint shear strength is not included in the RMEI classification.Nevertheless, the nature of the joint factor is the same given that it incorporates the natural conditions of joints.The RMEI classification considers the joint spacing factor as a less important factor (RF = 1) compared with the nature of the joint factor, which is weighted with RF = 2.The Q-system provides a more important rating to the K b factor, which is an indicator of the rock block size, compared with the K d factor, which represents the joint shear strength [26].The discordance of the relative importance of the parameters included in these two classification systems demonstrates how the field evaluation can be highly affected by analyst judgement.Consequently, the relative importance of the rock mass parameters is difficult to determine by using an accurate alternative.

Background of the Comparative Methods Base on Kirsten's Index
Van Schalkwyk assumed the erosion condition as a function of erosion depth (Table 1) [29].They also used recent data from Moore, Kirsten's index limitations, and scale effects to update their findings [38].Accordingly, the scour threshold lines were changed (Table 1).Annandale analysed the findings of Moore, Van Schalkwyk, and other scour data using Kirsten's index to offer a scour threshold related to unit stream power dissipation (Π UD ) (Table 1) [21,29,31].Table 1 illustrates that the erosion depth beyond two leads to scour occurrence, while the values below two are insignificant and might be the outcome of loose rock blocks being removed from the surface [21].In 2000, Kirsten adopted a similar approach for determining the scour threshold to show the scour and no-scour conditions by using data from the study of Dooge and Moore [18,30,31].A comparison of the scouring onset values for Kirsten's index versus Π UD is presented in Figure 4, and the corresponding equations are shown in Table 2.
Water 2021, 13, 3205 8 of 42 UD scour occurrence, while the values below two are insignificant and might be the outcome of loose rock blocks being removed from the surface [21].In 2000, Kirsten adopted a similar approach for determining the scour threshold to show the scour and no-scour conditions by using data from the study of Dooge and Moore [18,30,31].A comparison of the scouring onset values for Kirsten's index versus Π UD is presented in Figure 4, and the corresponding equations are shown in Table 2.  Pells classified the erosion condition in certain classes with minor modifications on the depth of eroded area according to the same concept of the study of Van Schalkwyk and added information on the extent of the eroded area (Table 3).Neither depth nor volume is a complete descriptor of the erosion severity.However, whether the conditions in the  Pells classified the erosion condition in certain classes with minor modifications on the depth of eroded area according to the same concept of the study of Van Schalkwyk and added information on the extent of the eroded area (Table 3).Neither depth nor volume is a complete descriptor of the erosion severity.However, whether the conditions in the first and second columns of Table 3 should be simultaneously or separately considered is ambiguous in this classification.Pells noted that neither depth nor volume is a complete descriptor.Thus, if the conditions should be simultaneously considered, then it might result in a problem, for example, whether class 1 or 2 should be considered for the maximum depth of 0.5 m and the general extent of 8 m 3 /100 m 2 .Pells plotted the calculated index values versus Π UD by using the eGSI index (Equation ( 6)) and the RMEI index (Equation ( 7)) to manually determine the selected erosion classes.These values were separated by scour threshold lines (Figure 5a,b) and were originally proposed by Pells and recently modified by Douglas [39].Pells plotted the calculated index values versus Π UD by using the eGSI index (Equation ( 6)) and the RMEI index (Equation ( 7)) to manually determine the selected erosion classes.These values were separated by scour threshold lines (Figure 5a,b) and were originally proposed by Pells and recently modified by Douglas [39].However, no finding has been obtained on the optimisation introduced into the original scour threshold lines.

Discussion of EIM and Pells' Methods
The easiness and wide pertinence to different stream conditions make the erodibility index method (EIM) especially alluring for practice.In any case, this approach is not without confinement.As the name of this approach suggests, it consolidates an empirical index to classify and identify the rock.Along these lines, the EIM cannot reveal various erosion mechanisms (i.e., fatigue failure, block removal, or brittle failure).The results of EIM are general and do not incorporate the erosion mechanisms.The geometry of rock is considered 2D, and the 3D nature of the discontinuity orientations is not accounted.However, this method can illustrate the discontinuity structure relative to the flow direction.However, no finding has been obtained on the optimisation introduced into the original scour threshold lines.

Discussion of EIM and Pells' Methods
The easiness and wide pertinence to different stream conditions make the erodibility index method (EIM) especially alluring for practice.In any case, this approach is not without confinement.As the name of this approach suggests, it consolidates an empirical index to classify and identify the rock.Along these lines, the EIM cannot reveal various erosion mechanisms (i.e., fatigue failure, block removal, or brittle failure).The results of EIM are general and do not incorporate the erosion mechanisms.The geometry of rock is considered 2D, and the 3D nature of the discontinuity orientations is not accounted.However, this method can illustrate the discontinuity structure relative to the flow direction.
The data analysed by Annandale consider several types of flow conditions (137 case studies), while Pells used the predominant channel flow in spillways.Under the channel flow conditions, block removal is considered to be the dominant mechanism for scour to occur; thus, the analysis of variables representing blocks (e.g., K b ) in those cases would be expected to show great importance.Under the jet conditions wherein fracturing of the rock can be more dominant, the other variables may be shown to be more relevant (e.g., M s ).
Rock scour is a highly complicated mechanism, and its evaluation should start by identifying the relevant rock mass parameters to assess this mechanism.For this purpose, in 2019, Boumaiza et al. developed a method on the basis of real data and previous empirical methods [35].They examined a set of rock mass parameters to determine those that are considered related parameters to evaluate the rock mass erosion.In the proposed approach, various parameters are assessed, such as E doa , NPES, J o , J s , K d , K b , V b , and UCS.Finally, K d , V b , J o , E doa , and NPES were chosen as the relevant parameters by sensitivity analysis.

Semi-Theoretical Methods for Computing Scour of the Plunge Pools
Several semi-theoretical equations can be used for computing the final scour depth in the plunge pool bottom in alluvial and weathered rock.Some semi-theoretical equations are identified in Table 4 [40][41][42][43][44][45][46][47][48][49][50].The reduced models of the Froude scale are used to produce most semi-theoretical equations based on the real data analysis and experiments.Some of these equations are proposed on the basis of the prototypes for ski-jump conditions.Several equations are also considered (Table 5).In 1985, Mason and Arumugam proposed Equation ( 8) for the intermediate outlet conditions of the free surface reservoir and dam structures on the basis of the several prototype cases and the models of numerous erodible beds [5].

Castillo and Carrillo (2016) [50]
The different parameters of Equation ( 8) are illustrated in Figure 6, where Y s represents the scouring depth below the bedrock, H n is the net energy head, g is the gravitational acceleration, Г is the experimental coefficient, D s stands for the sum of the scour depth and tailwater depth, q represents the specific flow rate, Y 0 is the tailwater depth, and d stands for the block and bed material size [50].

Analysts
Г v w x y z d Hartung (1959) [51] 1.4 0 0 0.64 0.36 0.32 d85 Chee and Padiyar (1969) [52] 2.126 0 0 0.67 0.18 0.063 dm Bisaz and Tschopp (1972) [53] 2.76 0 0 0.5 0.25 1 d90 Martins (1975) [44] 1.5 0 0 0.6 0.1 0 -Machado (1982) [54] 1.35 0 0 0.5 0.3145 0.0645 d90 The different parameters of Equation ( 8) are illustrated in Figure 6, where Ys represents the scouring depth below the bedrock, Hn is the net energy head, g is the gravitational acceleration, Г is the experimental coefficient, Ds stands for the sum of the scour depth and tailwater depth, q represents the specific flow rate, Y0 is the tailwater depth, and d stands for the block and bed material size [50].In Equation ( 8) and Figure 6, , , , , and  are experimental and empirical constants settled in Table 4 by different authors.Various coefficients of Table 4 contrasted with general equations. represents the vertical interval between the tailwater level and the reservoir, and t0 stands for the energy loss in the chute; h0 is the vertical distance between the outlet exit and the tailwater level, and Hn = H0 = HB − t0 stands for the net energy head at the outlet exit.Ui, i, and i are the jet velocity, thickness, and angle in the initial point, respectively.Uj, Bj, and j are the total jet velocity, thickness, and angle in the impact point, respectively.Table 6 shows the abbreviation list of the plunge pool semi-theoretical approaches.In Equation ( 8) and Figure 6, x, y, z, v, and w are experimental and empirical constants settled in Table 4 by different authors.Various coefficients of Table 4 contrasted with general equations.H represents the vertical interval between the tailwater level and the reservoir, and t 0 stands for the energy loss in the chute; h 0 is the vertical distance between the outlet exit and the tailwater level, and H n = H 0 = H B − t 0 stands for the net energy head at the outlet exit.U i , B i , and θ i are the jet velocity, thickness, and angle in the initial point, respectively.U j, B j , and θ j are the total jet velocity, thickness, and angle in the impact point, respectively.Table 6 shows the abbreviation list of the plunge pool semi-theoretical approaches.

Semi-Analytical Methods
Semi-analytical methods are derived by calibration of at least one parameter using experiments, or empirical calibration of coefficients using data.A jet plunging into a pool leads to the development of average and fluctuating dynamic pressures that can cause the breakup and removal of rock.When a jet plunges into a pool, it creates dynamic pressures below the point of jet impact.These pressures consist of average and fluctuating dynamic pressures.If free air is introduced into a plunge pool by the jet, then the average dynamic pressure decreases with the increase in the air concentration [3,[55][56][57].In 2019, Maleki and Fiorotto analysed the block stability and potential rock scour in the plunge pool bottom on the basis of new experimental study by focusing on the mean dynamic pressure and fluctuating pressure coefficients produced by rectangular water jets [58].In the fluvial erosion case, the comprehensive scour model (CSM) of Bollaert is one of the adopted semi-analytical methods [59].This model was proposed in 2002 for plunge pools with circular jets and later in 2010 developed for riverbed by Bollaert.Plunging jets occur in different types of hydraulic structures, including the characteristic location of the spillways, overtopping dams, valves, and dam structural gates [59].In such a manner, the CSM could be used for plunging jets.The CSM of Bollaert comprises three unique techniques, for example, comprehensive fracture mechanics (CFM), quasisteady impulsion (QSI), and dynamic impulsion (DI) (Figure 7).In the main strategy, the CSM explains the fracture propagation by utilising linear elastic fracture mechanics.In the subsequent technique, the second law of Newton (impulsion) and net uplift pressure are used to determine the single block uplift.Bollaert proposed the QSI technique to assess the eroded area in walls of plunge pools and determine the scour depth.The QSI technique is developed for protruding rock blocks on the basis of the different flow conditions and relative flow direction.The hydraulic parameters are determined along the critical area in the pool bottom and rock mass surface for every fracture mechanism.Figure 7 illustrates the CSM.
Figure 7 illustrates the main physical and mechanical processes used for determination of rock scouring (CSM).CFM comprises three modules: falling jet module, plunge pool module, and rock mass module.The latter executes the fracture mechanism.
The falling jet module (Figure 8) clarifies how the geometrical and hydraulic parameters of the plunging jets are changed from outlet exit to impact point of plunge pool surface.Issuance jet velocity (V i ), initial turbulence intensity (T u ), and issuance jet diameter (D i ) identify the jet characteristics at the outlet exit.T u represents the extent of the root-meansquare (RMS) estimation of the fluctuating velocity of the plunging jet to the mean axial velocity of the plunging jet.The plunging jet pathway relies upon air drag and ballistics.The falling jet module can determine the jet pathway length L from issuance to the impact point, the longitudinal distance of the plunging jet from issuance to the impact point, the jet velocity (V j ), and the jet diameter (D j ) at the impact point (Equations ( 9) and ( 10)) [60,61].The falling jet module (Figure 8) clarifies how the geometrical and hydraulic parameters of the plunging jets are changed from outlet exit to impact point of plunge pool surface.Issuance jet velocity (Vi), initial turbulence intensity (Tu), and issuance jet diameter (Di) identify the jet characteristics at the outlet exit.Tu represents the extent of the root-mean-square (RMS) estimation of the fluctuating velocity of the plunging jet to the mean axial velocity of the plunging jet.The plunging jet pathway relies upon air drag and ballistics.The falling jet module can determine the jet pathway length L from issuance to the impact point, the longitudinal distance of the plunging jet from issuance to the impact point, the jet velocity (Vj), and the jet diameter (Dj) at the impact point (Equations ( 9) and ( 10)) [60,61].The plunge pool module depicts the plunging jet whilst crossing the plunge pool.This module characterises and computes the pressures of water at the plunge pool bottom on the rock surface.The Y/Dj ratio plays an important role in defining coefficients, which is straightforwardly identified with jet diffusion.Y is fundamental and represents the plunge pool water depth.The fluctuating dynamic pressures C′pa and the mean dynamic pressure coefficient Cpa, which are determined by using Equations ( 11)-( 13), are The plunge pool module depicts the plunging jet whilst crossing the plunge pool.This module characterises and computes the pressures of water at the plunge pool bottom on the rock surface.The Y/D j ratio plays an important role in defining coefficients, which is straightforwardly identified with jet diffusion.Y is fundamental and represents the plunge pool water depth.The fluctuating dynamic pressures C pa and the mean dynamic pressure coefficient C pa , which are determined by using Equations ( 11)-( 13), are the important pressures straightforwardly estimated under the centreline of the jet impact point.These pressure coefficients were defined on the basis of the experimental data analysis of [55].
where α i is the air concentration in % and β is the volumetric air-to-water ratio.Table 7 illustrates T u for various outlet structural types, and the polynomial coefficients for various turbulence intensities.The input and output parameters of the plunge falling jet and pool modules are noted in Table 8. Figure 9a,b demonstrate the C pa and C pa as functions of the pool water depth to jet diameter (Y/D j ) on the basis of the curve fitting of Ervine results [55].
Table 7.The range of T u due to the type of outlet structure (left side) [61]; polynomial and regression coefficients for various T u (right side) [62].

CFM
The CFM method recognises the final scour depth in the pool bottom around the area of the turbulent jet impingement for the close-ended joints.This method presents two different approaches for brittle failure and fatigue failure mechanism.
For this purpose, the maximum dynamic pressure in the close-ended fissures P max is determined by using Equation (15) [59]: where C max p represents the dynamic pressure coefficient and is provided from multiplication of the amplification factor Г + with C pa , which is determined by using Table 9 and by superposition by C pa .Table 9. Conditions for computing the amplification factor (Г + ) [63].

Curve of Maximum Values
Curve of Minimum Values In 2016, Bollaert modified and simplified the P max and defined the new equation for P max (MPa) (Equation ( 16)) [64].
where ρ is the density of water with a dimension of kg/m 3 .
As noted in the past area, the minimum and maximum pressures of the cycles are used to determine ∆Pc.The f c follows the supposition of an ideal resonator framework and relies upon the air concentration (Equation ( 14)) in the joint α i and on the joint length L f that appeared in Figure 10 [59]. Figure 10a shows the L f in the closed-end rock joints, and Figure 10b  Crack propagation represents the stresses occurring at the joint tip because of the fatigue mechanism owing to the cyclic nature of pressures inside the rock joints.This concept is depicted by linear elastic fracture mechanics.The propagation of cracks can be fragile or time subordinate diffusion.The fragile or brittle failure occurs when an applied stress intensity (KI) surpasses the fracture toughness (KIc).Fatigue failure occurs in the contrary case.Stresses are described by KI (Equation ( 17)) [65]: Crack propagation represents the stresses occurring at the joint tip because of the fatigue mechanism owing to the cyclic nature of pressures inside the rock joints.This concept is depicted by linear elastic fracture mechanics.The propagation of cracks can be fragile or time subordinate diffusion.The fragile or brittle failure occurs when an applied stress intensity (K I ) surpasses the fracture toughness (K Ic ).Fatigue failure occurs in the contrary case.Stresses are described by K I (Equation ( 17)) [65]: where K I is the stress intensity (MPa/m 1/2 ), which is produced by dynamic pressures in the pool bottom and applied to the crack tip that may cause joint propagation; σ water is water stress (MPa) and determined as 80% of P max because of the varying values of P max inside the fracture; f represents the boundary correction factor and relies on the joint type and persistency.Three relevant joint geometries are presented in Figure 10.These geometries are a semi-elliptical crack (EL), single edge crack (SE), and centre cracked (CC).In each of these configurations, the boundary correction factors (f ) are defined by Equations ( 18)- (20).
where C is coefficient in Equation ( 18) graphically determined on the basis of the a B , a c , c W ratios and Ø. Bollaert (2016) simplified Equation ( 17) and proposed Equation ( 21) for determining the stress intensity [64].
As previously mentioned, stress intensity (K I ) is produced by dynamic pressures in the pool bottom and applied to the joint tip that may cause joint propagation.If the stress intensity (K I ) is higher than the fracture toughness K Ic , then crack propagation may occur.
The K Ic relies on the mineralogical type of rock, tensile strength (T), and UCS.The changes are conducted to illustrate the rate of loading and in situ stress.Along these lines, the in situ fracture toughness K I,ins is developed on the basis of data analysis and constrained by using Equations ( 22) and ( 23) and Table 10 [3].
where σ c indicates the confinement horizontal in situ stress and is determined by using vertical stress and K 0 (σ h /σ v ).UCS, σ c , and T are in MPa.If the applied stress intensity surpasses the in situ fracture toughness (K I ≥ K I,ins ), then brittle failure will occur; otherwise, failure mechanics and crack propagation are communicated as Equation ( 24) [66]: where C and m rely on the rock characteristics that are recognised on the basis of fatigue tests; dN represents the number of pressure cycles; and ∆K I defines the most and the least stress intensity differences.Variables m and C must be determined to handle timedependent crack propagation into the model [67].
Water 2021, 13, 3205 18 of 42  22) and ( 23) for the in situ fracture toughness value K I,ins [59].24) and proposed Equation ( 25) as the new equation for the fatigue process.Fracture toughness (K Ic ) was considered in the new equation [68]:

Type of
The flowchart of CFM is presented in Figure 11.The ultimate scour depth and time of scouring could be computed on the basis of the maximum dynamic pressure in the pool bottom (P max ), stress intensity (K I ), in situ fracture toughness (K I,ins ), and fatigue failure mechanism owing to this method.

DI Method
Bollaert proposed the DI method for the first time, which is used for predicting the ultimate scour depth in the plunge pools for open-ended joints and single rock blocks [59].This method explains the movement of rock blocks by their mass on the basis of the uplift pressures.DI is based on the maximum impulse (Imax) in an open-ended rock mass fracture (single rock block) that is provided by the time integral of net forces.Net forces consist of the submerged weight of the block, shear and interlocking forces, and forces

DI Method
Bollaert proposed the DI method for the first time, which is used for predicting the ultimate scour depth in the plunge pools for open-ended joints and single rock blocks [59].This method explains the movement of rock blocks by their mass on the basis of the uplift pressures.DI is based on the maximum impulse (I max ) in an open-ended rock mass fracture (single rock block) that is provided by the time integral of net forces.Net forces consist of the submerged weight of the block, shear and interlocking forces, and forces produced by pressures under and over the block [59,62].The final purpose in the DI method is to determine the uplift height (h up ) by using net uplift velocity (V ∆tpulse ), which is determined by net impulsion (I net ) (kN/s), and the block mass (m s ) (kg).The DI method is explained stepwise in the following paragraphs.The flowchart of DI is presented in Figure 12; Figure 13 shows the impulse dynamics on a block.
where  ∆ is the uplifting velocity (m/s), which is also named 'launch velocity' computed from Equation (36) and net impulsion (Inet) shown in Equation (28).
where Fu and Fo are considered the forces produced by pressures under and over the block (kN), and are produced because of the dynamic pressure in the pool bottom;  is the submerged weight of block (kN), which is determined by using Equation (29); and  is the shear force (kN).Shear force ( ) is considered zero because it acts in a direction parallel to a surface.This variable is zero for a cubic block.Time integral () is in seconds, taken over the ∆ period (∆ is the period in which positive differences exist in forces produced by pressures over and under the block), and is determined by Equation (34).The block mass  is in kilograms and determined by Equation (35).
where  is the block volume (m 3 ); ϒ is the unit weight of rock (kN/m 3 ); ϒ is the unit weight of water (kN/m 3 ); and Xb, Yb, and Zb are the block width, length, and height, respectively.
where V ∆tpulse is the uplifting velocity (m/s), which is also named 'launch velocity' computed from Equation ( 36) and net impulsion (I net ) shown in Equation (28).
where F u and F o are considered the forces produced by pressures under and over the block (kN), and are produced because of the dynamic pressure in the pool bottom; G b is the submerged weight of block (kN), which is determined by using Equation (29); and F sh is the shear force (kN).Shear force (F sh ) is considered zero because it acts in a direction parallel to a surface.This variable is zero for a cubic block.Time integral (dt) is in seconds, taken over the ∆t period (∆t is the period in which positive differences exist in forces produced by pressures over and under the block), and is determined by Equation (34).The block mass m s is in kilograms and determined by Equation (35).
where V b is the block volume (m 3 ); Υ r is the unit weight of rock (kN/m 3 ); Υ w is the unit weight of water (kN/m 3 ); and X b , Y b , and Z b are the block width, length, and height, respectively.The shear force F sh was dismissed in these relations.The nature of the rock joints is rough in real situations.The roughness of the rock joints plays an important role in determining the shear strength.The shear strength of rock must be determined when designing rock-engineering structures.In this regard, the shear resistance must be determined in the block uplift process.The block uplift velocity relies upon the lateral pressure differences and frictional resistance of the sidelong discontinuities.For this reason, in 2014, Pan et al. developed several analytical equations for computing shear resistance force (F sh ) for two unique modes [69].The first mode (Equation ( 30)) is used when the pressure fluctuations are equal in the opposite surface (top and bottom) of the block.The second mode (Equation ( 31)) is used when the pressure fluctuations at the top of a rock block are more prominent than that at the bottom of a rock block.
where µ stands for the frictional coefficient, Υ is the submerged unit weight (kg/m 3 ), and ∆h is the pressure head difference.Bollaert determined the difference between the forces under and over the block by using a nondimensional (C I ) coefficient (Equation ( 33)), which is produced by curve fitting on the basis of some experiments [59].Variable C I is the coefficient for determining the difference between pressures under and over the block.Thus, Equation ( 33) should be used for computing net uplift force (F u − F o ).
where Y is the water depth in the plunge pool and D j is the jet diameter at the impact.Both variables are in meters. where 2g is the incoming kinetic energy (m), V j is the jet speed at impact (m/s), and A is the area of the block surface over or under the block (m).
where T up is a nondimensional time coefficient and assumed by large-scale laboratory experiments, L f is the total length of the joint (m), and c is the mean wave celerity (m/s).
If X b , Y b , and Z b are considered as block width, length, and height, respectively, then the total length of the joint is where V b stands for the block volume (m 3 ), and ρ s is the density of rock (kg/m 3 ).Considering Equations ( 26)-( 29) and ( 32)-( 35), uplifting velocity (V ∆tpulse ) is computed similar to Equation (36): Finally, the uplift height is directly computed by using Equation (37).
After h up is determined, the ratio of uplift height to block height (h up /Z b ) must be characterised.The capability of rock scour must be evaluated by utilising Figure 14, which depends on field data analysis.In 2005, Bollaert and Schleiss determined that the block might be separated from the rock mass when h up /Z b is greater than 20% (h up /Z b > 0.20) [61].
After hup is determined, the ratio of uplift height to block height (hup/Zb) must be characterised.The capability of rock scour must be evaluated by utilising Figure 14, which depends on field data analysis.In 2005, Bollaert and Schleiss determined that the block might be separated from the rock mass when hup/Zb is greater than 20% (hup/Zb > 0.20) [61].The Bollaert portrayal for the DI method has two ambiguities, especially in enlisting (Fu − Fo).A comparable delineation was determined for DI and portraying CI coefficient in several articles of Bollaert [59,[61][62][63].
The CI, (Fu − Fo), and uplift height (hup) are difficult to determine considering this description.The final dimension of (Fu − Fo) is in m instead of kN by considering this description.The unit weight of water (ϒ ) and area of block surface (A) are not intended for converting the DI coefficient (CI) to force.
In 2004, Annandale corrected Bollaert's problem in converting CI to force and proposed a modified equation for determining (hup) [71].However, Annandale did not consider the nondimensional time coefficient ( ) in Equation (37).
Bollaert published a new article and proposed the simplified CSM (SCSM) because of the difficulty in understanding the CSM [64].In SCSM, the DI method is explained stepwise.Bollaert corrected the problem in previous publications.This time, he converted CI to pressure instead of force and did not consider the block surface area (A).We present Equation ( 38) for computing the hup by considering its explanation about the DI method and block surface area (A).In 2004, Annandale corrected Bollaert's problem in converting C I to force and proposed a modified equation for determining (h up ) [71].However, Annandale did not consider the nondimensional time coefficient (T up ) in Equation (37).
Bollaert published a new article and proposed the simplified CSM (SCSM) because of the difficulty in understanding the CSM [64].In SCSM, the DI method is explained stepwise.Bollaert corrected the problem in previous publications.This time, he converted C I to pressure instead of force and did not consider the block surface area (A).We present Equation (38) for computing the h up by considering its explanation about the DI method and block surface area (A).
where the block surface area (A) is in m 2 .
In 2006, Manso developed the CSM of Bollaert and studied the effects of the plunge pool geometry on the rock scour with a high velocity jet.Manso performed different experiments by considering the various plunge pool geometries (12 geometries) [72].This work was conducted to determine the effect of the plunge pool geometry on pressure produced by circular jets.In 2014 and 2016, Duarte developed the CSM focused on the jet aeration and considered various types of fluid on the basis of air content (β) to propose a new equation for determining block displacement (uplift height h up ) [73,74].
In 2009, Asadollahi proposed a semi-analytical methodology for assessing the single block stability in the plunge pools by using three-dimensional block stability (BS3D) [75].Asadollahi developed this methodology by numerical BS3D code.BS3D is used to compute the block uplift acceleration by using Federspiels' experimental studies [76][77][78].In this methodology, the uplift height can be determined on the basis of Newton's second law after the aforementioned acceleration is computed.Asadollahi assessed the erosion in the plunge pools on the basis of Bollaert DI technique.Moreover, Asadollahi modified Bollaert's criteria for deciding the occurrence of block removal and proposed the h up /Zb > 0.25 ratio based on data analysis of several case studies and using Martins' experimental data [79].For this reason, Asadollahi actualised a Fortran code dependent on the calculation presented by Tonon (2007) [80].BS3D considers all expansive displacement techniques for rock blocks subject to conventional powers.This methodology had admissible results compared with the observed block uplift.

QSI Method
In 2010, Bollaert proposed the QSI technique to assess the scour in walls of the plunge pools for parallel flow condition relative to the rock surface [81].Bollaert developed this methodology on the basis of Reinius' study, which was for identification of the effect of the pressure fluctuations on the rectangular rock block uplift and protruding block cases [82].The effective forces that are applied to the block and cause the rock block uplift or ejection follow: rock block submerged weight, quasi-steady uplift force produced because of the pressures underneath the rock block, and the turbulent uplift forces produced because of the pressure fluctuations (Figure 15b) [81].The uplift forces produced because of the pressure fluctuations are significant to cause the block uplift in protruding rock block cases.In 2015, George studied various flow condition analyses and concluded that the quasisteady system is more effective than the other lifting forces when considering the parallel flow condition.George also found that the forces produced by pressure fluctuations could be neglected from calculations.The flowchart of QSI is presented in Figure 16 [83].
Figure 15a shows the various areas that developed as a result of the falling jet in the plunge pools.Figure 15b illustrates the effective lifting forces that cause the rock block uplift, where e block represents the block protrusion, h block is the block height, L block denotes the block length, and β block is the angle between the block edges.In the QSI, only quasi-steady lift force is considered in the calculations.The quasisteady lift force is determined by using Equation ( 39) [81]: where C uplift denotes the net uplift pressure coefficients and computed by deducting the surface pressure coefficient (C surf ) from the pressure coefficient inside the joints (C joint ) (C joint − C surf ).Table 11 shows the potential net static uplift pressures for various arrangements of rock blocks [82].Figure 15a shows the various areas that developed as a result of the falling jet in the plunge pools.Figure 15b illustrates the effective lifting forces that cause the rock block uplift, where eblock represents the block protrusion, hblock is the block height, Lblock denotes the block length, and βblock is the angle between the block edges.
In the QSI, only quasi-steady lift force is considered in the calculations.The quasi-steady lift force is determined by using Equation (39) [81]: where Cuplift denotes the net uplift pressure coefficients and computed by deducting the surface pressure coefficient (Csurf) from the pressure coefficient inside the joints (Cjoint) (Cjoint − Csurf).Table 11 shows the potential net static uplift pressures for various arrangements of rock blocks [82].As referenced, L block is the block length (m) that appeared in Figure 15b, and V X,max represents the diffusive jet velocity in different locations.The different distances form wall jet issuance points and are computed by using Equation ( 40) [84].The jet velocity profile and jet thickness have an inverse relationship where the wall jet velocity decreased, and the wall jet velocity profile smoothed when the jet thickness is increased (Figure 15a).
V X,max = 3.5 V Z bottom X h downstream (40) where h downstream is the initial thickness of flow (Figure 15a) at the point of deflection and determined by using Equation (41), where q downstream is the discharge rate in the downstream side and q total is the total discharge rate of the falling jet.
The V Zbottom is the initial wall jet velocity (m/s).V Zbottom relies upon the angles that the jet makes with the water surface at the impacting point and water depth Z, and computed by utilising the model of Hartung and Häusler, which is presented in Equation ( 42) [85].
where Z core (m) represents the vertical distance that the core jet needs to be deflected inside plunge pools and commonly considered as 4-5 times the D j .The V i represents the initial plunging jet velocity with a dimension of m/s, and Z is the water depth.As referenced, Lblock is the block length (m) that appeared in Figure 15b, and VX,max   As referenced, Lblock is the block length (m) that appeared in Figure 15b, and VX,max represents the diffusive jet velocity in different locations.The different distances form wall jet issuance points and are computed by using Equation ( 40) [84].The jet velocity   Block shape is one of the important parameters in rock block removal.For this purpose, the common blocks shape was identified on the basis of the existing rock mass geomechanical data and field observations.In that capacity, a plate-like formed block will be ejected just by jet impact.Finally, the ultimate scour depth can be determined by utilising Bollaert's DI technique after the lifting forces are determined.
On the basis of Lesleighter's study on Paradise Dam in Australia, the DI method overestimates the scour depth, where the prediction of CFM is close to the actual case [87].
The basis of the semi-analytical methods is experimental.For this purpose, Table 13 illustrates some information about the setup of these experiments.Quantitative:

Speed at displacement
Numerous parameters affect the plunge pool scour.These parameters can be studied in two different aspects: hydraulic and rock mass.
The most critical parameter in the hydraulic aspect is uplift pressure, which depends on many other geometrical and hydraulic parameters; for example, tailwater depth, initial jet velocity, jet diameter, and plunging jet height affect pool bottom and total uplift pressures.Increasing the tailwater depth will decrease the uplift pressure; conversely, the uplift pressure will increase by raising the initial jet velocity, jet diameter, and plunging jet height.
Numerous geometrical and geomechanical parameters in the rock mass aspect affect the ultimate scour depth.These parameters include joint opening, block volume, block shape, joint roughness, and fracture toughness.This selection of parameters can be comprehensively studied as future challenges.

Scouring at the Multiblock System with Numerical Methods
Most previous methods used for computing the scour depth or to assess the hydraulic erodibility were developed for the analysis of a single rock block related to the geomechanical, geometrical, and hydraulic parameters.In most plunge pool cases, the stability of a single rock block is assessed as a representative block of a characteristic blocky layer in various hydraulic conditions.This assessment is continued for nether layers because the representative block started to be stable for removal.Limited methods have been proposed to assess the scouring for multiblock systems, which are presented in Section 4. Figure 17 presents the numerical approach.
Water 2021, 13, x FOR PEER REVIEW 29 of 44 Numerous parameters affect the plunge pool scour.These parameters can be studied in two different aspects: hydraulic and rock mass.
The most critical parameter in the hydraulic aspect is uplift pressure, which depends on many other geometrical and hydraulic parameters; for example, tailwater depth, initial jet velocity, jet diameter, and plunging jet height affect pool bottom and total uplift pressures.Increasing the tailwater depth will decrease the uplift pressure; conversely, the uplift pressure will increase by raising the initial jet velocity, jet diameter, and plunging jet height.
Numerous geometrical and geomechanical parameters in the rock mass aspect affect the ultimate scour depth.These parameters include joint opening, block volume, block shape, joint roughness, and fracture toughness.This selection of parameters can be comprehensively studied as future challenges.

Scouring at the Multiblock System with Numerical Methods
Most previous methods used for computing the scour depth or to assess the hydraulic erodibility were developed for the analysis of a single rock block related to the geomechanical, geometrical, and hydraulic parameters.In most plunge pool cases, the stability of a single rock block is assessed as a representative block of a characteristic blocky layer in various hydraulic conditions.This assessment is continued for nether layers because the representative block started to be stable for removal.Limited methods have been proposed to assess the scouring for multiblock systems, which are presented in Section 4. Figure 17 presents the numerical approach.Multi-block assessments are more important than single-block assessments because of the location dependent assessment of this system.This system can assess the scouring for various rock block shape and geometrical conditions.In 2009, Wibowo evaluated rock erosion by modelling in universal distinct element code (UDEC) software for dam spillways, based on the block theory of Goodman [88,89].Figure 18 shows the UDEC model.The red blocks in this model are more vulnerable to removal than those shown in blue.
Water 2021, 13, x FOR PEER REVIEW 30 of 44 Multi-block assessments are more important than single-block assessments because of the location dependent assessment of this system.This system can assess the scouring for various rock block shape and geometrical conditions.In 2009, Wibowo evaluated rock erosion by modelling in universal distinct element code (UDEC) software for dam spillways, based on the block theory of Goodman [88,89].Figure 18 shows the UDEC model.The red blocks in this model are more vulnerable to removal than those shown in blue.
In 2016, Castillo and Carrillo proposed a numerical simulation by using computational fluid dynamics (CFD) in FLOW-3D software to help in determining the bed erosion in plunge pools [50].Their methodology solves the Navier-Stokes equations and combines various turbulence models, a block transport model, a semi-theoretical model of the rock mass erosion, and a methodology for computing the free surface of the fluid.This methodology evaluates the hydraulic erodibility of the bed rock by using the combination of the aforementioned models by accurate characterization.In the final step of their study, they conducted a comparison between the pressures and the mean dynamic pressure coefficients obtained exactly at the stagnation area with the parametric approach and studied the local effects of the block's movement.
Li and Liu, in 2010, proposed another numerical method on the basis of the discrete fracture network (DFN) for evaluating scouring in the plunge pools.The DFN is developed based on a numerical analysis of the hydraulic erodibility in the downstream of dams and the Monte Carlo technique.Ejectable 2D rock blocks were determined by considering the geometry of the plunge pools and discontinuity behaviour.Li and Liu simulated the fluctuating pressures inside the joint and the distribution pressure fluctuating inside the fractured media.Li also analysed the rock block stability on the basis of the empirical equations as a function of rock joint hydraulic pressures [90].
In 2011, Dasgupta, proposed a methodology for evaluating the plunge pool scouring on the basis of three-dimensional computational fluid dynamic (CFD) programming for computing the representative erosive agent of water that is applied to the rock surface [91].Dasgupta used 2D UDEC to simulate the rock mass.The results of ANSYS FLUENT (plunge pool bottom pressures) were applied to the UDEC to assess the rock block ejection and brittle failure.Figure 19a highlights the 3D flow pattern of the plunging jet of a single valve.Figure 19b shows the primary results of the surface erosion.Figure 19b(1), Figure 19b(2), and Figure 19b(3) show the results of the block removal simulation, brittle failure of the rock block simulation, and combination of the two latter simulations, respectively.The fluctuation of the transient pressure is applied to the rock surface as the Figure 18.Two-dimensional ejectable blocks in an unlined spillway [88].
In 2016, Castillo and Carrillo proposed a numerical simulation by using computational fluid dynamics (CFD) in FLOW-3D software to help in determining the bed erosion in plunge pools [50].Their methodology solves the Navier-Stokes equations and combines various turbulence models, a block transport model, a semi-theoretical model of the rock mass erosion, and a methodology for computing the free surface of the fluid.This methodology evaluates the hydraulic erodibility of the bed rock by using the combination of the aforementioned models by accurate characterization.In the final step of their study, they conducted a comparison between the pressures and the mean dynamic pressure coefficients obtained exactly at the stagnation area with the parametric approach and studied the local effects of the block's movement.
Li and Liu, in 2010, proposed another numerical method on the basis of the discrete fracture network (DFN) for evaluating scouring in the plunge pools.The DFN is developed based on a numerical analysis of the hydraulic erodibility in the downstream of dams and the Monte Carlo technique.Ejectable 2D rock blocks were determined by considering the geometry of the plunge pools and discontinuity behaviour.Li and Liu simulated the fluctuating pressures inside the joint and the distribution pressure fluctuating inside the fractured media.Li also analysed the rock block stability on the basis of the empirical equations as a function of rock joint hydraulic pressures [90].
In 2011, Dasgupta, proposed a methodology for evaluating the plunge pool scouring on the basis of three-dimensional computational fluid dynamic (CFD) programming for computing the representative erosive agent of water that is applied to the rock surface [91].Dasgupta used 2D UDEC to simulate the rock mass.The results of ANSYS FLUENT (plunge pool bottom pressures) were applied to the UDEC to assess the rock block ejection and brittle failure.Figure 19a highlights the 3D flow pattern of the plunging jet of a single valve.Figure 19b shows the primary results of the surface erosion.Figure 19b(1), Figure 19b(2), and Figure 19b(3) show the results of the block removal simulation, brittle failure of the rock block simulation, and combination of the two latter simulations, respectively.The fluctuation of the transient pressure is applied to the rock surface as the dynamic pressure and the time-dependent pressures inside cracks and fractures.The rock block removal at the water-rock interface explains the scouring mechanism.The methodology of Dasgupta demonstrates promise to utilise numerical techniques for combining the three-dimensional rock mass geometry alongside complicated stream cases.
George presented a non-numerical (experimental-analytical) methodology by using block theory to evaluate the stability of the 3D rock blocks and potential hydraulic erodibility in dam spillways [83].The methodology of George and Sitar proposed to evaluate the hydraulic erodibility on the basis of multiblock systems [92,93].Scaled physical model testing and a prototype field experiment (the first of its kind) were performed to examine the role of the 3D geologic structure on the block erodibility.They assessed the possibility of ultimate block failure and the most probable failure mechanism by using Monte Carlo simulations and a first-order reliability method [93].This study was the first to consider noncubic/nonrectangular block geometries and showed how the geometry influenced a block response.The block theory method was shown to reasonably predict block stability under channel flow conditions.Figure 20 shows the spillway geometry and situation of the removable block inside a rock mass.The methodology of Dasgupta demonstrates promise to utilise numerical techniques for combining the three-dimensional rock mass geometry alongside complicated stream cases.
George presented a non-numerical (experimental-analytical) methodology by using block theory to evaluate the stability of the 3D rock blocks and potential hydraulic erodibility in dam spillways [83].The methodology of George and Sitar proposed to evaluate the hydraulic erodibility on the basis of multiblock systems [92,93].Scaled physical model testing and a prototype field experiment (the first of its kind) were performed to examine the role of the 3D geologic structure on the block erodibility.They assessed the possibility of ultimate block failure and the most probable failure mechanism by using Monte Carlo simulations and a first-order reliability method [93].This study was the first to consider noncubic/nonrectangular block geometries and showed how the geometry influenced a block response.The block theory method was shown to reasonably predict block stability under channel flow conditions.Figure 20 shows the spillway geometry and situation of the removable block inside a rock mass.
In 2018, Gardner developed the experimental-analytical methodology of George and proposed a new numerical methodology to evaluate rock scour in the jointed rock mass [83,94].Gardner performed the coupled fluid-solid simulations to assess the waterrock interaction behaviour and potential of the 3D polyhedral rock block erosion.This methodology was developed on the basis of the coupled discrete element method (DEM)lattice Boltzmann method (LBM) programme.In the first step, Gardner developed a new open-source DEM code, which is written using C++ language, to evaluate the mechanical behaviour and the kinematic response of the jointed rock mass.Then, Gardner developed a DEM code and coupled it with LBM to analyse the ejection and movement of 3D polyhedral rock blocks in water.In this methodology, various block shapes and geometries were analysed, and a comparison between its coupled solid-CFD approach with real dataset was performed.In 2018, Gardner developed the experimental-analytical methodology of George and proposed a new numerical methodology to evaluate rock scour in the jointed rock mass [83,94].Gardner performed the coupled fluid-solid simulations to assess the water-rock interaction behaviour and potential of the 3D polyhedral rock block erosion.This methodology was developed on the basis of the coupled discrete element method (DEM)-lattice Boltzmann method (LBM) programme.In the first step, Gardner developed a new open-source DEM code, which is written using C++ language, to evaluate the mechanical behaviour and the kinematic response of the jointed rock mass.Then, Gardner developed a DEM code and coupled it with LBM to analyse the ejection and movement of 3D polyhedral rock blocks in water.In this methodology, various block shapes and geometries were analysed, and a comparison between its coupled solid-CFD approach with real dataset was performed.

Critical Analysis and Application of Different Approaches on Real Cases
To compare various approaches of hydraulic erodibility assessment, we applied the existing methods to real data and presented their results in the present section.It is not possible to compare semi-theoretical approaches with semi-analytical and numerical methods, because of the different nature of the methods.
More detailed data and in situ tests are required to analyse and compare the semi-analytical methods, such as the CSM method and the existing dataset, making it difficult to compare this approach with other approaches.
There are no clear methodology and or software for the application of numerical methods for rock erosion prediction.The flow parameter should be determined by hydraulic software and implemented to geomechanical software.The principal question is, which flow parameter should be deduced from hydraulic software and how should it be implemented in the geomechanical software because of their existing limitations.The application and comparison have been performed on the methods used most often, such as Kirsten, Annandale, Van Schalkwyk, RMEI, and eGSI.For comparison purposes, the field data collected from various case studies conducted by Pells (2016) on unlined rocky spillways are used in this comparative critical analysis.Table A1 illustrates the result of the application of various semi-theoretical approaches on several case studies.These results are obtained by using the threshold lines proposed by these five approaches.
The efficiency of the five comparative scour thresholds is determined herein according to the number of case studies with poorly evaluated scour conditions.According to the results shown in Table A1 and Figure 21a, it can be seen that the minimum

Critical Analysis and Application of Different Approaches on Real Cases
To compare various approaches of hydraulic erodibility assessment, we applied the existing methods to real data and presented their results in the present section.It is not possible to compare semi-theoretical approaches with semi-analytical and numerical methods, because of the different nature of the methods.
More detailed data and in situ tests are required to analyse and compare the semianalytical methods, such as the CSM method and the existing dataset, making it difficult to compare this approach with other approaches.
There are no clear methodology and or software for the application of numerical methods for rock erosion prediction.The flow parameter should be determined by hydraulic software and implemented to geomechanical software.The principal question is, which flow parameter should be deduced from hydraulic software and how should it be implemented in the geomechanical software because of their existing limitations.The application and comparison have been performed on the methods used most often, such as Kirsten, Annandale, Van Schalkwyk, RMEI, and eGSI.For comparison purposes, the field data collected from various case studies conducted by Pells (2016) on unlined rocky spillways are used in this comparative critical analysis.Table A1 illustrates the result of the application of various semi-theoretical approaches on several case studies.These results are obtained by using the threshold lines proposed by these five approaches.
The efficiency of the five comparative scour thresholds is determined herein according to the number of case studies with poorly evaluated scour conditions.According to the results shown in Table A1 and Figure 21a, it can be seen that the minimum committed error is related to RMEI with 41%, and the maximum error is associated with Kirsten and Annandale methods with 81% error.This error is 81, 78, 81, 51 and 41%, respectively, for the Annandale, Van Schalkwyk, Kirsten, eGSI, and RMEI approaches.
Extensive' classification, the performance of the eGSI has less error, but for the other approaches, it is extremely high.
According to Figure 21c, it can be deduced that all the compared methods, except the RMEI method, underestimate the degree of rock mass erosion.The underestimation percentage for Kirsten, Van Schalkwyk, Annandale, RMEI, and eGSI approaches are, respectively, 100, 96.5, 90, 73 and 33%.Underestimation can have hazardous consequences for hydraulic structures, and overestimation can increase costs.

Summary and Discussion
Since 1930, several methods have been proposed to predict the hydraulic erodibility of materials.Several theoretical, semi-theoretical, numerical, and semi-analytical approaches have been developed for assessing the hydraulic erodibility of rocks.Most of these methods are developed for soils, granular materials, and diving jet cases.Some approaches have been proposed for computing and assessing the hydraulic erodibility of rocks in unlined spillways (open channels).Studies on erodibility in spillways and scour in plunge pools have two important aspects, namely, geomechanical and hydraulic aspects.Most semi-theoretical studies use the erodibility index (i.e., Kirsten index), except for Pells' method, in which RMEI and eGSI indices are used, to assess the geomechanical aspect of erodibility.
In semi-analytical methods, specifically, in the CSM of Bollaert, several factors, such as dynamic pressure, fatigue, fracture toughness, and stress intensity in the CFM method According to Figure 21b, it can be seen that the committed error for 'Negligible to Moderate' for all methods except eGSI is almost below 50%, and this shows that the performance of these methods for lesser erosions has less error.Conversely, for the 'Large to Extensive' classification, the performance of the eGSI has less error, but for the other approaches, it is extremely high.
According to Figure 21c, it can be deduced that all the compared methods, except the RMEI method, underestimate the degree of rock mass erosion.The underestimation percentage for Kirsten, Van Schalkwyk, Annandale, RMEI, and eGSI approaches are, respectively, 100, 96.5, 90, 73 and 33%.Underestimation can have hazardous consequences for hydraulic structures, and overestimation can increase costs.

Summary and Discussion
Since 1930, several methods have been proposed to predict the hydraulic erodibility of materials.Several theoretical, semi-theoretical, numerical, and semi-analytical approaches have been developed for assessing the hydraulic erodibility of rocks.Most of these methods are developed for soils, granular materials, and diving jet cases.Some approaches have been proposed for computing and assessing the hydraulic erodibility of rocks in unlined spillways (open channels).Studies on erodibility in spillways and scour in plunge pools have two important aspects, namely, geomechanical and hydraulic aspects.Most semitheoretical studies use the erodibility index (i.e., Kirsten index), except for Pells' method, in which RMEI and eGSI indices are used, to assess the geomechanical aspect of erodibility.
In semi-analytical methods, specifically, in the CSM of Bollaert, several factors, such as dynamic pressure, fatigue, fracture toughness, and stress intensity in the CFM method (for close-ended joints) and net impulsion in the DI method (for open-ended joints), are employed to compute the ultimate scour depth and time in such scouring.The QSI method of Bollaert's CSM computes scour depth in the plunge pool walls.The channel bottoms based on the total lift (quasi-steady lift (QSL) force) are applied to a protruding block, where the QSL force depends on the flow velocity and uplift pressure.Dasgupta proposed a numerical technique in 2011.They utilised 2D UDEC for modelling the rock mass to freely assess the block evacuation and fragile crack alongside the 3D CFD programming (ANSYS FLUENT) to decide the erosive limits.The erodibility file (for example, Kirsten index) is first created to survey excavatability, and the heaviness of the erodibility list boundaries is inapplicable for evaluating erodibility in hydraulic structures.In most cases, contradictions are observed between the erodibility predictions of existing methods and the observations made after dams are used for operation.An example is the erosion observed in the drainage channel of the Mokolo Dam in South Africa, where the erosion process generated a 30 m deep gully.In most semi-theoretical methods, limited data are used to assess erodibility.Results are limited and do not cover every type of geomechanical situation, and the hydraulic conditions are observed at the dam spillways around the world.In the aforementioned methods [18,21,24], the effects of spillway geometry, especially spillway surface roughness, are not comprehensively studied.The relationship between the geomechanical and hydraulic parameters is not elucidated.
Regarding the consideration of erosive force of water, the velocity of flowing water (V), shear stress applied to a rock surface (τ b ), unit stream power dissipation of water (Π UD ), stress intensity (K I ) and lifting force (F L ) have been used in several studies as the hazard parameters during erosion.A unique method for determining the erosive force of water for rock mass erodibility evaluation is lacking.The existing methods also have several limitations and can be used only in specific cases.For example, various equations exist to determine the unit stream power dissipation of water (Π UD ), which was developed as a parameter on the basis of the internal flow conditions.Stress intensity (K I ) was initially developed for metallurgical cases and used only to determine the possibility of crack propagation in intact rocks but not in rock masses.In addition, stress intensity (K I ), which is used in Bollaert's method, is based on the maximum pressure in pool bottoms, in which Bollaert's CFM method can compute pressure heads just under the centreline of jets but not the exact pressure of water applied to joint tips.Developments for unlined spillway structures are also lacking.In Bollaert's DI approach, which is on the basis of impulsion and Newton's second law, the geomechanical and geometrical parameters of rock masses are not comprehensively considered, and only block size is used.Shear force (F sh ) is deemed zero owing to the vertical fracture consideration.The single rock block models and Bollaert's CSM are restricted when the fractures directions are not symmetrical and vertical.This situation is common in volcanic and heterogeneous rock types.Twodimensional nonrectangular obstructs are as yet prohibitive in their capacity to have contact to a rock mass, and the procedure is intrinsically 3D.Table 14 illustrates the advantages and disadvantages of various approaches.Table 15 shows the detailed information references about erodibility phenomena.

Figure 1 .
Figure 1.Main mechanical processes of rock erosion and the occurrence timescale [3].

Figure 1 .
Figure 1.Main mechanical processes of rock erosion and the occurrence timescale [3].

Figure 2 .
Figure 2. Classifications of the scour prediction methods.

Figure 2 .
Figure 2. Classifications of the scour prediction methods.
x = horizontal direction of flow; y = depth of water; Z b = the vertical distance of the bottom surface from the reference datum; ρ = water density (kg/m 3 ); V = velocity of Water 2021, 13, 3205 5 of 42 flowing water; g = gravitational acceleration (m/s 2 ); Q = total discharge amount (m 3 /s);

Figure 3 .
Figure 3. (a) Variation of properties over a differential flow section in an open channel.(b) Cross-sectional view of the flow channel [19].

Figure 3 .
Figure 3. (a) Variation of properties over a differential flow section in an open channel.(b) Cross-sectional view of the flow channel [19].

Figure 4 .
Figure 4. (a) Original dataset as presented in the study of Van Schalkwyk [24]; (b) after inclusion of data from Moore [24]; (c) scour threshold line presented by Annandale; (d) comparison of scour threshold lines.

Figure 4 .
Figure 4. (a) Original dataset as presented in the study of Van Schalkwyk [24]; (b) after inclusion of data from Moore [24]; (c) scour threshold line presented by Annandale; (d) comparison of scour threshold lines.

Figure 5 .
Figure 5. (a) Erosion classes determined according to the eGSI index; (b) erosion classes determined according to the RMEI index [39].

Figure 5 .
Figure 5. (a) Erosion classes determined according to the eGSI index; (b) erosion classes determined according to the RMEI index [39].

Figure 7
Figure7illustrates the main physical and mechanical processes used for determination of rock scouring (CSM).CFM comprises three modules: falling jet module, plunge pool module, and rock mass module.The latter executes the fracture mechanism.The falling jet module (Figure8) clarifies how the geometrical and hydraulic parameters of the plunging jets are changed from outlet exit to impact point of plunge pool surface.Issuance jet velocity (Vi), initial turbulence intensity (Tu), and issuance jet diameter (Di) identify the jet characteristics at the outlet exit.Tu represents the extent of the root-mean-square (RMS) estimation of the fluctuating velocity of the plunging jet to the mean axial velocity of the plunging jet.The plunging jet pathway relies upon air drag and ballistics.The falling jet module can determine the jet pathway length L from issuance to the impact point, the longitudinal distance of the plunging jet from issuance to the impact point, the jet velocity (Vj), and the jet diameter (Dj) at the impact point (Equations (9) and (10))[60,61].

Figure 9 .Figure 9 .
Figure 9. (a) Diversity of the mean dynamic pressure coefficient (along the jet centre line) as a performance of the dimensionless ratio of [55]; (b) nondimensional root-mean-square dynamic pressure coefficient for jet velocities higher than 20 m/s [62].

44 Figure 10 .
Figure 10.Encountered rock mass layer circumstances for (a) irregularly and (b) totally jointed rocks; (c) proposed structure for the essential geometrical designs of the discontinuously jointed rock [59].

Figure 10 .
Figure 10.Encountered rock mass layer circumstances for (a) irregularly and (b) totally jointed rocks; (c) proposed structure for the essential geometrical designs of the discontinuously jointed rock [59].

Figure 13 .
Figure 13.Conceptualisation of impulse dynamics on a block, based on the Bollaert DI method [59].Kinetic and potential energies have been used in Equation (26) for determining the uplift height (h up ) (m) from Equation (27) and comparing with block height (Z b ) (m).

Figure 14 .
Figure 14.Proposed criteria to evaluate the rock scour potential by DI [63,70].The Bollaert portrayal for the DI method has two ambiguities, especially in enlisting (F u − F o ).A comparable delineation was determined for DI and portraying C I coefficient in several articles of Bollaert [59,61-63].The C I , (F u − F o ), and uplift height (h up ) are difficult to determine considering this description.The final dimension of (F u − F o ) is in m instead of kN by considering this description.The unit weight of water (Υ w ) and area of block surface (A) are not intended for converting the DI coefficient (C I ) to force.In 2004, Annandale corrected Bollaert's problem in converting C I to force and proposed a modified equation for determining (h up )[71].However, Annandale did not consider the nondimensional time coefficient (T up ) in Equation(37).Bollaert published a new article and proposed the simplified CSM (SCSM) because of the difficulty in understanding the CSM[64].In SCSM, the DI method is explained stepwise.Bollaert corrected the problem in previous publications.This time, he converted C I to pressure instead of force and did not consider the block surface area (A).We present

Figure 15 .
Figure 15.(a) Plane jet deflection on a flat bottom and wall jet velocity profiles; (b) effective forces on block uplift [81].Figure 15.(a) Plane jet deflection on a flat bottom and wall jet velocity profiles; (b) effective forces on block uplift [81].

Figure 15 .
Figure 15.(a) Plane jet deflection on a flat bottom and wall jet velocity profiles; (b) effective forces on block uplift [81].Figure 15.(a) Plane jet deflection on a flat bottom and wall jet velocity profiles; (b) effective forces on block uplift [81].

Figure 17 .
Figure 17.Flowchart of the numerical approach.Figure 17.Flowchart of the numerical approach.

Figure 17 .
Figure 17.Flowchart of the numerical approach.Figure 17.Flowchart of the numerical approach.

Water 2021 ,
13,  x FOR PEER REVIEW 31 of 44 dynamic pressure and the time-dependent pressures inside cracks and fractures.The rock block removal at the water-rock interface explains the scouring mechanism.

Figure 20 .
Figure 20.Schematic of the spillway and removable block geometry for reliability analysis [83].

Figure 20 .
Figure 20.Schematic of the spillway and removable block geometry for reliability analysis [83].

Figure 21 .
Figure 21.(a) The percentage of poorly estimated cases for each approach.(b) Committed errors of various methods according to different erosion classes.(c) Over-and underestimation rates of various semi-theoretical approaches.

Figure 21 .
Figure 21.(a) The percentage of poorly estimated cases for each approach.(b) Committed errors of various methods according to different erosion classes.(c) Over-and underestimation rates of various semi-theoretical approaches.

Table 1 .
Classification of the erosion extent.

Table 2 .
Scour threshold equations by different researchers.

Table 2 .
Scour threshold equations by different researchers.

Table 5 .
Coefficients of five simplified scour formulae.

Table 6 .
Symbol notation list of the plunge pool semi-theoretical method.

Table 8 .
[59]t and output parameters of the falling jet and plunge pool modules[59].

Table 10 .
Parameters of Equations (

Table 11 .
[82]uplift stagnation pressures for the various setups of block projection and joint points[82].

block ) h block /e block C joint C surf C uplift
1 Water 2021, 13, x FOR PEER REVIEW 27 of 45

Table 11 .
[82]uplift stagnation pressures for the various setups of block projection and joint points[82].

Table 11 .
[82]uplift stagnation pressures for the various setups of block projection and joint points[82].

Table 11 .
[82]uplift stagnation pressures for the various setups of block projection and joint points[82].

Table 11 .
[82]uplift stagnation pressures for the various setups of block projection and joint points[82].

Table 11 .
[82]uplift stagnation pressures for the various setups of block projection and joint points[82].

Table 11 .
[82]uplift stagnation pressures for the various setups of block projection and joint points[82].

Table 11 .
[82]uplift stagnation pressures for the various setups of block projection and joint points[82].

Table 11 .
[82]uplift stagnation pressures for the various setups of block projection and joint points[82].
As referenced, Lblock is the block length (m) that appeared in Figure15b, and VX,max represents the diffusive jet velocity in different locations.The different distances form

Table 11 .
[82]uplift stagnation pressures for the various setups of block projection and joint points[82].

Table 12
[86]strates the downstream and upstream digressed pieces of the absolute stream rate for various impingement angles δ, which are introduced by Reich[86].

Table 12 .
[86]ream and downstream strayed pieces of the absolute discharge rate for various points of impingement δ[86].

Table 13 .
Hydraulic geometrical characteristics of the experiments.

Table 14 .
Advantages and disadvantages of the various approaches.