1. Introduction
Plate-fin heat exchangers are widely used for heat dissipation in automotive engines because of their compact and lightweight structure, good heat transfer performance, and low production cost [
1]. Serrated staggered fins are common enhanced heat exchange surfaces in plate-fin heat exchangers. The principle is that the fins are periodically staggered at certain intervals from each other along the fluid flow direction, so that the boundary layer formed by the fluid near the fin surface enters the rear row of fins before it fully develops, making full use of the boundary layer separation effect. At the same time, the tail vortexes generated by the fluid on the upstream fins also have an excitation and enhancement effect on the heat transfer of the downstream fins [
2,
3]. In recent years, with the improvement of engine performance requirements of the extended-range hybrid vehicle engine, the requirements of the heat exchanger have been gradually increased, and many scholars have attempted to introduce new algorithms for heat exchanger optimization design research. At present, new algorithms, such as the genetic algorithm [
4], annealing simulation algorithm [
5], and model search algorithm [
6], have been successfully applied to heat exchanger optimization design research [
7]. However, these optimization algorithms have rarely been applied in engineering practice, and there is a lack of computational procedures for heat exchanger optimization design which can guide engineering application. At the same time, there is little research on plate-fin heat exchangers, which are widely used for automotive engine heat dissipation. Therefore, it is of great interest to conduct optimization design research and program development for plate-fin heat exchangers.
The optimal design of heat exchangers using traditional methods such as the logarithmic mean temperature difference (LMTD) method and the Efficiency-number of Heat-transfer Units (η-NTU) method is costly and time-consuming [
8]. With the rapid development of computational fluid dynamics (CFD) and computer technology, it has become possible to use computers to optimize the design of heat exchangers with high efficiency. Studies by many researchers have shown that simulations of various types of heat exchangers using CFD are reliable [
9,
10,
11,
12]. Therefore, CFD numerical simulations can minimize unnecessary tooling production costs, test work, and R&D time, and provide great convenience for efficient optimal design of heat exchangers.
Juan et al. [
13] analyzed the effects of structural parameters on thermal characteristics of transverse direction (TD) type serrated fin by numerical simulation methods and found that the vortexes generated near the TD fin region enhanced the turbulence intensity to reduce the thickness of boundary and improve the synergy between the velocity and temperature gradient. Minsung et al. [
14] used numerical simulation methods to investigate the aerothermal performance of a slanted-pin fin heat exchanger under high-speed-bypass stream condition, and found that the design of slanted-pin fin can make the heat exchanger lighter and smaller. However, in order to obtain the heat transfer and flow resistance characteristics of the heat exchanger as a whole, it is necessary to conduct numerical simulation analysis of the whole heat exchanger. Nevertheless, the fin scale is very small relative to the overall heat exchanger scale, and it is extremely inefficient to directly simulate numerically heat exchangers with complex structures, such as plate-fin type without the necessary simplifications. At present, porous media model is the most applied simplification method, which was first proposed by Patankar et al. [
15] for the simulation of steam generators. Zhang et al. [
16] studied the flow and heat transfer characteristics in finned-tube heat exchangers based on the porous media method, and analyzed the effect of air flow on the resistance characteristics and heat transfer performance, and the simulation results were verified by the measured data. Therefore, the simplification of the overall heat exchanger using the porous media method will help to achieve an efficient numerical simulation analysis of heat exchangers with complex structures, such as the plate-fin type.
Generally, numerical simulation results need to be combined with optimization algorithms to obtain the optimal solution. Many scholars have made many efforts to this end [
17,
18,
19]. The non-dominated sequencing genetic algorithm (NSGA-II) is widely used for optimization design due to its advantages such as fast operation and good convergence of the solution set. Alireza et al. [
20] adopted CFD and NSGA-II to carry out multi-objective optimization designs for the combustion chamber of CO steam boiler, and studied the multi-objective optimization problem with two input parameters and two objectives. Jing et al. [
21] used numerical simulation method and NSGA-II for a multi-objective optimization of mini U-channel cold plate with SiO
2 nanofluid to obtain the optimal performance. However, the Pareto frontier obtained by NSGA-II is an optimal solution set, and these studies ultimately use subjective analysis to select the optimal solution from the Pareto solution set. Technique for order preference by similarity to an ideal solution (TOPSIS) can effectively avoid the subjectivity of data, and can well depict the comprehensive impact of multiple impact indicators [
22]. Aminu et al. [
23] comprehensively evaluated the Pareto optimal solution set of NSGA-II using TOPSIS and obtained the optimal performance of concentrated Photovoltaic-Thermoelectric hybrid system.
In this study, in order to obtain the optimal performance of the plate-fin heat exchanger, the fin angle of the heat exchanger and the inlet flow rate are optimally designed. Firstly, a fluid-solid coupling model of the heat exchanger is established based on commercial CFD software Fluent, and the serrated staggered fins are simplified using porous media. By changing the fin angle, oil flow rate and water flow rate, 45 numerical simulation test cases are designed. Additionally, based on the numerical simulation results, support vector machine regression (SVR) is used to establish regression models of heat exchanger performance (heat transfer, oil pressure drop, and oil outlet temperature) with respect to fin angle and inlet flow rate, and these regression models are also taken as the optimized objective functions. Then, the Pareto optimal solution set of the objective function is obtained using NSGA-II. Finally, the optimal solution that leads to the best performance of the heat exchanger is preferred by the TOPSIS comprehensive evaluation method.
3. Multi-Objective Optimization
3.1. Optimization Objectives and Optimization Parameters
Heat exchanger as one of the heat dissipation components of the engine, heat transfer quantity and pressure drop are important indicators to evaluate its performance. The greater the heat transfer quantity, the better, and the smaller the pressure drops, the better. However, heat transfer quantity and pressure drop are contradictory to each other, and increasing heat transfer will inevitably lead to the increase of pressure drop, and vice versa. At the same time, during the operation of the heat exchanger, the closer the oil outlet temperature is to the water inlet temperature, the better. According to the simulation and η-NTU calculation results, the oil pressure drop is much larger than the water pressure drop during the working process of the heat exchanger. Therefore, the optimization objectives established in this paper are heat transfer quantity, oil pressure drop, and oil outlet temperature.
The law of heat transfer and flow in the heat exchanger is complex, and selecting the appropriate inlet flow rate on the oil side and the water side can give full play to the performance of the heat exchanger. In addition, changing the key parameters of the fins can effectively adjust the pressure drop and heat transfer quantity of the heat exchanger. There are many key geometric parameters of fins, but studies on fin angles are rare. Therefore, the optimized parameters for this study are determined as fin angle, oil flow rate, and water flow rate.
3.2. SVR Regression Model
After the data of 45 test cases are obtained by numerical simulation, it is necessary to fit them to obtain the regression models of heat transfer quantity, oil pressure drop and oil outlet temperature with respect to fin angle, oil flow rate and water flow rate. These regression models will be used as objective functions of multi-objective optimization.
Support vector regression (SVR) is an extension of support vector machine (SVM), which is widely used for regression of engineering problems. The given data set is set as D = {(x1, y1), (x2, y2), ..., (xn, yn)}, where xi Rn is the vector of input variables and yi Rn is the corresponding scalar output (target) value. The target of SVR is to minimize the “distance” to the farthest sample point of the fitted hyperplane, so as to accurately predict the target {yi} corresponding to a set of input samples {xi}.
SVR creates a “spacing band” on both sides of the hyperplane with a spacing of
ε (tolerance bias), and does not calculate the loss for all samples that fall into the spacing band. For the linear problem, SVR constructs the following linear model [
31,
32]:
Thus, its optimization objective is:
SVR allows for the presence of samples outside the spacing band, but the losses should be as small as possible so that the optimization objective of SVR can be formalized as:
where
is the empirical error term of SVR, consisting of a loss function
LC (
f (
xi),
yi) and a factor
C. The loss function is expressed as:
and the factor
C represents the weights. To explain the error beyond the limit
ε, the relaxation variables
ξ and
ξ* are introduced. At this point, all sample data meet the condition:
The relaxation variable transforms the SVR problem into a dual optimization problem, and its optimization objectives are:
The optimization problem with constraints can be transformed into an unconstrained optimization problem by Lagrange multiplier method. With the introduction of Lagrange multipliers
and
, the linear model can be rewritten as:
For nonlinear problems, the nonlinear problem can be mapped to a linear problem by boosting the sample dimension. However, when the data dimension itself is large, boosting can make the computational effort increase dramatically. To overcome the contradiction between high-dimensional feature space and computational complexity, SVR will define appropriate kernel functions [
33]. With the help of kernel functions, the results of the inner product of samples in the high-dimensional space can be computed directly in the low-dimensional space, thus greatly reducing the computational effort. The choice of kernel function requires that Mercer’s theorem is satisfied, i.e., the kernel function is semi-positive definite for any Gram matrix in the sample space [
34]. The Gaussian radial basis function (RBF), also called Radial Basis Function, can map the data to infinite dimensions with the expression [
35,
36]:
where
i,
j = 1…
m,
σ is the width of the RBF. After the SVR is mapped into the kernel function, its model can be rewritten as:
where
(
x) is the mapping function that maps
x to a higher dimensional space.
Since the RBF kernel has good universality and depends on only one parameter σ, the RBF kernel function is chosen to build the SVR regression model in this study. Meanwhile, the data of 45 simulation cases are fitted with the fin angle, oil flow rate, and water flow rate as input parameters, and the regression models of heat transfer quantity, oil pressure drop, and oil outlet temperature are finally obtained. These models will be used as the objective function for multi-objective optimization.
3.3. Multi-Objective Optimization Based on NSGA-II Algorithm
There are two main types of multi-objective optimization algorithms: ordinary gradient methods and gradient-free direct methods. The first type of method relies on the quality of the initial guess, which is easy to fall into local extremes and is only applicable to continuous smooth functions. The gradient-free direct method is more suitable for the study of nonlinear phenomena. Among them, genetic algorithms are the most widely used [
37,
38]. Such algorithms are insensitive to the discontinuity of the objective function, do not get trapped in local optima, and are suitable for parallel processing.
NSGA-II is a kind of genetic algorithm that can efficiently order the nondominated solutions while providing a set of Pareto optimal solutions well distributed along the Pareto frontier and considering an elite strategy of accelerating convergence [
39]. The algorithm is widely used to minimize or maximize two or more objective functions under given constraints and boundary conditions. The result of its optimization represents the set of solutions with the best compromise between the objective functions.
NSGA-II generates a random population in the initial state, and then the population individuals undergo crossover (parents produce offspring) and mutation (small random changes in offspring). The algorithm then sorts individuals based on non-dominance rank and crowding degree, and selects higher quality individuals to form the next generation. The population is driven towards the optimal Pareto frontier while maintaining the diversity of the population [
40]. The algorithm runs until a predefined number of generations is reached.
In this study, NSGA-II was used to optimize three conflicting objective functions (heat transfer quantity, oil pressure drop, oil outlet temperature). The population size, crossover probability, mutation probability, and maximum number of generations are set to 2000, 0.9, 0.2, and 5000, respectively. The boundary conditions are:
- (1)
30° ≤ Fin angle ≤ 90°;
- (2)
L/min ≤ Oil flow rate ≤ 15 L/min;
- (3)
L/min ≤ Water flow rate ≤ 15 L/min.
3.4. Multi-Objective Decision Making Based on TOPSIS Algorithm
The optimal solution set with 2000 solutions is eventually obtained using NSGA-II, and the multi-objective decision making method can prefer the optimal solution from it. TOPSIS is one of the most commonly used multi-objective decision making methods for selecting the best compromise solution between the incommensurable and conflicting objective functions [
22]. The optimal solution based on TOPSIS decision making is closest to the positive ideal solution and farthest from the negative ideal solution. The steps of TOPSIS multi-objective decision making are as follows:
Step 1: Create a decision matrix (aij)m × n, where m is the decision point and n is the number of objective functions.
Step 2: Normalize the decision matrix using the Euclidean method:
Step 3: Develop a weighted normalized decision matrix:
where
.
Step 4: Determine positive ideal solutions (
X+) and negative ideal solutions (
X−):
where
J+ is the indicator of positive standard,
J− is the indicator of negative standard, and
and
are the maximum and minimum values of each column, respectively.
Step 5: Calculate the distance of each objective alternative from the positive ideal solution and the negative ideal solution:
Step 6: Calculate the relative proximity of each objective alternative to the ideal state:
The solution of TOPSIS is the point with the largest value of on the Pareto optimal solution set. If the solution of TOPSIS is not reasonable, the weights (wj) of the objective function can be reassigned and the results are recalculated. In this study, the weights of all optimization objectives are set to be equal, i.e., the weights of heat transfer quantity, oil pressure drop, and oil outlet temperature are all 1/3. The TOPSIS algorithm is used to make decisions among 2000 optimal solutions, and the optimal result is preferentially selected.
In summary, the overall optimization process for the plate-fin heat exchanger is shown in
Figure 5. Firstly, the data of numerical simulation was obtained using the fluid-solid coupling model of the plate-fin heat exchanger, and then the regression models of heat transfer quantity, oil pressure drop, and oil outlet temperature with respect to the fin angle, oil flow rate, and water flow rate are obtained by SVR. Additionally, these regression models are used as the objective function for multi-objective optimization. NSGA-II is utilized to optimize on the objective function to obtain the Pareto optimal solution set. Finally, TOPSIS is applied to decide the optimal solution in the solution set.