Development of Bubble Characteristics on Chute Spillway Bottom

Chute aerators introduce a large air discharge through air supply ducts to prevent cavitation erosion on spillways. There is not much information on the microcosmic air bubble characteristics near the chute bottom. This study was focused on examining the bottom air-water flow properties by performing a series of model tests that eliminated the upper aeration and illustrated the potential for bubble variation processes on the chute bottom. In comparison with the strong air detrainment in the impact zone, the bottom air bubble frequency decreased slightly. Observations showed that range of probability of the bubble chord length tended to decrease sharply in the impact zone and by a lesser extent in the equilibrium zone. A distinct mechanism to control the bubble size distribution, depending on bubble diameter, was proposed. For bubbles larger than about 1–2 mm, the bubble size distribution followed a—5/3 power-law scaling with diameter. Using the relationship between the local dissipation rate and bubble size, the bottom dissipation rate was found to increase along the chute bottom, and the corresponding Hinze scale showed a good agreement with the observations.


Introduction
Chute aerators, regarded as a cost-effective method, have been widely used owing to their proven success as a countermeasure in the Grand Coulee dam since 1960 [1].Previous studies have provided a significant amount of information on the macroscopic air-water flow properties [2][3][4][5][6][7][8][9][10][11][12][13], such as the air concentration, air entrainment and detrainment.Rutschmann and Hager [2], combining both model and prototype data, proposed two approaches to calculate the air entrainment coefficient β = q a /q, where q a = air discharge and q = water discharge.The air entrainment coefficient was also investigated by, among others, Shi [14], Pinto [15,16], and Low [17].Chanson [3] investigated the characteristics of chute aerator flow on the distributions of various bubble parameters (void fraction, velocity, air discharge, turbulent diffusivity, etc.) for different inflow Froude numbers.He also provided a better understanding of the air entrainment processes and found a strong detrainment process occurred in the impact zone.Kramer et al. [5] divided the chute aerator into five zones: inflow region; air detrainment region; minimum air concentration region; air entrainment region; and uniform mixture flow region, and systematically investigated the air detrainment gradients in the far-field of chute aerators.Pfister and Hager [6][7][8] measured the air concentration and divided the chute aerator into three zones: jet zone, reattachment and spray zone, and far-field zone.The authors also discussed the air entrainment and proposed two comprehensive methods to calculate the average and bottom air concentration development using deflectors and offsets along the chute.Avoiding the upper aeration effect, Bai et al. [9] divided the chute aerator into three zones and found that an intensive air detrainment occurred in the impact zone and the detrained air escaped into the cavity zone.The previous experimental investigations are summarized in Table 1.Bai [9,10] 4.0-9.00.15 3.3-7.4Model (W = 0.25 m) Shi [14] 14.0 0.058 18.6 Model (W = 0.20 m) Pinto [15,16] 19. 5 Many studies concerning the minimum or critical air concentration to avoid cavitation erosion have been conducted [18,19].Peterka [18] found that when the average air concentration was between 0.01 and 0.06, no cavitation erosion was observed.Up to now, there is no reliable design guideline for the distance required between two aerators to produce sufficient bottom air concentration [5].Based on prototype observations on Russian dams, Semenkovand Lantyaev [20] provided an average bottom air concentration decay of 0.40% to 0.80% per chute meter.
However, air concentration governs the comprehensive performance of the bubbles, including their chord length and frequency form.In particular, Robinson [21], Xu et al. [22], and Brujan [23] adopted advanced techniques to investigate the mechanism of erosion damage and found that the interaction between the air bubbles and cavitation bubbles was the mechanism that assisted in the prevention of erosion damage.These studies indicate that the bubble size and frequency are significant parameters for the mechanism of erosion damage.However, negligible research has been conducted on the spillway bottom air microcosmic characteristics such as the bottom bubble frequency and chord size.Such investigations are difficult because of the complex nature of the flow and strong impaction of the bottom.In the present study, the development of the bottom bubble characteristics was measured and analyzed.

Hydraulics Model
The experiments were conducted in a rectangular chute that was 0.25-m wide with a variable bottom slope 5.71 • ≤ α ≤ 18.2 • (Figure 1).The measured section was 5-m long and consisted of a smooth convergent nozzle with a width of 0.25 m and height of 0.15 m.The approach flow depth was h 0 = 0.15 m and the emergence angle was 0 • ≤ θ 0 ≤ 14.1 • .The water discharge was supplied by a large tank (3.1-m wide, 3.2-m long, and 9.2-m high) and was measured using a rectangular sharp-crested weir with an accuracy of approximately 1%.In order to gain sufficient velocity, the nozzle was located in the bottom of the tank.A range of approach flow velocities was set from 4.0 m/s to 9.0 m/s.Observations indicated that the approach of the flow was stable.Hence, the Froude number F 0 (V 0 /(gh 0 ) 0.5 ) was within the range of 3.3 to 7.4 (Table 2), where g was the gravity.x was the streamwise coordinate along the chute bottom, h s was the offset height, L was the cavity length, and L m and L D were the specific lengths [9].All measurements were performed on the channel centre line.The air-water flow properties were measured with a phase-detection needle probe.The working principle of a conductivity probe was based on the difference in the voltage indices at the platinum tip between the air and water phases.The response time of this sensor was less than 10 μs.The signals from the conductivity probe were recorded at a scan rate of 100 kHz per for a 40 s scan period.The accuracy of the air concentration was ΔC/C = 3%, and the bubble frequency was Δf/f = 1%, where C was the air concentration and f was the bubble frequency.The length of the impact zone was very short.A few studies in the impact zone have been performed with a space of 0.03 m along the chute bottom.

Air Concentration and Bubble Frequency Profiles
A typical photo of the offset-aerator flow is shown in Figure 2, where four zones were introduced.The typical distributions of the air concentration and bubble frequency in the different zones are depicted in Figure 2   All measurements were performed on the channel centre line.The air-water flow properties were measured with a phase-detection needle probe.The working principle of a conductivity probe was based on the difference in the voltage indices at the platinum tip between the air and water phases.The response time of this sensor was less than 10 µs.The signals from the conductivity probe were recorded at a scan rate of 100 kHz per for a 40 s scan period.The accuracy of the air concentration was ∆C/C = 3%, and the bubble frequency was ∆f /f = 1%, where C was the air concentration and f was the bubble frequency.The length of the impact zone was very short.A few studies in the impact zone have been performed with a space of 0.03 m along the chute bottom.

Air Concentration and Bubble Frequency Profiles
A typical photo of the offset-aerator flow is shown in Figure 2, where four zones were introduced.The typical distributions of the air concentration and bubble frequency in the different zones are depicted in Figure 2 as a function of the dimensionless distance perpendicular to the bottom, i.e., z/h 0 .Owing to a sufficient depth of the approach flow, the unaerated water prevented the upper aeration effect.bottom, i.e., z/h0.Owing to a sufficient depth of the approach flow, the unaerated water prevented the upper aeration effect.In the cavity zone (0 < x/L < 1), air entrainment occurs in the form of air bubbles entrapped into the lower flow and the air concentration decreased from the lower surface to the vertical direction.At x/L = 1, the lower jet impacted on the chute bottom.Further downstream (x/L > 1), the distributions of air concentration were different from the cavity zone.Figure 2a shows that the air concentration distributions downstream of the cavity zone, exhibiting an air concentration peak in the turbulent shear region, whereas the air concentration increases from the chute bottom to its maximum value, and then decreases in the vertical direction.
Note that: (1) for x > L, the maximum value air concentration Cm decreased sharply in the impact zone, while in the equilibrium zone, the peak air concentration Cm decreased slightly.(2) Experimental results indicated the existence of a maximum value fm in bubble frequency (Figure 2b), but its location differed from the locus of the maximum void fraction as noticed by Bai et al. [10].In the cavity, the location of fm approached C = 0.50 with the development of the low jet aeration.The location of fm was approximately the same as that of the section maximum air concentration Cm at a section in the impact zone, but lower than that in the equilibrium zone.
For more than 25 experiments conducted in this model, only one set of experiments is shown herein to illustrate the basic chute aerator flow property [9].However, equations presented hereafter are based on all experiments.In the cavity zone (0 < x/L < 1), air entrainment occurs in the form of air bubbles entrapped into the lower flow and the air concentration decreased from the lower surface to the vertical direction.At x/L = 1, the lower jet impacted on the chute bottom.Further downstream (x/L > 1), the distributions of air concentration were different from the cavity zone.Figure 2a shows that the air concentration distributions downstream of the cavity zone, exhibiting an air concentration peak in the turbulent shear region, whereas the air concentration increases from the chute bottom to its maximum value, and then decreases in the vertical direction.
Note that: (1) for x > L, the maximum value air concentration C m decreased sharply in the impact zone, while in the equilibrium zone, the peak air concentration C m decreased slightly.(2) Experimental results indicated the existence of a maximum value f m in bubble frequency (Figure 2b), but its location differed from the locus of the maximum void fraction as noticed by Bai et al. [10].In the cavity, the location of f m approached C = 0.50 with the development of the low jet aeration.The location of f m was approximately the same as that of the section maximum air concentration C m at a section in the impact zone, but lower than that in the equilibrium zone.
For more than 25 experiments conducted in this model, only one set of experiments is shown herein to illustrate the basic chute aerator flow property [9].However, equations presented hereafter are based on all experiments.

Bottom Air Bubble Characteristics
The entrained air from the chute aerator aims to prevent erosion damage.Recent investigations have indicated that the bubble size and frequency are the significant factors preventing erosion damage [22,23].The air bubbles generated in the cavity zone travelled through the impact zone to continue downstream.For the chute aerator flows, a large amount of air entrained in the cavity zone was only partially transported to the downstream flow.Under the effect of rollers, large bubbles broke into small bubbles before being transported downstream zone simultaneously in the impact zone.[7] is shown in Figure 3, which displays a similar trend downstream of the chute aerator.However, the present values are generally larger than those of Pfister and Hager [7].At a sufficient water depth, the detrained air bubbles cannot escape into the atmosphere through the upper surface; it can only escape into the cavity zone, which may cause the level of C b to be larger than that of previous studies.In addition, for a shorter length of the impact zone, the fewer measurement sections may cause an insufficient sharp decrease.(4) It can be concluded that the air transportation regularity is different in the impact and equilibrium zones.
The entrained air from the chute aerator aims to prevent erosion damage.Recent investigations have indicated that the bubble size and frequency are the significant factors preventing erosion damage [22,23].The air bubbles generated in the cavity zone travelled through the impact zone to continue downstream.For the chute aerator flows, a large amount of air entrained in the cavity zone was only partially transported to the downstream flow.Under the effect of rollers, large bubbles broke into small bubbles before being transported downstream zone simultaneously in the impact zone.[7] is shown in Figure 3, which displays a similar trend downstream of the chute aerator.However, the present values are generally larger than those of Pfister and Hager [7].At a sufficient water depth, the detrained air bubbles cannot escape into the atmosphere through the upper surface; it can only escape into the cavity zone, which may cause the level of Cb to be larger than that of previous studies.In addition, for a shorter length of the impact zone, the fewer measurement sections may cause an insufficient sharp decrease.(4) It can be concluded that the air transportation regularity is different in the impact and equilibrium zones.From the experimental data analysis, bottom air concentration Cb can be expressed by the following equation (Figure 3):

Bottom Air Concentration and Bubble Frequency
( ) where R is the correlation coefficient, (Cb)L is the bottom air concentration at x = L, (Cb)Lm is the bottom air concentration at x = Lm.
The experimental data of bottom bubble frequency fbh0/V0 along the chute are displayed in Figure 4.There are a few reports on the chute bottom bubble frequency of the lower jet downstream of the chute aerator.Note the following: (1) fbh0/V0 decreased in both the impact and equilibrium zones, whereas decay rate of fbh0/V0 differed slightly.For F0 = 7.4, fbh0/V0 = 11.2 at x/L = 1.06, and it decreased to 10.88 at x/L = 1.17 in the impact zone (Figure 4d).In contrast, the bottom bubble frequency decreased slightly more in the equilibrium zone, with fbh0/V0 = 9.44 at x/L = 1.31 decreasing to 6.22 at x/L = 1.91.(2) In contrast, fbh0/V0 decreased slightly in the impact zone, and there was a large variation in air bubble concentration Cb. (3) Furthermore, fbh0/V0 increased with increasing F0.From the experimental data analysis, bottom air concentration C b can be expressed by the following equation (Figure 3): where R is the correlation coefficient, (C b ) L is the bottom air concentration at x = L, (C b ) Lm is the bottom air concentration at x = L m .The experimental data of bottom bubble frequency f b h 0 /V 0 along the chute are displayed in Figure 4.There are a few reports on the chute bottom bubble frequency of the lower jet downstream of the chute aerator.Note the following: (1) f b h 0 /V 0 decreased in both the impact and equilibrium zones, whereas decay rate of f b h 0 /V 0 differed slightly.For F 0 = 7.4, f b h 0 /V 0 = 11.2 at x/L = 1.06, and it decreased to 10.88 at x/L = 1.17 in the impact zone (Figure 4d).In contrast, the bottom bubble frequency decreased slightly more in the equilibrium zone, with f b h 0 /V 0 = 9.44 at x/L = 1.31 decreasing to 6.22 at x/L = 1.91.(2) In contrast, f b h 0 /V 0 decreased slightly in the impact zone, and there was a large variation in air bubble concentration C b .(3) Furthermore, f b h 0 /V 0 increased with increasing F 0 .From the experimental data analysis, bottom air concentration Cb can be expressed by the following equation (Figure 3): ( ) where R is the correlation coefficient, (Cb)L is the bottom air concentration at x = L, (Cb)Lm is the bottom air concentration at x = Lm.The experimental data of bottom bubble frequency fbh0/V0 along the chute are displayed in Figure 4.There are a few reports on the chute bottom bubble frequency of the lower jet downstream of the chute aerator.Note the following: (1) fbh0/V0 decreased in both the impact and equilibrium zones, whereas decay rate of fbh0/V0 differed slightly.For F0 = 7.4, fbh0/V0 = 11.2 at x/L = 1.06, and it decreased to 10.88 at x/L = 1.17 in the impact zone (Figure 4d).In contrast, the bottom bubble frequency decreased slightly more in the equilibrium zone, with fbh0/V0 = 9.44 at x/L = 1.31 decreasing to 6.22 at x/L = 1.91.(2) In contrast, fbh0/V0 decreased slightly in the impact zone, and there was a large variation in air bubble concentration Cb. (3) Furthermore, fbh0/V0 increased with increasing F0.Overall, the present data could be described by where R was the correlation coefficient.

Bottom Bubble Chord Length Distributions
The chord length distribution of the air bubbles shows microscopic variations passing a measurement point in the stream-wise direction.For all the flow conditions, the experimental results exhibited a broad range of the chord length, varying from less than 0.1 mm to more than 30.0 mm.In this study, the chord length interval was 0.1 mm.The chord lengths are not the bubble diameters, but are the characteristic stream-wise bubble sizes.
Garrett et al. [24] and Deane et al. [25] indicated that the bubble chord length distribution Na depended on the dissipation rate ε and bubble diameter dab, and the scaling law followed as where A was a constant, Na = Nt/t, Nt was the number of air bubbles detected in the scan period t (t = 40 s).
The distributions of the bottom bubble chord length in the impact and equilibrium zones are depicted in Figure 5.In the impact zone, the air concentration was within the range of 0.10 < C < 0.90, such that the intermediate region neither corresponded to the high-air concentration region nor to the bubbly flow region.In the equilibrium zone, the maximum air concentration was less than 0.3, so that the zone corresponded to the bubbly flow region.Note the following: (1) the distributions of bottom bubble chord length were skewed with a large proportion of small bubbles.Compared with the impact region, the range of bubble chord length distributions decreased.In the impact zone, the largest probability of bubble chord lengths was concentrated in the range 0.1 ≤ dab ≤ 30 mm, while in the equilibrium zone, the highest probability of bubble chord lengths was concentrated in the range 0.1 ≤ dab ≤ 10 mm.(2) The curve of distribution of the bottom bubble chord sizes was also sharper and narrower along the chute.The power-law scaling of bubble density on diameter is about β = −5/3 for larger than about 1-2 mm bubbles in the impact and equilibrium zones.(3) In the impact zone, β showed a slight increase from −1.2 at x/L = 1.06 to −5/3 at x/L = 1.15 due to the air detrainment.In the equilibrium zone, β showed a slightly increase from -5/3 at x/L = 1.21 to −2.4 at x/L = 1.64 due to the buoyancy effect.Overall, the present data could be described by where R was the correlation coefficient.

Bottom Bubble Chord Length Distributions
The chord length distribution of the air bubbles shows microscopic variations passing a measurement point in the stream-wise direction.For all the flow conditions, the experimental results exhibited a broad range of the chord length, varying from less than 0.1 mm to more than 30.0 mm.In this study, the chord length interval was 0.1 mm.The chord lengths are not the bubble diameters, but are the characteristic stream-wise bubble sizes.
Garrett et al. [24] and Deane et al. [25] indicated that the bubble chord length distribution N a depended on the dissipation rate ε and bubble diameter d ab , and the scaling law followed as where A was a constant, N a = N t /t, N t was the number of air bubbles detected in the scan period t (t = 40 s).
The distributions of the bottom bubble chord length in the impact and equilibrium zones are depicted in Figure 5.In the impact zone, the air concentration was within the range of 0.10 < C < 0.90, such that the intermediate region neither corresponded to the high-air concentration region nor to the bubbly flow region.In the equilibrium zone, the maximum air concentration was less than 0.3, so that the zone corresponded to the bubbly flow region.Note the following: (1) the distributions of bottom bubble chord length were skewed with a large proportion of small bubbles.Compared with the impact region, the range of bubble chord length distributions decreased.In the impact zone, the largest probability of bubble chord lengths was concentrated in the range 0.1 ≤ d ab ≤ 30 mm, while in the equilibrium zone, the highest probability of bubble chord lengths was concentrated in the range 0.1 ≤ d ab ≤ 10 mm.(2) The curve of distribution of the bottom bubble chord sizes was also sharper and narrower along the chute.The power-law scaling of bubble density on diameter is about β = −5/3 for larger than about 1-2 mm bubbles in the impact and equilibrium zones.(3) In the impact zone, β showed a slight increase from −1.2 at x/L = 1.06 to −5/3 at x/L = 1.15 due to the air detrainment.In the equilibrium zone, β showed a slightly increase from -5/3 at x/L = 1.21 to −2.4 at x/L = 1.64 due to the buoyancy effect.

Bottom Turbulent Dissipation Rate
Xu et al. [22,23] and Brujan et al. proved that the size of air bubbles near the bottom was a significant parameter for the collapse process of cavitation bubbles.There was little information to investigate the air bubble variations on the chute bottom.The equilibrium zone belonged to the bubbly flow region and the dynamic behavior of the flow was considered to be a steady state air-water two-phase flow.
Hinze [26] indicated that any fluid gas bubbles or droplets were likely to fragment if the pressure forced on its surface exceeds the restoring forces associated with surface tension.
where Wc = the critical Weber numbers, u = the turbulent velocity field on the scale of the bubble.If the fluctuating velocity field was assumed to be described by Kolmogorov's inertial subrange, so that u 2 = 2ε 2/3 dab 2/3 , then only bubbles with diameters larger than the Hinze scale would fragment [17][18][19].
where dc = the Hinze scale, ε = turbulent dissipation rate, which was estimated rather than measured directly.Recent experiments suggest that Wc lied in the range 0.585-4.7,and we have used the value Wc = 1.0 (Table 3).For this to exceed a critical value, Wc required that d > dc.

Bottom Turbulent Dissipation Rate
Xu et al. [22,23] and Brujan et al. proved that the size of air bubbles near the bottom was a significant parameter for the collapse process of cavitation bubbles.There was little information to investigate the air bubble variations on the chute bottom.The equilibrium zone belonged to the bubbly flow region and the dynamic behavior of the flow was considered to be a steady state air-water two-phase flow.
Hinze [26] indicated that any fluid gas bubbles or droplets were likely to fragment if the pressure forced on its surface exceeds the restoring forces associated with surface tension.
where W c = the critical Weber numbers, u = the turbulent velocity field on the scale of the bubble.If the fluctuating velocity field was assumed to be described by Kolmogorov's inertial subrange, so that u 2 = 2ε 2/3 d ab 2/3 , then only bubbles with diameters larger than the Hinze scale would fragment [17][18][19].
where d c = the Hinze scale, ε = turbulent dissipation rate, which was estimated rather than measured directly.Recent experiments suggest that W c lied in the range 0.585-4.7,and we have used the value W c = 1.0 (Table 3).For this to exceed a critical value, W c required that d > d c .
Water 2018, 10, 1129 9 of 11 If the Hinze scale was known, the turbulent energy dissipation rate ε on the chute bottom could be estimated from Equation (7).Hinze [26] used the definition that 95% of the air is contained in bubbles with a diameter less than d c and the calculated turbulent energy dissipation rate ε for Deane's [18] data was much smaller.Garrett et al. [24] used the results of Martínez-Bazán et al. [30,31] and established a one-to-one relationship between the local dissipation rate and bubble size.
The bottom dissipation rate ε calculated from Equation ( 7) is shown in Figure 6 and tends to increase in the measured distance.From Equation ( 6), a change in dissipation rate ε from 0.10 to 0.75 m 2 s −3 corresponded to a change in Hinze scale from 5.5 to 2.4 mm.Thus, the estimate of 2.4 mm for the Hinze scale was in good agreement with the observed change in power-law scaling of the bubble distribution at around 1-2 mm (Figure 5b).Because bubbles smaller than the Hinze scale are stabilized by surface tension forces, a reasonable estimate for this scale is the diameter of the smallest bubbles observed to fragment.In the present study, the bubble chord lengths were smaller than the bubble diameters, due to the measurement technique utilizing conductivity probes.A challenge remains to observe the bubble diameter in a model facility, or even a prototype facility.If the Hinze scale was known, the turbulent energy dissipation rate ε on the chute bottom could be estimated from Equation (7).Hinze [26] used the definition that 95% of the air is contained in bubbles with a diameter less than dc and the calculated turbulent energy dissipation rate ε for Deane's [18] data was much smaller.Garrett et al. [24] used the results of Martínez-Bazán et al. [30,31] and established a one-to-one relationship between the local dissipation rate and bubble size.
The bottom dissipation rate ε calculated from Equation ( 7) is shown in Figure 6 and tends to increase in the measured distance.From Equation (6), a change in dissipation rate ε from 0.10 to 0.75 m 2 s −3 corresponded to a change in Hinze scale from 5.5 to 2.4 mm.Thus, the estimate of 2.4 mm for the Hinze scale was in good agreement with the observed change in power-law scaling of the bubble distribution at around 1-2 mm (Figure 5b).Because bubbles smaller than the Hinze scale are stabilized by surface tension forces, a reasonable estimate for this scale is the diameter of the smallest bubbles observed to fragment.In the present study, the bubble chord lengths were smaller than the bubble diameters, due to the measurement technique utilizing conductivity probes.A challenge remains to observe the bubble diameter in a model facility, or even a prototype facility.

Conclusions
This study presented new information on the bottom air bubble characteristics in the downstream of chute offset-aerators, which were measured with an accurate measurement technique utilizing conductivity probes.For eliminating the upper aeration, the bottom concentration, bubble frequency, bubble chord length, and bottom dissipation rate were analyzed along the chute longitudinal direction.The present study complemented the existing studies on the bubble properties of chute aerator flows.
The bottom air concentration was found to decrease sharply in the impact zone and to decrease slightly in the equilibrium zone.The bottom air bubble frequency was found to decrease slightly along the chute bottom.A formula to predict the bottom bubble frequency in the impact and equilibrium zones was proposed.
Observations showed that the range of probability of bottom bubble chord lengths tended to decrease sharply in the impact zone and to a lesser extent in the equilibrium zone.The bubbles larger than the Hinze scale were subject to fragmentation and showed −5/3 power-law scaling with diameter.Using the relationship between the local dissipation rate and bubble size [24], the bottom dissipation rate was found to increase along the chute bottom.The Hinze scale estimated from Equation ( 6) showed a good agreement with the observed in power-law scaling of the bubble distribution.

Figure 1 .
Figure 1.Definition sketch with relevant parameters (a) and side view of chute aerator flow with high speed camera (b).

Figure 1 .
Figure 1.Definition sketch with relevant parameters (a) and side view of chute aerator flow with high speed camera (b).

4. 1 .
Bottom Air Concentration and Bubble Frequency Bottom air concentration C b was defined experimentally using the values measured closest to the model chute bottom (z = 1-1.5 mm).The experimental data of bottom air concentration C b is presented in Figure 3.Note the following: (1) the decay rate of C b was different in the impact (L < x < L m ) and equilibrium zones (L m ≤ x < L D ); C b decreased sharply in the impact zone and decreased slightly in the equilibrium zone.In the impact zone, bottom air concentration C b = 0.64 at x/L = 1.06 and it decreased to 0.18 at x/L = 1.16.In contrast, the bottom air concentration decreased slightly in the equilibrium zone, with concentration C b = 0.145 at x/L = 1.31 decreasing to 0.064 at x/L = 1.91.(2) With the increase in F 0 , the values of C b increase in the present experiments, in agreement with the results of Pfister and Hager [7].(3) The comparison of the measured values of C b with the values from the formula proposed by Pfister and Hager Bottom air concentration Cb was defined experimentally using the values measured closest to the model chute bottom (z = 1-1.5 mm).The experimental data of bottom air concentration Cb is presented in Figure 3.Note the following: (1) the decay rate of Cb was different in the impact (L < x < Lm) and equilibrium zones (Lm ≤ x < LD); Cb decreased sharply in the impact zone and decreased slightly in the equilibrium zone.In the impact zone, bottom air concentration Cb = 0.64 at x/L = 1.06 and it decreased to 0.18 at x/L = 1.16.In contrast, the bottom air concentration decreased slightly in the equilibrium zone, with concentration Cb = 0.145 at x/L = 1.31 decreasing to 0.064 at x/L = 1.91.(2) With the increase in F0, the values of Cb increase in the present experiments, in agreement with the results of Pfister and Hager [7].(3) The comparison of the measured values of Cb with the values from the formula proposed by Pfister and Hager

Figure 6 .
Figure 6.Calculated values of turbulent energy dissipation rate ε on the chute bottom.Figure 6. Calculated values of turbulent energy dissipation rate ε on the chute bottom.

Figure 6 .
Figure 6.Calculated values of turbulent energy dissipation rate ε on the chute bottom.Figure 6. Calculated values of turbulent energy dissipation rate ε on the chute bottom.

Table 2 .
Experimental conditions of chute aerator flow.

Table 2 .
Experimental conditions of chute aerator flow.

Table 3 .
Critical Weber numbers for the splitting of air bubbles in water flows.

Table 3 .
Critical Weber numbers for the splitting of air bubbles in water flows.