1. Introduction
Recently, with the increase in the usage of disposable masks because of the COVID (Corona virus disease)-19 pandemic, used masks are being commonly discarded in toilets. If the cloth wastes such as a disposable mask or large waste such as a baby diaper is disposed into a toilet, the flow path of the pump that transports wastewater is blocked, and as a result, the function of the wastewater transport system is lost. Therefore, the demand for special pumps for wastewater transportation is increasing, and it has gained attention as an industry with the potential for future growth.
As an example of the special pumps, grinder and vortex pumps are widely used for transporting wastewater. However, these special pumps have low efficiency and high maintenance costs, contributing to large operating costs. Several studies explored the treatment of sewage containing solid waste to solve these problems [
1,
2,
3,
4]. Lu et al. [
1] studied the hydraulic performance and pressure fluctuation characteristics of a grinder pump. Through a numerical analysis, the hydraulic performances of the pump when the flow path is in a clogging state and in a normal operating state were compared, and steady and unsteady Reynolds-averaged Navier-Stokes (RANS) analyses was performed. It was found that, as the degree of clogging of the grinder cutter increases, the total head of the pump declines parabolically, with the best efficiency point shifting to the low flow rate region and the high efficiency area narrowing. Schivley and Dussourd analyzed and designed a vortex pump, using a one-dimensional analytical model [
2]. They compared the calculated performance parameters with those measured by using the laboratory model, and they computed the overall hydraulic characteristics of the pump and compared these characteristics with those of many test pumps. They improved the theoretical formula used in the preliminary design of a vortex pump. Ohba et al. [
3] theoretically analyzed the flow characteristics inside a vortex pump and secured the theoretical reliability by comparing the predicted values with experimental results. Their theoretical formula could predict not only the hydraulic performance of the vortex pump, but also the velocity component inside the pump. The obtained results were in good agreement with the experimentally measured values. Steinmann et al. [
4] analyzed the internal flow of a vortex pump through numerical analysis and experiments to investigate the unsteady cavitating flow of a vortex pump. The Rayleigh–Plesset cavitation model was used to simulate cavitation under the overload condition of the vortex pump, and an acrylic glass window was installed in the experimental apparatus, to observe this phenomenon. Under the best efficiency point condition, the numerically derived head and shaft power of the pump were about 6% higher than the experimentally measured values. Conversely, cavitation under overload was observed more in experimental results than in numerical analysis.
In addition to the grinder and vortex pumps mentioned above, single-channel pumps designed to transport relatively larger solid wastes were actively studied [
5,
6,
7]. The single-channel impeller with a single free annulus passage can smoothly transfer sewage-containing solid wastes. However, this impeller has an unsymmetrical annulus flow passage, and it is difficult to stabilize the unsteady flow-induced vibration, due to the interaction between the rotating impeller and stationary volute [
5]. Shi and Tsukamoto [
6] numerically analyzed the pressure fluctuation due to the impeller–diffuser interaction in a diffuser pump. They confirmed that the flow characteristics due to this interaction can be analyzed through an unsteady flow analysis. Feng et al. [
7] analyzed the unsteady flow characteristics between the impeller and diffuser of a radial pump by unsteady RANS (URANS) analysis and laser Doppler velocimetry (LDV); they identified two types of rotor–stator interaction effects. One is the downstream effect induced by the impeller, which has an unsteady flow characteristic because of the highly distorted flow field and the wake of the impeller. The other is the upstream effect induced by the stator, which causes unsteady pressure and velocity fluctuations. Such a single-channel pump has the advantage of being able to transport relatively large waste solids, but it has the disadvantage that the fluid-induced vibration is greater than that in the existing special pumps (e.g., grinder and vortex pumps) because of its asymmetric structural characteristics.
The special pumps for wastewater treatment introduced so far clearly have advantages and disadvantages, depending on their type. Grinder and vortex pumps have low vibration during operation, but the size of transportable solid matters is relatively small, and the maintenance costs are relatively high. On the other hand, a single-channel pump can easily transport large solid matters, but in some cases, relatively severe vibration occurs. To solve these problems, the authors intended to design a two-vane pump in this study. The impeller of the two-vane pump is composed of two blades that are symmetrical in the rotational axis. Therefore, it is structurally simple compared to the grinder and vortex pumps, and the relatively large flow path of the impeller has a small number of blades, so waste solids can be transferred smoothly. In addition, due to the symmetrical impeller geometry, the fluid-induced vibration is relatively less than that of a single-channel pump.
In the past, design optimization using numerical analysis has been widely used for fluid-based machinery [
8,
9,
10,
11]. For example, Lee et al. [
10] performed an optimization to improve the efficiency of a low-speed axial flow fan, using a gradient-based search algorithm. Lee and Kim [
11] optimized axial compressor blades to improve the efficiency, using numerical optimization techniques such as conjugate direction methods and the golden section method, combined with a three-dimensional (3D) thin-layer Navier-Stokes solver. Recently, with the rapid advances in computing resources, optimal design is being actively researched based on machine-learning techniques [
12,
13].
In this study, a two-vane pump was designed to develop a series of special pumps for transporting wastewater. Compared to a single-channel pump, this pump can transport relatively smaller waste solids but experiences less vibration during operation because of the axial symmetry of the impeller. Considering the characteristics of this pump, a three-objective optimization design based on machine learning was performed to maximize the size of the waste solids that can be transported while simultaneously improving the hydraulic performances of the pump. For the three-objective optimization, the geometric design variables, i.e., inlet and outlet blade angles and cross-sectional area of the volute, were chosen, and the volume of the waste solids and the head and efficiency of the pump were considered as objective functions. Their relationship was predicted by using an artificial neural network (ANN) [
14] based on machine learning. Then, the genetic algorithm (GA) [
15] was used to find the optimal solutions, and the Pareto-optimal front surface [
16] was derived as the final optimization result.
3. Optimization Techniques
The three-objective optimization problem was defined as follows:
Maximize: F(x) = [F1(x), F2(x), F3(x)]
Design variable bound: LB ≤ x ≤ UB, x∈R,
where
F(
x) is the vector of real-valued objective functions;
x is the vector of the design variables;
LB and
UB denote the vectors of the lower and upper bounds, respectively; and
R is a real number [
24].
Figure 4 shows the procedure of a typical optimization design. Once the operating conditions of the design target are determined, the type of turbomachinery and airfoil (or hydrofoil) are selected. Then a preliminary design is performed through one-dimensional mean-line analysis, and the initial blade topology is derived. Subsequently, the optimization design is performed for improving the performance of the turbomachinery.
In this study, the optimization design was based on machine learning. First, the objective functions and constraints were defined according to the design goal. Subsequently, the design variables and their ranges were chosen. Thereafter, a database of the correlation between the design variables and the objective functions was established within the design space, using the design of experiment (DOE). In the next step, a predictive model was constructed by using ANN [
14] to correlate the design variables and objective functions of the two-vane pump. In this step, the predictive model was trained by using machine learning. This process is described in detail in
Section 3.2. The GA [
15] was used to find the optimal design solution, considering the correlation of each objective function in the constructed prediction model. The GA is a well-known stochastic searching algorithm based on the mechanism of natural selection, genetics, and evolution. It evaluates various points in the design space and can be applied to find a global solution to any given problem. This algorithm proceeds as shown in
Figure 5. Finally, the Pareto-optimal front surface was derived by using MATLAB R2018b [
25], and the optimization procedure was terminated.
3.1. Design Variables and Objective Functions
In order to maximize the hydraulic performances and size of waste solids, the efficiency (
η) and total head (
Ht) of the pump and the volume of waste solids (
Vs) were selected as the objective functions:
where
P,
ρ,
g,
Q,
τ, and
ω are the total pressure, density of the working fluid, acceleration due to gravity, flow rate, torque, and angular velocity, respectively. Further,
Ds in Equation (6) is defined as the distance between the leading edge of the blade and the other blade at the impeller inlet, as shown in
Figure 6.
Figure 7 shows the design variables considered in this work. Their ranges are listed in
Table 1. The inlet and outlet blade angles and the cross-sectional area of the volute were chosen for the optimization. They are defined as shown in
Figure 7. As mentioned earlier, the blade angle distribution is the same in the span direction from the hub to the shroud and is defined by using the fourth-order Bézier curve [
26] in the streamwise direction, as shown in
Figure 7a. In order to adjust the inlet (
β1) and outlet (
β2) blade angles, two control points (CP
1 and CP
2) were fixed. The distribution of the cross-sectional area of the volute was changed linearly in the circumferential direction, as shown in
Figure 7b.
3.2. Surrogate Modeling Based on Machine Learning
In the present optimization design, supervised machine learning was adopted, considering the reasonable computation costs. This optimization technique is a very efficient approach for optimizing a system without analytical representation, fitting a surrogate model.
In order to construct the input data, the central composite design was used as the DOE. In all, 15 experimental points were arranged for three design variables, and the objective function values were derived through RANS analysis at each design point. These values are listed in
Table 2. In this table, the values of the design variables and objective functions are normalized by dividing them by the corresponding reference value.
A radial basis neural network (RBNN) model [
27], which is a type of ANN, was used in this optimization study. The RBNN model has a hidden layer of radial neurons and an output layer of linear neurons, as shown in
Figure 8. The hidden layer uses a series of radial primitives to nonlinearly modify the input space to the intermediate space. The output of the hidden layers then executes a linear combiner to produce the desired targets [
27]:
where
wj is the weight, and
φj is radial basis function, which uses the Gaussian function. Several parameters are needed to construct a surrogate model: the input weight and the center and width of a unit for the Gaussian function. In the present RBNN model, the input weights were chosen by machine learning. The input vector with the worst performance was chosen as the center of a new hidden-layer Gaussian function [
25]. Then the RBNN only needed to determine the width of the Gaussian function (spread constant). The network training was performed by adjusting the cross-validation error (CV) by changing the spread constant (SC), as shown in
Figure 9. SC
1, SC
2, and SC
3 correspond to the total head, efficiency, and size of waste solids, respectively. The SC values for each objective function were chosen by a K-fold CV test [
28].
The data sample
was partitioned into K disjoint subsets (K-fold CV), as shown in
Figure 10. Of these, (K-1) folds were used to train the RBNN network, and the last fold (the Kth set) was used for evaluation. This process was repeated K times, using a different fold for evaluation each time. The network training was performed by adjusting the CV error by changing SC. The CV error at a particular SC value was calculated as follows:
where
is the prediction error for the Kth set. The predicted values
were determined by using the constructed RBNN model from the sample points in the (K-1) subsets. In the present study, K was set as 15. According to the results of the K-fold CV test, the final SC values SC
1, SC
2, and SC
3 were set as 0.3, 1.9, and 9.9, respectively.
5. Discussion
The size of waste solids is directly related to the blade inlet angle. As the blade inlet angle increases, the size of waste solids increases. The flow passage areas of the AODs in the meridional direction are compared in
Figure 14. Their values are normalized by using the maximum area of AOD 1 in the meridional direction. It was confirmed that the inlet passage area of AOD 1 with the greatest increase in the blade inlet angle was the largest, and the inlet passage area of AOD 5 with the reduced inlet blade angle was the smallest among AODs.
AOD 3, which satisfies the design target, was selected for further analysis of the optimum design. The 3D geometry for AOD 3 is shown in
Figure 15. Through optimization, an inlet blade angle of about 1.289 times the reference design and an outlet blade angle of about 0.128 times the reference design were realized. That is, considering the change in the flow path inside the impeller, the blade wrap angle was reduced compared to the reference design. Reduction in the blade wrap angle led to a shorter flow path, thus reducing the friction loss. In the case of the volute, the cross-sectional area of the optimum design was less than that in the reference design. Furthermore, the cross-sectional area continuously decreased from the tongue to the outlet, and then it decreased by approximately 18.9% at the outlet, as shown in
Figure 16.
Figure 17 shows the velocity distribution at 50% span of the blade. In the optimum design, the stagnation point is formed at the leading edge of the blade, whereas in the reference design, it is located relatively downstream. This phenomenon occurs because the inlet angle of the blade is not designed to fit the operating condition, and it can increase the incidence angle to cause flow separation and reduce the hydraulic performance of the pump. In addition, the stagnation point is formed on the suction side (SS) of the blade, and the flow proceeds non-uniformly. In particular, a very high-velocity region is distributed at the leading edge (LE) of the pressure side (PS), and severe flow separation occurs at the SS, resulting in a blockage inside the passage.
The distribution of the streamlines and the inlet velocity component at 50% span are analyzed and shown in
Figure 18. As described in
Figure 17, the inlet flow collides at the SS of the blade in the case of the reference design; conversely, in the optimum design, the inlet flow collides precisely with the LE of the blade. Through the optimization, the blade inlet angle (
βb) of the optimum design is greater than that in the reference design, and, therefore, the inlet flow proceeds smoothly, as shown in
Figure 17b.
The rotational velocity (
U) of the blade, the axial velocity (
Ca), and the relative velocity (
W) of the flow at the inlet of the impeller can be represented in the velocity triangle diagram. When the three velocity components are known, the flow angle (
β) can be derived, and the incidence angle (
I) can be calculated by comparison with the blade angle (
βb). The incidence angle distribution from the hub to the tip is shown in
Figure 19. The incidence angle of the reference design gradually decreases from the hub to about 80% span and then increases again to the tip span. The tendency of the optimum model is similar to that of the reference design, but the incidence angle is small overall in all spans, compared with the reference design. The largest incidence angle in the reference design is 24° at the tip region, whereas that in the optimum design is 18° at 15% span.
In rotating-fluid machines, the flow is driven by the centrifugal force. Hence, the incidence angle at the tip region is smaller than that at the hub; this improves the hydraulic performance. In particular, since the pump used in this study is a special pump for transporting wastewater, the blade angle distribution from the hub to the tip should be maintained. From this point of view, it is judged that the incidence angle distribution of the optimum design shown in
Figure 19 is ideal.
Figure 20 shows the distribution of streamlines and vortices on the iso-surface of the velocity invariant (>5 × 10
5 s
−2), to visualize the flow characteristics inside the impeller. As shown in
Figure 17 and
Figure 18, in the case of the reference design, the inlet flow angle is significantly different from the blade inlet angle, and an apparent inter-blade vortex is observed inside the passage, as shown in
Figure 20a. On the other hand, the inlet flow angle of the optimum design agrees well with the blade inlet angle, and the flow proceeds smoothly into the passage, as shown in
Figure 20b.
In order to analyze the blade loading, the pressure distribution on the blade in the optimum and reference design are compared, as shown in
Figure 21. The pressure values are normalized by using the maximum pressure value of the reference design. Similar to
Figure 17 and
Figure 18,
Figure 21 shows that the static pressure of the reference design at LE is relatively less than that in the optimal design, owing to the rapid increase in the velocity in the PS of the blade LE. In addition, the pressure on the PS and SS of the reference design shows a larger overall region than in the optimum design; in other words, the blade loading of the reference design is larger.
Figure 22 shows the internal velocity distribution from the volute tongue to the outlet. In the reference design, a relatively high-velocity region is distributed near the tongue, but the velocity decreases toward the outlet. On the other hand, the optimum design shows a relatively uniform velocity distribution over the entire area of the volute. In the reference design, the cross-sectional area of the volute is designed to be larger than necessary, and as the flow diffuses near the outlet, it is judged that the static pressure recovery occurs more than in the optimum design.
Figure 23 shows the iso-surface contours for a low velocity of 1.0 m/s, to identify the low-velocity region inside the pump. As shown in
Figure 22, a relatively wide low-velocity area is formed at the outlet of the reference design. Since the two-vane pump considered in this work is used for transporting wastewater, the solids move inside the impeller and volute. Now, the presence of low-velocity region can cause solids to stagnate. Therefore, the pump should be designed so that no such low-velocity region is created. From this point of view, the optimum design is an ideal design with a significant reduction in the area of the low-velocity region.
6. Conclusions
In the present study, the shapes of a blade and volute were optimized to improve the hydraulic performance of a two-vane pump for transporting wastewater via 3D RANS analysis based on supervised machine-learning optimization. Considering the waste solids, the blade inlet and outlet angles and the cross-sectional area of the volute were selected as the design variables. Three objective functions were considered to maximize the hydraulic performances, i.e., the total head and efficiency of the pump, while maximizing the size of the waste solids.
In this analysis, the RBNN predictive model was trained by using machine learning, and the GA was used to find global optimal solutions. The network training was performed by adjusting the CV error by changing the SC. Through training, an accurate and reliable RBNN surrogate model was constructed, and five arbitrary optimum design points were derived. They were compared with the numerical results to verify the predicted accuracy. The machine-learned surrogate model showed very accurate predictions compared with the numerical results and had relative errors less than 5%.
Through the optimization, the pump was made compact, and the efficiency was improved by about 14% compared to the reference design, while satisfying the design goal of the total head (>10 m). As the cross-sectional area of the optimum volute was reduced, excessive diffusion at the outlet was reduced, and, therefore, the velocity distribution inside the volute was more uniform compared with the reference design. In addition, as the inlet angle of the optimum blade increased, the water flowing from the impeller inlet proceeded smoothly into the flow path, and the flow separation and formation of the inter-blade vortex were significantly reduced.
In future work, the hydraulic performances of the optimum design will be investigated in more detail experimentally, and the characteristics of the flow, containing solid wastes, will be also analyzed.