Establishing Surrogate Model to Predict the Optimal Thermodynamic and Economic Performance of a Packed Bed Humidiﬁer via Multi-Objective Optimization

: In this paper, to optimize the thermodynamic and economic performance of a packed bed humidiﬁer, a multi-objective optimization combined response surface method with a genetic algorithm is employed. The critical parameters, including geometric and thermodynamic parameters, are designated as the impact factors, and the objective functions contain unit humidiﬁcation capacity of volume and unit humidiﬁcation capacity of cost in a Box–Behnken design. The results of the analysis of variance demonstrated that the quadratic regression models of objectives are reliable and robust. It is found that the liquid–gas ratio, the interaction of the liquid–gas ratio, and inlet water temperature are simultaneously the strongest inﬂuence factors for the thermodynamic and economic indicators among the independent and interactive parameters. In addition, the optimal parameter group is found out through a genetic algorithm, and the actual optimal results are obtained as 0.11 kgs − 1 m − 3 for thermodynamic performance and 15.86 kg$ − 1 for economic performance. Furthermore, it is shown that the thermodynamic performance improves by 56% and the economic performance increases by 6.55%, compared with optimum experimental design points. During the optimization design process, the computational time to ﬁnd the optimal values reduces from 69,000 s with previous mathematical models to 10 s with established regression models. Additionally, a series of Pareto-optimal points for possible best thermodynamic and economic performance give the reference for the designers of packed bed humidiﬁers.


Introduction
Facing the rapid growth of the global economy and the improvement of the people's living standards, the demand for water and energy is becoming more and more intense, which promotes the agenda of energy structure optimization tending to be more costeffective and energy-efficient. In addition to chemical engineering and office buildings, humidifiers are also widely used in people's daily lives for regulating air humidity. Among them, the packed bed humidifier has tremendous applications thanks to its high energy effectiveness and large contact area [1]. Packings such as wooden slates [2], corrugated aluminum sheets [3], polypropylene [4], and ceramic foam [5] have shown great humidification performance. Designers expect to achieve maximum humidification effectiveness with minimum costs. However, the relationship between humidification capacity and the input cost of the humidifier is proverbially a positive correlation in the physical process. Therefore, it is significant to optimize the overall performance of a packed bed humidifier referred to thermodynamics and economics. process based on the genetic mechanism of Darwin's bio evolution theory. Moreover, it is verified that GA is a more complete and sounder model than particle swarm algorithm. A more accurate correlation for the temperature prediction of solar panels was established based on GA, while the independent factors and interaction factors between solar radiation, ambient temperature, wind velocity, and ambient relative humidity were considered [17]. It was found that the prediction accuracy of the solar panel temperature was improved by approximately 2-4 times. For further improvement, combination optimizations between RSM and GA [19] and RSM and the particle swarm algorithm [20] were developed to improve the experimental design. It was demonstrated that most performance parameters after multi-objective optimization were better than those without optimization. The effects of pitch, diameter, pipe length, and inlet velocity on the separation efficiency and pressure drop were investigated based on the RSM and GA [21]. The optimization results showed that the optimization degree of the separation efficiency had a certain limitation at the optimal points, while the pressure drop can be reduced by approximately 50%. Recently, Shahverdian [22] proposed a dynamic multi-objective optimization method in order to improve the power output and simultaneously reduce the cooling water consumption of a photovoltaic module. It was found that the temperature difference between the average and maximum value in a year evaluated by the novel method decreased, respectively, by 79.79% and 54.53%, compared with no-cooling and a constant water flow of 0.1 LPM. These achievements obtained by multi-objective optimization provide a robust foundation for optimizing the overall performance of the packed bed humidifier.
From the literature survey, it was found that the improvement of thermodynamic performance of the packed bed humidifier is generally at the expense of input cost, while the performance optimization is subject to the univariate optimization (single-objective optimization). In addition, the overall optimization of thermodynamic and economic performance for the packed bed humidifier is missing. Moreover, a previous physical model can become too complicated and time-consuming in the optimization process. To simplify the optimization process, save time cost, and improve overall performance simultaneously, in this study, the combination optimization of RSM and GA is employed to assess and optimize the thermodynamic and economic performance of a packed bed humidifier filled with corrugated wire mesh. Surrogate models for quickly searching the optimal thermodynamic and economic performances are established. Furthermore, the influences of independent and interactive factors on designated objective functions are addressed, and an explicit influence degree of the important parameters is given based on the analysis of variance. Finally, a series of Pareto-optimal points for possible optimal thermodynamic and economic performance is given with optimal parameter combination, which can provide significant references for the designer of packed bed humidifiers.

Mechanism Description
The humidification process is restricted by the thermodynamic state of the humid air. Its humidification ability will be lost once the humid air reaches saturation. In this case, in order for the physical process to continue, the temperature of the humid air should be raised to bring it back to the unsaturated state. The schematic diagram of the humidification process is shown in Figure 1. The packed bed humidifier is mainly composed of cylindrical shell and packing, and the packing can provide a relatively large direct contact area for the heat and mass transfer process. In this study, a counter-flow humidification process with corrugated wire mesh packing, Mellapak 700 Y, is selected as the optimization object. Commonly, hot water is sprayed directly into the packing layer to form a uniform liquid film on the packing surface under the action of gravity, while the relatively cold humid air enters from the bottom of the humidifier. Owing to the direct contact of two-phase fluid, the thermodynamic process of simultaneous heat and mass exchange is intense and intricate in the packing zone. In addition to the temperature difference, the water molecules in the water stream driven by the pressure difference are constantly captured by the humid air, such that the enthalpy of the humid air and water shift along the process curve, as shown in Figure 1b. Finally, the humid air stream with high temperature and humidity is discharged from the top of the humidifier, while the water stream with low temperature is exhausted from the bottom. tion object. Commonly, hot water is sprayed directly into the packing layer to form a uniform liquid film on the packing surface under the action of gravity, while the relatively cold humid air enters from the bottom of the humidifier. Owing to the direct contact of two-phase fluid, the thermodynamic process of simultaneous heat and mass exchange is intense and intricate in the packing zone. In addition to the temperature difference, the water molecules in the water stream driven by the pressure difference are constantly captured by the humid air, such that the enthalpy of the humid air and water shift along the process curve, as shown in Figure 1b. Finally, the humid air stream with high temperature and humidity is discharged from the top of the humidifier, while the water stream with low temperature is exhausted from the bottom.

Governing Equations
In order to accurately calculate the thermodynamic characteristics of the heat and mass transfer process in the packing zone, the finite difference method [23] is applied from the perspective of the control volume, as shown in Figure 2. By evenly dividing the whole packing zone into abundant control volumes, the differential governing equations can be first established based on the typical control volume, and then the whole heat and mass transfer process is solved by an iterative solver.

Governing Equations
In order to accurately calculate the thermodynamic characteristics of the heat and mass transfer process in the packing zone, the finite difference method [23] is applied from the perspective of the control volume, as shown in Figure 2. By evenly dividing the whole packing zone into abundant control volumes, the differential governing equations can be first established based on the typical control volume, and then the whole heat and mass transfer process is solved by an iterative solver. tion object. Commonly, hot water is sprayed directly into the packing layer to form a uniform liquid film on the packing surface under the action of gravity, while the relatively cold humid air enters from the bottom of the humidifier. Owing to the direct contact of two-phase fluid, the thermodynamic process of simultaneous heat and mass exchange is intense and intricate in the packing zone. In addition to the temperature difference, the water molecules in the water stream driven by the pressure difference are constantly captured by the humid air, such that the enthalpy of the humid air and water shift along the process curve, as shown in Figure 1b. Finally, the humid air stream with high temperature and humidity is discharged from the top of the humidifier, while the water stream with low temperature is exhausted from the bottom.

Governing Equations
In order to accurately calculate the thermodynamic characteristics of the heat and mass transfer process in the packing zone, the finite difference method [23] is applied from the perspective of the control volume, as shown in Figure 2. By evenly dividing the whole packing zone into abundant control volumes, the differential governing equations can be first established based on the typical control volume, and then the whole heat and mass transfer process is solved by an iterative solver.   Many calculation models were proposed for solving the thermodynamic process of heat and mass transfer, based on the finite difference method, while it was verified that the Poppe model is preferred at counter-flow conditions [24]. In terms of the Poppe assumptions, the following mathematical models can be derived in detail. It is assumed that (1) the humidification process is in a steady-state; (2) the energy loss to the environment, kinetic, and potential energy items are neglected; (3) the distribution of two-phase fluids in the packing zone is even; and (4) flow resistances in pipes, water pump, and air-blower are ignored. For the whole control volume, the mass balance equation for these two streams can be written as: where Z is the total packing height, ω is the absolute humidity of humid air, and the energy balance equation between water stream and the air stream is obtained as: The enthalpy, h, of humid air consists of the enthalpy of dry air (subscript da), and water vapor (subscript w). Its expression for unsaturated and saturated cases is written as: Further, the specific heat capacity, c pa , of each node of the water is set to a constant value with the above assumptions. Therefore, the following equation can be written at constant pressure: Substituting Equation (4) into Equation (2) and divide by control volume height, dZ, on both sides, and Equation (5) can be found upon rearrangement: At the gas-liquid interface, where the simultaneous heat and mass exchange occurs, the mass balance equation of the air side can be given by: Considering the heat loss on the water side due to evaporation in the energy balance equation, the energy balance equation to calculate the sensible heat and latent heat captured by the humid air stream from the water side is expressed as: where subscript s means saturation state, h d is the mass transfer coefficient. Its empirical correlation for gas absorption system is employed in this work, which can be referred to [25] as follows: where where f is the corresponding geometric coefficient, 5.23 for d p > 0.015 m and 2 for d p < 0.015 m, and D a is the diffusion coefficient from the water to the air. In order to calculate the heat transfer coefficient, h c , Lewis factor, which represents the indication of the relative rates of heat transfer and mass transfer, is utilized according to  [26], and the empirical correlation [27] used in this study is shown as: ω+0.622 (12) where ω sw is the saturated humidity ratio calculated at the local water temperature, and the contact area of the control volume, dA, can be calculated by: Then, the following differential equations for unsaturated and saturated cases are found upon rearrangement: where m r is the liquid-gas ratio, m w /m da kg dry air. During the humidification process, it is highlighted that excess water vapor will be carried by the humid air, which will be referred to as the supersaturated cases. Considering the all-state thermodynamic of humid air, the governing equations at supersaturated conditions should also be built. In addition to the enthalpy of dry air and basic water vapor, the enthalpy of supersaturated air, h ss , should include the enthalpy of excess water vapor, expressed by: Additionally, the relevant differential equations established for the supersaturated air should also be accordingly replaced by:

Unit Humidification Capacity of Volume
In order to characterize the thermodynamic performance of the packed bed humidifier, considering the effect of packing size, a novel thermodynamic evaluation criterion, unit humidification capacity of volume (UHCV), is established based on humid air in our previous work [8]: where V p is the required packing volume, calculated by:

Unit Humidification Capacity of Cost
Here, once referred to the economic performance, the input cost should contain the investment cost and operation cost. The investment cost mainly includes the cost of airblower, water pump, and packing, while the operation cost means the power consumption of electrical components. For this case, a novel economic criterion, unit humidification capacity of cost (UHCC), is defined [8] as follow: where C b and C wp are the cost of air-blower and water pump, depending on the required power, and the packing cost, C p , is associated with the price of packing 316 L. The cost of each component can be referred to the price shown in Table 1; where C e is the operation cost, calculated according to the running time of 10 years with 7200 h per year, and the electrovalence of P e = 0.085$ kW −1 h −1 : The detailed calculation formulas of N e , N b , and N wp can be referred to in [5].

Multi-Objective Optimization Design
In this study, multi-objective optimization is employed to obtain the potential optimal thermodynamic and economic performance of the packed bed humidifier in a combination of a response surface method and genetic algorithm. In this process, the important parameters affecting the given objectives should be specified and combined based on the RSM first. The experimental design method selected currently is the Box-Behnken design (BBD) method, thanks to its efficient and convenient design scheme. As shown in Table 2, thermodynamic parameters (liquid-gas ratio, inlet air temperature, inlet relative humidity, and inlet water temperature) and geometric parameter (specific surface area of the packing) are chosen as the factor items, and the response items contain the unit humidification capacity of volume and the unit humidification capacity of cost. Based on the established regression models of responses, multi-objective optimization can be accomplished through the GA program, following the natural evolution process. In the genetic mechanism, genes are derived from the range of the factors in Table 2, and response UCHV, as well as UHCC, are simultaneously expected to be maximum objectives.

Results
According to the differential governing Equations (12)-(14) and (18)-(20) which were verified to be robust and reliable for the packed bed humidifier based on the experimental results [8], an integrated solver can be coded in the MATLAB software, as shown in Figure 3. The typical conditions of the humidification process are presented in Table 3. For the experimental design, the first step of response surface analysis is to combine design variables reasonably and generate parameter groups. The design results are used to fit the regression equations of the objective functions. After the experimental design, the analysis of variance (ANOVA) is carried out to assess the influence degree of the important factors and their interaction on the responses. Then, the genetic optimization process can be promoted to search for the optimal parameter combination and thermodynamic and economic performance of the packing bed humidifier, according to the fitting regression equations.

Regression Models of Objective Functions
Based on the BBD method, 46 groups of experimental designs are first combined, and then the relevant design results are obtained by the established mathematical models, as exhibited in Table 4. It is found that the maximum UHCV appears as 0.11 kgs −1 m −3 at group 26, while the value of UHCC is 11.15 kg$ −1 in this case. Moreover, the maximum value of

ANOVA Results and Model Diagnosis
Firstly, the accuracy of the established models for the objective functions should be evaluated through the ANOVA, and the relevant assessment indicators are presented in Table 5 in detail. The influences of mentioned independent and interactive parameters on the designative responses are revealed based on the ANOVA results, as presented in Tables 6 and 7, respectively.   The quadratic regression models obtained by the BBD method can be further diagnosed by the comparisons between predicted values and actual values, as shown in Figure 4. The quadratic regression models obtained by the BBD method can be further diagnosed by the comparisons between predicted values and actual values, as shown in Figure  4.

Thermodynamic Performance Results
According to the ANOVA, the important parameters affecting the value of the UHCV contain parameters A, B, D, E, AD, AE CE, A 2 , and B 2 . It is found that parameter B (inlet air temperature) is the weakest influence factor among the four important independent parameters, and A (liquid-gas ratio) has the strongest effect on UHCV. Moreover, the interaction parameter AD (liquid-gas ratio and inlet water temperature) has the greatest impact on the UHCV compared with the other interaction parameters AE and CE. Therefore, in order to analyze the interactive influences of parameter AD on the UHCV of the packing bed humidifier, the contour and 3D surface plots of parameter AD are presented in Figure 5 based on the RSM.

Thermodynamic Performance Results
According to the ANOVA, the important parameters affecting the value of the UHCV contain parameters A, B, D, E, AD, AE CE, A 2 , and B 2 . It is found that parameter B (inlet air temperature) is the weakest influence factor among the four important independent parameters, and A (liquid-gas ratio) has the strongest effect on UHCV. Moreover, the interaction parameter AD (liquid-gas ratio and inlet water temperature) has the greatest impact on the UHCV compared with the other interaction parameters AE and CE. Therefore, in order to analyze the interactive influences of parameter AD on the UHCV of the packing bed humidifier, the contour and 3D surface plots of parameter AD are presented in Figure 5 based on the RSM.

Economic Performance Results
According to the ANOVA, parameters A, B, C, D, E, AB, AD, A 2 , and B 2 are significant factors influencing the UHCC according to ANOVA results. It is illustrated that parameter A (liquid-gas ratio) has the strongest influence on UHCC than the other important single parameters, and parameter C (inlet relative humidity) is the weakest one. Among the important interaction factors, the influence of AD (liquid-gas ratio and inlet water temperature) is brightly greater than that of AB (liquid-gas ratio and inlet air temperature). Therefore, in order to analyze the interactive influences of parameter AD on the UHCC of the packing bed humidifier, the contour and 3D surface plots of parameter AD are presented in Figure 6.

Economic Performance Results
According to the ANOVA, parameters A, B, C, D, E, AB, AD, A 2 , and B 2 are significant factors influencing the UHCC according to ANOVA results. It is illustrated that parameter A (liquid-gas ratio) has the strongest influence on UHCC than the other important single parameters, and parameter C (inlet relative humidity) is the weakest one. Among the important interaction factors, the influence of AD (liquid-gas ratio and inlet water temperature) is brightly greater than that of AB (liquid-gas ratio and inlet air temperature). Therefore, in order to analyze the interactive influences of parameter AD on the UHCC of the packing bed humidifier, the contour and 3D surface plots of parameter AD are presented in Figure 6.

Economic Performance Results
According to the ANOVA, parameters A, B, C, D, E, AB, AD, A 2 , and B 2 are significant factors influencing the UHCC according to ANOVA results. It is illustrated that parameter A (liquid-gas ratio) has the strongest influence on UHCC than the other important single parameters, and parameter C (inlet relative humidity) is the weakest one. Among the important interaction factors, the influence of AD (liquid-gas ratio and inlet water temperature) is brightly greater than that of AB (liquid-gas ratio and inlet air temperature). Therefore, in order to analyze the interactive influences of parameter AD on the UHCC of the packing bed humidifier, the contour and 3D surface plots of parameter AD are presented in Figure 6.

Multi-Objective Optimization
Through the natural evolution process, a series of Pareto-optimal points can be found in Figure 7 compared with the results of the BBD method. In this genetic programming, the size of the population is 100, and the number of generations is set as 600. In addition, non-dominant inheritance within the population allows chromosomes to undergo selection, crossover, and mutation so that the optimal solution of the design parameter group can be found out by means of seeking a better fitness in a non-dominated population as presented in Figure 8.

Multi-Objective Optimization
Through the natural evolution process, a series of Pareto-optimal points can be found in Figure 7 compared with the results of the BBD method. In this genetic programming, the size of the population is 100, and the number of generations is set as 600.
In addition, non-dominant inheritance within the population allows chromosomes to undergo selection, crossover, and mutation so that the optimal solution of the design parameter group can be found out by means of seeking a better fitness in a non-dominated population as presented in Figure 8.

Discussion
According to the experimental design, it is found from Table 4 that it is extremely difficult for the thermodynamic and economic performance of the packed bed humidifier to be simultaneously optimal as expected at the same condition, based on the simplified orthogonal combination. Furthermore, it is expected that there is still great potential for the UHCV and UHCC to have a higher value through multi-objective optimization.
Based on the ANOVA results, it is found from Table 5 that the F-value of 2387.81 for UHCV and 414.45 for UHCC illustrate that the models are extremely significant. In addition, it is demonstrated that the p-value of both UHCV and UHCC is less than 0.0001, which shows that the regression models are effective, among which the UHCV model has better accuracy. Moreover, according to the values of R-Squared and Adj R-Squared, it is clear each of the models has a great fitting quality. It is also found that the value of Pred R-Squared is in reasonable agreement with the Adj R-Squared value both for UHCV and UHCC, and the differences are both less than 0.2, which confirmed this fact again. Furthermore, the Adeq Precision value of 214.256 for UHCV and 84.571 for UHCC are both greater than 4, which are satisfactory and an adequate signal to navigate the experimental design.
From the results in Tables 6 and 7, it is said that a p-value less than 0.05 indicates the relevant model term is significant. Therefore, it is summarized that parameters A, B, D, E, AD, AE, CE, A 2 , and B 2 have significant influences on UHCV, while parameters A, B, C, D, E, AB, AD, A 2 , and B 2 are significant model terms for UHCC. Here, it is found that in addition to the independent parameters, the interactive parameters such as AB, AD, AE, and CE also have significant impacts on the thermodynamic and economic performance. This is the reason why multi-objective optimization on the overall performance of the packed bed humidifier is developed in this study.
Through model diagnosis, it is shown from Figure 4 that most design points inhabit near the diagonal, and the relationship between predicted values and actual values is linear, which confirms that the proposed models can be employed to predict the thermodynamic and economic performance of the packed bed humidifier. Figure 5 exhibits the case that the other independent parameters are constant in the middle value of the variable range. The relationship between liquid-gas ratio and UHCV is a quadratic parabola correlation, while there is a linear relationship between inlet water temperature and UHCV, as shown in Figure 5a. Therefore, a significant interaction influence appears between these two parameters. It is found that higher values of the liquid-gas ratio and the inlet water temperature will lead to a larger UHCV evidently. This is because either increasing the inlet water temperature or increasing the liquid-gas ratio will be beneficial to the heat load of the hot side. Then, the moisture absorption capacity of humid air will rise sharply with the increase in its temperature, which will cause a higher value of the UHCV, as shown in Figure 5b. It is observed that a maximum value of 0.11 kgs −1 m −3 for the UHCV appears at the conditions of the liquid-gas ratio of 5 and inlet water temperature of 90 • C.

Economic Performance Analysis
It is highlighted that Figure 6 is plotted when the other parameters are set as the middle value. From Figure 6a, it is found that there is a slight linear relationship between the inlet water temperature and response UHCC, while the influence relationship between the liquid-gas ratio and UHCC is a quadratic parabola. The difference is that the impact trend of these two variables on UHCC is the opposite. It can be seen that the increase in inlet water temperature can enhance the value of the UHCC, while a higher value of the liquid-gas ratio will rapidly reduce the value of the UHCC. It is found that a maximum UHCC of 15.51 kg$ −1 emerges at the conditions of the liquid-gas ratio of 1 and inlet water temperature of 90 • C. That is because increasing the inlet water temperature will not only improve the moisture absorption capacity of the humid air but also increase the energy consumption of the electric heater. Therefore, the comprehensive improvement effect is not obvious. However, the decrease in the liquid-gas ratio will greatly reduce the power consumption of the water pump and electric heater. Simultaneously, since the spray water temperature is unchanged, the moisture absorption of the humid air is relatively small compared to the cost increase. As a result, the interaction of parameter AD shows a sloping surface on the UHCC from the front view in Figure 6b.

Multi-Objective Optimization Analysis
As analyzed in Sections 4.1 and 4.2, it is summarized that it is laborious to obtain the desirable thermodynamic and economic performance simultaneously for the packed bed humidifier based on BBD. To overcome this problem, a multi-objective optimization is implemented to predict the best performance for the packed bed humidifier employing the Non-dominated Sorting Genetic Algorithm-II. After natural selection, as shown in Figure 7, it can be found that the optimal points obtained are all located in the upper right of the experimental values, which indicates that the overall performance of the packing bed humidifier can be substantially improved.
It is found from Figure 8 that population B (inlet air temperature), C (inlet relative humidity), D (inlet water temperature), and E (specific surface area) have a favorite habitat around 16 • C, 0.6, 90 • C, and 700 m 2 m −3 , respectively, while population A (liquid-gas ratio) is randomly distributed in the specific range. Additionally, a higher liquid-gas ratio is beneficial to thermodynamic performance, while a lower value of liquid-gas ratio can improve the economic performance in combination with the optimal parameter group.

Verification of Optimal Solution
It was shown that the optimal performance obtained through multi-objective optimization is an optimization interval. Here, in order to verify the accuracy of the optimal values, two Pareto-optimal points recorded as point 1 (0.032 kgs −1 m −3 for UHCV, 15.99 kg$ −1 for UHCC) and point 2 (0.11 kgs −1 m −3 for UHCV, 11.74 kg$ −1 for UHCC) are selected as the comparison objects as shown in Figure 7. The corresponding optimal parameter groups are, respectively, 1, 11.62 • C, 0.6, 90 • C, and 700 m 2 m −3 , and 5, 16.71 • C, 0.6, 90 • C, and 700 m 2 m −3 , successively. Through introducing the optimal groups into the established mathematical model, the actual performances are obtained as 0.031 kgs −1 m −3 for UHCV, 15.86 kg$ −1 for UHCC corresponding to point 1, and 0.11 kgs −1 m −3 for UHCV, 11.60 kg$ −1 for UHCC corresponding to point 2, respectively. It was found that the errors between predicted and actual values are, respectively, +3.22% and +0.82% at point 1 based on the actual value, while +0.63% and +1.21% are at point 2, which demonstrate that the predicted optimal results are desirable and reliable. Furthermore, through multi-objective optimization, the value of UHCV increases by 56% compared with the result of group 22 with the same UHCC, and the value of UHCC raises by 6.55% compared with the value of group 26 with the same UHCV, respectively. Furthermore, the computational time to find the optimization results using a mathematical model and regression model is extremely different. As shown in Figure 9, the quickest computational time for optimization results employing the original mathematical model is 69,000 s, through the simplified orthogonal method, while the computational time to search the optimization results drops to only 10 s for the novel regression model through genetic algorithm. Therefore, the established regression model can give the optimization results more quickly, which justifies the usefulness and robustness of the surrogate model. specific design conditions given in Table 2. Maybe the coming genetic algorithm-III can improve this limitation.

Conclusions
To calculate the thermal characteristics of the combined heat and mass transfer process in the packed bed humidifier, an integrated mathematical model and the relevant objective functions of thermodynamic and economic performance are first established. Then, the quadratic regression models of objective functions are fitted based on the experimental design results through the response surface method. The degrees of the important influences of independent and interactive parameters on the objectives are summarized by means of the analysis of variance. Through multi-objective optimization, some conclusions can be written as follows: 1. Feasible regression models of unit humidification capacity of volume and unit humidification capacity of cost are established based on an experimental design scheme consisting of 46 groups, while the thermodynamic model is more accurate than that of economic; 2. The liquid-gas ratio has the strongest influence among the independent parameters on thermodynamic and economic performance, simultaneously, and the interaction of the liquid-gas ratio and the inlet water temperature has the greatest impact on thermodynamic and economic performance compared with other interaction factors; 3. The optimal parameter groups are found out through genetic programming: the liquid-gas ratio is randomly distributed in the design range, inlet air temperature is around 16 °C, inlet relative humidity is 0.6, inlet water temperature is 90 °C, and specific surface area is 700 m 2 m −3 , respectively. Furthermore, the Pareto-optimal points obtained are all superior to the values of experimental design; 4. After multi-objective optimization, the thermodynamic performance improves by 56%, and the economic performance elevates by 6.55%, compared with the values of basic experimental design, respectively; 5. The computational time to find the optimization results using an established regression model can be reduced by 99.98%, compared with the original mathematical model.  Although the proposed method is time-saving, efficient, and accurate, it also has some problems, such as premature convergence and nonstandard coding, especially the specific design conditions given in Table 2. Maybe the coming genetic algorithm-III can improve this limitation.

Conclusions
To calculate the thermal characteristics of the combined heat and mass transfer process in the packed bed humidifier, an integrated mathematical model and the relevant objective functions of thermodynamic and economic performance are first established. Then, the quadratic regression models of objective functions are fitted based on the experimental design results through the response surface method. The degrees of the important influences of independent and interactive parameters on the objectives are summarized by means of the analysis of variance. Through multi-objective optimization, some conclusions can be written as follows:

1.
Feasible regression models of unit humidification capacity of volume and unit humidification capacity of cost are established based on an experimental design scheme consisting of 46 groups, while the thermodynamic model is more accurate than that of economic; 2.
The liquid-gas ratio has the strongest influence among the independent parameters on thermodynamic and economic performance, simultaneously, and the interaction of the liquid-gas ratio and the inlet water temperature has the greatest impact on thermodynamic and economic performance compared with other interaction factors; 3.
The optimal parameter groups are found out through genetic programming: the liquid-gas ratio is randomly distributed in the design range, inlet air temperature is around 16 • C, inlet relative humidity is 0.6, inlet water temperature is 90 • C, and specific surface area is 700 m 2 m −3 , respectively. Furthermore, the Pareto-optimal points obtained are all superior to the values of experimental design; 4.
After multi-objective optimization, the thermodynamic performance improves by 56%, and the economic performance elevates by 6.55%, compared with the values of basic experimental design, respectively; 5.
The computational time to find the optimization results using an established regression model can be reduced by 99.98%, compared with the original mathematical model.

Conflicts of Interest:
The authors declare no conflict of interest.