Abstract

To explore the rectification parameters of the diversion piers optimized for the forebay of the pump station with a lateral angle of 45°, the orthogonal experiment and computational fluid dynamics methods are used to analyze the flow characteristics of the diversion piers under different parameter combinations. The flow pattern in the forebay of the side water inlet is improved. The rectification effect of the diversion piers under 16 schemes is analyzed, considering the length, width, radian, and relative height of the diversion piers. Combined with numerical simulation, a better rectification scheme is provided, and finally, a reasonable range of values for the rectification parameters of the forebay diversion pier of the side 45° bend angle pump station is obtained.

1. Introduction

The forebay of the pumping station is an important water diversion structure, and its hydraulic characteristics directly affect the performance of the pump. The forebay can be divided into frontal water intake and lateral water intake [14]. Some pumping stations can only adopt the form of lateral water intake because the layout of pumping stations is often restricted by objective factors, such as water flow characteristics and geographic location. If the water flow from the lateral inflow is oblique to the axis of the inlet pool, then backflow or vortex easily occurs. A number of studies have been carried out to analyze the flow pattern in the forebay and pump. At present, scholars mainly carried out a combination of model test and numerical simulation [58]. At the same time, substantial practical engineering problems have been solved using mathematical methods [912]. With the development of computer technology, the accuracy of numerical simulation is gradually enhanced and has become the main research method [1316]. Sweeney [17] used narrow water distribution holes and splitter plate rectification measures to improve the bad flow pattern of the forebay of the lateral diversion pump station and obtained a stable forebay inlet flow pattern. The water inlet conditions of the pump are satisfied, and the hydraulic characteristics are significantly improved. Spence [18] used numerical simulation and model tests to optimize a series of parameters to explore the possible adverse effects of centrifugal pump operation. Various rectification measures are adopted to improve the flow pattern, which is beneficial to better manage the operation of the pump station. Zhou et al. [19] used the pump station model test to analyze the causes of the bad flow pattern in the forebay and proposed corresponding measures to improve the flow pattern. The first domestic research on the vortex at the inlet of the pump inlet channel was based on the Froude model to simulate the prototype flow pattern by increasing the flow velocity by approximately 0.5 times. The results are relatively satisfactory. Cheng et al. [20] quantitatively analyzed the influence of the position, height, angle, and length of the Y-shaped diversion pier on the flow pattern of the pump station’s forebay by using Fluent software and obtained the optimal parameters of the Y-shaped diversion pier. Orthogonal experiment was used to analyze multifactor and multilevel experiments. Compared with the full-factor method, the number of experiments using orthogonal experiment was greatly reduced. Using the standardized orthogonal table to design the test plan and carry out the range and variance analysis of the test results, the influence law of each factor on the index can be obtained. Then, the optimized production conditions or process conditions can be obtained. Zhou et al. [21] used the orthogonal test method, selected four-factor three-level orthogonal table, analyzed the orthogonal test data, put forward an optimized design model pump scheme, and obtained the influence law of each geometric parameter on the index. Shi et al. [2224] used a five-factor four-level orthogonal table to design the test plan for the relevant control dimensions of the pump device. Computational fluid dynamics (CFD) numerical simulation was used to calculate each plan, analyze the influence of each parameter on the hydraulic performance of the flow channel, and then verify it with a model test. Finally, the priority order of each control parameter to the index and the optimal plan were obtained. In the present study, numerical simulation is carried out for the forebay of the pumping station with a 45° lateral bending angle. The orthogonal test method is used to design the geometric parameters of the diversion pier, and the optimized parameter combination is obtained. Then, the CFX numerical calculation is used for analysis and verification. The improvement law of the parameters of the diversion pier on the flow pattern in the forebay of the lateral inlet is obtained.

2. Study Area and Method

2.1. Study Area

A pumping station hub with lateral water intake is considered an example. Three pump units are installed in the pump room; each pump has a design flow of 4.0 m³/s, and the flow channel inlet size is 2.36 m (width) × 1.8 m (length). The intake level of the pumping station is 1.9 m, the elevation of the bottom of the drainage canal is 0.0 m, and the elevation of the forebay bottom is −1.0 m. The main functions of the project are flood control and drainage. The left side of the pumping station is a self-draining culvert. When the pumping station is in operation, the self-draining culvert control gate is closed. The total width of the outlet of the forebay is B, and the length from the center of the inlet of the forebay to the center of the outlet is L. Figure 1 shows the plan layout of the intake building of the pumping station, where the lateral bend angle is 45°.

2.2. Study Method

CFX is used to numerically simulate the model, and the overall calculation domain is divided into unstructured grids. The hydraulic loss is used as the evaluation index. After grid independence verification, when the total grid number is 3.58 million, the hydraulic loss is relatively small. Therefore, the influence of the number of grids on the numerical calculation can be ignored.

The inlet boundary condition adopts the quality inlet, the inlet flow channel outlet is selected as the calculated outlet boundary, the standard atmosphere is set to 1, and the “steel cover assumption” is selected for the surface of the water body, that is, symmetry surface boundary treatment and nonslip walls.

3. Orthogonal Optimization Design

3.1. Selection of Test Factors and Levels

Diversion piers are used to improve the flow pattern of the forebay. The geometric parameters of the diversion piers have a certain influence on the inlet channel and the water inlet conditions. The axial uniformity and weighted average angle of the inlet section of the inlet flow passage are used as the quantitative indicators of the test, and the geometric parameters of the diversion pier greatly affect the indicators, which need to be integrated for orthogonal test design. The results of the orthogonal test are analyzed, the optimal parameter combination of the diversion pier is determined, and then the influence of various geometric parameters of the diversion pier on the hydraulic performance of the forebay is studied. The key to this research is to reasonably select the range of geometric parameters of the diversion pier.

For the forebay with different lateral bends, the length of the diversion piers is mainly selected according to the range of its bad flow patterns. The maximum horizontal width of the diversion pier is 1.0 m because the distance between 1# and 2# inlet channels is 1.0 m. The arc-shaped diversion pier can play a good diversion function, and the range of the diversion pier’s radian should be selected through repeated trials. The relative height of the diversion pier is used to study the influence of the height of the diversion pier on the hydraulic characteristics of the flow pattern in the forebay, considering that different water levels may occur in the forebay. The schematic of the geometric parameters of the diversion pier is shown in Figure 2, where l is the length of the diversion pier, b is the width of the diversion pier, α is the radian of the diversion pier, and h is the relative height of the diversion pier. The expression is as follows: h = H1/H2. In the formula, H1 is the actual height of the diversion pier, m; H2 is the water depth in the forebay, m.

3.2. Orthogonal Table

In this study, the parameters of the orthogonal experiment are the length, width, radian, and relative height of the diversion pier. The factors are determined according to the design scope. The length of the diversion pier is based on the length of the swing zone in the forebay in the original design. The level of the test factor is determined according to the possible variation range of each factor parameter in the actual project. When the lateral angle is 45°, the orthogonal design of the geometric parameters of the diversion pier in the forebay selects four factors, and each factor selects four levels. The specific parameters of the factor level are shown in Tables 1 and 2. The factor codes A, B, C, and D in the table represent the length (m), width (m), radian (°), and relative height of the diversion pier, respectively.

According to the selected factors and the number of levels, the orthogonal table L16 (45) is selected for orthogonal design, and one column is the error column.

3.3. Orthogonal Test Results of the Geometric Parameters of the Diversion Pier

The preliminary research analysis shows the position of the bad flow pattern in the forebay, the reason for its formation, and the approximate range of the recirculation zone when the lateral angle is 45°. The bad flow pattern mainly affects the 1# inlet flow passage and causes harm to the inlet conditions of the 1# inlet flow passage, thereby threatening the normal stability of the pump unit.

Therefore, the engineering measures of adding diversion piers between the 1# and 2# inlet channels are considered to improve the flow pattern of the forebay, and the influence of the geometric parameters of the diversion piers on the hydraulic performance of the forebay is studied. The four geometric parameters and levels of the selected diversion piers are designed according to the orthogonal table L16 (45), and each plan is numerically simulated by CFD. Moreover, the axis of the inlet section of the 1# inlet channel under each plan is obtained. The uniformity and weighted average angle of flow velocity, test plan, and calculation results are shown in Table 2. Intuitive analysis and variance analysis of the orthogonal test results are carried out, and the influence trend of various factors on the improvement of the inlet conditions of the inlet channel is analyzed. Then, the parameter combination of the optimization of the geometric parameters of the diversion pier is obtained. Finally, the optimized parameter combination obtained is analyzed and verified through numerical calculation to obtain the improvement law of the flow characteristics of the forebay by the geometric parameters of the diversion pier.

3.4. Visual Analysis

Range refers to the difference between the maximum value and the minimum value of the sum of each level index in the same factor, that is, in any column. The value reflects the degree of influence of this factor on the index, that is, the degree of influence of the parameter on the uniformity of the inlet flow velocity of the 1# inlet flow channel. The larger value indicates greater influence on the uniformity of the flow velocity at the inlet of the inlet flow channel. Table 2 shows that the difference between the parameters and the uniformity of axial flow velocity is D> B > A > C. The order of the factors is as follows: the relative height of the diversion pier, the width of the diversion pier, the length of the diversion pier, and the radian of the diversion pier. The research shows that, among the four factors that affect the rectification effect of the diversion pier, the relative height of the diversion pier is the key factor, the length of the diversion pier, the width of the diversion pier are general factors, and the radian of the diversion pier is the secondary factor.

A trend chart is drawn to analyze the specific influence trend of the changes in various factors on improving the uniformity of the flow velocity of the inlet section of the inlet flow channel. The abscissa is the level of each factor, and the ordinate is the average of the sum of the results of the corresponding factors. The relationship among the ki value, the factors, and the uniformity of the flow velocity of the inlet section of the inlet flow channel is shown in Figure 3.

Figure 3 shows that, within the range of the number of factors studied, various factors have different effects on the uniformity of the axial flow velocity of the inlet section of the 1# inlet channel. When the length of the diversion pier is 8 m, the index is the smallest, and the index increases below 8 m and above 8 m. The 20° radian of the diversion pier leads to the smallest index, and the index is lower than 20° and higher than 20°. As the width of the diversion pier increases, the index gradually increases, forming a trend that the index increases when the width of the diversion pier continues to increase. When the relative height of the diversion pier changes from 0.4 to 1.0, the index also increases continuously. This finding indicates that the index reaches the maximum value when the height of the diversion pier is flush with the water surface.

3.5. Analysis of Variance

Variance analysis should be performed on the test results to estimate the size of the error in the test process and analyze whether the influence of various factors on the uniformity of the axial flow velocity at the inlet section of the 1# inlet flow passage is significant. The orthogonal table leaves an error column to calculate the experimental error. The sum of squared deviations excludes the difference between the factor levels and only represents the size of the experimental error because this column is blank. The results of variance calculation are shown in Table 3. The results show that the variance of the error term is smaller than the variance of each factor. This finding indicates that each factor affects the test index, and no factor is included in the error term. In addition, if the variance and error of the radian of the diversion pier are small, then the change in this factor level has no significant impact on the test results. In addition, the variance and error of other factors are large. Thus, the change in the factor level affects the test results, which have a significant impact. The significant degree of the specific influence of each factor can be judged on the basis of the comparison of the variance ratio and the critical value of F.

Table 3 shows that the change in the length of the diversion pier has a certain influence on the index, that is, the uniformity of the axial flow velocity of the inlet section of the 1# inlet channel. The change in the width of the diversion pier significantly affects the index. The change in the radian level of the diversion pier has no significant impact on the index, and the parameters can be arbitrarily selected within the test range according to the specific situation. The change in the relative height of the diversion pier has a very significant impact on the index.

3.6. Determination of Optimal Solution

For the uniformity of axial flow velocity, the larger value is better. General calculation and analysis are performed on the orthogonal test results and the trend graph of the relationship between various factors and indicators, and then, a better test plan is obtained. The maximum average value ki of the axial flow velocity uniformity of the inlet section of the inlet channel of 1# is selected under the level i of the diversion pier. Figure 3 shows that, among the four levels considered by the length of the diversion pier, the ki value is the largest when the length is 7 m. Among the four levels considered by the width of the diversion pier, the ki value is the largest when the width is 1.0 m. Among the four levels considered by the radian of the diversion pier, the ki value is the largest when the radian is 25°. Among the four levels considered by the relative height of the diversion pier, the ki value is the largest when the relative height is 1.0 m. Therefore, the optimal plan is as follows: the length of the diversion pier is 7 m; the width of the diversion pier is 1.0 m; the radian of the diversion pier is 25°; the relative height of the diversion pier is 1.0. Therefore, the design scheme is plan A1B4C4D4.

4. Analysis of Numerical Simulation Results Based on Orthogonal Experimental Design

4.1. Feature Section Selection

Two characteristic sections of the horizontal section of the front pond are selected for flow field analysis, namely, two horizontal longitudinal sections of the front pond surface layer and bottom layer. Figure 4 shows the horizontal section, namely, the surface layer a-a and the bottom layer b-b to study the flow pattern and axial velocity distribution of the forebay surface layer and bottom layer, respectively. Among them, the surface layer a-a section is 0.05 m away from the water surface, and the coordinates of the surface layer and bottom layer are Za = 1.85 m and Zb = 0.05 m, respectively.

4.2. Selection of the Flow Pattern Scheme for the Horizontal Profile of the Forebay

The optimal combination of diversion pier parameters is obtained in accordance with the calculation of orthogonal experiment. According to this parameter setting, several schemes are numerically simulated again, and the data of the scheme with the highest axial flow velocity uniformity among the 16 schemes are compared. The 16 schemes in the orthogonal experimental design are named L1L16. The program parameter settings are shown in Table 4.

A representative scheme is selected for numerical simulation analysis to study the flow pattern of the water in the forebay. The accuracy of the orthogonal experiment results is analyzed and verified to further confirm the effect of the orthogonal experiment, and the selected typical schemes are shown in Table 5.

4.3. Analysis of Numerical Simulation Results of Each Typical Scheme

Figure 5(a) and 5(b) show the horizontal cross-sectional flow field diagrams of the surface layer and bottom layer of Scheme 4. The flow field diagram shows that the flow pattern at the bottom of the forebay is the same as that of the surface layer and the flow pattern is smooth. The diversion piers under the scheme can effectively improve the flow pattern of the side forebay. Figures 5(c) and 5(d) are the horizontal cross-sectional flow field diagrams of the surface layer and bottom layer of Scheme 7, and a low-velocity zone is formed on the left side wall of the forebay. A large-scale recirculation zone appears in the low-velocity zone of the bottom layer, and the recirculation center position is approximately at X −6 m. From the bottom layer to the surface layer, the recirculation center position gradually moves in the downstream direction. The range of the swirling zone gradually decreases, the surface layer forms a small-scale recirculation zone, and the recirculation center position is at X of −2 m. The main flow line is skewed due to the large radian of the diversion pier, thereby affecting the water inlet conditions of the 1# inlet channel. Figures 5(e) and 5(f) show the horizontal cross-section flow field diagrams of the surface layer and bottom layer of Scheme 11. The water flow velocity in the middle of the forebay is relatively large with a tendency to decrease on both sides. Moreover, a large-scale recirculation zone appears on the bottom layer. The center position of the recirculation is at X −5 m. From the bottom layer to the surface layer, the center of the recirculation gradually moves in the direction of the downstream flow, and the range of the swirling zone gradually decreases to form a small-scale recirculation zone on the surface layer. Comparing plan 11 with plan 7, the return center position is at X −2 m, and the water flow in front of the 1# inlet channel of plan 11 is smooth. Figures 5(g) and 5(h) show the horizontal cross-sectional flow field diagrams of the surface layer and bottom layer of Scheme 15. A large-scale recirculation zone is found in the bottom layer, and the recirculation center position is approximately at X −6 m. Scheme 15 has a slightly smaller bottom recirculation zone than scheme 11. From the bottom layer to the surface layer, the center position of the recirculation gradually moves in the downstream direction, and the scope of the recirculation zone gradually decreases, forming a small-scale recirculation zone on the surface layer. The center position of the surface layer recirculation is at −1.5 m. The recirculation zone formed at the bottom layer is far from the inlet of the 1# inlet channel. Thus, it slightly affects the flow pattern of the water in front of the 1# inlet channel.

The axial flow field diagrams of the surface layer and bottom layer horizontal section of each typical scheme are compared and analyzed, as shown in Figure 5. In the L4 scheme of Figures 5(a) and 5(b), the surface layer and bottom layer have the smoothest flow lines and the best flow pattern, which is the optimal scheme. For the excellent scheme directly observed in the orthogonal experiment, the design scheme is observed when the index is the largest, that is, the fourth scheme A1B4C4D4. The maximum flow rate uniformity is 90.75%, which is 7.29% higher than the original scheme. This conclusion shows that the results of orthogonal experimental design and numerical simulation analysis are consistent.

5. Conclusion

The rectification measures for the forebay of the pump station with a 45° lateral bending angle are studied. The orthogonal test method is used to design the various geometric parameters of the diversion pier, and the optimized parameter combination is obtained by analyzing the orthogonal test results. Then, CFX numerical calculation is used to verify the analysis. Conclusions are as follows:(1)Adding diversion piers between the 1# and 2# inlet channels can effectively improve the flow pattern of the forebay. Through the orthogonal experiment, the different geometric parameters of the diversion pier lead to a large difference in the uniformity of the flow velocity at the inlet section of the 1# inlet channel.(2)Through the analysis, the range indicates that the factors from primary to secondary are DBAC. That is, among the four factors that affect the rectification effect of the diversion pier, the relative height of the diversion pier is the key factor, the length of the diversion pier and the width of the diversion pier are the general factors, and the radian of the diversion pier is the secondary factor.(3)Various factors have different effects on the uniformity of the axial flow velocity at the inlet section of the inlet flow channel. The length of the diversion pier at 8 m leads to the smallest index, and the indexes below 8 m and above 8 m are large. Similarly, the diversion pier radian of 20° leads to the smallest index, and the indexes below 20° and higher than 20° are large. As the width of the diversion pier increases, the index gradually increases, forming a trend that when the width of the diversion pier continues to increase, the index also increases. Similarly, when the relative height of the diversion pier is from 0.4 to 1.0, the index also increases, that is, when the height of the diversion pier is flushed with the water surface, the index reaches the maximum value.

Data Availability

The curve data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.