Characteristics of Positive Surges in a Rectangular Channel

A positive surge is an unsteady open channel flow motion characterized by an increase of flow depth. In previous experimental studies, a positive surge was typically induced by either a sudden increase of discharge in a channel or by the rapid closure of a downstream sluice gate, thus leading to a steep initial profile. However, in many instances, the evolution of a positive surge is of a progressive manner (e.g., in the downstream navigation canal during the emptying operation of lock chambers). In the present work, the inception and development of a positive surge induced by a progressive increase of discharge was investigated in a rectangular channel with a smooth bed. Both undular and breaking surges were studied. The results demonstrate that the maximum wave height at the first wave crest of an undular surge is in very close agreement with the McCowan theory. Additionally, the wave amplitude essentially shows a linearly increasing trend with an increasing surge Froude number up to Fr0 = 1.26 to 1.28, whereas it tends to suggest a power law reduction for larger surge Froude numbers. Moreover, the dispersion of undular surges is consistent with the linear wave theory only for surge Froude numbers close to unity. Overall, the present study demonstrates the unique features of positive surges induced by a progressive increase of discharge.


Introduction
A positive surge is an unsteady open channel flow motion characterized by an increase of flow depth [1].Positive surges are widely observed in man-made and natural channels.They may be induced by the emptying of lock chambers during lock operations [2] or by sluice gates installed along water supply canals for irrigation and power purposes [3].For some flow conditions, the front of a positive surge is followed by a train of undulations (i.e., an undular surge) and even behaves as a nearly vertical water wall (i.e., a breaking surge).The passage of undular and breaking surges is associated with significant variations in hydrodynamics (i.e., flow velocity, water level and Reynolds stresses), and therefore they can be expected to affect the mixing and advection processes of sediments in channels [4,5].In the case of navigation channels, these surges also potentially influence tow's maneuverability, thereby affecting the navigation safety and functional efficiency of locks [6].
The first systematic study of positive waves was performed by Favre [7], who investigated the characteristics of undular surges in a rectangular flume.In the experimental studies of Benet and Cunge [8], and Treske [6], the focus was on the properties of both undular and breaking surges in trapezoidal channels.New experiments were carried out by Soares Frazão and Zech [9] to investigate the characteristics of undular surges and to validate the proposed numerical scheme.More recent experimental work of the same kind has been performed by Chanson and co-workers [10][11][12][13][14][15].Koch and Chanson [10] focused on the free-surface and turbulent flow properties of both undular and breaking surges.Hydraulic studies by Chanson [11] and Gualtieri and Chanson [12] were mostly concerned with the effect of bed roughness on the wave height attenuation and dispersion of an undular surge.An attempt was made by Chanson [13] to compare the free-surface profile of the undulations with both the sinusoidal and cnoidal wave functions.In the studies by Gualtieri and Chanson [14] and Leng and Chanson [15], the free-surface properties were compared to the analytical solutions of Lemoine [16] and Anderson [17], based on the linear wave theory and the Boussinesq equations, respectively.In some other studies, the upstream propagation of positive surges in a sloping rectangular channel was also investigated [18].
Undular and breaking surges in previous experimental studies have typically been generated in two different ways.One way is to induce a sudden increase in discharge in a channel, while another way is to generate undular and breaking surges by the rapid closure of a downstream sluice gate, thus leading to steep initial water profiles.However, in many instances, these surges occurred in a progressive manner.For example, in the case of the emptying operation of lock chambers, a positive surge with a rather smooth water profile is first generated in the downstream channel; subsequently, the front of the positive surge steepens and undulations slowly grow at the front head during propagation; and finally, an undular surge with some breaking or a breaking surge may form, depending on the flow conditions.
This work is focused on the inception and development of a positive surge induced by the progressive opening of an upstream plate gate.The free-surface characteristics were systematically investigated during the downstream propagation of positive surges.

Experimental Setup and Experiment
The experimental arrangement is schematically shown in Figure 1.The experimental facility consisted of a rectangular upstream reservoir, a rectangular flume and a pipeline connecting the two parts together.Both the length and width of the upstream reservoir were 1 m.The horizontal channel was made of a smooth poly (vinyl chloride) (PVC) bed and Plexiglass sidewalls, 38 m long and 0.3 m wide, with a height of 0.35 m.The flume was connected to a pipeline of the same width with a curved entrance at the upstream end of the channel bed.The downstream end of the flume was closed by a vertical rigid wall.The pipeline system was equipped with a plate gate and its shaft.The width and height of the gate were 0.3 m and 0.1 m, respectively.The gate was driven vertically by a stepper motor, which could be used to exactly reproduce any desirable opening function.The total traveling height of the gate was equal to its height.the characteristics of undular surges and to validate the proposed numerical scheme.More recent experimental work of the same kind has been performed by Chanson and co-workers [10][11][12][13][14][15].Koch and Chanson [10] focused on the free-surface and turbulent flow properties of both undular and breaking surges.Hydraulic studies by Chanson [11] and Gualtieri and Chanson [12] were mostly concerned with the effect of bed roughness on the wave height attenuation and dispersion of an undular surge.An attempt was made by Chanson [13] to compare the free-surface profile of the undulations with both the sinusoidal and cnoidal wave functions.In the studies by Gualtieri and Chanson [14] and Leng and Chanson [15], the free-surface properties were compared to the analytical solutions of Lemoine [16] and Anderson [17], based on the linear wave theory and the Boussinesq equations, respectively.In some other studies, the upstream propagation of positive surges in a sloping rectangular channel was also investigated [18].Undular and breaking surges in previous experimental studies have typically been generated in two different ways.One way is to induce a sudden increase in discharge in a channel, while another way is to generate undular and breaking surges by the rapid closure of a downstream sluice gate, thus leading to steep initial water profiles.However, in many instances, these surges occurred in a progressive manner.For example, in the case of the emptying operation of lock chambers, a positive surge with a rather smooth water profile is first generated in the downstream channel; subsequently, the front of the positive surge steepens and undulations slowly grow at the front head during propagation; and finally, an undular surge with some breaking or a breaking surge may form, depending on the flow conditions.This work is focused on the inception and development of a positive surge induced by the progressive opening of an upstream plate gate.The free-surface characteristics were systematically investigated during the downstream propagation of positive surges.

Experimental Setup and Experiment
The experimental arrangement is schematically shown in Figure 1.The experimental facility consisted of a rectangular upstream reservoir, a rectangular flume and a pipeline connecting the two parts together.Both the length and width of the upstream reservoir were 1 m.The horizontal channel was made of a smooth poly (vinyl chloride) (PVC) bed and Plexiglass sidewalls, 38 m long and 0.3 m wide, with a height of 0.35 m.The flume was connected to a pipeline of the same width with a curved entrance at the upstream end of the channel bed.The downstream end of the flume was closed by a vertical rigid wall.The pipeline system was equipped with a plate gate and its shaft.The width and height of the gate were 0.3 m and 0.1 m, respectively.The gate was driven vertically by a stepper motor, which could be used to exactly reproduce any desirable opening function.The total traveling height of the gate was equal to its height.The initial flow depths were recorded with point gauges in the upstream reservoir and the channel.The unsteady water depth in the channel was measured using a set of displacement meters spaced along the channel centerline, and a sampling frequency of 200 Hz was used for all the measurements.For all experiments, flow measurements were started 60 s before the gate opening, The initial flow depths were recorded with point gauges in the upstream reservoir and the channel.The unsteady water depth in the channel was measured using a set of displacement meters spaced along the channel centerline, and a sampling frequency of 200 Hz was used for all the measurements.For all experiments, flow measurements were started 60 s before the gate opening, and were stopped when the surge reached the downstream rigid wall to avoid wave reflection interference.
It was assumed that the plate gate was raised at a constant speed until complete gate opening.The drop height H was defined as the initial difference in elevation between the water surfaces of the upstream reservoir and the channel.Hydraulic considerations indicated that the unsteady water depth h in the channel should essentially depend on the initial water depth in channel h 0 , the drop height H, and the amount of time needed to fully open the gate t v .In the present study, the initial water depth in the channel was set to 8 cm by referring to the studies of Treske [6] and Koch and Chanson [10].A summary of the experimental measurements is shown in Table 1.Overall, a large number of experimental runs were carried out in four series, denoted as A, B, C, and D, corresponding to four different times taken to fully open the gate of 20, 25, 30, and 40 s with four different drop heights.The experimental conditions were selected to generate both undular and breaking surges, taking into account the influence of undulation shape deviations on the characteristics of surges with close Froude numbers, which is defined in the system of coordinates in translation with the surge and is given by where S = surge front celerity positive downstream and g = gravity acceleration.

Basic Flow Patterns
A positive surge was generated in the flume by raising the plate gate at a constant speed.During the downstream propagation of the positive surge, different wave types may occur at the front head.Here, the overall developing processes observed at the front of positive surges have been consolidated into four major flow stages, and are described below.It should be noted that the number of flow stages for each experimental run was variable, depending on the drop height and the time taken to fully open the gate.

•
In Stage 1, the front of the positive wave is rather smooth.

•
Stage 2 corresponds to the generating process of the undular surges, and can be divided into two sub-stages: Stage2a and Stage2b.In Stage 2a, the front head of the positive wave gradually becomes "rough" and develops into a series of cascading steps at the end of this stage; in Stage 2b, the cascading steps develop into a train of well-formed undulations (i.e., non-breaking undular surges) (Figure 2).Additionally, some sidewall shock waves can be observed in this stage that develop upstream of the first wave crest and intersect next to the first crest.This stage ends when some small waves break at the first wave crest (i.e., breaking undular surges) (Figure 3).

•
In Stage 3, the wave amplitude decreases and the free-surface undulations become flatter.This stage ends with the disappearance of fluctuating characteristics of water depth at the surge front.
The tests indicate that transitions between different flow stages can be characterized solely with the surge Froude number (Table 2).Moreover, the surge Froude numbers corresponding to these transitions are close to those of Chanson [11,13] and Leng and Chanson [15] for both smooth and rough beds, although these studies focused on the positive surges induced by the rapid closure of a downstream sluice gate.The tests indicate that transitions between different flow stages can be characterized solely with the surge Froude number (Table 2).Moreover, the surge Froude numbers corresponding to these transitions are close to those of Chanson [11,13] and Leng and Chanson [15] for both smooth and rough beds, although these studies focused on the positive surges induced by the rapid closure of a downstream sluice gate.The tests indicate that transitions between different flow stages can be characterized solely with the surge Froude number (Table 2).Moreover, the surge Froude numbers corresponding to these transitions are close to those of Chanson [11,13] and Leng and Chanson [15] for both smooth and rough beds, although these studies focused on the positive surges induced by the rapid closure of a downstream sluice gate.The tests indicate that transitions between different flow stages can be characterized solely with the surge Froude number (Table 2).Moreover, the surge Froude numbers corresponding to these transitions are close to those of Chanson [11,13] and Leng and Chanson [15] for both smooth and rough beds, although these studies focused on the positive surges induced by the rapid closure of a downstream sluice gate.In a rectangular horizontal channel, when neglecting friction loss, the combined application of the continuity and momentum equations across the surge front yields where h conj and h 0 are the conjugate and initial depths, respectively.Equation ( 2) is the famous Bélanger equation.In the present study, the conjugate depth h conj was defined as the average of the first wave crest and the trough depths for the undular surges, while it was the water depth immediately behind the surge front for breaking surges.The ratio of conjugate depths, h conj /h 0 , is plotted against the surge Froude number, Fr 0 , in Figure 5 for both undular and breaking surges.Although the experimental trend is comparable to that produced by Equation ( 2), the entire dataset shows generally higher values of h conj /h 0 for 1.07 < Fr 0 < 1.57.This may be partially attributed to the uncertainty of estimating the conjugate depth; that is, in the undular surges, the conjugate depth was calculated based on a symmetrical undulation profile which is inconsistent with the experimental data, while the measurement of the conjugate depth was affected adversely by the large free-surface fluctuations behind the surge front in the breaking surges.It is worth noting that the data in flow Stage 2b are more scattered than in flow Stages 3 and 4 for a given surge Froude number, which might imply larger deviations of wave shapes for non-breaking surges than for breaking undular and breaking surges.
Water 2018, 10, x FOR PEER REVIEW 5 of 12

Ratio of Conjugate Depths
In a rectangular horizontal channel, when neglecting friction loss, the combined application of the continuity and momentum equations across the surge front yields ( ) where hconj and h0 are the conjugate and initial depths, respectively.Equation ( 2) is the famous Bélanger equation.In the present study, the conjugate depth hconj was defined as the average of the first wave crest and the trough depths for the undular surges, while it was the water depth immediately behind the surge front for breaking surges.The ratio of conjugate depths, hconj/h0, is plotted against the surge Froude number, Fr0, in Figure 5 for both undular and breaking surges.
Although the experimental trend is comparable to that produced by Equation ( 2), the entire dataset shows generally higher values of hconj/h0 for 1.07 < Fr0 < 1.57.This may be partially attributed to the uncertainty of estimating the conjugate depth; that is, in the undular surges, the conjugate depth was calculated based on a symmetrical undulation profile which is inconsistent with the experimental data, while the measurement of the conjugate depth was affected adversely by the large free-surface fluctuations behind the surge front in the breaking surges.It is worth noting that the data in flow Stage 2b are more scattered than in flow Stages 3 and 4 for a given surge Froude number, which might imply larger deviations of wave shapes for non-breaking surges than for breaking undular and breaking surges.

Free-Surface Undulation Characteristics (1) Wave height of the first wave crest
The wave height of the first wave crest (h1c − h0) is plotted in a dimensionless form against Fr0 in Figure 6.For non-breaking undular surges in Stage 2b, the values of (h1c − h0)/h0 show a trend of linear increase with Fr0 in general.For breaking undular surges in Stage 3, a sharp decrease in wave height is observed shortly after the appearance of some wave breaking at the first wave crest, and subsequently, (h1c − h0)/h0 exhibits the same increasing trend as the data in flow Stage 2b, but at a much smaller rate.It is worth noting that in flow Stage 2b, (h1c − h0)/h0 is close to the solitary wave theory for surge Froude numbers Fr0 less than 1.1, while the data show consistently higher values of (h1c − h0)/h0 for larger surge Froude numbers.

Free-Surface Undulation Characteristics (1) Wave height of the first wave crest
The wave height of the first wave crest (h 1c − h 0 ) is plotted in a dimensionless form against Fr 0 in Figure 6.For non-breaking undular surges in Stage 2b, the values of (h 1c − h 0 )/h 0 show a trend of linear increase with Fr 0 in general.For breaking undular surges in Stage 3, a sharp decrease in wave height is observed shortly after the appearance of some wave breaking at the first wave crest, and subsequently, (h 1c − h 0 )/h 0 exhibits the same increasing trend as the data in flow Stage 2b, but at a much smaller rate.It is worth noting that in flow Stage 2b, (h 1c − h 0 )/h 0 is close to the solitary wave theory for surge Froude numbers Fr 0 less than 1.1, while the data show consistently higher values of (h 1c − h 0 )/h 0 for larger surge Froude numbers.
The wave height data of the first wave crest are further compared with previous experimental studies in Figure 7.It is found that present data are consistent with the trend observed by Peregrine [19], Koch and Chanson [10] and Chanson [11].In Figure 7, the data are also compared with the McCowan theory, which shows the maximum wave heights attained by solitary waves.It can be seen that there is an excellent agreement between this paper's results and McCowan's theory [20].(2) Wave amplitude, length and steepness The wave amplitude, length and steepness are shown in Figure 8 in dimensionless forms (i.e., aw/h0, Lw/h0, and aw/Lw).Both wave amplitude and wave length were calculated according to the definitions by Koch and Chanson [10]; that is, the wave amplitude, aw, was half of the difference between the water depth at the first wave crest and at the first wave trough, while the wave length, Lw, was defined between the first and second crests.The data are compared with the linear wave theory solution of Lemoine, the Boussinesq equation solution of Anderson, earlier experimental studies and field observations [21][22][23].
It can be seen from Figure 8a that the wave amplitude increases monotonically with an increasing surge Froude number to a local maximum for non-breaking undular surges in Stage 2b.However, an abrupt decrease in wave amplitude is observed immediately after the appearance of The wave height data of the first wave crest are further compared with previous experimental studies in Figure 7.It is found that present data are consistent with the trend observed by Peregrine [19], Koch and Chanson [10] and Chanson [11].In Figure 7, the data are also compared with the McCowan theory, which shows the maximum wave heights attained by solitary waves.It can be seen that there is an excellent agreement between this paper's results and McCowan's theory [20].
Water 2018, 10, x FOR PEER REVIEW 6 of 12 The wave height data of the first wave crest are further compared with previous experimental studies in Figure 7.It is found that present data are consistent with the trend observed by Peregrine [19], Koch and Chanson [10] and Chanson [11].In Figure 7, the data are also compared with the McCowan theory, which shows the maximum wave heights attained by solitary waves.It can be seen that there is an excellent agreement between this paper's results and McCowan's theory [20].(2) Wave amplitude, length and steepness The wave amplitude, length and steepness are shown in Figure 8 in dimensionless forms (i.e., aw/h0, Lw/h0, and aw/Lw).Both wave amplitude and wave length were calculated according to the definitions by Koch and Chanson [10]; that is, the wave amplitude, aw, was half of the difference between the water depth at the first wave crest and at the first wave trough, while the wave length, Lw, was defined between the first and second crests.The data are compared with the linear wave theory solution of Lemoine, the Boussinesq equation solution of Anderson, earlier experimental studies and field observations [21][22][23].
It can be seen from Figure 8a that the wave amplitude increases monotonically with an increasing surge Froude number to a local maximum for non-breaking undular surges in Stage 2b.However, an abrupt decrease in wave amplitude is observed immediately after the appearance of (2) Wave amplitude, length and steepness The wave amplitude, length and steepness are shown in Figure 8 in dimensionless forms (i.e., a w /h 0 , L w /h 0 , and a w /L w ).Both wave amplitude and wave length were calculated according to the definitions by Koch and Chanson [10]; that is, the wave amplitude, a w , was half of the difference between the water depth at the first wave crest and at the first wave trough, while the wave length, L w , was defined between the first and second crests.The data are compared with the linear wave theory solution of Lemoine, the Boussinesq equation solution of Anderson, earlier experimental studies and field observations [21][22][23].

numbers.
Figure 8c presents the dimensionless wave steepness data.It can be seen that the wave steepness exhibits a similar trend to wave amplitude-aw/Lw shows an increasing and decreasing trend for nonbreaking and breaking undular surges, respectively.Additionally exhibited in this figure are the data sets obtained in undular hydraulic jumps [24,25].For surge Froude numbers less than 1.29, the present data set is in excellent agreement with undular jump data.(

3) Dispersion characteristics
The undular surges investigated in this study have the property of net mass transport, which is in conjunction with the finite amplitude water waves [26].For a finite amplitude wave, the linear wave theory yields a dispersion relationship between the wave length Lw, wave period T and water depth h: It can be seen from Figure 8a that the wave amplitude increases monotonically with an increasing surge Froude number to a local maximum for non-breaking undular surges in Stage 2b.However, Water 2018, 10, 1473 8 of 12 an abrupt decrease in wave amplitude is observed immediately after the appearance of wave breaking at the first wave crest for breaking undular surges in Stage 3; subsequently, a gradual declining trend was held before the disappearance of free-surface undulations.The relationship between a w /h 0 and Fr 0 for undular surges is well described by the following equations a w h 0 = 7.50Fr 0 −14.79 for 1.26-1.28< Fr 0 < 1.45-1.50 Overall, the present data exhibit a qualitative agreement with all previous studies, as indicated in Figure 8a.Moreover, although neither the linear wave theory nor the Boussinesq equation solution can accurately predict the wave amplitude before surge breaking, the data tend to get close to Lemoine's solution for surge Froude numbers in the range of 1.2 to 1.28.
The dimensionless wave length data are shown in Figure 8b.Although the experimental trend is close to that predicted by the theories of both Lemoine and Anderson and earlier experimental studies, the data show consistently lower wave lengths, especially for very low or high Froude numbers.
Figure 8c presents the dimensionless wave steepness data.It can be seen that the wave steepness exhibits a similar trend to wave amplitude-a w /L w shows an increasing and decreasing trend for non-breaking and breaking undular surges, respectively.Additionally exhibited in this figure are the data sets obtained in undular hydraulic jumps [24,25].For surge Froude numbers less than 1.29, the present data set is in excellent agreement with undular jump data. (

3) Dispersion characteristics
The undular surges investigated in this study have the property of net mass transport, which is in conjunction with the finite amplitude water waves [26].For a finite amplitude wave, the linear wave theory yields a dispersion relationship between the wave length L w , wave period T and water depth h: However, the undular surge is typical of a translation wave, i.e., water particles underneath the water surface do not move backwards, indicating marked differences between undular surges and finite amplitude waves.Therefore, a dimensionless dispersion parameter D can be defined for undular surges as Plots of D versus Fr 0 for undular surges are presented in Figure 9.It is shown that a definite correlation exists between the two dimensionless variables, with essentially no influence of undular surge types.In the experimental range 1.07 < Fr 0 < 1.50, the following empirical equation can be proposed by fitting both data sets (i.e., for non-breaking and breaking undular surges) It is clearly demonstrated that for a surge Froude number Fr 0 less than 1.1, the dispersion parameter D is close to 1, and hence the dispersion relationship based on the linear wave theory is assumed to hold.

Energy and Momentum Fluxes Properties
A positive surge is seen as a quasi-steady flow situation by an observer travelling at the surge speed S. In a rectangular, horizontal channel, if both the rate of energy dissipation and friction loss are negligible, the equations of conservation of energy and momentum can be rewritten as [26] where E is the specific energy; M is the specific momentum; h is the flow depth; E* and M* correspond to the dimensionless specific energy and specific momentum, respectively; and h c is the critical depth in the system of coordinates in translation with the undular surge front, and is defined as This pair of equations may be viewed as the parametric representation of M* = f (E*), with h/h c as a parameter.The function E* − M* has two branches (the red dash line in Figure 10), intersecting at (1.5, 1.5).The lower branch of the curve E* − M* corresponds to a supercritical flow while the upper branch represents a subcritical flow.Figure 10 shows a comparison of M* = f (E*) and the entire dataset in this study.It is expected that the initial flow data are located on the supercritical branch while the corresponding conjugate flow data are on the subcritical branch.It can also be concluded from Figure 10, that at the subcritical branch of the curve E* − M*, some overlaps exist between regions corresponding to different flow stages, indicating marked differences between travelling positive surges and stationary hydraulic jumps [27].
Water 2018, 10, x FOR PEER REVIEW 9 of 12 where E is the specific energy; M is the specific momentum; h is the flow depth; E* and M* correspond to the dimensionless specific energy and specific momentum, respectively; and hc is the critical depth in the system of coordinates in translation with the undular surge front, and is defined as    Figure 11 shows a comparison of the dimensionless specific energy in the initial flow (E 0 *) and the corresponding conjugate flow (E conj *).It can be found that the values of E conj * are basically larger than that of E 0 * for non-breaking undular surges in Stage 2b.This is probably due to the pressure distribution beneath an undular surge that is less than hydrostatic beneath the wave crest and greater than hydrostatic beneath the wave trough.With the ensuing development of undular surges, E conj * approaches E 0 * for breaking undular surges in Stage 3, and is essentially less than E 0 * for breaking surges in Stage 4. This may be primarily attributed to energy dissipation due to wave breaking, and hence a smaller specific energy in the conjugated flow than in the initial flow is obtained in breaking surges.
In Figure 12, a comparison of the dimensionless specific momentum in the initial flow (M 0 *) and the corresponding conjugate flow (M conj *) is presented.The values of M conj * are consistently larger than that of M 0 * for the entire present test range.The main reason for this may be the deviation from the hydrostatic and the complicated velocity distribution in undular and breaking surges.

Conclusions
The free-surface properties of a positive surge induced by the progressive opening of an upstream plate gate were investigated experimentally in a rectangular channel with a smooth bed.

Conclusions
The free-surface properties of a positive surge induced by the progressive opening of an upstream plate gate were investigated experimentally in a rectangular channel with a smooth bed.

Conclusions
The free-surface properties of a positive surge induced by the progressive opening of an upstream plate gate were investigated experimentally in a rectangular channel with a smooth bed.

Figure 1 .
Figure 1.Sketch of the experimental facility.

Figure 1 .
Figure 1.Sketch of the experimental facility.

Figure 4 .
Figure 4. Breaking surges in Stage 4: (a) Lateral view; (b) looking upstream at the incoming wave.

Figure 4 .
Figure 4. Breaking surges in Stage 4: (a) Lateral view; (b) looking upstream at the incoming wave.

Figure 4 .
Figure 4. Breaking surges in Stage 4: (a) Lateral view; (b) looking upstream at the incoming wave.

Figure 4 .
Figure 4. Breaking surges in Stage 4: (a) Lateral view; (b) looking upstream at the incoming wave.

Figure 5 .
Figure 5. Variation of h conj /h 0 with Fr 0 for both undular and breaking surges.

Figure 7 .
Figure 7. Wave height at the first wave crest for undular surges.

Figure 7 .
Figure 7. Wave height at the first wave crest for undular surges.

Figure 7 .
Figure 7. Wave height at the first wave crest for undular surges.

Figure 8 .
Figure 8. Characteristics of undular surges as functions of Fr 0 : (a) Dimensionless wave amplitude a w /h 0 ; (b) dimensionless wave length L w /h 0 ; (c) wave steepness a w /L w .
equations may be viewed as the parametric representation of M* = f(E*), with h/hc as a parameter.The function E* − M* has two branches (the red dash line in Figure10), intersecting at (1.5, 1.5).The lower branch of the curve E* − M* corresponds to a supercritical flow while the upper branch represents a subcritical flow.Figure10shows a comparison of M* = f(E*) and the entire dataset in this study.It is expected that the initial flow data are located on the supercritical branch while the corresponding conjugate flow data are on the subcritical branch.It can also be concluded from Figure10, that at the subcritical branch of the curve E* − M*, some overlaps exist between regions corresponding to different flow stages, indicating marked differences between travelling positive surges and stationary hydraulic jumps[27].

Figure 11
Figure11shows a comparison of the dimensionless specific energy in the initial flow (E0*) and the corresponding conjugate flow (Econj*).It can be found that the values of Econj* are basically larger than that of E0* for non-breaking undular surges in Stage 2b.This is probably due to the pressure distribution beneath an undular surge that is less than hydrostatic beneath the wave crest and greater than hydrostatic beneath the wave trough.With the ensuing development of undular surges, Econj* approaches E0* for breaking undular surges in Stage 3, and is essentially less than E0* for breaking surges in Stage 4. This may be primarily attributed to energy dissipation due to wave breaking, and hence a smaller specific energy in the conjugated flow than in the initial flow is obtained in breaking surges.In Figure12, a comparison of the dimensionless specific momentum in the initial flow (M0*) and the corresponding conjugate flow (Mconj*) is presented.The values of Mconj* are consistently larger

Figure 9 .
Figure 9. Variation of D with Fr 0 for undular surges.

Water 2018 ,
10, x FOR PEER REVIEW 10 of 12than that of M0* for the entire present test range.The main reason for this may be the deviation from the hydrostatic and the complicated velocity distribution in undular and breaking surges.

Figure 10 .
Figure 10.Dimensionless relationship between the momentum and energy fluxes for both undular and breaking surges.

Figure 10 .
Figure 10.Dimensionless relationship between the momentum and energy fluxes for both undular and breaking surges.

Water 2018 ,
10, x FOR PEER REVIEW 10 of 12than that of M0* for the entire present test range.The main reason for this may be the deviation from the hydrostatic and the complicated velocity distribution in undular and breaking surges.

Figure 10 .
Figure 10.Dimensionless relationship between the momentum and energy fluxes for both undular and breaking surges.

Figure 11 .
Figure 11.Comparison of E 0 * and E conj * for both undular and breaking surges.

Water 2018 ,
10, x FOR PEER REVIEW 10 of 12than that of M0* for the entire present test range.The main reason for this may be the deviation from the hydrostatic and the complicated velocity distribution in undular and breaking surges.

Figure 10 .
Figure 10.Dimensionless relationship between the momentum and energy fluxes for both undular and breaking surges.

Figure 12 .
Figure 12.Comparison of M 0 * and M conj * for both undular and breaking surges.

Table 2 .
Surge Froude number values for flow stage transitions.

Table 2 .
Surge Froude number values for flow stage transitions.

Table 2 .
Surge Froude number values for flow stage transitions.

Table 2 .
Surge Froude number values for flow stage transitions.