Next Article in Journal
Fed-Batch Sucrose Crystallization Model for the B Massecuite Vacuum Pan, Solution by Deterministic and Heuristic Methods
Next Article in Special Issue
Optimization Design of a Two-Vane Pump for Wastewater Treatment Using Machine-Learning-Based Surrogate Modeling
Previous Article in Journal
Effect of Regenerated Cellulose Fibers Derived from Black Oat on Functional Properties of PVA-Based Biocomposite Film
Previous Article in Special Issue
Multi-Condition Optimization of Cavitation Performance on a Double-Suction Centrifugal Pump Based on ANN and NSGA-II
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Applications of an Improved Aerodynamic Optimization Method on a Low Reynolds Number Cascade

School of Mechanical Engineering, Shanghai Jiao Tong University, Shanghai 200240, China
*
Author to whom correspondence should be addressed.
Processes 2020, 8(9), 1150; https://doi.org/10.3390/pr8091150
Submission received: 18 August 2020 / Revised: 5 September 2020 / Accepted: 9 September 2020 / Published: 14 September 2020

Abstract

:
The effect of cascade aerodynamic optimization on turbomachinery design is very significant. However, for most traditional cascade optimization methods, aerodynamic parameters are considered as boundary conditions and rarely directly used as the optimization variables to realize optimization. Given this problem, this paper proposes an improved cascade aerodynamic optimization method in which an incidence angle and nine geometric parameters are used to parameterize the cascade and one modified optimization algorithm is adopted to find the cascade with the optimal aerodynamic performance. The improved parameterization approach is based on the Non-Uniform Rational B-Splines (NURBS) method, the camber line superposing thickness distribution molding (CLSTDM) method, and the plane cascade design method. To rapidly and effectively find the cascade with the largest average lift-drag ratio within a certain range of incidence angles, modified particle swarm optimization combined with the modified very fast simulated annealing algorithm (PSO-MVFSA) is adopted. To verify the feasibility of the method, a cascade with NACA4412 and a practical cascade are optimized. It is found that the average lift-drag ratios of two optimal performance cascades are respectively increased by 13.38% and 15.21% in comparison to those of two original cascades. Meanwhile, through optimizing the practical cascade of the Blade D500, under different volume flow rates, the pressure coefficient of the optimized cascade is increased by an average of more than 6.12% compared to that of the prototype, and the average efficiency is increased by 11.15%. Therefore, this improved aerodynamic optimization method is reliable and feasible for the performance improvement of cascades with a low Reynolds number.

1. Introduction

Blade design is of great importance to the efficiency and properties of turbomachinery. Achieving the aerodynamic design of a blade is a very complex and arduous task due to complicated flow phenomena and the interactions among various parameters [1]. A small geometric change of one blade can lead to a deterioration of the aerodynamic performance of the whole machine [2]. In general, the cascade design method is one of the most popular design methods employed for turbomachinery [3,4]. In this method, the blade is stacked by the sections on the different radiuses and shown in Figure 1a. Additionally, the section is projected onto the plane to form the cascade, as shown in Figure 1b. Therefore, the blade performance is affected by the cascade design. With the rapid development of the computer technique, more and more parameterization methods and optimization algorithms have been proposed and used in the design process of a cascade. This means that the time required to design a new cascade is becoming shorter and the aerodynamic performance of the new cascade is becoming better.
A good parameterization method can not only use fewer design variables to describe an airfoil accurately, but also rapidly re-construct one airfoil in the optimization process [5,6]. The parameterization methods are usually divided into two categories: The constructive method and the deformative method. Each method has been continuously improved by many researchers [7]. The deformative parameterization method is the simpler of the two methods. In this method, a standard airfoil is deformed to generate one new airfoil, in order to satisfy a certain condition. The Hicks-Henne function [8], radial basis function [9], Bezier function [10], B-Spline function [11], and Non-Uniform Rational B-Splines (NURBS) function [12,13,14] are usually used as deformative functions to generate a new airfoil based on the original airfoil. In particular, NURBS [12] is the most popular function due to its ability of local control and its conics description over the curve. However, these deformative functions are used to generate one new airfoil based on the point coordinates of an airfoil. To relate the airfoil shape to the airfoil geometric feature parameters, the camber line superposing thickness distribution molding (CLSTDM) method [15] was proposed. In this method, several airfoil geometric parameters are used to parameterize the half-thickness distribution curve and the mean camber curve through two polynomials. Then, these two curves are coupled to form a whole airfoil. The feasibility of this method has been proved by several works [15,16,17]. In the CLSTDM method [15,16], the blade contours described by many coordinate points can be transformed into functions controlled by several parameterized variables. Moreover, it is convenient for designers to use this method to parameterize one blade based on their experience. Nowadays, it is widely used in the optimization of turbomachinery. However, the aerodynamic parameters are not considered in the parameterization and still act as the boundary condition of flow field computing.
To rapidly find the airfoil with the optimal aerodynamic property, many intelligent optimization technologies have been used in the process. Among all of the published optimization algorithms, the genetic algorithm (GA) [18,19,20,21] and the simulated annealing (SA) algorithm [10,11,22], as two traditional intelligent optimization algorithms, have been widely used in airfoil optimization. They aim to find the airfoil with the optimal performance precisely. However, much time is required to complete the search, which leads to a poor computational efficiency [23,24]. A new intelligent algorithm, known as particle swarm optimization (PSO) and proposed by Kenney and Eberhart [25], can be used to solve this problem. PSO is a population-based, self-adaptive searching optimization method. The principle of this method is based on animal social behaviors, such as birds’ migration. However, it is easy to become trapped into a local extreme value or converged to precocity by the standard PSO. Therefore, researchers began looking for improved methods to solve this problem. Shi [26] used a linearly decreasing inertia weight to balance the global and local searching character. However, the local searching capability of this method was weak. Simultaneously, it is hard to predict the maximum iteration number. Clerc [27] set a constriction factor determined by two learning factors to cancel the boundary limits of velocity, and to balance the global and local searching capability. Hu [28] adopted a stochastic inertia weight to replace the linearly decreasing inertia weight. This method could accelerate the convergence velocity and avoid being trapped into a local best solution. However, these methods still have the risk of local convergence.
Most approaches to parameterization of the cascade have been used with only the help of geometric feature parameters, and the aerodynamic parameters were not referred to. The efficiency of optimization was thus negatively affected. Therefore, in this paper, considering the cascade aerodynamic characters, one aerodynamic parameterization approach for a low Reynolds number cascade is proposed. To find the airfoil with the optimal performance during a certain range of incidence angles, a modified PSO-MVFSA algorithm is studied. Furthermore, two cases, such as the cascade with NACA4412 and the blade of FAN D500, are selected to verify the feasibility of the improved parameterization and optimization method.

2. Aerodynamic Parameterization Method

2.1. CLSTDM-NURBS Method

In this paper, for the CLSTDM method, the pressure side and suction side are obtained through the camber curve superposing the half-thickness distribution curve. This method is used to describe an airfoil and is shown in Figure 2. Due to its good ability to design a complex geometry, the two-order NURBS function defined by Equation (1) is used to describe the camber curve and the half-thickness distribution curve.
{ C ( u ) = i = 0 n B 2 , i ( u ) ω i Q i B 2 , i ω i ( u ) B 2 , i ( u ) = 2 ! i ! ( 2 i ) ! u i ( 1 u ) 2 i ,
where u is the knot vector, n is the order of NURBS ( n = 2 ), i is the mark of the control points ( i = 0 , 1 , 2 ), C ( u ) is the coordinate of the point of the fitting curve parameterized by the NURBS function, B 2 , i ( u ) is the Bernstain function, ω i represents the weight coefficients ( ω 0 = 1 ,   ω 2 = 1 , ω 1 = ω ), and Q i is the control point. To solve the two-order NURBS function, the De Boor algorithm [29], which provides a fast and numerically stable way of finding a point on a B-spline curve with the given u in the domain, is adopted and programed.
For camber parameterization, the camber curve is divided into two two-order NURBS curves (solid line), and these two NURBS curves are respectively controlled by several geometric control points ( Q 0 , Q 1 , Q 2 , Q 3 , Q 4 ) shown in Figure 3. Five geometric feature parameters of an airfoil, such as the chord line L = 1.0 , the leading edge angle χ 1 , the trailing edge angle χ 2 , and the coordinate of the maximum camber point ( B x ,   B y ) , are adopted to derive the coordinates of the geometric control points of NURBS. To ensure the continuous property of two NURBS curves at the location of the maximum camber point Q 2 , two lines Q 1 Q 2 ¯ and Q 2 Q 3 ¯ need to be collinear. Moreover, all of the geometric control points of NURBS are derived as follows in Equation (2).
[ Q 0 Q 1 Q 2 Q 3 Q 4 ] = [ ( 0 , 0 ) ( B y / tan ( χ 1 ) , B x ) ( B x , B y ) ( L B y / tan ( χ 2 ) , B y ) ( L , 0 ) ]
For parameterization of the half-thickness distribution curve, the curve is divided into three parts, including the leading edge half-thickness, the middle half-thickness, and the trailing edge half-thickness, as shown in Figure 4. The leading edge (LE) part is described by one two-order NURBS curve, the middle part by two two-order curves, and the trailing edge (TE) part by one two-order NURBS curve. Four NURBS curves are controlled by nine geometric control points ( P 0 , P 1 , P 2 , P 3 , P 4 , P 5 , P 6 , P 7 , P 8 ), which are derived by seven geometric feature parameters of the airfoil, such as the leading edge radius R 1 , the trailing edge radius R 2 , the chord line L = 1.0 , the thickness gradient angles α 1 , α 2 , and the coordinate of the maximum half-thickness point ( T x , T y ) . As is the case for the camber curve, the half-thickness distribution curve also requires continuity. Therefore, at the three geometric control points P 2 , P 4 , P 6 , it is necessary to maintain collinearity for the relative lines. Utilizing geometric principles, these control point coordinates are shown in Equation (3).
[ P 0 P 1 P 2 P 3 P 4 P 5 P 6 P 7 P 8 P 9 ] = [ ( 0 , 0 ) ( 0 , R 1 cos ( α 1 ) R 1 tan ( α 1 ) ( 1 sin ( α 1 ) ) ) tan ( α 1 ) ( R 1 ( 1 sin ( α 1 ) ) , R 1 cos ( α 1 ) ) ( ( T y R 1 cos ( α 1 ) ) / tan ( α 1 ) + R 1 ( 1 sin ( α 1 ) ) , T y ) ( T x , T y ) ( ( T y R 2 cos ( α 2 ) ) / tan ( α 2 ) + L R 2 ( 1 sin ( α 2 ) ) , T y ) ( L R 2 ( 1 sin ( α 2 ) ) , R 2 cos ( α 2 ) ) ( L , R 2 cos ( α 2 ) R 2 tan ( α 2 ) ( 1 sin ( α 2 ) ) ) ( L , 0 ) ]

2.2. Improved Aerodynamic Parameterization

Unlike conventional parameterization approaches in which the airfoil is only parameterized with the geometric feature parameters [30,31,32], an improved aerodynamic parametrization approach combining the plane cascade design method and the CLSTDM-NURBS is proposed, in which the incidence angle, i , and nine geometric parameters are used as control variables.
For one specific cascade, it can be assumed for the airflow condition that the inlet flow angle β 1 is constant and along the tangential line. In this case, a change of the incidence angle can cause a change of the geometric inlet angle of the cascade, and the variation of the incidence angle Δ i is equal to that of the geometric inlet angle Δ β 1 A . From the definition of the geometric inlet angle, it is clear that the slope of the mean camber curve at the leading point is changed with the changing of the geometric inlet angle. This means that the variation of the incidence angle Δ i has an indirect influence on the modified value of the leading edge angle χ 1 . The relationship is shown in Equation (4).
χ 1 = χ 1 Δ i
Simultaneously, the deviation angle δ can also be affected by a change of the incidence angle. In order to determine the deviation angle, one semi-rational formula of the deviation angle was used, based on the plane cascade experiments by Howell [33,34], which was suitable for a low Reynolds number cascade and only applicable in the application conditions ( τ = 0.7 ~ 2.0 , T y = 0.05 ~ 0.12 , and B x = 0.4 ~ 0.5 ). The definition formulas of the incidence angle i and the deviation angle δ have been proposed in the cascade design process [3,35].
τ = t L
δ = β 2 A β 2
In this paper, one special equation combining a semi-empirical formula [33,34] is shown in Equation (7), where φ is equal to 0.5 for the rotor blade. Therefore, if the aerodynamic and geometric parameters i , τ , β 1 , β 2 , and B x are given, the geometric outlet angle β 2 A can be calculated by Equation (7). Additionally, the revised trailing edge angle χ 2 can be obtained by Equation (8). However, in Equation (7), it is found that the cascade solidity τ and the outlet flow angle β 2 cannot be determined.
( β 2 A β 2 ) = [ 0.23 ( 2 B x ) 2 0.002 β 2 + 0.18 ] ( β 2 A i β 1 ) ( 1 τ ) φ
χ 2 = χ 2 Δ β 2 A
In order to obtain the abovementioned two parameters, the diffusion factor D related to the cascade solidity τ , the inlet flow angle β 1 , and the outlet flow angle β 2 , is introduced as a constraint. This coefficient can be used to control the aerodynamic load of the cascade. The definition of the diffusion factor is shown in Equation (9).
D = ( 1 sin β 1 sin β 2 ) + sin β 1 2 τ ( c t g β 1 c t g β 2 )
Therefore, during the process of parameterization, the leading edge angle and trailing edge angle shown in Figure 3 are independent variables and no longer depend on the incidence angle. In this way, the control variables of cascade parameterization are determined, such as B x , B y , T x , T y , R 1 , R 2 , α 1 , α 2 , and i .

3. The Improved Airfoil Aerodynamic Optimization Method

3.1. Modified PSO-MVFSA

With the extensive use of standard particle swarm optimization (PSO), some drawbacks have been found, such as the local extreme minimum and the precocity. To solve these problems, many works have been published [26,27,28]. Some only adjusted the change of inertia weight to avoid being trapped into a local optimal solution. Moreover, some used two learning factors to balance the local searching and global searching. In this paper, the control variables of cascade parameterization are grouped as a particle. Additionally, the best particle position X i , j corresponding to the optimal cascade is obtained by the modified PSO, as shown in Equation (10) [27].
{ V i , j 1 = ϕ · [ ω V i , j 0 + c 1 γ 1 ( P i , j X i , j 0 ) + c 2 γ 2 ( P g , j X i , j 0 ) ] ϕ = 2 | 2 ( c 1 + c 2 ) ( c 1 + c 2 ) 2 4 ( c 1 + c 2 ) | ω = μ min + ( μ max μ min ) γ 4 + σ · γ 3 X i , j 1 = X i , j 0 + V i , j 1 ,
where V is the particle velocity, subscript i , j is the particle sequence, ω is the stochastic inertia weight [28], ϕ is the constriction factor [27], c 1 , c 2 represents two learning factors, γ 1 , γ 2 represents the random number uniformly distributed in ( 0 , 1 ) , P i , j is the best position in its flight history, P g , j is the best position in the particle swarms, μ max is the maximum stochastic inertia weight, μ min is the minimum stochastic inertia weight, σ is the variance of the stochastic inertia weight, and γ 3 , γ 4 is the random number of the standard normal distribution.
However, it is very difficult to judge whether the results are the local optimal solutions or the global solutions when only using the modified PSO. Due to the probability of the simulated annealing (SA) algorithm jumping out of the local value, it can be used to solve the problem. Furthermore, considering that the low searching velocity of the standard SA algorithm can lead to a greater time consumption, modified very fast simulated annealing (MVFSA) [36] is adopted to search for the global optimal minimum. Since the Cauchy distribution depending on temperature was better than the Gaussian distribution, the coefficient of variable perturbation y i j based on the Cauchy distribution has been redefined by Equation (11). In order to improve the efficiency further, the random probability of the acceptance, P r , has been redefined based on the Boltzmann–Gibbs distribution, as shown in Equation (12).
y i j = T max exp ( c k 1 / N 1 ) · sgn ( u 0.5 ) [ ( 1 + 1 T ) | 2 u 1 | 1 ] ,
P r = [ 1 ( 1 h ) ( E ( M i j ) E ( M i j ) ) / T ] 1 / ( 1 h ) ,
where u is a random number ranging from 0 to 1 , T max is the initial simulated high temperature, N 1 is the syllogism coefficient, k is the number of the marked annealing stage, c is a given constant value, M i j is the disturbed variable, M i j is the undisturbed variable, and h is a real number (set h = 0.5 ). Furthermore, E ( M i j ) is adopted to evaluate the current energy of the particle.

3.2. Verification of Modified PSO-MVFSA

Due to its local optimums and premature results, the Rastrigin function is usually used as a test function. The definition of the function is shown in Equation (13).
f ( x ) = j = 0 D x j 2 10 · cos ( 2 π x j ) + 10       x j [ 20 , 20 ]
In theory, when the independent variable x j is equal to 0 , the global minimum optimum of the function is f ( x j ) = 0 . In Figure 5, three different solutions of the function are respectively shown. When the independent variable range is [ 20 , 20 ] , the function looks smooth and monotonic, as shown in Figure 5a. With the increase of the resolution, many small peaks and valleys are observed in Figure 5b,c. This means that the function has many widespread local minima.
To prove the feasibility and accuracy of the modified PSO-MVFSA, two other PSO-based optimization algorithms, as two comparison algorithms, are adopted to search for the optimal minimum of the Rastrigin function. As the function is two-dimensional, each particle consists of two variables for three kinds of PSO algorithms. The initial particle number, the total iteration number of PSO, and the total iteration number of MVFSA are set to 30, 80, and 40, respectively. Some other coefficients are listed in Table 1, such as the positive value α , the Markov chain length J , and so on. By means of these searching optimization methods, the optimal results are obtained and shown in Table 2. It can be seen that the optimal function value obtained by PG-PSO is smaller than that of the standard PSO. However, the computing time is longer than that of the standard PSO. For the modified PSO-MVFSA, not only its time consumption, but also its ability for minimum searching, are optimal in comparison with the other two algorithms.

3.3. The Improved Airfoil Aerodynamic Optimization

3.3.1. Fitness Function

For improved cascade aerodynamic optimization, nine control variables are used to parameterize the cascade and grouped as a particle. In all flights, the selection of the particle with the best position through the optimum target is very important. The lift-drag ratio of the airfoil is related to its capability regarding the power output and aerodynamic loss [16]. Therefore, the optimization proposition is set up as shown in Equation (14).
{ f ( S 1 , S 2 S 10 ) = i = 1 n ( C L / C D ) i / n s . t . χ 1 > arctan ( B y / B x )    χ 2 > arctan ( B y / ( 1 B x ) ) α 1 > arctan ( T y / T x )     α 2 > arctan ( T y / ( 1 T x ) ) 0.4 D 0.6 ,
where S i is the control variable, C L / C D is the lift-drag ratio, ( B x , B y ) is the coordinate of the maximum camber point, ( T x , T y ) is the coordinate of the maximum half-thickness point, α 1 , α 2 represent the LE and TE thickness gradient angles, χ 1 , χ 2 represent the LE and TE angles, and n is the number of incidence angles. In this paper, the average lift-drag ratio of the airfoil is evaluated under all airflow incidence angles. Additionally, the fitness function can be obtained as shown in Equation (15).
F = c 1 f ( S 1 , S 2 S 10 ) - ( C L / C D ) r e f ( C L / C D ) r e f + c 2 α 1 arctan ( B y / B x ) arctan ( B y / B x ) c 3 α 2 + arctan ( B y / ( 1 B x ) ) arctan ( B y / ( 1 B x ) )      + c 4 X 1 arctan ( T y / T x ) arctan ( T y / T x ) + c 5 X 2 arctan ( T y / ( 1 T x ) ) arctan ( T y / ( 1 T x ) ) + c 6 ( D i D )
where c i is the weighting factor and subscript r e f presents the original airfoil.

3.3.2. Aerodynamic Optimization Process

During the evaluation of airfoil aerodynamic performance two-dimensional computational fluid dynamics (2D CFD) code which utilized the standard k ω turbulence model [37,38] to solve the Reynolds Average Navier-Stokes (RANS ) equations was used, and it was only used to simulate the low Reynolds number airfoils [23,39]. Due to the fact that the lift-drag ratio of one airfoil could be obtained quickly by this code, this 2D CFD code was integrated into MATLAB as a fast solver to evaluate the fitness value of each particle. The whole aerodynamic optimization flowchart is shown in Figure 6. Moreover, the detailed process is introduced in the following steps:
(a)
Inputting coordinate points of one airfoil, and setting coefficients of the modified PSO-MVFSA;
(b)
Selecting one angle from the sets of incidence angles and nine geometric variables as an initial particle;
(c)
Conducting perturbation of nine geometric variables of the initial particle to generate the initial particle swarm with one incidence angle, based on the super Latin square method; constructing airfoil swarm by the improved parameterization method; and evaluating the airfoil fitness value by Equation (15);
(d)
Finding the best previous position of the particle and the best position of the swarm, re-calculating the velocity and the position of each particle by adopting Equation (10), and re-calculating the fitness values of the new swarms;
(e)
If Δ E > 0 ( Δ E = E ( M i j ) E ( M i j ) is the error function), the data related to the particle are unchanged; if not, replace the data and particles with new data and particles;
(f)
Repeating steps (d) and (e) until the total iteration of the modified PSO is reached;
(g)
Obtaining the swarm particles ( S ) i that satisfy [ F o p t F ( ( S ) i ) ] / F o p t < 10 % ;
(h)
Putting the swarm particles ( S ) i into the optimization process of MVFSA, conducting perturbation, and re-evaluating the fitness by Equation (10);
(i)
If Δ E > 0 , the corresponding particle is preserved; if not, the corresponding particle as a basic particle is re-disturbed and re-evaluated;
(j)
If Δ E > 0 , the re-disturbed particle is retained; if not, the particle is accepted with the acceptance probability equation;
(k)
Repeating steps (g) (h) (i) until the total iteration of MVFSA is reached;
(l)
Outputting the particle with the largest function F from the preserved particle swarm at different incidence angles respectively;
(m)
Selecting the best of the optimal particles by MVFSA as the final optimal particle.

4. Applications in Cascade Optimization

4.1. Optimization of a Cascade with an NACA4412 Profile

NACA4412 is adopted as the original basic profile of the cascade and evolved into an advanced airfoil by this improved method. The initial particle number, the total iteration number of the modified PSO, and the total iteration number of MVFSA are set to 30, 100, and 50, respectively. Each cascade is described by the improved aerodynamic parameterization method. The average lift-drag ratio of the airfoil of the cascade is calculated by the flow solver. Moreover, the constant airflow boundary condition is set so that the Reynolds number is as high as 1 × 10 5 , the Mach number is 0.1 , and the solidity of the cascade is 1.5 . Moreover, some other coefficients of the modified PSO-MVFSA are shown in Table 3.
In actual airflow conditions, the incidence angle of the cascade is not constant and limited to a certain range. Therefore, the aim of multi-point optimization is not to obtain the best aerodynamic performance under one constant incidence angle, but to reach the best whole aerodynamic performance in the whole working range, referred to in Equation (12). In this paper, the incidence angle is uniformly distributed from 0° to 15° by a step of 1.0°. The aerodynamic performance of each cascade parameterized by one incidence angle and eight geometric variables is calculated by the CFD simulation. Based on Equation (13), the fitness value corresponded to each cascade is figured out. Then, the global minimum value is found by the user of the improved PSO-MVFSA method. NACA4412 airfoil and the optimal airfoil are shown in Figure 7, in which the blue curve is the optimal airfoil and the red curve denotes the original airfoil. From this figure, it can be observed that the suction side and pressure side of the optimized airfoil are changed in terms of geometry in comparison with that of NACA4412.
To evaluate the aerodynamic performance of the cascade, the pressure coefficient C p is adopted, which is defined by Equation (16).
C p = ( p p 1 2 ρ u 2 ) = p p p 0 p ,
where p , p 0 , p are the current pressure, the stagnation pressure, and the inlet airflow pressure, respectively; ρ is the density; and u is the inlet airflow velocity.
The increasing curvature of the suction side close to LE can lead to the intensifying of the velocity increasing and the pressure decreasing. Then, the pressure coefficient of the suction side of the optimized C p near LE is larger than that of the original cascade. However, it is due to the geometric change of the pressure side close to LE that the decreasing of the velocity and the increasing of the pressure are alleviated. These phenomena can be observed in Figure 8. In the figure, the detailed pressure distributions on the suction and pressure side are shown, which respectively correspond to four incidence angle conditions of i = 0 ° , i = 3 ° , i = 6 ° , and i = 12 ° . It is also clear that the optimized airfoil surface pressure distribution presented by the blue curve is better than that of the original airfoil described by the red curve under each condition. Therefore, the pressure differences between the pressure side and suction side of the optimized airfoil are larger than those of the original airfoil. The cascade performances under the incidence angles ranging from 0° to 15° are shown in Figure 9. The average lift-drag ratio of the optimized cascade is increased by 13.38% in comparison with that of the cascade with NACA4412.
For further analysis, CFD simulation software CFX was used to calculate the cascade performance. The velocity contours of the original and optimized cascades at incidence angles of i = 0 ° and i = 12 ° are shown in Figure 10. It is clear that the velocity difference between the suction side and pressure side of the optimized cascade is bigger than that of NACA4412. Therefore, the diffusion factor D of the optimized cascade can be increased. It can also be observed that the airflow separation point on the suction side of the optimized cascade moves backward along the airfoil in comparison with that of NACA4412 and the area of the wake becomes smaller than that of NACA4412. Therefore, it can be considered that the aerodynamic load of the cascade is increased and the aerodynamic loss is controlled with the help of the improved aerodynamic optimization algorithm.

4.2. Blade D500 Optimization

Blade D500 is used in an axial-flow fan for an evaporator system. It is a kind of low Reynolds number 3D blade. In this work, in order to verify the feasibility of the improved method applied to the practical blade design, the cascade at the 50% radius of Blades D500 was selected and optimized. Two performance parameters, including the pressure coefficient C p defined by Equation (16) and the efficiency η defined by Equation (17), were used to evaluate the aerodynamic performance of Blade D500. In order to obtain the aerodynamic performance parameters, CFX was used to calculate the airflow field of Blade D500.
η = Q v · P N ,
where Q v is the volume flow rate, P is the static pressure increase, and N the aerodynamic power.

4.2.1. Validation of the CFD Simulation Based on Experiments

The experiments of the axial fan with Blade D500 shown in Figure 11 were conducted in the Key Lab for Power Machinery and Engineering of SJTU. In order to match the practical situation, the rotor was mounted on the guard grill, which was connected with a short bell month at the inlet of the tube. The diameter of the tube was increased to 600 mm to avoid destroying the flow field at the rotor outlet and simulate the real situation in the evaporator system as much as possible. The pressure probes were mounted at points A and B, by which the flow rate and the pressure increase could be calculated. Under eight volume flow rates, the pressure increases were calculated.
Under the same conditions as those used in the experiments, CFX was adopted to simulate the aerodynamic performance of the axial fan. The grids consisted of the structured hexahedral meshes shown in Figure 12 that were generated by TurboGrid. The total number of nodes was 1,260,626 for the full channel after studying the grid independence, and the minimum value of the mesh quality was 0.3. Due to the excellent ability to simulate the flow with a fiercely adverse pressure gradient of the standard k ω model, it was selected as the turbulent model to simulate the flow field. The mass flow and static pressure were set at the inlet and outlet, respectively.
After a series of experiments and CFD simulations of the axial fan with Blade D500, the pressure coefficients were obtained and are shown in Figure 13. Based on Figure 13, it could be found that the pressure coefficients C s p calculated by CFD are slightly higher than those calculated by experiments. However, the trends of the two curves are very similar. In Figure 14, the efficiency is compared for values obtained by calculations and experiments. It can be observed that the relative error between the two series of average efficiencies is under 8.35%. Considering the abovementioned results, the CFD simulation can be used to feasibly estimate the aerodynamic performance of an axial fan.

4.2.2. Optimization of Blade D500

For aerodynamic optimization of the cascade at the 50% radius of Blade D500, the airflow condition was set as Re = 1 × 10 5 , M a = 0.1 , i = 2 ° 10 ° , D = 0.4 0.6 . After a series of optimization iterations, the airfoil of the optimal cascade was obtained and is shown in Figure 15, as well as the airfoil of the original cascade. It can be observed that the airfoil of the optimized cascade is different from that of the original cascade. These geometric changes of the optimal cascade can lead to an increase of the curvature of the suction side, while the pressure side is changed a little. Additionally, aerodynamic performance comparisons are shown in Figure 16. It could be found that the pressure differences of the optimized cascade are larger than those of the original airfoil along the chord line within the range of the incidence angle. After a series of simulation calculations, it could be found that the average lift-drag ratio of the optimized cascade is increased by 15.21% in comparison with that of the original cascade. Therefore, it can be considered that the aerodynamic performance of the optimized cascade is better than that of the original one.
To evaluate the feasibility of the improved aerodynamic optimization method, the optimized 3D Blade D500 and the original 3D Blade D500 were simulated by CFD under the same airflow condition. Moreover, the simulation results of two blades were compared. The streamlines of the suction surfaces of the two blades are shown in Figure 17. From the figure, it is clear that the streamline between the middle section and the hub of the original blade is not very good, and a large turbulence loss must have occurred in this region. In contrast, the streamline in the same region of the optimized blade becomes very smooth. The pressure coefficients C s p of the two fans are shown in Figure 18. In this figure, it can be seen that the pressure coefficient of the optimized fan is better than that of the original fan, and the average pressure coefficient is increased by 6.12%. Meanwhile, from Figure 19, it can be observed that the efficiency of the optimized fan is higher than that of the original one, and averagely increased by 11.15%. Therefore, from the above analysis, it can be considered that the aerodynamic performance of optimized Blade D500 is improved and the improved aerodynamic optimization method is feasible.

5. Conclusions

In this paper, an improved optimization method especially applied to a cascade with a low Reynolds number is proposed. In this method, the incidence angle and eight geometric parameters are used as the control variables altogether. Meanwhile, the modified PSO-MVFSA algorithm is adopted to optimize the objective cascades, and two cascade cases, such as NACA4412 and Blade D500, are selected to testify the method. Some valuable results are as follows:
(a)
Since the aerodynamic parameter, such as the incidence angle, is considered as one of the control variables, the relationship between the geometry of the airfoil and the aerodynamic performance of the cascade is learnt. Therefore, during the whole optimization process, an improvement of the aerodynamic performance can give rise to a direct modification of the geometry so that the optimization becomes more targeted and more efficient;
(b)
In this study, particular effort was devoted to designing a fitness function which is suitable for optimizing a cascade with a low Reynolds number. Furthermore, the combination of PSO and MVFSA succeeded in increasing the optimization efficiency and avoided the local optimal to reach a global solution, which was verified by the Rastrigin function and two cascade cases;
(c)
Based on the analysis of the results from the two cascade cases, such as NACA4412 and Blade D500, it was demonstrated that the average lift-drag coefficient of the optimized cascade was improved, whilst the drag coefficient was kept at a low level. Therefore, it can be considered that the modified PSO-MVFSA can be adopted as an efficient and robust optimizer to solve the problems of multi-variable optimization confronted in cascade design.

Author Contributions

S.Z. and B.Y., investigation, methodology, software, and writing—original draft preparation. H.X. and M.S., validation. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

This research was supported by the National Science and Technology Major Project of China (2017-II-0006-0019).

Conflicts of Interest

The authors declare that they have no conflict of interest with regards to this work. They declare that they do not have any commercial or associative interests that represent a conflict of interest in connection with the work submitted.

Abbreviations

NURBSNon-Uniform Rational B-Splines
CLSTDMCamber Line Superposing Thickness Distribution Molding
PSOParticle Swarm Optimization
PSO-MVFSAParticle Swarm Optimization-Modified Very Fast Simulated Annealing

Nomenclature

B x Vertical coordinate of the maximum camber point
B y Horizontal coordinate of the maximum camber point
C p Pressure coefficient
C L Lift coefficient
C D Drag coefficient
D Diffusion factor
L Aerodynamic chord length
M Variable
N 1 Syllogism coefficient
N Shaft power
R 1 LE radius
R 2 TE radius
P Control points of the thickness distribution curve
Q Control points of the camber curve
Q v Volume flow rate
T Temperature
χ 1 LE angle
χ 2 TE angle
T x Vertical coordinate of the maximum half thickness
T y Horizontal coordinate of the maximum half thickness
P r Acceptance probability
α 1 Thickness gradient angle of LE
α 2 Thickness gradient angle of TE
β 1 A Geometric inlet angle
β 2 A Geometric outlet angle
β 1 Inlet flow angle
β 2 Outlet flow angle
y i j Coefficient of variable perturbation
ω Weight
i Incidence angle
μ Stochastic inertia weight
λ Constriction factor
τ Cascade solidity
γ Random number in ( 0 , 1 )
η Pressure efficiency
c Learning factor
( x , y ) Cartesian coordinates

References

  1. Li, Z.; Zheng, X. Review of design optimization methods for turbomachinery aerodynamics. Prog. Aerosp. Sci. 2017, 93, 1–23. [Google Scholar] [CrossRef]
  2. Garzon, V.E.; Darmofal, D.L. Impact of geometric variability on axial compressor performance. J. Turbomach. 2003, 125, 692–703. [Google Scholar] [CrossRef] [Green Version]
  3. Eck, B. Ventilatoren; Springer: Berlin/Heidelberg, Germany, 1991. [Google Scholar]
  4. Batyaev, E.A.; Kurzin, V.B. Method of optimization of blade shapes in aerodynamic design of the fan cascade. J. Appl. Mech. Tech. Phys. 2002, 43, 701–705. [Google Scholar] [CrossRef]
  5. Wu, H.-Y.; Yang, S.; Liu, F.; Tsai, H.-M. Comparisons of three geometric representations of airfoils for aerodynamic optimization. In Proceedings of the 16th AIAA Computational Fluid Dynamics Conference, Orlando, FL, USA, 23–26 June 2003. [Google Scholar]
  6. Samareh, J.A. Survey of shape parameterization techniques for high-fidelity multidisciplinary shape optimization. AIAA J. 2001, 39, 877–884. [Google Scholar] [CrossRef]
  7. Mahmood, S.M.H.; Turner, G.M.; Siddappaji, K. Flow characteristics of an optimized axial compressor rotor using smooth design parameters. In Proceedings of the ASME Turbo Expo 2016, GT2016, Seoul, Korea, 13–17 June 2016. [Google Scholar] [CrossRef]
  8. Hicks, R.M.; Henne, P.A. Wing design by numerical optimization. J. Aircr. 1978, 15, 407–412. [Google Scholar] [CrossRef]
  9. Buhmann, M.D. Radial basis functions. New Algorithms Macromol. Simul. 2004, 37, 31–65. [Google Scholar] [CrossRef]
  10. Liu, J.-L. A novel Taguchi-simulated annealing method and its application to airfoil design optimization. In Proceedings of the 17th AIAA Computational Fluid Dynamics Conference, Toronto, ON, Canada, 6–9 June 2005. [Google Scholar] [CrossRef]
  11. Chao, S.-M.; Whang, A.J.-W.; Chou, C.-H.; Su, W.-S.; Hsieh, T.-H. Optimization of a total internal reflection lens by using a hybrid Taguchi-simulated annealing algorithm. Opt. Rev. 2014, 21, 153–161. [Google Scholar] [CrossRef]
  12. Ghaly, W.S.; Mengistu, T.T. Optimal geometric representation of turbomachinery cascades using NURBS. Inverse Probl. Eng. 2003, 11, 359–373. [Google Scholar] [CrossRef]
  13. Sonoda, T.; Yamaguchi, Y.; Arima, T.; Olhofer, M.; Sendhoff, B.; Schreiber, H.-A. Advanced high turning compressor airfoils for low reynolds number condition: Part 1—Design and optimization. In Proceedings of the ASME Turbo Expo 2003, Collocated with the 2003 International Joint Power Generation Conference, Atlanta, GA, USA, 16–19 June 2003; Volume 6, pp. 437–450. [Google Scholar] [CrossRef]
  14. Schreiber, H.-A.; Steinert, W.; Sonoda, T.; Arima, T. Erratum: “Advanced High-Turning Compressor Airfoils for Low Reynolds Number Condition—Part II: Experimental and Numerical Analysis” [Journal of Turbomachinery, 2004, 126(4), pp. 482–492]. J. Turbomach. 2005, 127, 646. [Google Scholar] [CrossRef] [Green Version]
  15. Chen, N.; Zhang, H.; Ning, F.; Xu, Y.; Huang, W. An effective turbine blade parameterization and aerodynamic optimization procedure using an improved response surface method. In Proceedings of the ASME Turbo Expo 2006: Power for Land, Sea, and Air, Barcelona, Spain, 8–11 May 2006; Volume 6, pp. 1169–1180. [Google Scholar] [CrossRef]
  16. Chen, N.; Zhang, H.; Xu, Y.; Huang, W. Blade parameterization and aerodynamic design optimization for a 3D transonic compressor rotor. J. Therm. Sci. 2007, 16, 105–114. [Google Scholar] [CrossRef]
  17. Mazaheri, K.; Khayatzadeh, P.; Nezhad, S.T. Airfoil Shape Optimization Using Improved Simple Genetic Algorithm (ISGA). In Proceedings of the 5th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics, Sun City, South Africa, 1–4 July 2007. [Google Scholar]
  18. Deb, K.; Pratap, A.; Agarwal, S.; Meyarivan, T. A fast and elitist multiobjective genetic algorithm: NSGA-II. IEEE Trans. Evol. Comput. 2002, 6, 182–197. [Google Scholar] [CrossRef] [Green Version]
  19. Goldberg, D.E.; Korb, B.; Deb, K. Messy genetic algorithms: Motivation, analysis, and first results. Complex Syst. 1989, 3, 493–530. [Google Scholar]
  20. Goldberg, D.E. Genetic Algorithms in Search Optimization and Machine Learning; Addison-Wesley: Boston, MA, USA, 1989. [Google Scholar]
  21. De Oliveira, A.; Lorena, L. A constructive genetic algorithm for gate matrix layout problems. IEEE Trans. Comput. Des. Integr. Circuits Syst. 2002, 21, 969–974. [Google Scholar] [CrossRef] [Green Version]
  22. Gopalaramasubramaniyan, G.; Kumar, V.S. Parameter optimization of the sheet metal shearing process in the manufacturing of leaf spring assembly using the grey Taguchi method and simulated annealing algorithm. Adv. Mater. Res. 2011, 314, 2458–2463. [Google Scholar] [CrossRef]
  23. Yang, B.; Xu, Q.; He, L.; Zhao, L.H.; Gu, C.G.; Ren, P. A novel global optimization algorithm and its application to airfoil optimization. J. Turbomach. 2015, 137, 041011. [Google Scholar] [CrossRef]
  24. Chandrasekar, K.; Ramana, N. Performance comparison of GA, DE, PSO and SA approaches in enhancement of total transfer capability using FACTS devices. J. Electr. Eng. Technol. 2012, 7, 493–500. [Google Scholar] [CrossRef] [Green Version]
  25. Kennedy, J.; Eberhart, R. Particle swarm optimization. In Proceedings of the ICNN’95—International Conference on Neural Networks, Perth, Australia, 27 November–1 December 2002; Institute of Electrical and Electronics Engineers (IEEE): Cracow, Poland, 2002; Volume 4, pp. 1942–1948. [Google Scholar]
  26. Shi, Y.; Eberhart, R.C. Empirical study of particle swarm optimization. In Proceedings of the 1999 Congress on Evolutionary Computation-CEC99 (Cat. No. 99TH8406), Washington, DC, USA, 6–9 July 1999; Institute of Electrical and Electronics Engineers (IEEE): Cracow, Poland, 2003; pp. 1945–1950. [Google Scholar]
  27. Clerc, M.; Kennedy, J. The particle swarm—Explosion, stability, and convergence in a multidimensional complex space. IEEE Trans. Evol. Comput. 2002, 6, 58–73. [Google Scholar] [CrossRef] [Green Version]
  28. Zhang, J.; Fu, F.-W.; Chang, Z.-L. An upper bound on the size of m-ary t-sEC/AUED codes. J. China Univ. Posts Telecommun. 2006, 13, 95–97. [Google Scholar] [CrossRef]
  29. Smith, P.W. A practical guide to splines (Carl de Boor). SIAM Rev. 1980, 22, 520–521. [Google Scholar] [CrossRef]
  30. Lee, S.W.; Kwon, O.J. Robust airfoil shape optimization using design for six sigma. J. Aircr. 2006, 43, 843–846. [Google Scholar] [CrossRef]
  31. Vavalle, A.; Qin, N. Iterative response surface based optimization scheme for transonic airfoil design. J. Aircr. 2007, 44, 365–376. [Google Scholar] [CrossRef]
  32. Grasso, F. Usage of numerical optimization in wind turbine airfoil design. J. Aircr. 2011, 48, 248–255. [Google Scholar] [CrossRef] [Green Version]
  33. Howell, A.R. Fluid dynamics of axial compressors. Proc. Inst. Mech. Eng. 1945, 153, 441–452. [Google Scholar] [CrossRef]
  34. Howell, A.R. Flow in Cascades. In Aerodynamics of Turbines and Compressors (HAS-1); Princeton University Press: London, UK, 1964. [Google Scholar]
  35. Üçer, A.S.; Stow, P.; Hirsch, C. Fluid Mechanics, Thermodynamics of Turbomachinery; Pergamon Press: Oxford, UK, 1978. [Google Scholar]
  36. Ingber, L.; Rosen, B. Genetic algorithms and very fast simulated reannealing: A comparison. Math. Comput. Model. 1992, 16, 87–100. [Google Scholar] [CrossRef] [Green Version]
  37. Wilcox, D.C. Comparison of two-equation turbulence models for boundary layers with pressure gradient. AIAA J. 1993, 31, 1414–1421. [Google Scholar] [CrossRef]
  38. Menter, F.R. Performance of popular turbulence model for attached and separated adverse pressure gradient flows. AIAA J. 1992, 30, 2066–2072. [Google Scholar] [CrossRef]
  39. Yang, B. A New Blade Design Scheme for Reversible Axial Flow Fan & Research on the Combined Cascades. Ph.D. Thesis, Shanghai Jiao Tong University, Shanghai, China, 2001. [Google Scholar]
Figure 1. An explanation of blades and a cascade.
Figure 1. An explanation of blades and a cascade.
Processes 08 01150 g001
Figure 2. Airfoil parameterization.
Figure 2. Airfoil parameterization.
Processes 08 01150 g002
Figure 3. Camber curve parameterization.
Figure 3. Camber curve parameterization.
Processes 08 01150 g003
Figure 4. Half-thickness distribution curve parameterization.
Figure 4. Half-thickness distribution curve parameterization.
Processes 08 01150 g004
Figure 5. Rastrigin function: (a) x i [ 20 , 20 ] ; (b) x i [ 5 , 5 ] ; and (c) x i [ 0.5 , 0.5 ] .
Figure 5. Rastrigin function: (a) x i [ 20 , 20 ] ; (b) x i [ 5 , 5 ] ; and (c) x i [ 0.5 , 0.5 ] .
Processes 08 01150 g005
Figure 6. Flow chart of the whole aerodynamic optimization process.
Figure 6. Flow chart of the whole aerodynamic optimization process.
Processes 08 01150 g006
Figure 7. Geometry comparison: NACA4412 airfoil and optimal airfoil.
Figure 7. Geometry comparison: NACA4412 airfoil and optimal airfoil.
Processes 08 01150 g007
Figure 8. Surface pressure distribution comparison in multi-point optimization.
Figure 8. Surface pressure distribution comparison in multi-point optimization.
Processes 08 01150 g008
Figure 9. Airfoil performance of the average lift-drag ratio (cal., Re = 1.0 × 10 5 ,   M a = 0.1 ).
Figure 9. Airfoil performance of the average lift-drag ratio (cal., Re = 1.0 × 10 5 ,   M a = 0.1 ).
Processes 08 01150 g009
Figure 10. Velocity contours.
Figure 10. Velocity contours.
Processes 08 01150 g010
Figure 11. The experiment of Blade D500.
Figure 11. The experiment of Blade D500.
Processes 08 01150 g011
Figure 12. Grid of the axial fan.
Figure 12. Grid of the axial fan.
Processes 08 01150 g012
Figure 13. Static pressure coefficients calculated by CFD and the experiment.
Figure 13. Static pressure coefficients calculated by CFD and the experiment.
Processes 08 01150 g013
Figure 14. Average efficiencies calculated by CFD and the experiment.
Figure 14. Average efficiencies calculated by CFD and the experiment.
Processes 08 01150 g014
Figure 15. Geometric comparisons of the optimized and original cascades.
Figure 15. Geometric comparisons of the optimized and original cascades.
Processes 08 01150 g015
Figure 16. Surface pressure distributions.
Figure 16. Surface pressure distributions.
Processes 08 01150 g016
Figure 17. Surface streamline comparisons of Blade D500.
Figure 17. Surface streamline comparisons of Blade D500.
Processes 08 01150 g017
Figure 18. Pressure coefficient comparisons of Blade D500.
Figure 18. Pressure coefficient comparisons of Blade D500.
Processes 08 01150 g018
Figure 19. Efficiency comparisons of D500.
Figure 19. Efficiency comparisons of D500.
Processes 08 01150 g019
Table 1. Coefficients of three algorithms based on particle swarm optimization (PSO).
Table 1. Coefficients of three algorithms based on particle swarm optimization (PSO).
Maximum Stochastic Inertia Weight,
μ max
Minimum Stochastic Inertia Weight,
μ min
Variance Stochastic Inertia Weight,
σ
Learning Factor,
c 1
Learning Factor,
c 2
0.80.40.22.251.85
The Initial Simulated High Temperature,
T max
The Final Cooling Simulated Temperature,
T min
The Syllogism Coefficient,
N 1
The Positive Value,
α
Markov Chain Length,
J
300.000190.815
Table 2. Results of three algorithms based on PSO.
Table 2. Results of three algorithms based on PSO.
AlgorithmsComputing Time
T (Seconds)
Optimal Particle
x o p t
Optimal Function Value
f ( x o p t )
Standard PSO [22]1.08942(0.001225, −0.000958)0.000481
PG-PSO [21]1.14177(−0.001495, 0.0000077)0.000443
PSO -MVFSA0.898254(0.000177, −0.000476)0.000051
Table 3. Coefficients of modified PSO-modified very fast simulated annealing (MVFSA).
Table 3. Coefficients of modified PSO-modified very fast simulated annealing (MVFSA).
Maximum Stochastic
Inertia Weight,
μ max
Minimum Stochastic
Inertia Weight,
μ min
Variance Stochastic Inertia Weight,
σ
Learning
Factor,
c 1
Learning
Factor,
c 2
0.80.40.22.251.85
The Initial Simulated High Temperature,
T max
The Final Simulated Cooling Temperature,
T min
The Syllogism Coefficient,
N 1
The Positive Value,
α
Markov Chain Length,
J
400.000190.820

Share and Cite

MDPI and ACS Style

Zhang, S.; Yang, B.; Xie, H.; Song, M. Applications of an Improved Aerodynamic Optimization Method on a Low Reynolds Number Cascade. Processes 2020, 8, 1150. https://doi.org/10.3390/pr8091150

AMA Style

Zhang S, Yang B, Xie H, Song M. Applications of an Improved Aerodynamic Optimization Method on a Low Reynolds Number Cascade. Processes. 2020; 8(9):1150. https://doi.org/10.3390/pr8091150

Chicago/Turabian Style

Zhang, Shuyi, Bo Yang, Hong Xie, and Moru Song. 2020. "Applications of an Improved Aerodynamic Optimization Method on a Low Reynolds Number Cascade" Processes 8, no. 9: 1150. https://doi.org/10.3390/pr8091150

APA Style

Zhang, S., Yang, B., Xie, H., & Song, M. (2020). Applications of an Improved Aerodynamic Optimization Method on a Low Reynolds Number Cascade. Processes, 8(9), 1150. https://doi.org/10.3390/pr8091150

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop