Research Article  Open Access
Aerodynamic Optimization Design of a Multistage Centrifugal Steam Turbine and Its OffDesign Performance Analysis
Abstract
Centrifugal turbine which has less land occupation, simple structure, and high aerodynamic efficiency is suitable to be used as small to medium size steam turbines or waste heat recovery plant. In this paper, onedimensional design of a multistage centrifugal steam turbine was performed by using inhouse onedimensional aerodynamic design program. In addition, threedimensional numerical simulation was also performed in order to analyze design and offdesign aerodynamic performance of the proposed centrifugal steam turbine. The results exhibit reasonable flow field and smooth streamline; the aerodynamic performance of the designed turbine meets our initial expectations. These results indicate that the onedimensional aerodynamic design program is reliable and effective. The offdesign aerodynamic performance of centrifugal steam turbine was analyzed, and the results show that the mass flow increases with the decrease of the pressure ratio at a constant speed, until the critical mass flow is reached. The efficiency curve with the pressure ratio has an optimum efficiency point. And the pressure ratio of the optimum efficiency agrees well with that of the onedimensional design. The shaft power decreases as the pressure ratio increases at a constant speed. Overall, the centrifugal turbine has a wide range and good offdesign aerodynamic performance.
1. Introduction
With the continuous reduction of fossil energy and the enhancement of people’s environment awareness, there has been increasing attention to the high efficiency utilization of energy, where the turbine is energy conversion component; if the turbine is improved effectively, it can help to improve the efficiency utilization of energy.
Centrifugal turbine is a new type of turbine engine which has many advantages. In the centrifugal turbine, the gas flows outward the center and the channel area of the flow path increase naturally as the fluid volume flow increases, during the expansion process. That meets the principle of aerodynamic and geometric matching. In addition, centrifugal turbine can be used to achieve multistage design [1] more easily than centripetal turbine. Thus it can avoid the limitation of supersonic flow. This is more useful for design conditions, especially for offdesign performance conditions.
For turbine as the key component, many scholars in this area have done some relevant researches, but there is not much research on centrifugal turbine. Domestically, mainly Jing and Peng had studied the pneumatic analysis of a rocket centrifugal turbine prototype starter and designed the modification turbine [2]; YinGe et al. [3] and Xin et al. [4, 5] had researched the singlestage centrifugal turbine design and its offdesign performance. Abroad, the preliminary hydrodynamic design of a small centrifugal turbine for the ORC was studied by Casati and others of the Delft University of Technology. They introduced the optimization based on the intermediate streamline method for evaluating turbine design and performance [6–10].
In this paper, onedimensional design method of centrifugal steam turbine is proposed by drawing on the conventional turbine [11]. And numerical simulation and optimization of multistage centrifugal turbine were studied by referring to the inlet and outlet thermal parameters of a small axial steam turbine.
2. Aerodynamic Design of a Steam Centrifugal Turbine
2.1. OneDimensional Design Method
A onedimensional design of centrifugal turbine has been developed, which is based on the onedimensional design method of conventional turbine. The onedimensional thermodynamic calculation program was developed by FORTRAN. The main design principles of centrifugal turbine onedimensional design program are as follows:(1)The expansion flow is assumed to be adiabatic, steady, and onedimensional in the cascade channel.(2)The properties of working fluids are obtained by calling Refpro 9.0, which is applicable to a variety of working fluid.(3)In order to simplify the design of the blade, the blade is designed to be of straight and constant height.(4)Stator and rotor velocity coefficients and are based on previous experience [12].(5)Each stage of the rotor absolute flow angle is 90 degrees.(6)The program is mainly composed of continuity and energy equations to achieve onedimensional design.
Figure 1 presents the design process of the centrifugal onedimensional aerodynamic program. The thermodynamic parameters of the centrifugal turbine, inlet stagnation temperature , inlet stagnation pressure , outlet pressure , mass flow rate , and rotation speed , are based on original conventional turbine. The other parameters, impeller diameter ratio (the influence of diameter ratio on the centrifugal turbine wheel efficiency is referred to in [3]; the optimum wheel efficiency calculated at the optimum degree of reaction and speed ratio increases as the diameter ratio decreases; ), each stage outlet flow angle of stator , the radial gap ( or , because the passage area along the flow path remains constant in axial flow turbine or the flow area reduces in radial flow turbine, so the radial gap can be large in traditional turbine, while, in centrifugal turbine, the flow area increases with the working fluid expansion. If the radial gap is too large, it will make the working fluid compressive in the gap, which is bad for turbine working. So the radial gap for centrifugal turbine should be small.), and stage number , are previously estimated. Iterative and screening methods are used to search the maximum wheel efficiency of the centrifugal turbine.
Figure 2 shows the overall stages diagram. The superscript represents the stagnant state. The second numbers of subscript are as follows: 1 represents the stator and 2 represents rotor. Taking the first stage as an example, the introduction of design process is as follows:(1)In the stator, the steam expands from state 0 to state 1. Lines 01 represent the actual expansion and lines 01 is the ideal expansion. In this process, the pressure energies are transferred to kinetic energies. Then the thermal parameters, velocity and geometry parameters can be calculated by using (2)Then, in the rotor, the steam continues to expand from state 1 to state 2. Lines 12 represent the actual expansion and lines 12 is the ideal expansion. The fluid kinetic energies are transferred to mechanical energies, which make the turbine outputs shaft power. The thermal parameters and the velocity triangles can be calculated by using(3)The rotor actual outlet density is compared with estimated density . If ( is a minimum), the assumption density is valid. Then the rotor outlet point is determined. If it is not satisfied, then reassume and repeat the above steps until the conditions are met.(4)Similar to the axial flow turbine, the relevant adiabatic efficiency is given by (3), which is referred to in [13].
The wheel efficiency is calculated at different speed ratio to select the optimum speed ratio and other corresponding parameters at maximum efficiency. Geometry and thermodynamic parameters are identified at the optimum speed ratio. The abovementioned steps are also adopted to design other stages.
2.2. OneDimensional Design Results
The onedimensional design program of the centrifugal turbine is used to design the multistage centrifugal turbine. The initial design thermal parameters are derived from a small axial turbine and summarized in Table 1.

Because there is no stage number information about the original axial turbine from the product information, the different stage number was tested to design the centrifugal turbine by using the onedimensional design program. The results show that if the stage number is 1 and 2, there is supersonic in the centrifugal turbine. As the stage number increased, there is no supersonic in the centrifugal turbine, but the size of the centrifugal turbine becomes more and more large. When the stage number is 3, and the diameter ratio is 1.12, the centrifugal turbine is subsonic, and it can satisfy the enthalpy drop as well as high efficiency. So the centrifugal turbine is designed with three stages.
Then, the same blade height and the absolute airflow angle of each stage are 90 degrees as the criteria, impeller diameter ratio is estimated to be 1.12 for each stage, each stage outlet flow angle of stator is assumed to be 12 degrees, and radial clearance between stator and rotor is set to be 2 mm. The optimal efficiency and reasonable structure are the objectives. The design parameters are calculated by using the abovementioned onedimensional calculation program, with iterative and screening methods. In this paper, the designed threestage centrifugal turbine just meets the required enthalpy drop and has a very high wheel efficiency. The main aerodynamic parameters of each stage are shown in Table 2, and the main aerodynamic parameters of overall centrifugal turbine are shown in Table 3. The geometry parameters are shown in Table 4, and the speed triangle data is shown in Table 5. The speed triangle schematic diagram is shown as Figure 3.




3. Airfoil Design and Optimization
The working fluid flows are outward the center in the centrifugal turbine. The passage area along the flow path increases naturally with the expansion of the working fluid, which makes the variety of specific volume matches the change of flow passage area. It can be seen that the centrifugal turbine is much superior to the conventional turbine from the structure. In the centrifugal turbine, the blade plane flow channel is expanding outward with the increase of the radius. It is obviously inappropriate to use the airfoil of conventional turbine at this time. So designing a suitable airfoil for centrifugal turbine is needed.
3.1. Parametric Expression of Airfoil
Angle and thickness design method is adopted to design airfoil, by using the inlet and outlet geometry angle, blade height, and leading and trailing edge diameters, which are based on onedimensional design results in Tables 3 and 4. On the BladeGen platform, twodimensional structure of the stator and rotor is identified. Stator and rotor of the first stage are an example, as shown in Figures 4(a) and 4(b). The blade surface includes four patches, namely, the leading edge and the trailing edge and the suction side and the pressure side. The mean camber line is controlled by cubic Bezier curve. The leading and trailing edges are both ellipse. Blade number is determined initially by referring to the relative pitch and expelling coefficient and the data is provided in [12]. The final numbers of stator and rotor for each stage are both 65.
(a) Initial profile of the firststage stator
(b) Initial profile of the firststage rotor
3.2. Blade Optimization for ThreeStage Centrifugal Steam Turbine
In the optimization process of blade, the blade number, the blade inlet and outlet geometry angles, and the leading and trailing edge thicknesses are the fixed parameters, and two control points of the tangent angle and two control points of the blade thickness are taken as the optimization parameters. Workbench is as the optimization software, and NLPQL algorithm is used as optimization method.
In the stator optimization process, the minimum loss coefficient is the object and stator back pressure is the constraint condition. Stators of threestage centrifugal turbine are optimized, respectively. The optimization mathematical model of stator is expressed as The design variables , , respectively, are the horizontal and vertical coordinates of control points of the tangent angle. and are the horizontal and vertical coordinates of control points of the blade thickness. A total of eight variables are involved in optimization (the stator of the first stage is an example, as shown in Figure 5). represents the pressure at the trailing edge of the stator, and represents the stator back pressure of onedimensional design ( represents stage number).
(a) Tangential angular distribution of medial camber line and control points selection of rotor blade
(b) The thickness distribution and control points selection of rotor blade
Then the rotor is added behind the optimized stator. The maximum shaft power is the object and rotor back pressure is the constraint condition. The optimization mathematical model of rotor is expressed as The same as the stator, the design variables and are the horizontal and vertical coordinate points for controlling the tangent angle of the rotor mean camber line, respectively. and are the horizontal and vertical coordinate points for controlling the thickness of the blade. A total of eight variables are involved in optimization (the rotor of the first stage is an example, as shown in Figure 6). represents the pressure at the trailing edge of the rotor, and represents the rotor back pressure of onedimensional design.
(a) Tangential angular distribution of medial camber line and control points selection of stator blade
(b) The thickness distribution and control points selection of stator blade
The optimization strategy and process are shown in Figures 7 and 8, respectively.
Figure 8 shows the blade optimization process. All of those were simulated automatically on the workbench. Firstly, the value range of optimization parameters (, , , and ), constraint condition, and object are set artificially. Then, Design Exploration selects a set of data in the value range of optimization parameters, automatically. The blade geometry and parameterization are obtained by geometry software. After that, the blade mesh is generated by TurboGrid, and centrifugal turbine is simulated by CFX. Finally, if the results satisfy the constraint and object, the parameter values of blade control points are the optimization values. If not, the Design Exploration will select another set of data in the value range of optimization parameters, automatically. Then repeat the above steps until the conditions are met. Since the optimized stages deviate from the design condition, when they are calculated together, the leading and trailing edge thicknesses of stator and rotor are needed to be changed slightly for local adjustment, then the centrifugal turbine can obtain a better performance.
Since the SST model is used in turbomachinery, in most cases, it requires to be very small (). So the layer mesh quality requirement of SST model is higher than epsilon model. If the SST uses the same mesh generated for epsilon, the layer mesh quality for SST is not good and cannot satisfy the requirement. The global size factors were changed to increase the mesh number, and the mesh quality was improved. So the grid number of SST is more than epsilon. It is about 8627,000, while the grid number of epsilon model is 4010,000. If the multistage centrifugal turbine simulation turbulence model was SST, a lot of computing memory will be needed; meanwhile a lot of calculating time will be cost. We found that the calculation results of SST and epsilon models are almost the same, as shown in Table 6, so there is no difference between SST model and epsilon model to simulate the multistage centrifugal turbine. In order to relax the demand on computer memory and to raise the efficiency when the computational grid number is very large, the epsilon model is used as the turbulence model to simulate the centrifugal turbine.

The optimized centrifugal turbine is simulated by CFX, based on NavierStokes equations. Stators and rotors use structured grid generated in TurboGrid, and the total grid number is about 4010,000; the blade global size factors of mesh are 1.2, 1.2, 1.25, 1.25, 1.25, and 1.2 for each stator and rotor, respectively. The computational model mesh used for simulation is shown in Figure 9. The boundary conditions are inlet total pressure, total temperature, back pressure, and adiabatic wall. The turbine rotation speed is 6500 rpm. A noslip boundary condition is applied at all the solid walls. The statorrotor interface is Frozen Rotor. Separate periodic conditions are used for rotor and stator regions. epsilon is used as the turbulence model, and the fluid is Water Ideal Gas.
(a) Computational domain for the simulation of stages
(b) Mesh distribution along the height and tangential direction of blade
The optimal parameters of stator and rotor are almost consistent with onedimensional aerodynamic design. The optimized stator and rotor of whole stages are both straight blades and the blade diagram is referred to in Figure 10.
(a) Optimized stators and rotors
(b) 3D model of stator and rotor blades
4. Results and Discussion
4.1. Design Condition Aerodynamic Performance
The numerical simulation results of the optimized centrifugal turbine are basically consistent with the onedimensional aerodynamic design, which shows that the onedimensional design program is reliable and effective. The overall performance data of centrifugal turbine is shown in Table 7. Comparison of performance data at each stage is shown in Table 8. Seen from Table 8, the performance data has some deviations at each stage. This is because the optimized stages influence each other during the matching process. So it results in some deviations in the performance data. And these deviations are gradually accumulated in the process of fluid flow. As the stage increases, the deviations are more. But, from the overall performance data seen in Table 7, the overall performance data is within 4% deviation. This is due to the fact that there are positive and negative deviations at each stage, which can offset each other. So the overall performance data of onedimensional design and simulation results is not quite different. From the total shaft power and wheel efficiency, the simulation results of overall centrifugal turbine are better than onedimensional design values. The threestage centrifugal turbine meets the design requirements and has good performance.


Figures 11–13 show the distribution of pressure, Mach number, and velocity streamline at the middle plane of the threestage centrifugal turbine, respectively. From the diagrams, the pressure distribution is uniform; the main flow in the impeller is the pressure flow. The flow field in the cascade channel is reasonable, and the streamline is smooth and the aerodynamic performance is in line with the expectation. As can be seen from Table 8, the stator and rotor outflow angles of the first and second stage have deviation from the onedimensional design value. From the analysis of the flow field diagram, there is a flow separation in the small area near the trailing edge; the fluid expands at the chamfered part of the stator outlet and the gap between the rotor and stator, which results in deflection of the flow angle. The flow loss caused by the separation and the deviation of the inflow angle caused by the flow deflection are not yet considered in the onedimensional design procedure. This is also the reason why the simulation results deviate from the onedimensional design values.
4.2. OffDesign Condition Aerodynamic Performance
Figures 14, 15, and 16 show the relationship between mass flow , efficiency , and shaft power at different rotation speeds with variable pressure ratio . The results indicate that as decreases, the mass flow increases until the critical flow is reached and then the mass flow remains constant, when the speed is constant. This is because the flow process appeared supersonic, and the flow reached its maximum. And the flow passage had the blocking phenomenon, but the turbine was still able to work. There is an optimal efficiency point on the efficiency with pressure ratio curve. At both sides of the optimal efficiency point, the efficiency decreases with the increase or decrease of . This is due to the change of the back pressure, and the enthalpy drop is changed, the speed ratio is deviated from the optimum value, the flow angel is deflected, and the blade surface occurs the flow separation. Those lead to the flow loss increase and efficiency decrease. For shaft power, the power decreases as π increases. Because of increase in back pressure, the enthalpy drop of the whole stage is reduced. Then the work capacity of working fluid is reduced. So the shaft power decreases with the increasing π. When the rotation speed changes, the optimum efficiency point moves in the direction of pressure ratio decreasing and the optimum efficiency value decreases with the increase of the rotation speed. But the efficiency trend with the pressure ratio is consistent at different speeds. For the mass flow, the flow characteristics at different rotation speeds are basically the same. It can be seen that the mass flow variation at different rotation speeds is the same and the change has little effect on mass flow. At different rotation speeds, the trends of shaft power curves are basically the same. But the curve of shaft power with the pressure ratio is steeper with the increase of rotation speed.
5. Conclusions
The threestage centrifugal turbine is simulated by CFX with Water Ideal Gas as working fluid, drawing on the conventional turbine aerodynamic design method and offdesign condition performance research method. Conclusions are as follows.
The numerical simulation results of the whole stages are basically consistent with the onedimensional design results, and the aerodynamic performance meets the expected requirements, which indicates the reliability and effectiveness of the centrifugal turbine design. Because each stage of centrifugal turbine is optimized, the simulated overall efficiency is 3.08% higher than the onedimensional design efficiency and the shaft power is 3.75% more than the onedimensional design shaft power. At the design condition, the streamline of threestage centrifugal steam turbine is smooth at cascade flow channel, the pressure distribution is uniform, and the flow field is reasonable.
At the offdesign condition, when the speed is constant, the pressure ratio π reaches the critical ratio and the mass flow reaches the maximum. If π continues to decrease, the maximum flow value remains the same, and the flow passage has the blocking phenomenon. But the turbine is still able to work. At the same π, the speed decreases; then the mass flow increases; but the impact is quite small.
At a constant speed, the efficiency is highest at optimum pressure ratio π. The pressure ratio of the optimum efficiency agrees well with that of the onedimensional design. When the speed decreases, the efficiency curve moves to the place where π increases, and the corresponding maximum efficiency increases.
The shaft power curve trend with pressure ratio is similar at different rotation speeds. But as the rotation speed increases, the power curve drops faster when the pressure ratio increases.
Overall, it can be seen that the centrifugal turbine has a wide range and good offdesign performance, from the numerical simulation results of the threestage centrifugal turbine and the analysis of the offdesign conditions.
Nomenclature
:  Temperature, K 
:  Pressure, Pa 
:  Enthalpy, kJ/kg 
:  Entropy, kJ/(K·kg) 
:  Rotational speed, r/min 
:  Mass flow, kg/s 
:  Diameter ratio 
:  Diameter, m 
:  Radius, m 
:  Absolute velocity, m/s 
:  Relative velocity, m/s 
:  Circumferential velocity, m/s 
:  Height, m 
:  Stage number 
:  Shaft power, kW 
:  Mach number. 
:  Density, kg/m^{3} 
:  Radial symmetry cone angle, ° 
:  Radial clearance between stator and rotor, m 
:  Absolute angle, ° 
:  Relative angle, ° 
:  Stator flow coefficient 
:  Rotor flow coefficient 
:  Speed ratio 
:  Wheel efficiency 
:  Reaction degree 
:  Pressure ratio. 
0:  First stage 
in:  Inlet 
out:  Outlet 
:  Stator 
:  Stage number. 
:  Stagnation condition. 
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by National Natural science Foundation of China no. 51536006 and Shanghai Science and Technology Committee with Grant no. 17060502300 and funded by Opening Project of Shanghai Key Laboratory of Multiphase Flow and Heat Transfer in Power Engineering no. 13DZ2260900.
References
 E. Coomes, R. Dodge, and D. Wilson, “Design of a high powerdensity ljüngstrom turbine using potassium as a working fluid,” in Proceedings of the 21st Intersociety Energy Conversion Engineering Conference, vol. 3, San Diego, CA, USA, 1986. View at: Google Scholar
 L. Jing and S. Peng, “Modification design and flowfield analysis of a centrifugal turbine stage with discrete nozzles,” Journal of Aerospace Power, vol. 23, no. 6, pp. 1047–1053, 2008. View at: Google Scholar
 L. YinGe, T. Xin, and L. XianQiao, “The thermal design and analysis of centrifugal turbine,” Journal of Engineering Thermophysics, vol. 37, no. 10, pp. 2103–2010, 2016. View at: Google Scholar
 T. Xin, L. YinGe, L. XianQiao et al., “Research on offdesign characteristics of centrifugal turbine,” Journal of Engineering Thermophysics, vol. 37, no. 6, pp. 1201–1208, 2016. View at: Google Scholar
 T. Xin, L. Hui, S. YanPing et al., “Design of the auxiliary components of the centrifugal turbine,” Journal of Engineering Thermophysics, vol. 0253231X, pp. 04–0000, 2017. View at: Google Scholar
 M. Pini, G. Persico, E. Casati et al., “Preliminary design of a centrifugal turbine for organic rankine cycle applications,” Engineering for Gas Turbines & Power, vol. 135, no. 4, 2013. View at: Google Scholar
 E. Casati, S. Vitale, M. Pini, G. Persico et al., “Centrifugal turbines for miniorganic rankine cycle power systems,” Journal of Engineering for Gas Turbines & Power, vol. 136, no. 136, Article ID 122607, 11 pages, 2014. View at: Google Scholar
 G. Persico, M. Pini, V. Dossena, and P. Gaetani, “Aerodynamic design and analysis of centrifugal turbine cascades,” in Proceedings of the ASME Turbo Expo: Turbine Technical Conference and Exposition, GT '13, vol. 6C, p. 12, San Antonio, Texas, Tex, USA, 2013. View at: Publisher Site  Google Scholar
 G. Persico, M. Pini, V. Dossena et al., “Aerodynamics of centrifugal turbine cascades,” Journal of Engineering for Gas Turbines & Power, vol. 137, no. 11, Article ID 112602, 2015. View at: Google Scholar
 M. Pini, A. Spinelli, G. Persico, and S. Rebay, “Consistent lookup table interpolation method for realgas flow simulations,” Computers & Fluids, vol. 107, pp. 178–188, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 L. YanSheng and L. GuiLin, Radial Turbine and Centrifugal Compressor, Machinery Industry Press, Bei Jing, China, 1987.
 W. XinJun and L. Liang, Principle of Steam Turbine, XI’an Jiaotong University Press, Xi'an, Shaanxi, China, 2014.
 S. L. Dixon and C. A. Hall, Fluid Mechanics and Thermodynamics of Turbomachinery, Elsevier Inc Press, 2010.
Copyright
Copyright © 2017 Hui Li and DianGui Huang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.