1. Introduction
Industrial dust and industrial soot emissions are the main pollution sources affecting the quality of the atmospheric environment [
1]. As effective facilities for controlling dust, especially fine particles, bag filters have the characteristics of high efficiency and good economy and are widely used for dust treatment in coal, electric power, steel, cement, and waste incineration industries [
2,
3]. The pulse-jet cleaning process is an important part of the bag filter workflow, and its performance is directly related to the effect of energy saving and emission reduction. By adopting advanced methods, optimized cleaning parameters can be obtained, thus improving the overall operating efficiency of bag filters.
In order to optimize the pulse-jet cleaning parameters of bag filters, relevant scholars obtained better cleaning parameters by constructing a digital simulation and optimization model for bag filters, which can improve cleaning performance. Kang et al. [
2] adopted a numerical simulation for a slit-nozzle optimization design to remove the caked-on dust of the dead zone. Park et al. [
4] relied on numerical simulation to study the influence of filter bag length on the distribution of filtration velocity and proposed a scheme to improve the uniformity of filtration velocity. Fan et al. [
5] conducted a numerical simulation of the bag filter dust-cleaning process based on fluid dynamics to study the effects of factors such as pulse-jet pressure, nozzle diameter, pulse-jet height, and filter bag length on the flow field distribution inside the bag, established a quadratic polynomial prediction model with four influencing factors, and modified the models to provide theoretical guidance for the optimization of the bag filter dust-cleaning system.
Complex and computationally expensive simulations and physical experiments are often required to quantify the economic and engineering performance of complex products such as bag filters. The surrogate model approach has been widely recognized to avoid the computational burden of directly adopting the finite element model for parameter optimization.
The surrogate model refers to the mapping between the design variables and the optimization objective based on the existing design point data with low computational effort and a short calculation period, while the calculation results are similar to the numerical analysis or physical test results. The focus of surrogate modeling is the construction of approximate models to evaluate system performance and obtain the relationship between the inputs and outputs, so as to obtain the influence of each specific variable on the optimization target, which can bring a more refined design experience to the designers [
6]. Common surrogate model methods include the polynomial response surface method (PRSM), radial basis function (RBF), and Kriging.
Xu et al. [
7] applied the PRSM surrogate model to the improvement of bentonite in shield construction. The results displayed the optimal parameters for strongly weathered granite strata at different shield advance speeds, thus providing a technical guarantee for efficient construction. Thakre et al. [
8] used the PRSM surrogate model in the parameter sensitivity analysis of the selective laser sintering process, which provided conditions for realizing the robust performance of additive manufacturing processes.
Sun et al. [
9] introduced the global optimization technique based on the RBF surrogate model into the correction and optimization of key parameters of honeycomb sandwich plates in spacecraft, which reduced the dynamics analysis error. Xu et al. [
10] used the RBF surrogate model to optimize the parameters of the annular jet pump to improve its hydraulic performance in submarine trenching and dredging.
Wang et al. [
11] employed the Kriging surrogate model to optimize the design of drum brake stability, and it was found that the introduction of Kriging could greatly improve the solution efficiency. Li et al. [
12] adopted the Kriging model for the design parameter optimization of a flat-type permanent magnet linear synchronous motor for improving average thrust while reducing thrust fluctuation.
Different surrogate models have their own applicable scenarios [
13]. The PRSM is suitable for solving low-order problems in low-dimensional spaces, Kriging is appropriate for solving complex nonlinear problems in low-dimensional spaces and slight nonlinear problems in high-dimensional spaces, and RBF is fit for solving high-order nonlinear problems in high-dimensional spaces. Moreover, the optimization results of a single surrogate model for different problems are highly variable and unpredictable, which brings some risks to engineering optimization.
To address the above issues, a hybrid surrogate model approach based on multiple surrogate models can reduce the risk of using a single surrogate model for optimization while improving the accuracy of the model. Li et al. [
14] combined three surrogate models to construct a hybrid surrogate model for the structural design of permanent magnet drives, and the optimization design effect was significantly improved. Xie et al. [
15] applied the hybrid surrogate model to the comprehensive dynamic performance optimization of rail vehicles to obtain more reasonable suspension parameters, resulting in significant improvements in both the smoothness and stability of the vehicle operation. He et al. [
16] applied the hybrid surrogate model to the shape and structural optimization design of the multi-bubble pressure cabin in an underwater vehicle, which enabled the pressure cabin to have a stronger load capacity and a significantly lower weight.
Compared with different single surrogate models, optimization approaches based on the hybrid surrogate model can combine the superior characteristics of sub-surrogate models to make the model more accurate, thus improving the accuracy of prediction and optimization results. The pulse-jet cleaning process of the bag filter is a complex, nonlinear, dynamic process. In this paper, a more suitable hybrid surrogate model approach will be selected to optimize the pulse-jet cleaning-related parameters to improve the performance of the bag filter.
2. Analysis of the Bag Filter Dust-Cleaning Optimization Problem
The dust-cleaning system of the pulse-jet bag filter is mainly composed of the compressed air cylinder, the clean air chamber, the tube sheet, several filter bags, and the corresponding pulse-jet valves, injection pipes, nozzles, etc. The dust-cleaning process starts with the initiation of the pulse-jet valves, and the pulse-jet time is only roughly 0.1 s to 0.5 s. The schematic diagram of the pulse-jet cleaning state is shown in
Figure 1 [
17]. When the pulse-jet valve is activated, the high-speed pulse airflow immediately enters the filter bag from the compressed air cylinder through the injection pipe and nozzle. At the same time, it will also cause a low-pressure area at the entrance of the filter bag, thus inducing a secondary airflow into the bag filter, while the secondary airflow rate is several times the initial pulse-jet compressed air [
18]. As a large amount of airflow enters the filter bag, the pressure inside the filter bag rises sharply.
The strong airflow impact will make the filter bag vibrate, while the pressure difference between the inside and the outside of the filter bag will make the bag, along with its surface dust cake, move outward in a radial direction. Under the action of tension, the radial velocity of the filter bag will gradually drop to zero, and then the filter bag will shrink inward in a radial direction [
19]. Because the dust cake attached to the outer wall of the filter bag is not subject to the tension, the dust cake will overcome the adhesion with the filter bag under the action of inertia and fall into the dust hopper under the action of gravity, thus completing the dust-cleaning process.
In the process of pulse-jet cleaning, there are many factors that affect the dust-cleaning performance [
18] such as the air cylinder volume, pulse-jet pressure, pulse-jet height, pulse-jet time, structure of nozzles, filter material type, filter bag diameter and length, clean air chamber volume, and so on.
An optimization model should be established for the pulse-jet system, and then the multiple relevant design parameters can be optimized simultaneously to obtain the optimum dust-cleaning performance. In the study of the optimization of dust-cleaning performance, four parameters which have a significant impact on the dust-cleaning performance, as well as being easily accessible in engineering, were selected. They are the pulse-jet pressure (
P), filter bag length (
L), filter bag diameter (
D), and nozzle diameter (
d). According to the actual conditions of engineering, the design intervals of the variables are shown in
Table 1.
In previous studies on the optimization of bag filter dust-cleaning performance, dust detachment from the filter bag was mainly evaluated by the sidewall peak pressure. In the process of dust cleaning, the sidewall pressure of the filter bag is unevenly distributed in the axial direction, and the pressure is relatively small in the bottom part of the filter bag, so the effect of dust cleaning is also poor. In other words, when the conditions at the bottom of the filter bag meet the cleaning requirements, the entire bag’s cleaning performance will also be guaranteed. When considering the force required, the dust cake and the filter bag should be considered as a whole. In practical engineering applications, it is difficult to accurately obtain the outer pressure of the filter bag due to the complex characteristics of the dust cake. Therefore, the average pressure on the inner side wall of the filter bag at 1 m from the bag bottom is selected as the optimization target, which is recorded as
F. That is, the maximum
F is the optimization goal of the cleaning performance. The multi-parameter optimization model for pulse-jet cleaning of the bag filter can be expressed as:
3. Hybrid Surrogate Model for Multi-Parameter Optimization
The surrogate model transforms the black box problem into a smooth and continuous explicit mathematical problem, which greatly reduces the computational complexity while satisfying the accuracy requirements. Different surrogate models have different characteristics and apply to different optimization problems. If they can be combined to build on their strengths and avoid their weaknesses, the optimization results may be better. Therefore, the hybrid surrogate model method is proposed to expect better optimization results.
Set as follows,
is the n-dimensional input variable,
is the output response, and the corresponding response data set is
for the training sample data set
of size
. The functional relationship between
and
can be expressed as:
where
is the surrogate model output function, and
is the error when approximating
with
.
The essence of the hybrid surrogate model is a weighted linear superposition of different single surrogate models called sub-surrogate models [
20], which can be expressed as:
where
is the output of the hybrid surrogate model,
is the number of sub-surrogate models,
is the output of the
sub-surrogate model, and
is the weight of
.
Using two sub-surrogate models cannot provide enough sub-surrogate model hybrid modes, and using four or more sub-surrogate models will greatly increase the amount of calculation. In order to make better use of the advantages of each sub-surrogate model and reduce the influence of its defects on the hybrid surrogate model, there are three sub-surrogate models selected by the proposed hybrid surrogate model method, which are the PRSM model, the RBF model, and the Kriging model.
3.1. Sub-Surrogate Model
3.1.1. PRSM
PRSM was first proposed by the mathematicians Box and Wilson in 1951. The process of PRSM is to construct a mapping relationship between the design variables and the target response through least squares regression based on the given input and output values. It can be expressed as:
where
is the
polynomial,
is the corresponding polynomial coefficient to be solved,
is the polynomial matrix, and
is the coefficient matrix.
The most common usage of PRSM in engineering applications is the second-order PRSM [
21]. Its general formula can be expressed as:
where
is the estimated value of the response, and
is the
component of the n-dimensional variable
.
,
,
, and
are the coefficients of the objective function to be solved, respectively, which can be arranged in a certain order to form the column vector
.
The polynomial response surface method uses polynomials with explicit expressions to replace complex simulation analysis models, which are easy to construct. The models constructed by PRSM have good continuity and derivability to reduce the influence of numerical noise. At the same time, the influence of each variable on the output target response can be judged by the coefficient of each component [
22]. The engineering practice shows that the PRSM can effectively reduce the number of simulation analyses. However, if there are too many design variables, the number of sample points will increase accordingly, and the computational burden will increase significantly.
3.1.2. RBF
RBF is a method that uses discrete multivariate data to fit unknown functions. The basic principle is to linearly weight the model constructed with a typical radial function as the basis function, transforming the multi-dimensional problem into a one-dimensional problem with the Euclidean distance between the known sample points and the unknown point to be measured as the independent variable [
23]. The function value corresponding to the unknown sample point
based on RBF can be expressed as:
where
is the Euclidean distance, and
is the weight coefficient of the radial function for the
sample point
.
The response vector [
24] can be expressed as:
where
is the radial function and
is the Euclidean distance. When the sample points do not coincide and
is positive definite [
25], the above equation has a unique solution
. After that, the function value at the unknown sample point can be predicted.
Among these common radial functions, the multiquadric (MQ) function
shows good effects in terms of accuracy, stability, and computational efficiency [
26]. The RBF is an interpolation-type surrogate model, and the approximation accuracy is strongly influenced by the shape coefficient
. The optimal value of
depends on the distribution of sample points and approximate targets.
The RBF has the advantages of simple structure, good flexibility, high computational efficiency, and fast convergence, making it suitable for solving problems with multiple variables. However, the model is sensitive to numerical noise, so it is slightly insufficient for fitting strongly non-linear responses [
27].
3.1.3. Kriging
Kriging is an unbiased estimation model with minimum estimation variance [
28]. This model adds a random process on the basis of the global regression model, which not only considers the influence of the distance relationship between the sample points on the output, but also the influence of the spatial distribution and location relationship between the sample points on the output [
29]. It can be generally expressed as:
where
is the global regression model to reflect the overall trend of the system response in the design space,
is the polynomial function vector,
is the regression coefficient vector, and
is a stochastic process with a mean of zero, a variance of
, and a non-zero covariance, which reflects the approximation of local deviation. The covariance satisfies:
where
is a spatial correlation function, indicating the spatial correlation between sample points
and
, which plays a decisive role in the accuracy of the simulation [
30].
where
is the kernel function of the correlation function
;
is the difference between the
component of the sample points
and
, namely
;
is the dimension of the sample points, that is, the number of design variables; and
is the model parameter of the kernel function in the
direction.
The GAUSS model has the characteristics of being smooth and differentiable everywhere [
31]; therefore, it is usually adopted in engineering applications to construct the correlation function model.
The Kriging model needs to minimize the mean square error of prediction
under the condition of unbiased estimation.
The Kriging model is based on the information of known sample points, fully considering the spatial correlation characteristics of variables. Moreover, the model has both local and global statistical properties. These characteristics make Kriging advantageous in solving problems of high nonlinearity to achieve desirable results [
32].
3.2. Method for Determining the Weighting Factor
There are two key aspects of the modeling process of the hybrid surrogate model: one is to choose suitable single surrogate models as its sub-surrogate models, and the other is to determine the hybrid strategy, that is, to determine the weight coefficients through an effective calculation method.
The inverse proportional averaging method regards the accuracy indexes of each sub-surrogate model as irrelevant and has no systematic deviation [
33]. The proportion of each sub-surrogate model in the hybrid surrogate model will be calculated according to this method.
where
is the approximate capability evaluation value of the
sub-surrogate model. In this paper, the root mean square error (RMSE) is picked as
.
3.3. Hybrid Surrogate Model Optimization Algorithm
The optimization algorithm proposed in this paper is a multi-parameter optimization strategy based on a hybrid surrogate model. The initial sample set is obtained through simulation experiments, and the hybrid surrogate model will be constructed. The accuracy of the hybrid surrogate model should be inspected in combination with the actual problem, and if it meets the requirements, the global optimal solution of the model will be obtained by the genetic algorithm; otherwise, the hybrid agent model needs to be updated by selecting a suitable addition strategy until the model accuracy meets the requirements. The process of optimization based on the hybrid surrogate model is displayed in
Figure 2.
4. Optimization of Dust-Cleaning Performance of Bag Filter
4.1. Simulation and Modeling
This study relies on the bag filter laboratory, as shown in
Figure 3a, and the digital model of the bag filter, as shown in
Figure 3b, is constructed based on physical equipment. In order to facilitate the study, a single filter bag is taken as the object, and some simplifications are applied. The finite element analysis model is shown in
Figure 3c.
In this paper, the main settings of the numerical simulation are as follows:
The pressure inlet boundary condition is adopted at the inlet of the injection pipe. Considering the change of outlet pressure during the opening and closing process of the pulse-jet valve, the user-defined function (UDF) is adopted to define the function of pulse-jet pressure changing with time. According to the physical experiment, the curve of pulse-jet pressure changing with time is shown in
Figure 4.
The boundary condition of the filter bag is set to the porous-jump medium.
The pulse-jet cleaning process satisfies the law of conservation of energy, the law of conservation of mass, and the law of conservation of momentum.
Too few sample points make it difficult to guarantee the surrogate model’s accuracy, while too many sample points will increase the computation cost. To ensure the space-filling of sample points and improve the computational efficiency, the Latin hypercube sampling method [
34] based on the maximizing minimum distance criterion is employed to perform the experimental design, and 70 sets of sample points are obtained to construct the surrogate models. To confirm the correctness of the proxy model, 30 sets of validation sample points are obtained using the same method to test the accuracy. The distribution of modeling sample points and validation sample points is shown in
Figure 5.
According to the experimental design conditions of 100 sets of sample points, the ANSYS is used to simulate the pulse-jet dust-cleaning process of the bag filter, as shown in
Figure 6.
The surrogate models will be constructed based on the data of the 70 sample points. Meanwhile, considering the large range of data variation in each dimension, it is necessary to normalize the input variables to reduce their impact on the accuracy of the surrogate models [
35].
where
is the normalized design variable,
is the initial design variable, and
and
are the minimum and maximum values of the initial design variable. It can be obtained:
where
,
,
, and
are the normalized values of
P,
L,
D, and
d, respectively.
According to the normalization and simulation results of the modeling sample points, the sub-surrogate models can be constructed.
- 2.
RBF sub-surrogate model
In this paper, the MQ function is employed as the kernel function of the RBF sub-surrogate model, which can be expressed as:
where
,
,
, and
denote the four-dimensional components of the
modeling sample point, respectively. The solution can be obtained:
,
.
- 3.
Kriging sub-surrogate model
In this paper, the quadratic function
is used as the global approximation function in the Kriging sub-surrogate model.
4.2. Accuracy Evaluation
In order to ensure the validity of sub-surrogate models, error analysis is required. Commonly used indicators for evaluating the model accuracy are root mean square error (RMSE), coefficient of certainty (
) [
36], and relative average absolute error (RAAE) [
37].
According to the formulas, the accuracy indexes of three sub-surrogate models are shown in
Table 2. It can be seen that the accuracy of the RBF sub-surrogate model is better than the other two.
4.3. Construction of Hybrid Surrogate Model
The basic form of the hybrid surrogate model for pulse-jet cleaning performance can be expressed as:
Three sub-surrogate models can construct multiple hybrid surrogate models. When the number of sub-surrogate models is 2, there are three hybrid surrogate models. When the number of sub-surrogate models is 3, there is one hybrid surrogate model.
According to formula 13, the weights of the second-order PRSM model and the RBF model in the hybrid surrogate model RH-HSM are
= 0.477680 and
= 0.522320.
- 2.
Hybrid surrogate model constructed by the second-order PRSM and Kriging (PK-HSM)
According to formula 13, the weights of the second-order PRSM model and the Kriging model in the hybrid surrogate model PK-HSM are
= 0.488495 and
= 0.511505.
- 3.
Hybrid surrogate model constructed by the RBF and Kriging (RK-HSM)
According to formula 13, the weights of the RBF model and the Kriging model in the hybrid surrogate model RK-HSM are
= 0.510827 and
= 0.489173.
- 4.
Hybrid surrogate model constructed by three sub-surrogate models (PRK-HSM)
According to formula (13), the weights of the second-order PRSM model, the RBF model, and the Kriging model in the hybrid surrogate model PRK-HSM are
= 0.318415,
= 0.348172, and
= 0.333413.
The accuracy indexes of four hybrid surrogate models are shown in
Table 3. It can be seen from
Table 3 that the accuracy of the hybrid surrogate model RK-HSM is better than the other three hybrid surrogate models.
The accuracy of the hybrid surrogate model is affected by the accuracy of its sub-surrogate model. Comparing the accuracy of three sub-surrogate models and four hybrid surrogate models, it can be seen that hybrid surrogate models are better, among which the hybrid surrogate model RK-HSM performs best. The model accuracy indexes RMSE, , and RAAE of RK-HSM are improved by 9.91%, 4.41%, and 15.60%, respectively, compared with the sub-surrogate model PRSM. With RMSE as the overall evaluation metric, the performance of the seven surrogate models is ranked as RK-HSM > RBF > PR-HSM > PRK-HSM > Kriging > PK-HSM > Second-order PRSM; with as the overall evaluation metric, the ranking is similarly RK-HSM > RBF > PR-HSM > PRK-HSM > Kriging > PK-HSM > Second-order PRSM; and, with RAAE as the overall evaluation metric, the ranking is RK-HSM > PR-HSM > PRK-HSM > RBF > Kriging > PK-HSM > Second-order PRSM. The accuracy metrics of the sub-surrogate model RBF and Kriging are better than some constructed hybrid agent models, which is related to the construction principle of the surrogate model.
Applying the best-performing hybrid surrogate model RK-HSM in this study, the values of the average pressure
on the inner side wall of the filter bag at 1 m from the bag bottom during the pulse-jet dust-cleaning process under different parameter combinations can be obtained, as shown in
Figure 7, which can be used to study the influence of different parameters on the pulse-jet dust-cleaning process. The parameter conditions corresponding to the images in
Figure 7 are shown in
Table 4.
As can be seen from
Figure 7, the output response
is comprehensively affected by the multiple parameters. The change rules of the average pressure on the inner side wall of the filter bag at 1 m from the bag bottom
can be roughly understood through these images. However, as these response surfaces are interlaced, it is difficult to describe the variation completely and accurately. In order to obtain the optimal parameter combination which results in the maximum
F, the genetic algorithm is chosen in this paper.
4.4. Genetic Algorithm for Seeking Optimal Parameters
A genetic algorithm [
38] is a bionic algorithm conceived to search for the optimal solution based on the principle of biological evolution. As a global search heuristic algorithm, the genetic algorithm is suitable for solving complex optimization problems with strong robustness [
39]. In this paper, combining the characteristics of the genetic algorithm and the common experience of parameter setting, the operational parameters of the genetic algorithm are set as shown in
Table 5.
Under the initial experience setting condition (
= 250 kPa,
= 8 m,
= 145 mm, and
= 16 mm), the simulation result of
F is −2526.6148 Pa. For the seven surrogate models constructed above, including three sub-surrogate models and four hybrid surrogate models, the genetic algorithm is applied to seek the optimal combination of the four variable parameters (
P,
L,
D, and
d) involved in the pulse-jet dust-cleaning process to maximize
F. The results are shown in
Table 6.
As can be seen from
Table 6, the adoption of the genetic algorithm based on surrogate models can increase the simulation result of average pressure
F in the pulse-jet dust-cleaning process by more than 30% compared with the result under the initial empirical parameter condition (
= 250 kPa,
= 8 m,
= 145 mm, and
= 16 mm). In addition, the comprehensive performance of the hybrid surrogate models is better than that of the sub-surrogate models. Among these surrogate models, the hybrid surrogate model RK-HSM performs more outstandingly, increasing the
F by 52.30% to −1205.1605 Pa, which significantly improves the dust-cleaning effect of the bag filter. In summary, the hybrid surrogate model RK-HSM performs superiorly in both model accuracy and genetic algorithm seeking, that is, the hybrid surrogate model RK-HSM is more suitable for solving the optimization of bag filter dust-cleaning performance.
5. Conclusions
This paper proposes a hybrid surrogate model construction method based on the second-order PRSM, the RBF, and the Kriging as sub-surrogate models. Based on different hybrid modes, four hybrid surrogate models are constructed, namely PR-HSM, PK-HSM, RK-HSM, and PRK-HSM. On this basis, combined with the genetic algorithm-parameter-seeking strategy, the hybrid surrogate model optimization methods are provided for optimizing the pulse-jet cleaning process of industrial bag filters.
For the multi-parameter optimization of the bag filter pulse-jet cleaning process, the accuracy indexes RMSE, , and RAAE of the hybrid surrogate model RK-HSM are improved by 9.91%, 4.41%, and 15.60%, respectively, compared with the sub-surrogate model PRSM, which greatly enhances the reliability and practicality of the surrogate model method. What’s more, the accuracy metrics of the sub-surrogate models RBF and Kriging are better than some constructed hybrid agent models, which is related to the construction principle of the surrogate model. Therefore, the choice of the surrogate model should be appropriate to the problem to be solved, rather than aiming at the complexity of the method.
Aiming at maximizing the average pressure F on the inner side wall of the filter bag at 1 m from the bag bottom, the genetic algorithm is used to search for the optimal parameter combination, which includes the pulse-jet pressure (P), filter bag length (L), filter bag diameter (D), and nozzle diameter (d). The simulation results show that among the seven surrogate models, including three sub-surrogate models and four hybrid surrogate models, the hybrid surrogate model RK-HSM performs more prominently, increasing the F by 52.30% to −1205.1605 Pa, compared with the results under the initial empirical parameter condition. The application of the optimization method based on the hybrid surrogate model RK-HSM can greatly improve the dust-cleaning effect.