1. Introduction
Thanks to the advantages of fast energy replenishment, low emission, smooth operation, and high energy conversion efficiency, vehicular hydrogen fuel cell systems are being developed in various countries [
1,
2,
3], despite their main drawback of high costs affecting their practical application. A hydrogen-based proton exchange membrane (PEM) FCS consists of a fuel cell stack, air supply system, hydrogen storage, thermal management system, and water management system [
4], as shown in
Figure 1. An electric motor drives an air compressor to pressurize and humidify the filtered air that is then provided to the fuel cell reactor where hydrogen and oxygen (in the pressurized) air undergo a chemical reaction that outputs electrical energy to supply the motor and air compressor. The exclusive compressors used in FCSs form one of the core components of the air supply system [
5] and account for up to 16.89% of the entire cost—the second-highest total cost of the FCS next to the cost of fuel stack [
4,
6]. Importantly, the aerodynamic performance of the compressors directly determines the comprehensive performance of the FCS [
7], which in turn has an impact on the overall performance and cost of the vehicle.
The most commonly used compressors in vehicular hydrogen fuel cells (HFC) are the screw, scroll, slide, roots, and radial compressors. They each have unique advantages [
8,
9]. For the screw, radial, and roots compressors, a large flow rate, for instance, more than 90 g/s, is easy to achieve. Of these, high pressure can be easily realized for the screw and radial-type compressors [
10]. The scroll, screw, and slide compressors are all in the type of positive displacement [
10,
11,
12,
13,
14,
15], which can boost the pressure through the shrinkage of inner chambers. In contrast, radial compressors inhale the gas and elevate the kinetic energy through a rotating impeller at high speed to boost the pressure of air. Eaton [
16] improved the R340 TVS series of roots compressors by a P-Series Roots positive displacement design and shifted the peak efficiency suit to an 80 kW PEM stack. Yuehua Li [
17] highlights screw and radial compressors as a good choice for FCSs due to their relatively low weight and high pressure. Of these, radial compressors are recognized for their high reliability and efficiency, especially those equipped with air bearings, which promotes pressure, rotating speed, and efficiency significantly.
Table 1 [
4,
13] is a compilation of data obtained from the US Department of Energy (DOE) comparing the performance of different types of compressors used in fuel cell vehicles (FCVs). Here, it is revealed that radial compressors have obvious advantages among them. As a result, the use of high-speed radial compressors as vehicular gas supply devices has become an international trend.
Since the application of fuel cells in vehicles is still in the exploratory phase from theory to practice, it is of great significance to improve the aerodynamic performance of radial compressors for FCSs, which has emerged as a crucial topic in this field.
The radial compressor required for FCVs has characteristics of a small flow rate and high total pressure ratio (up to 2.0–4.0) in a single stage, which has a narrow working area and tends to work along the surge boundary. By properly increasing the gas supply flow rate and gas supply pressure [
3,
4,
18], both the power density and efficiency of the FCS can be improved, as well as reduce the cost of the whole vehicle. Excessive gas supply pressures, however, can increase the power consumption and lower the output efficiency of the FCS. According to current research [
19], the FCSs consumption can be further reduced via isentropic efficiency promotion, which can be realized through the optimal design of the radial impeller or blade.
In the past 40 years, in the three-dimensional (3D) field, optimization of the blade design has focused mainly on two aspects: the optimization algorithm and the parametric optimization method. In terms of optimization algorithms, it has gone through the process from local optimization [
20] to global optimization. Omidi et al. [
21] used a hybrid method, comprising a genetic algorithm and a simulation package to realize the global optimization, simulate the function of a radial com impeller pressor, and evaluate the effects of losses in the impeller. In terms of the parametric optimization, blade parameterization methods are mainly divided into two categories: one is to fix the shape of each section and only change the position of the stacking line; the other is to reshape the geometry of each section, i.e., deform the pressure and suction surfaces, or directly change the mid-arc and thickness distribution [
22,
23]. Hehn et al. [
24] used three independently designed camber curves instead of ruled surfaces to optimize the blade geometry and aerodynamic analysis, which required 156 optimization variables. The combined geometric and aerodynamic analysis revealed that a forward-swept leading edge and a concave suction side at the tip of the leading edge are effective design features for reducing the shock wave strength and blade shape of the optimized compressor impeller and enable favorable impeller outlet flow. Hildebrandt et al. [
25] used 545 geometry parameters with the help of 3D-CFD and statistical correlations based on the linear Pearson and the ranked Spearman coefficients to evaluate the numerical aerodynamic analysis for the pressure slope optimization of a radial compressor impeller. Li et al. [
26] presented a 3-D multi-objective aerodynamic optimization method by integrating a self-adaptive, multi-objective differential evolutionary algorithm, 3-D blade parameterization method, and RANS. Control points on blade hub and shroud contours were selected as design variables. The total pressure ratio and isentropic efficiency were increased by 1.26% and 3.06%, respectively. Liu et al. [
27] redesigned a transonic radial compressor impeller by a multi-point, multi-objective optimization method. Camber curves of the blades and profiles at the tip and root sections of both main and splitter blades were parametrized. A genetic algorithm was used as the optimization method. The overall performances of baseline and optimum impellers were compared. Improvements in the total pressure ratio (by 5.3%) and isentropic efficiency (by 1.9%) were captured. Khalil Ekradi and Ali Madadi [
28] present a procedure for three-dimensional optimization of a transonic radial compressor impeller with splitter blades by integrating a 3D blade parameterization method, a genetic algorithm (GA), an artificial neural network, and a CFD solver. The isentropic efficiency is increased by 0.97% at the design point, and the total pressure ratio and mass flow rate are increased by 0.74% and 0.65%, respectively. These methods can alter the shape of the blade surface and achieve good design optimization results for the radial impeller. Still, the author believes they cannot be considered surface optimization methods due to the control variables based on curves rather than surfaces. In such methods, surface generation is dependent on the skinning process built into the mesh generation software, causing many design optimization variable problems and low optimization efficiency, or the surface is an inflexibility ruled blade.
Burguburu [
29] first proposed a semi-blade Bezier surface parameterization method that works with fewer control parameters and can attain high efficiency, surface smoothness, and better intuition, thus providing a new parameterization direction in the design optimization of axial flow compressor impeller. Cheng [
30] applied this method to optimize a single row rotor and single-stage transonic axial compressor impeller and achieved ideal optimization results. Huang et al. [
31] developed and proposed a full-blade surface parameterization method of an axial flow compressor impeller to compensate for the low degree of freedom of the semi-blade surface parameterization reformulation. The essence of these surface parameterization methods is to modify the original blade surface by the superposition of Bezier surfaces, which has a good construction convenience. However, the spatial morphology of radial compressor impeller blades with non-ruled surfaces is more distorted (greater inclination of the impeller inlet) compared with the shape of axial compressor blades. Furthermore, the flow inside of a turbomachine is viscous and compressible. These characters together with the complicated geometry of blades complicate the flow study. Therefore, the application of traditional parameterization methods based on surface superposition on radial compressor impeller blades is difficult.
In this paper, an improved MDOF surface parameterization method for the radial impeller is proposed and an effective automated aerodynamic optimization system for radial impeller blades is developed. The main blade and a splitter blade of the radial impeller blade are parametrically reshaped using the MDOF surface parameterization method, and the global optimization is carried out by using MIGA and the 3D CFD solver to explore the optimization performance after using the new approaches for the vehicle-mounted radial impeller, which is completely different from traditional parameterization methods.
2. MDOF Surface Parameterization Method for Radial Impeller Blades
Although computing power has developed rapidly in recent years, the current computing power is still not sufficient for many engineering applications in optimization. The “curse-of-dimensionality” problem is very difficult in the design optimization process that has to be solved due to the lack of computational power. In the design optimization of blades using the traditional parametric approach [
22,
23,
24,
25], the number of design variables increases geometrically with the number of blade rows, and the design space increases exponentially, making it difficult to obtain an optimized solution within a limited engineering time frame.
In recent years, Huang [
31] proposed a full-blade surface parameterization method, which considers the pressure and suction surface of the original blade as a whole, and then the variable value of each point on a Bezier surface is superimposed on the circumferential direction of the corresponding point of the original blade to form a new blade. Although this method can achieve the purpose of dimensionality reduction, it is not flexible to modify blade shape, and the blade optimization can often lead to blade thinning or even deformation to weak-strength blades. In particular, the surface superposition direction of the full-blade surface parameterization is circumferential, and the blade inclination at the inlet of the radial impeller is relatively large, so the modified impeller is easily unable to intersect with the casing line, thus making mesh generation and flow field calculation impossible.
In order to remedy this deficiency, this paper proposes the MDOF surface parameterization method, shown in
Figure 2, which is a parametric mapping and dimensionality reduction method that includes the following characteristics: firstly, it does well in maintaining the low-dimensional characteristics of the surface parameterization; secondly, it improves the flexibility and smoothness of the pressure and suction surface reshaping; thirdly, it maintains the blade thickness and mechanical strength in the process of optimization; finally, it superimposes Bezier surfaces on the suction and pressure surfaces of the original blade along their normal direction, respectively, and utilize the change in optimization variables in a one-dimensional direction to realize the 3D geometric deformation of blade profiles. The high order continuity of the Bezier surface [
32] can ensure that the surface smoothness of the optimized blade is not lower than that of the original blade, thus reducing the flow loss caused by the surface roughness and improving the quality and efficiency during blade manufacturing and processing [
33,
34,
35,
36].
In this paper, the fixed geometric parameters include the hub geometry data, inlet diameter, outlet diameter, blade axial length, blade radial height, and blade thickness. The variable geometry is the full 3D deformation of the main blade and splitter blade by the MDOF method. The process of the MDOF surface parameterization method is shown in
Figure 3, which can be described by the following steps:
Determination of the leading edge (L.E.) and trailing edge (T.E.) points of the blade by the monotonicity of transverse coordinates of the blade;
Encrypting the points of each cross-section by means of lateral interpolation;
Parameterizing the chord length of each section of the original blade. Since the Bezier surface is a unit mapping surface in the computational domain, to make each point of the original blade correspond to the Bezier surface it is necessary to parameterize the chord length of each section of the original blade. The parameterization method is expressed by Equations (1) and (2):
where
and
are the horizontal and vertical coordinates of the chord length parameterized, respectively,
I ∈ (1, Np) and Np refers to the number of points of each radial section,
j ∈ (1, Ns) and Ns refers to the total number of radial sections,
lc refers to the length of the
cth segment of the chord length of the
jth section in the radial direction,
Lj is the sum of the chord lengths of the
jth section in the radial direction,
lr refers to the length of the
rth segment of the radial length of the
ith section in chord direction, and
Li is the sum of the radial lengths of the segments of the
ith section in chord direction.
Generating the Bezier surface where the Bezier surface is defined by Equations (3)–(7).
In Equation (3),
is the coordinate value of each point on the Bezier surface where
,
,
),
are the control vertexes of the Bezier surface, for which there are a total of (
m + 1) × (
n + 1) control vertexes as variables,
and
are Bernstein basis functions calculated from Equations (4) and (5) (where
v and
u are two independent variables of the Bezier surface varying in the range [0,1]),
is calculated by Equation (6), and
is calculated by Equation (7);
Setting the variable value of the control vertexes of the Bezier surface and then calculating the variable value of each point on the Bezier surface and the original blade surface;
Calculating the unit normal vector of each point on the original blade surface;
Superimposing the optimized variable,
SZ, on the normal direction of the pressure surface, with the magnitude and direction change in the suction surface consistent with the pressure surface to finally generate a new blade as shown in
Figure 2.
Figure 4 shows the distribution of control vertexes for the MDOF surface parameterization method. The suction and pressure surfaces have the same number of control vertexes and distribution positions, i.e.,
with
m + 1 points in the
direction and
n + 1 points in the
direction. The number of distribution points in the
and
directions should be moderate. Excessive points will lead to an increase in the dimensionality of the optimization variables and mutual interference of different control vertexes on the aerodynamic performance due to the global support characteristic of the Bezier surface, which will reduce the optimization efficiency. In contrast, few points will result in insufficient optimization space and lead to poor optimization results.
The position of the red points of the leading edge of the pressure and suction surfaces are fixed to ensure first-order continuity of the leading-edge connection. The pressure surface control vertexes (green) are free moving points and the suction surface control vertexes (blue) follow the moving points. In order to ensure the mechanical strength of the blade and prevent it from being thinner or of wrong geometry, the following constraint is implemented in the optimization process: the variable value of the points on the suction surface must be consistent with the pressure surface, which enables a significant reduction in the number of optimization variables whilst maintaining relatively good flexibility in the modification.
Compared with traditional parameterization methods [
22,
23,
24,
25,
26,
27,
28,
29,
30,
31], which, for a single blade, can require hundreds of variables and full 3D optimization surfaces that are difficult to construct, the number of variables in the MDOF surface parameterization method is effectively reduced to 20, which successfully achieves dimensionality reduction in the parameterization method and easy construction of the surface. At the same time, the method can achieve better flexibility of reshaping, it takes into account the changes in pressure surface and suction surface, maintains the smoothness of the blade profile, increases the design optimization space, and provides a parametric basis for the flexible design of radial compressor impeller blades.