Next Article in Journal
An Experimental Study on Flexural Behaviors of Reinforced Concrete Member Replaced Heavyweight Waste Glass as Fine Aggregate under Cyclic Loading
Previous Article in Journal
Generation of Complete SAR Geometric Distortion Maps Based on DEM and Neighbor Gradient Algorithm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Detailed Design and Aerodynamic Performance Analysis of a Radial-Inflow Turbine

by
Amir Karimi Noughabi
1 and
Shervin Sammak
2,*
1
Department of Aerospace Engineering, Amirkabir University of Technology, 424 Hafez Ave, Tehran 15875, Iran
2
Center for Research Computing, University of Pittsburgh, 4420 Bayard St, Pittsburgh, PA 15213, USA
*
Author to whom correspondence should be addressed.
Submission received: 16 October 2018 / Accepted: 16 October 2018 / Published: 9 November 2018
(This article belongs to the Section Energy Science and Technology)

Abstract

:
Radial-inflow turbines (RI turbines) are critical components of air cycle machines (ACM), utilized in environment conditioning systems (ECS). In this work, a process is presented for the full design of a radial-inflow turbine, including volute, nozzle blade row and rotor blade row. The preliminary and detailed aerodynamic design approaches are explained in detail. A mean-line (1D) analysis method coupled with the preliminary design is employed to evaluate the aerodynamic performance of the RI turbine. Then, ANSYS CFX is used to perform unsteady 3D simulations of the designed RI turbine at design and off-design conditions to predict the aerodynamic performance of the turbine for future optimization.

1. Introduction

Compact air cycle machines (ACM) have been widely employed in many industrial applications, such as environment conditioning systems (ECS) and air cycle refrigeration. The coefficient of performance (COP) of the cooling cycles changes dramatically by efficiency and pressure ratio of the ACM [1,2]. The radial-inflow (RI) turbine is an indispensable member of the ACM, which is regularly used in compact air cooling systems because of its reliability, simplicity and robust construction, fast response, various working fluids (such as organic working fluids), low emission and low manufacturing cost [3,4,5]. Thus, the performance of the RI turbine is essential for these types of systems.
The aerodynamic design of the radial-inflow turbines is efficiently accomplished in well-known steps. First, a preliminary design is performed to calculate the basic stage component dimensions. An effective preliminary design provides important advantages in terms of quality and productivity for RI turbines and costs about 50% of the engineering time in the design process of a radial-inflow turbine [6]. Aungier [7] explains that the second step in the design of radial inflow turbine is the detailed design of the rotor and nozzle blades, according to the geometry dimensions obtained by the preliminary design (first step). In this step, the nozzle blade profile and the pressure and suction sides of the rotor blades are accurately determined. These detailed geometries can be ready for manufacturing. Rohlik [8] and Glassman [9] presented an analytical method to compute the geometrical dimensions of an RI turbine for maximum efficiency based on specific velocity. Rohlik considered boundary layer loss (profile loss) in nozzle and rotor and tip clearance loss in the design process. He studied the variations in turbine efficiency related to the rotor tip-radius ratio, the exit flow angle at nozzle and the ratio of nozzle blade height to rotor-inlet diameter. Glassman [9] provided a computer program to design RI turbines based on power, mass flow rate, inlet temperature and pressure. The design variables included a stator-exit angle, rotor-exit-tip to rotor-inlet radius ratio, rotor-exit-hub to tip radius ratio and distribution of rotor-exit tangential velocity. The program was able to compute the basic dimensions of the turbine stage, temperature and pressure at inlet and outlet sections of each component. Rodgers [10] developed an advanced program to design a small, cooled radial turbine, which operated at high inlet turbine temperature. This program was able to determine the rotor blade camber line in addition to the capabilities of the program developed by Glassman. Whitfield [11] presented a non-dimensional approach based on non-dimensional power to minimize Mach number at the rotor inlet and outlet so that the passage losses are reduced. The non-dimensional turbine design can then be transformed to the dimensional design by applying inlet stagnation conditions (total temperature and pressure) and mass flow rate. Ventura et al. [12] also provided an automated mean-line (1D) approach to the preliminary design of RI turbines. Authors used a brute-force search algorithm to select dimensions based on non-dimensional performance and geometry characteristics. Ebaid et al. [13] employed an optimization algorithm during preliminary design to remove selective open parameters such as rotor axial length. They also presented a detailed approach for nozzleless volute design. However, all aforementioned works did not describe an approach to achieve the second step in radial-inflow design. Although commercial software such as ANSYS-Vista RTD and RITAL do not provide any information about the rotor detailed design, stator geometry, stator profile and volute dimensional specifications, designers usually use them for the preliminary design [7,14,15]. In other similar works with preliminary design of the RI turbine, we cannot find a detailed design approach for all components of a typical RI turbine (Volute, stator and rotor) [16,17,18,19,20].
In this study, we redesign the radial-inflow turbine applied in the air cycle machines at Amikabir University of Technology to achieve higher efficiency such that the coefficient of performance of ACM is improved. To achieve this goal, first, a preliminary design approach has been developed based on the various available methods. Second, a method is applied to the detailed aerodynamic design of the turbine rotor blades based on Aungier [7] and then the detailed design of volute is established. After finalizing the turbine geometry, unsteady 3D simulations using ANSYS CFX are performed to compute the aerodynamic performances of the turbine, including the mass flow rate map, efficiency map, temperature drop and streamlines through the rotor and nozzle.

2. Design Process

The design process will be accomplished based on the predefined parameters such as total inlet temperature, total inlet pressure, mass flow rate, total-to-static pressure ratio, and specific speed (or number of rotor revolution). Equation (1) provides a relation between specific speed and rotor rotational speed. In this equation, Q 5 is the volume flow rate at rotor exit and Δ H i d is the isentropic expansion from the stage inlet total conditions to rotor exit static conditions. The value of these parameters for the RI turbine installed on the ACM at the moment are presented at Table 1.
n s = ω Q 5 ( Δ H i d ) 0.75
Aungier [7] provided the relations presented in Equations (2) and (3) to compute other basic design parameters such as total-to-static velocity ratio, ν t s , and total-to-static efficiency, η t s . The present program is able to take both of these parameters from the designer and compute them by Equations (2) and (3).
ν t s = 0.737 n s 0.2
η t s = 0.87 1.07 ( n s 0.55 ) 2 0.5 ( n s 0.55 ) 3
As a general rule, the preferred range of values for the specific speed is 0.45 to 0.75. However, the other design constraints such as Mach number, speed, mechanical integrity and size can impose the specific speed. In this work, the specific speed was determined by rotational speed due to mechanical integrity and sizing adaptation between the designed turbine and the previous turbine.
Figure 1 illustrates the flowchart of the design process. In preliminary design, a mean-line analysis method is used to compute the performance of the turbine designed by the preliminary design code based on the input value. In an iterating process, the final dimensions of the RI turbine are defined to achieve the design input parameters. The turbine rotor is finalized in the preliminary design method, and then exported into the detailed design code. In this step, the distribution of the blade thickness and angle at camber surface are computed and the suction and pressure surfaces of rotor blades are generated to be employed in a more accurate turbine analysis, such as quasi-3D method and CFD. In the final step of the design process, ANSYS-CFX software is applied to simulate flow in the designed turbine. This workflow enables the opportunity to modify the turbine model and therefore improve turbine performance.

3. Preliminary Design

The preliminary design procedure addresses the basic stage components defining the inlet volute, nozzle row, and rotor. A vaneless passage is automatically considered before and after nozzle row. The geometric and aerodynamic parameters are computed in the inlet and exit of all components shown with subscripts 1 to 5 following the station illustrated in Figure 2.

3.1. Rotor Preliminary Design

The velocity diagram is illustrated in Figure 3. In this figure, “C” and “W” are absolute and relative velocity, respectively. The blade and flow angle are computed related to the tangential direction.
In rotor design, it is assumed that the inlet blade angle is constrained to be 90°. Since the turbine is designed for optimum performance, the exit absolute tangential velocity is considered to be zero. On the other hand, it is a reasonable assumption that the absolute Mach number at rotor exit has a relative low value in a well-designed turbine. By these assumptions, the rotor tip radius and rotor inlet tangential velocity can be computed by the following equations and using Equations (4)–(6).
U 4 = ν t s C 0 t s     w h e r e     C 0 t s = 2 Δ H i d
r 4 = U 4 / ω  
C θ 4 = U 4 η t s / ( 2 ν t s 2 )
Rohlik [4] recommended rotor inlet absolute flow angle as a function specific speed presented in Equation (7). Glassman [9] suggested Equation (8) to estimate the number of rotor blades based on Jamieson expression.
α 4 = 10.8 + 14.2 n s 2
N R = π 30 ( 110 α 2 ) tan α 2
Aungier [7] provided the useful default relations to compute the hub radius, the ratio of meridional velocity at rotor exit to rotor inlet meridional velocity and the rotor axial length (Equations (9)–(11)). He also suggested the default equations to estimate the rotor blade thickness at the inlet and discharge stations. However, in this work, the rotor blade thickness was considered constant at all points of the rotor blade to simplify the final machining process. Though the hub radius, the blade thickness and the axial length are influenced by aerodynamic and mechanical considerations, it is more effective to use the default value for these parameters in the preliminary design. These parameters can be optimized according to aerodynamic and structural analysis.
r h 5 = 0.185 r 4
C m 5 / C m 4 = 1 + 5 ( b 4 / r 4 ) 2
Δ z R = 1.5 ( r s 5 r h 5 )
A well-designed rotor should satisfy a rotor exit relative Mach number less than unity. By this criteria, the mean blade pitch at the rotor exit and the rotor mean throat width could be computed by Equations (12) and (13).
s 5 = 2 π r 5 / N R
o 5 = s 5 C m 5 / W 5
In some other research, the hub and shroud contours were assumed to be elliptical and circular arcs, respectively [9,12]. However, in this work, the hub contour is considered the largest circular arc, which is compatible with the rotor dimensions to minimize passage curvature effects. Since the value of the rotor axial length was less than the value of the rotor tip ratio, a linear segment is added at the rotor inlet (Figure 4). The rotor shroud contour is defined by a power law relation to obtain reasonable passage area distribution, which is shown in Equations (14) and (15) [7].
r s = r s 5 + ( r 4 r s 5 ) ξ n
ξ = ( z s z 5 ) / ( Δ z R b 4 )
The value of n should be selected in the range of 2 through 9 so that the mean passage area should be average of the rotor area at inlet and exit stations. In this study, this parameter is considered to be equal to 7.
Two different types of rotor blade usually are used in RI turbines: Straight-line element (three dimensional) and radial element blades. Straight-line element blades are constructed by connecting corresponding blade camber-line. These blades are easy to design and manufacture and usually provide more efficient blades than radial element blades. Here the rotor blades are designed based on straight-line element presented in [7,21].

3.2. Nozzle Preliminary Design

The nozzle blade profile is usually designed and optimized according to flow velocity data [7,22]. However, it is necessary to define an initial profile for the nozzle blades. Though it is applicable to use NACA 4-digit and 5-digit airfoils to design and modify the nozzle blades in an iterative process, we have employed a general blade thickness distribution on a parabolic-arc camberline according to [23]:
x 2 + c 2 a b x y + ( c 2 a ) 2 4 b 2 y 2 c x c 2 4 a c 4 b y = 0
Equation (16) includes a general parabolic-arc camberline. Since a symmetric nozzle blade will be easier to design and manufacture, Equation (16) can be simplified to a straight line (x = 0) by assuming a = 0.5 c . The blade thickness distribution is adopted according to:
t =   t r e f + [ t m a x t r e f ] ξ e
t r e f = t L E + [ t T E t L E ] ( x / d )
where
{ ξ = x / d           x d ξ = ( c x ) / ( c d )     x > d
e = 0.4 d / c [ 0.95 ( 1 x c ) ( 1 ξ ) + 0.05 ]
Default values presented in Table 2 were considered for the constants in the Equations (17)–(20) according to Aungier [7].
Watanabe et al. [24] recommend a minimum ratio of the nozzle exit radius to the rotor inlet radius given by Equation (21). The tangential velocity, meridional velocity and other flow data at the nozzle exit can then be computed by conservation of angular velocity, mass balance and equation of state. In addition, the nozzle passage width (blade height) is assumed to be constant.
r 3 / r 4 = 1 + 2 b 4 sin α 4 / r 4
The nozzle throat width is defined by the sine rule for subsonic and choked flow exit according to:
{ o s 3 = sin α 3   M < 1 o s 3 = ρ 3 C m 3 ( ρ * a * )                                                 M = 1
To set the nozzle setting angle, γ 3 , an iterative process will be used based on the nozzle throat width, the number of nozzle blade, and exit blade pitch to chord ratio. Simpson et al. [25] and Glassman [9] suggested selecting 0.75 and 0.8 for the nozzle exit blade pitch to chord ratio, respectively. The number of nozzle blade should be the lowest value, which satisfies Equation (23) [7]:
4 s 3 sin ( β 2 α 3 ) c sin β 2 [ 1 + r 3 sin α 3 / ( r 2 sin α 2 ) ] 1
where α 2 and β 2 are absolute flow angle and blade angle at the nozzle inlet station, respectively.

3.3. Volute Design

The volute sizing depends on both the aerodynamic performances and the overall size of the stage in terms of maximum axial and radial dimensions. Since there are no dimensional restrictions for the RI turbine in the current work, an external volute was chosen to accelerate the flow to a higher discharge. The cross section of the volute is usually selected according to convenience rather than fluid dynamics principles so that the cross section shape was considered circular. Based on the design process of RI turbines presented by Aungier [7,21] and Ebaid et al. [13], a vaneless passage was considered between the nozzle inlet and the volute exit stations with the inlet-to-discharge radius ratio of 1.05 ( r 2 = 1.05 r 2 ). The volute cross section is illustrated in Figure 5.
The θ inlet volute area A 1 , and mean radius, r 1 , can be computed by applying the conservation of angular velocity ( r 1 C 1 = r 2 C θ 2 ), mass balance ( m ˙ = ρ 1 C 1 A 1 ) and thermodynamic equations in an iterative process. The area and mean radius at the other azimuth angle, θ , (Figure 5) are defined by:
A θ = ( 3 4 π + 1 ) ( r θ r 2 ) 2 = ( ρ 2 ρ θ ) θ b 2 r θ tan ( α 2 )

3.4. Mean-Line Analysis

A mean-line analysis was employed to estimate the RI turbine aerodynamic performance. The boundary layer loss (profile loss) is considered for all components of the RI turbine according to Equation (25). The detailed information to compute profile loss is presented in [7].
Y P = Δ P t P t e x , i d P e x = 2 Θ + Δ 2 ( 1 Δ ) 2
P t e x , i d and P e x are the ideal total pressure and actual static pressure at the discharge station of each turbine component. Total normalized boundary layer momentum thickness, Θ , and boundary layer thickness, Δ , are defined by Aungier. A summary of other losses considered for each component at mean-line analysis is shown in Table 3.

3.5. Preliminary Design Outputs

Table 4 shows the final geometric dimensions of the turbine stage computed by an iterative process between the preliminary design code and the mean-line analysis code. These parameters will then be used in rotor detailed design code to finalize the radial-inflow turbine geometry for CFD simulation. The final design dimensions for the volute are shown in Table 5. The reference point for the azimuth angle, θ, was illustrated in Figure 5.

4. Detailed Design of Rotor

The general gaspath design system is employed in the rotor blades detailed aerodynamic design (Figure 6). The Bezier polynomial curve [26] is used to evaluate the hub and shroud curve imported by the preliminary design system. This curve is well suited to the gaspath design process to generate a smooth curve for the blade surfaces. The suction and pressure sides of the rotor blade at the hub and shroud can be computed by Equations (26)–(28). As is shown in the Figure 6, the cross product of the tangent vector to blade camberline, S , and the vector, B , directed along the quasi-normal line determines the vector T , which defines the local normal direction to camberline surface. By determining the hub and shroud contours at the suction and pressure sides, the intermediate surface coordinates are defined by connecting the hub and shroud surface points along the quasi-normal line.
[ T x T y T z ] = S × B = det ( [ i j k sin θ c sin φ sin β c + cos θ c cos β c cos θ c sin φ sin β c sin θ c cos β c cos φ sin β c ( x c s x x h ) / L ( y c s y c h ) / L ( z c s z c h ) / L ] )
{ x = r c sin θ c ± 1 2 t r T x y = r c cos θ c ± 1 2 t r T y z = z c ± 1 2 t r T z
cot β c = r c d θ c d m
θc and βc are blade and polar angle with respect to the tangent at the hub and shroud camberlines, respectively. The angle φ (Figure 6) is the meridional surface slope angle with respect to z. It should be suitable to maintain a constant rotor tip radius for all meridional surfaces, according to Equation (29) and polar to Cartesian coordinate transformation.
{ θ 4 = θ c , 4 ± tan 1 [ t r , 4 2 r c , 4 ( T x cos θ c , 4 T y sin θ c ) ] x 4 = r 4 sin θ 4 y 4 = r 4 cos θ 4
The final geometry of the designed radial-inflow turbine for the air cycle machine is shown in Figure 7, Figure 8 and Figure 9, which show the rotor blade angle and the polar angle in the hub and shroud section.

5. Numerical Analysis

The CFD solver ANSYS-CFX has been used to solve the three-dimensional compressible unsteady Reynolds-averaged Navier-Stokes (URANS) equation. The high resolution discretization was employed for momentum and turbulence equations. The standard k-ε has been chosen to model the turbulent flow. The pressure inlet and pressure outlet boundary conditions were set for the volute inlet and the rotor outlet. A frozen rotor interface was set between the volute and the nozzle and a mixing plane condition was set at the interface between the stationary (nozzle) and rotational (rotor) frames. Therefore, the circumferential variation of the flow has no important effects on the stage efficiency and mass flow rate. The stage interface employs time-averaging to estimate the parameters in the interface between the rotor and nozzle stages. The physical time step is considered equal to 1 / ω , where ω is the revolution speed in rad/sec. Due to the difference in the number of blades in the rotor and nozzle stages, a pitch angle change, which is equal to the rotor-nozzle blade ratio, is considered for the flow.

5.1. Grid Generation

A structured O-H three-dimensional mesh was generated by ANSYS Turbogrid software for the nozzle and rotor computational domains (Figure 10). The grid independency shown in Figure 11 was performed to choose the appropriate number of cells for the meshes of rotor and nozzle blades. The number of cells in the O-type boundary layer and spanwise direction for both nozzle and rotor has been changed from 8 to 12 and from 30 to 60, respectively. According to Figure 11, it is obvious that the suitable number of cells for the nozzle and rotor blade domain should be more than 300,000. It was shown that the dimensionless wall distance,   y + , in the range of 30 to 200 had no obvious effects on prediction of the RI turbine aerodynamic performances so that the mesh was generated with this range of y + in this work.

5.2. CFD Validation

The validation of the current 3D viscous simulation was made against the experimental analysis of the RI turbine conducted by Hiett and Johnston [27]. Figure 12 shows a comparison of the stage performance obtained by numerical analysis and experiments. A good agreement is observed between numerical and experimental results; the deviation of stage efficiency at design point is negligible and the maximum efficiency deviation at off-design conditions is below 10%.

5.3. CFD Results

3D numerical simulations have been performed for design and off-design points. The mass flow rate map and the efficiency (total-to-total) map are illustrated in Figure 13 and Figure 14. The maximum total-to-total efficiency is about 87.8% at a total-to-total velocity ratio of 0.73. The total-total efficiency at design point is about 1.3% less than the maximum efficiency. However, the total-to-total efficiency is more than 80% for a wide range of velocity ratio ( 0.5 < v t t < 0.9 ) even at 80% of nominal speed. The corresponding points in Figure 13 and Figure 14 show that the new design of the turbine will have a good aerodynamic performance at a low mass flow rate and off-design rotational speed. It shows that we could achieve a good design for a small radial-inflow turbine employed in an air cycle machine. In Figure 15, the total-to-total efficiency of the designed turbine was compared with that of current AUT turbine at the design speed. As illustrated in the figure, a 15% shift in efficiency is achieved at the design conditions.
Figure 16 and Figure 17 present the evolution of total-to-static temperature ratio with pressure ratio and mass flow rate at various rotational speeds. The temperature ratio is important for the turbine in an air conditioning system due to the possibility of dehumidification and icing related to temperature drops. At the design point, the static temperature drops to about −25 °C, while the static and total pressure is slightly more than 1 atm, making dehumidification possible. Since the low exit temperature for an RI turbine employed in air conditioning systems is unavoidable, the systems use an auxiliary passage to transfer hot air to the turbine exit. Although the possibility of the dehumidification is reduced by a reduction of mass flow rate and rotational speed, more studies and simulations with wet air are required.
The meridional view of static pressure, static temperature and Mach number through stator and rotor flow passage have been presented in Figure 18. The exit station of stator has been slightly choked at design point with a maximum Mach number 1.05. The Mach and pressure contours presented in Figure 18 show that the flow expands gradually without any separation at hub. The static temperature at the exit station of the rotor change was from 273 K (0 °C) at the shroud to 235 K (−38 °C) at the hub. Since the mass flow rate at the design point is almost equal to the maximum mass flow rate, the range of rotor exit temperature will be more positive by the reduction in cooling demand.
Figure 19 illustrates the streamline and Mach number contour at mid span (50% blade height) for design point conditions. According to the figure, a small local recirculation is noticeable along the first 10% of the suction surface of the rotor blade. The recirculation bubble size increases by reduction of rotational speed and mass flow rate (at off-design conditions). In Figure 20, the streamline and Mach contour at mid-span for the 65% of nominal rotational speed and PRts = 2.0 is presented. Even at that off-design conditions, the local recirculation is only at first 20% of the suction surface. The effects of mass flow rate on the rotor leading edge recirculation are illustrated in Figure 21. As is shown in Figure 21, the recirculation is shifted from the suction surface to the pressure surface by a reduction of mass flow rate to 0.5 design value. The size of the recirculation is small enough (about 25% blade length) to have a noticeable negative effect on turbine efficiency. It shows that the new design of the RI turbine has an acceptable aerodynamic performance for a broad range of working conditions. Although circulation can be slightly improved by applying geometric modifications to the rotor blade leading edge (such as by using the elliptic leading edge), manufacturing costs will increase.

6. Conclusions

This paper presented a full design process for the radial-inflow turbine (RI turbine) used in an air cycle machine (ACM). The design process included a preliminary and detailed aerodynamic design for the volute, nozzle blade row, and rotor blade row. The preliminary design was coupled with a mean-line (1D) analysis to evaluate aerodynamic performance. An iteration process between the preliminary design and mean-line analysis was employed to aerodynamically achieve a high-quality turbine to be imported in the detailed design process. In detailed design, the suction and pressure surfaces of the rotor and the section surfaces of the volute were determined.
Unsteady 3D simulation using ANSYS CFX was performed to obtain the aerodynamic performance of the radial-inflow turbine designed by the presented process. The total-to-total efficiency was computed to be more than 80% for a wide range of velocity ratio ( 0.5 <   v t t < 0.9 ); even at 80% of the nominal speed, the maximum value was about 87.7%. The averaged rotor exit temperature was calculated at about −25 °C at design conditions; it is possible that dehumidification and even icing were observed at the rotor exit. Simulations showed that there is circulation at the rotor blade leading edge, located at the first 10% of the rotor suction surface for the design conditions and enlarged by the reduction of mass flow rate and rotational speed. However, at the critical mass flow rate, m ˙ = 0.5 m ˙ d , the size of the recirculation is small enough (about 25% blade length) to have a noticeable negative effect on turbine efficiency. These properties, along with the efficiency map of the turbine, shows that the design process can create a high-aerodynamic performance radial-inflow turbine, which operates at a high rotational speed in the air cycle machine.

Author Contributions

A.K.N. designed, conducted the simulation, and wrote the paper. S.S. provided guidance, suggestions for the proposed framework, and wrote the paper.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflicts of interest.

Nomenclature

A Area ( m 2 )
a Sound speed, the location of maximum camber along the chord
b Blade width passage (mm), maximum camber
C Absolute Velocity (m/sec)
C 0 Spouting Velocity (m/sec)
C p ,   C v Specific heats at constant pressure and volume (J/(kg·K))
c Nozzle chord
d Location of maximum thickness along the chord
D Diameter (mm)
H Total enthalpy (J/kg)
m ˙ Mass flow rate (kg/sec)
N Number of rotor/nozzle and splitter blades
n s Specific Speed
Q Volume flow rate
o Throat width (mm)
P Pressure (Pa)
P R Pressure Ratio
r Radius (mm)
R e Reynolds number
s Blade pitch (mm)
T Temperature (K)
T R Temperature ratio
t Blade thickness (mm)
U Blade local velocity (m/sec)
W Relative Velocity (m/sec)
x X-coordinate of nozzle airfoil
Y Loss coefficient
y Y-coordinate of nozzle airfoil
Z Axial dimension
Greek Symbols
α Absolute flow angle (deg.)
β Blade angle (deg.)
γ Specific heat ratio
Δ Total dimensionless boundary layer displacement thickness
η Efficiency
Θ Total dimensionless boundary layer momentum thickness
θ Boundary layer momentum thickness (mm), polar angle (deg.)
v Velocity ratio
ρ Density (kg/sec)
ϕ Mean surface angle with axial direction (deg.)
ω Angular Velocity (rad/sec)
Subscripts
1 , 2 , 3 , 4 , 5 State point in the turbine
a v e Average of a value
B L Blade loading
C L Clearance
H S Hub-to-shroud
i n s Incidence
i d Isentropic
h Hub
m Meridional component
m a x Maximum value
p Profile loss
R Rotor
r e f Reference value
s Shroud
t Total condition
t s Total-to-static
t t Total-to-total
θ Tangential component/circumferentially
* Parameter corresponding to sonic flow
Superscript
* Optimum

References

  1. Conceição, S.T.; Zaparoli, E.L.; Turcio, W.H. Thermodynamic Study of Aircraft Air Conditioning Air Cycle Machine: 3-wheel × 4-wheel. SAE Tech. Paper 2007. [Google Scholar]
  2. Santos, A.P.; Andrade, C.R.; Zaparoli, E.L. A Thermodynamic Study of Air Cycle Machine for Aeronautical Applications. Int. J. Thermodyn. 2014, 17, 117–225. [Google Scholar] [CrossRef]
  3. Simpson, A.T.; Spence, S.W.; Watterson, J.K. A Comparison of the Flow Structures and Losses Within Vaned and Vaneless Stators for Radial Turbines. J. Turbomach. 2009, 131, 031010. [Google Scholar] [CrossRef]
  4. Rohlik, H. Analytical determination of radial inflow turbine design geometry for maximum efficiency. In National Aeronautics and Space Administration; NASA Lewis Research Center: Cleveland, OH, USA, 1969; p. 4384. [Google Scholar]
  5. Feng, Z.; Deng, Q.; Li, J. Aerothermodynamic Design and Numerical Simulation of Radial Inflow Turbine Impeller for a 100kW Microturbine. Am. Soc. Mech. Eng. 2005, 1, 873–880. [Google Scholar]
  6. Zheng, Y.; Hu, D.; Cao, Y.; Dai, Y. Preliminary design and off-design performance analysis of an Organic Rankine Cycle radial-inflow turbine based on mathematic method and CFD method. Appl. Therm. Eng. 2017, 112, 25–37. [Google Scholar] [CrossRef]
  7. Aungier, R.H. Turbine Aerodynamics: Axial-Flow and Radial-Flow Turbine Design and Analysis; ASME: New York, NY, USA, 2006; pp. 10016–15990. [Google Scholar]
  8. Rohlik, H.E. Turbine Design and Application: Radial-Inflow Turbine; NASA: Washington, DC, USA, 1975; p. 392. [Google Scholar]
  9. Glassman, A.J. Computer Program for Design Analysis of Radial-Inflow Turbines; NASA-TN-D-8164; NASA: Washington, DC, USA, 1976; p. 66.
  10. Rodgers, C. Advanced Radial Inflow Turbine Rotor Program: Design and Dynamic Testing; NASA-CR-135080, SOLAR-ER-2519; NASA: Washington, DC, USA, 1976. [Google Scholar]
  11. Whitfield, A. The Preliminary Design of Radial Inflow Turbines. J. Turbomach. 1990, 112, 50. [Google Scholar] [CrossRef]
  12. Ventura, C.A.; Jacobs, P.A.; Rowlands, A.S.; Petrie-Repar, P.; Sauret, E. Preliminary Design and Performance Estimation of Radial Inflow Turbines: An Automated Approach. J. Fluids Eng. 2012, 134, 31102. [Google Scholar] [CrossRef]
  13. Ebaid, M.S.Y.; Bhinder, F.S.; Khdairi, G.H. A Unified Approach for Designing a Radial Flow Gas Turbine. J. Turbomach. 2003, 125, 598. [Google Scholar] [CrossRef]
  14. Al Jubori, A.; Daabo, A.; Al-Dadah, R.K.; Mahmoud, S.; Ennil, A.B. Development of micro-scale axial and radial turbines for low-temperature heat source driven organic Rankine cycle. Energy Convers. Manag. 2016, 130, 141–155. [Google Scholar] [CrossRef]
  15. Sauret, E.; Gu, Y. Three-dimensional off-design numerical analysis of an organic Rankine cycle radial-inflow turbine. Appl. Energy 2014, 135, 202–211. [Google Scholar] [CrossRef]
  16. Fiaschi, D.; Innocenti, G.; Manfrida, G.; Maraschiello, F. Design of micro radial turboexpanders for ORC power cycles: From 0D to 3D. Appl. Therm. Eng. 2016, 99, 402–410. [Google Scholar] [CrossRef] [Green Version]
  17. Li, Y.; Ren, X.D. Investigation of the organic Rankine cycle (ORC) system and the radial-inflow turbine design. Appl. Therm. Eng. 2016, 96, 547–554. [Google Scholar] [CrossRef]
  18. Costall, A.W.; Hernandez, A.G.; Newton, P.J.; Martinez-Botas, R.F. Design methodology for radial turbo expanders in mobile organic Rankine cycle applications. Appl. Energy 2014, 157, 729–743. [Google Scholar] [CrossRef]
  19. Rahbar, K.; Mahmoud, S.; Al-Dadah, R.K.; Moazami, N.; Mirhadizadeh, S.A. Development and experimental study of a small-scale compressed air radial inflow turbine for distributed power generation. Appl. Therm. Eng. 2017, 116, 549–583. [Google Scholar] [CrossRef]
  20. Zhou, A.; Song, J.; Li, X.; Ren, X.; Gu, C. Aerodynamic design and numerical analysis of a radial inflow turbine for the supercritical carbon dioxide Brayton cycle. Appl. Therm. Eng. 2018, 132, 245–255. [Google Scholar] [CrossRef]
  21. Aungier, R.H. Centrifugal Compressors: A Strategy for Aerodynamic Design and Analysis; ASME Press: New York, NY, USA, 2000. [Google Scholar]
  22. Jones, A.C. Design and Test of a Small, High Pressure Ratio Radial Turbine. J. Turbomach. 1996, 118, 362. [Google Scholar] [CrossRef]
  23. Aungier, R.H. Axial-Flow Compressors: A Strategy for Aerodynamic Design and Analysis; ASME Press: New York, NY, USA, 2003; pp. 10016–15990. [Google Scholar]
  24. Watanabe, I.; Ariga, I.; Mashimo, T. Effect of Dimensional Parameters of Impellers on Performance Characteristics of a Radial-Inflow Turbine. J. Eng. Power 1971, 93, 81. [Google Scholar] [CrossRef]
  25. Simpson, A.; Spence, S.; Watterson, J. Numerical and Experimental Study of the Performance Effects of Varying Vaneless Space and Vane Solidity in Radial Inflow Turbine Stators. In Proceedings of the ASME Turbo Expo 2008: Power for Land, Sea and Air GT2008, Berlin, Germany, 9–13 June 2008; Volume 135, pp. 1495–1505. [Google Scholar]
  26. Forrest, A.R. Interactive interpolation and approximation by Bézier polynomials. Comput. Des. 1990, 22, 527–537. [Google Scholar] [CrossRef] [Green Version]
  27. Hiett, G.F.; Johnston, I.H. Experiments Concerning the Aerodynamic Performance of Inward Flow Radial Turbines. Proc. Inst. Mech. Eng. Conf. Proc. 2006, 178, 28–42. [Google Scholar] [CrossRef]
Figure 1. Design process flowchart.
Figure 1. Design process flowchart.
Applsci 08 02207 g001
Figure 2. Design station for calculating geometric and aerodynamic parameters.
Figure 2. Design station for calculating geometric and aerodynamic parameters.
Applsci 08 02207 g002
Figure 3. Velocity diagram in a rotor blade.
Figure 3. Velocity diagram in a rotor blade.
Applsci 08 02207 g003
Figure 4. Rotor hub contour.
Figure 4. Rotor hub contour.
Applsci 08 02207 g004
Figure 5. An external circular volute section.
Figure 5. An external circular volute section.
Applsci 08 02207 g005
Figure 6. Gas path design system used for the rotor blade detailed design.
Figure 6. Gas path design system used for the rotor blade detailed design.
Applsci 08 02207 g006
Figure 7. The final geometry of the designed RI turbine.
Figure 7. The final geometry of the designed RI turbine.
Applsci 08 02207 g007
Figure 8. Rotor blade angle in the hub and shroud section.
Figure 8. Rotor blade angle in the hub and shroud section.
Applsci 08 02207 g008
Figure 9. Rotor polar angle in the hub and shroud section.
Figure 9. Rotor polar angle in the hub and shroud section.
Applsci 08 02207 g009
Figure 10. 3D view of the mesh generated for the newly designed turbine.
Figure 10. 3D view of the mesh generated for the newly designed turbine.
Applsci 08 02207 g010
Figure 11. Grid independency for the rotor and nozzle blade rows.
Figure 11. Grid independency for the rotor and nozzle blade rows.
Applsci 08 02207 g011
Figure 12. Comparison between experiments and CFD for the Ricardo RI turbine [27].
Figure 12. Comparison between experiments and CFD for the Ricardo RI turbine [27].
Applsci 08 02207 g012
Figure 13. Total-to-total efficiency map for the new design of RI turbine.
Figure 13. Total-to-total efficiency map for the new design of RI turbine.
Applsci 08 02207 g013
Figure 14. Mass flow rate map for the new design of RI turbine.
Figure 14. Mass flow rate map for the new design of RI turbine.
Applsci 08 02207 g014
Figure 15. Comparison between the efficiency of the designed turbine and the current version.
Figure 15. Comparison between the efficiency of the designed turbine and the current version.
Applsci 08 02207 g015
Figure 16. Evolution of total-to-static temperature with total-to-static pressure for the new design of RI turbine.
Figure 16. Evolution of total-to-static temperature with total-to-static pressure for the new design of RI turbine.
Applsci 08 02207 g016
Figure 17. Evolution of total-to-static temperature with mass flow rate for the new design of RI turbine.
Figure 17. Evolution of total-to-static temperature with mass flow rate for the new design of RI turbine.
Applsci 08 02207 g017
Figure 18. Meridional view of static pressure, static temperature and Mach number at the design point.
Figure 18. Meridional view of static pressure, static temperature and Mach number at the design point.
Applsci 08 02207 g018
Figure 19. The streamline and Mach number contour at mid span for design point conditions.
Figure 19. The streamline and Mach number contour at mid span for design point conditions.
Applsci 08 02207 g019
Figure 20. The streamline and Mach number contour at mid span for 65% of nominal rotational speed and PRts = 2.0.
Figure 20. The streamline and Mach number contour at mid span for 65% of nominal rotational speed and PRts = 2.0.
Applsci 08 02207 g020
Figure 21. The streamline for different mass flow rate at mid span.
Figure 21. The streamline for different mass flow rate at mid span.
Applsci 08 02207 g021
Table 1. Input parameters to design the RI turbine for the ACM.
Table 1. Input parameters to design the RI turbine for the ACM.
ParameterUnitValue
Total inlet temperatureK364
Total inlet pressurekPa541.822
mass flow rateKg/sec0.217
total-to-static pressure ratio-4.85
number of rotor revolution-75,000
Table 2. Constant parameters for the thickness distribution of the nozzle blade profile [7].
Table 2. Constant parameters for the thickness distribution of the nozzle blade profile [7].
ParameterValue
  t L E / c 0.025
  t T E / c 0.012
  t m a x / c 0.06
  d / c 0.4
Table 3. Loss models used in the mean-line analysis.
Table 3. Loss models used in the mean-line analysis.
Types of LossEmpirical Correlation
Volute
Circumferential distortion loss   Y θ =   [ r   1 C   1 r   3 C   θ 3 C   3 ] 2
Nozzle
Incidence loss coefficient   Y i n c = s i n 2 ( α 1 α * ) P t 1 P 1 P t 3 P 3
Rotor
Blade loading loss   Y B L = 1 24 [ Δ W W 3 ] 2
Hub-to-shroud loading loss   Y H S = 1 24 [ k m b 3 W 2 W 3 sin α 3 ] 2
Incidence loss   Y i n c = s i n 2 ( α 1 α 1 * ) P t 1 ´ P 1 P t 3 ´ P 3
Clearance loss Y C L = m ˙ C L Δ P / [ m ˙ ( P t 3 P 3 ) ]
Table 4. Output dimensions from preliminary design.
Table 4. Output dimensions from preliminary design.
ParameterUnitValue
Volute
Volute inlet area, A 1 mm2414.8
Volute   inlet   radius ,   r 1 mm69.0
Volute discharge radius, r 2 mm57.9
Nozzle
Nozzle inlet radius, r 2 mm55.1
Nozzle discharge radius, r 3 mm48.8
Nozzle exit throat width, o 3 mm3.6
Nozzle passage width, b 2 = b 3 mm2.9
Nozzle inlet blade angle, β 2 deg.28.6
Nozzle exit blade angle, β 3 deg.7.0
Number of nozzle blade, N N -20
Rotor
Rotor inlet radius, r 4 mm40.6
Rotor exit hub radius, r h 5 mm7.5
Rotor exit hub radius, r s 5 mm21.8
Rotor throat width, o s 5 mm4.7
Rotor inlet blade angle, β 4 deg.90.0
Rotor exit blade angle, β 5 deg.45.9
Rotor blade thickness, t r mm1.0
Number of rotor blades, N R -13
Table 5. Design dimensions of the turbine volute.
Table 5. Design dimensions of the turbine volute.
No.Azimuth Angle [deg.]Radius of Center [mm]Radius of Section [mm]
1057.901.00
23060.933.03
36062.234.33
49063.245.34
512064.116.21
615064.886.98
718065.587.68
821066.238.34
924066.858.95
1027067.439.53
1130067.9810.09
1233068.5110.62
1336069.0211.13

Share and Cite

MDPI and ACS Style

Karimi Noughabi, A.; Sammak, S. Detailed Design and Aerodynamic Performance Analysis of a Radial-Inflow Turbine. Appl. Sci. 2018, 8, 2207. https://0-doi-org.brum.beds.ac.uk/10.3390/app8112207

AMA Style

Karimi Noughabi A, Sammak S. Detailed Design and Aerodynamic Performance Analysis of a Radial-Inflow Turbine. Applied Sciences. 2018; 8(11):2207. https://0-doi-org.brum.beds.ac.uk/10.3390/app8112207

Chicago/Turabian Style

Karimi Noughabi, Amir, and Shervin Sammak. 2018. "Detailed Design and Aerodynamic Performance Analysis of a Radial-Inflow Turbine" Applied Sciences 8, no. 11: 2207. https://0-doi-org.brum.beds.ac.uk/10.3390/app8112207

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