Geometrical and functional specifications of analysed turbo-machinery.
The aero-thermodynamic design and performance of a compressor need to conquer many vital challenges like it is a gas-driven turbo-machinery component, involvement of extensive iterative process for the convergence of the design, enormous design complexity due to three-dimensional flow phenomena, and multiflow physics embedded within a dynamic state-of-the-art. In this chapter, a strong attempt is made to address the above-cited technical issues to achieve an optimized design and performance of a centrifugal compressor with backward swept blade profile producing total pressure ratio of 5.4 with an ingested mass flow rate of 5.73 kg/s. A mean-line design methodology was implemented to configure sizing of the compressor. An optimum grid size was well validated by carrying out computational analysis with three different mesh sizes within the same framework. Finally, a detailed three-dimensional numerical simulation was performed using Reynolds-averaged Navier-Stokes equations based on finite volume discretization method (RANS-FVDM) scheme. Consequently, the polytropic efficiency, total-to-total efficiency, stagnation pressure ratio at a fixed rotational speed, and the overall design and aero-thermodynamic performance of the centrifugal compressor are validated.
- numerical simulations
- Reynolds-averaged Navier-Stokes equations
- finite volume discretization method
- compressor design
Centrifugal compressors find usage over wide range of propulsion applications and are regarded as one among the key air-breathing propulsive engine components. The cognitive research and development of compressors is directed toward achieving a higher pressure ratio, higher efficiency, and reduced structural weight of compressor and the engine as well. Various compressor stages achieve gradual increase in the stagnation-to-flow pressure contributed by flow diffusion. Energy is added in the rotor blade section, increasing the total pressure and absolute component of flow velocity. Stator blade row diffuses the flow, thus reducing absolute velocity component and elevating static pressure. Blade topology requires adaptation of a cautious design procedure to achieve the designated pressure rise while minimizing aero-thermodynamic losses in order to run and achieve design pressure ratios and design efficiencies.
In this chapter, a strong attempt is made to enumerate the detailed procedure of the centrifugal compressor in Section 2, concepts and basics of Numerical Schemes in Section 3. Further, a case study is discussed in detail in Section 4, the corresponding results and discussions are presented in Section 5. Final conclusions are summarized in Section 6.
2. Centrifugal compressor and the basic steps involved in its design
The design of the centrifugal compressor is very difficult due to the reasons like it is a gas-driven turbo-machinery component, involvement of extensive iterative process for the convergence of the design, enormous design complexity due to three-dimensional flow phenomena, and multiflow physics embedded within a dynamic state-of-the-art. Hence, in this article, a detailed procedure for the design is presented based on the design methodology as explained in the literature [1, 2, 3]. The basic steps involved in the design of the centrifugal compressor are discussed in detail in the following sections.
2.1. Impeller and velocity triangles
Impeller is the main design task during the phase of the compressor design. There are various strategic geometric/design features to be identified and discussed with respect to impeller, its inlet and outlet design. In principal, aerodynamic losses occurring in majority of turbo-machines arise primarily due to the boundary layer growth, its separation on blade profile, and passage surfaces termed as profile losses widely configured under primary losses.
The velocity triangles of the impeller of a centrifugal compressor are diagrammatically shown in Figure 1.
The impeller design is the basic and major vital part of the design of the compressor. The following conditions for the design may be considered as limiting conditions
At the inlet of the impeller, very high relative velocity of the shroud should be avoided.
Separation of the flow in the passages of the impeller should be minimized so that the losses are minimized.
To meet the above limiting condition, the de Haller impeller ratio should be .
The stability of the diffuser strongly depends on and for both vaneless and vaned diffusers, . Further, it should be taken from to for better design.
For the backswept impellers, the ratio of the relative velocities should be increased for betterment in the design.
Little bit of the increase in the stagnation enthalpy will be attributed to that of the velocity, because of which the polytropic efficiency of impellers with nil back sweep will be less than that of backswept impellers.
Work coefficient for radial impellers with zero back sweep is around 0.9 and 0.5 for backswept impellers.
To with stand the stress levels, the maximum rotational speed of the backswept impeller is maintained less than that of the radial impeller with same tip diameter.
The above two conditions lead to decrease in the maximum possible work per stage for a backswept impellers.
Hence, if the maximum work per stage is not required, blade swept angles up to can be used.
2.2. Specific work
The specific work can be calculated using the following correlation, which is based on energy and Euler turbine equations
2.3. Slip factor
In general for the design of the compressor, the inlet swirl is considered as zero and hence, . Furthermore, the Wiesner’s correlation for slip factor is given by
where number of blades is denoted by Z and represents the angle between the radial direction and tangent to the rotor blade at the periphery.
should be less than such that the downstream diffuser does not get prone to stall when and are reduced by keeping very large .
2.4. Number of blades
The number of blades and the exit blade angle are dependent on each other and can be calculated using loading coefficient. Loading coefficient is the ratio between the outlet tangential flow velocity and blade speed, which is given by the correlation
For the sake of easy reference, even number of blades are chosen so that half of the blades can be considered as splitter blades. In general for a better design, 20 blades and blade angle at outlet are taken.
2.5. Blade peripheral speed
The peripheral speed of the impeller is calculated using the following correlations
where increase in the enthalpy is given by
The isentropic law can be applied, which is given by
The polytropic efficiency can be calculated by choosing isentropic efficiency as 0.83 and optimum specific speed as 0.8 for the better design and using the following correlation
Other blade parameters can be calculated using the regular trigonometric correlations from the velocity triangles of the blade like the following
2.6. Density at the inlet of the blade
The Mach number at the inlet can be calculated using the following correlation
Static density can be calculated from the following correlations
2.7. Rotational speed
The rotational speed can be calculated using the following correlation
2.8. Blade axial width at outlet
The outlet blade axial width is calculated with the assumptions that the effect of the thickness and boundary layer of the blade are neglected. The following correlations are used for the calculations
where,, the rotor diameter can be calculated using the correlation
2.9. Flow separation condition
Flow separation should be minimized in the design. To reduce the separation of the flow, the back sweep angle should be higher. Furthermore, the separation in the rotor with subsonic inlet relative Mach number and normal Reynold’s number can be minimized by using a lower limit of the ratio of the outlet to inlet velocity ratio as 0.8.
3. Numerical schemes
In any gas turbine, the flow of the fluid in the compressor is always more unsteady and turbulent, which further makes the design more complicated. For all turbulent flows, the governing equations are the unsteady Navier-Stokes (N-S) equations. But those equations are very much difficult to solve. The following are the governing equations in tensor notation for understanding the basic nature of the equations as presented by .
Energy equation is given by
The main disadvantage in solving N-S equations is computing. This is due to the reason that to represent even very small scale of velocity and pressure fluctuations using N-S equations, temporal resolution is required as the N-S equations are spatial fineness equations. Furthermore, the accuracy and resolution of the scheme reduce due to the increase in accumulative rounding off value errors because of the increase in grid points to achieve fine meshes.
Hence, more accurate schemes are required to achieve the solutions of turbulence models with highest resolution. One such numerical scheme, which can effectively solve the turbulence problems, is Reynolds-averaged Navier-Stokes (RANS) equations.
3.1. Reynolds-averaged Navier-Stokes (RANS) equations
The transport equations should be modified by introducing the components with averaged and fluctuating components for solving the turbulent models. RANS are the equations of the motion of the fluid flow with time-averaged equations. If flow turbulence properties are known, then suitable approximations can be made and high-resolution solutions to the N-S equations can be achieved by solving the RANS equations. The RANS equations in tensor notation are described below 
where is the mean rate of strain tensor.
The problems of physical engineering, which involve high turbulence when modeled based on RANS equations, are also called statistical turbulence models because the method strongly involves statistical averaging procedure. Consequently, the computational effort to solve RANS equations is very less when compared to the other schemes.
3.2. Governing equations for turbulence models
As discussed in the earlier articles, due to the involvement of more unknown parameters in the modeling of turbulent problems, the accuracy of the model becomes very much difficult and leads to erroneous results. Hence, it is required to develop a numerical procedure to close or converge the system of equations. One such two equations category turbulence model is K-epsilon model .
3.2.1. K-epsilon models
In this model, the velocity of the turbulent flow and the length scales are independently calculated using two different equations. The basic model is the standard k-epsilon model .
220.127.116.11. Standard k-model
The velocity of the turbulence can be calculated by forming a model for the corresponding kinetic energy using the following Eq. (4)
The length scale is represented by ε, which is the rate of dissipation and can be calculated from the Eq. (4)
where represents the turbulence kinetic energy due to the mean velocity gradients, represents the turbulence kinetic energy due to buoyancy, and represents the contribution of the fluctuating dilatation in compressible turbulence to the overall dissipation rate.
3.3. Mean: Line analysis
The difficulty in the design of the compressor is because it involves two vital phases of the design. One being the 1D design of the compressor and the other is the deep numerical analysis of the design. The difficulty in the first phase is overcome by using mean line theory. In principle, the mean line theory follows the preliminary design that is carried out by neglecting the air flow variations in radial direction and the location of the mean blade radius is considered for the analysis.
3.4. Finite volume discretization method(FVDM)
Any fluid flow problem can be solved with high resolution when the basic concept of continuum is considered during the numerical solution methodology. Finite volume discretization method (FVDM/FVM) is one such methodology developed by McDonald, MacCormack, and Paullay during 1970s. The basic structure of the FVM is as follows:
The physical problem is discretized into control volumes.
Each control volume is formulated with balance equations involving integrals.
Numerical integration is applied to approximate the integral equations.
Nodal value interpolations are used to approximate the derivatives and function values.
The solution is achieved by assembling the algebraic systems of the discrete control volumes.
4. Case study: Design of the centrifugal compressor
The case study presented in the literature  is considered for the present illustration for the design and numerical analysis of the centrifugal compressor. Main objective of the current work is to design a centrifugal compressor capable of delivering a pressure ratio of 5.4 with a mass flow ingestion rate of 5.73 kg/sec. The compressor is targeted to achieve 82% total-total efficiency with design constraints of 280 mm for the impeller and 340 mm for the diffuser in terms of diameter. Also, the compressor has to generate flow parameters relative to functional downstream components.
4.1. Centrifugal compressor specifications
The input parameters of the centrifugal compressor are tabulated in Table 1.
|S. No||Name of the parameter||Value|
|1||Mass flow rate||5.73 kg/s|
|4||Maximum impeller diameter||280 mm|
|5||Maximum diffuser diameter||340 mm|
|6||Rotational speed||38,000 rpm|
|7||Inlet pressure||101,325 Pa|
|8||Inlet temperature||288.15 K|
Preliminary design of the impeller and vaned diffuser is carried out by adapting a one-dimensional model approach. Comprehensive design of the centrifugal compressor stage was generated using ANSYS BladeGen module.
4.2. Impeller design
Impeller was designed using guidelines from various sources and in specific from . The number of impeller blade profiles was configured based upon a choice iterative method aiming for passage flow devoid of heavy separations, and it was streamlined and set at 19 blade profiles. A rotational speed of 38,000 rpm was set for the impeller.
4.2.1. Impeller inlet
Hub-to-tip ratio is a key owing to the secondary losses, which occurs in the flow regions near the end walls. The presence of any undesirable circulatory or cross-flows develops on account of rapid and steep flow turning through the blade channel accounting for annulus wall boundary layers. Therefore, impeller hub-to-tip ratio was set at 0.35 within the range of 0.3–0.6 prescribed in .
Mach number at impeller inlet was set at 0.65. Impeller inlet blade angles were setup by PCA Engineer’s Vista CCD tool analytical calculations. The inducer leading edge angles are 35, 56, and 63°, respectively, with incidences of 11.8, 3.7, and 0.1° at hub, mean, and shroud, respectively. Leading edge of the impeller was defined by an elliptic ratio of 6.
4.2.2. Impeller outlet
Impeller backswept angle was set at 0° to minimize impeller diameter, and a lean angle of 30° was also incorporated into the design. Impeller exducer height at 13.7032 mm and impeller diameter of 257.954 mm were set by PCA Engineer’s Vista CCD software impeller trailing edge that was defined as a square cutoff.
Data generated by Vista CCD tool were used to generate a 3D computer-aided design model of impeller. Inlet portion of 50 mm and horizontal was designed satisfying the diameter constraint of 280 mm.
4.3. Design of vaned diffuser
In order to obtain higher pressure ratio in a radial diffuser, the diffusion process has to be achieved across a relatively shorter radial distance. This requires the application of vanes, which provide greater guidance to flow inside diffusing passages. The vaned diffuser was designed by observing various flow parameters reflected at impeller exit after performing numerical simulations.
To circumvent flow separation, divergence of diffuser blade passages in vaned diffuser ring can be kept small by incorporating a large number of vanes. However, this can lead to higher friction losses. Thus, an optimum number of diffuser vanes must be employed and ensure flow passage divergence not to exceed 12°. Thus, final diffuser design contains 30 blades. The diffuser vane leading edge was at a radius of 136 mm, and the trailing edge of the diffuser vanes was set at a radius of 166 mm. Blade inlets and outlet angle were set at 64°. The leading and trailing edge were defined by an elliptic ratio of 6 and a radius of 0.25 mm.
Another method to prevent very steep velocity gradients at diffuser entry is by providing a small vane-less space (0.05d2–0.1d2) between impeller exit and diffuser entry. Therefore, a vane-less space of 7.023 mm was allotted.
The design methodology adapted was mainly focused aiming at decrease in Mach number and flow angle at diffuser exit by satisfying a diameter constraint of 340 mm.
4.4. Meshing strategy
Adapted meshing strategy for both the impeller and diffuser fluid domains was achieved using Ansys-Turbo grid module. A total node count of 7e-01 million was setup for the CFD solver, as it allows generation of refined quality hexahedral meshes required for the blade passages in turbo-machinery.
4.5. Numerical setup
In the present context, ANSYS tool is setup with pressure-based solver simulating a steady, three-dimensional viscous flow fields using complete set of Navier-Stokes code solving for Reynolds-averaged Navier-Stokes equations based on finite volume discretization method . A high-resolution scheme is used to solve for continuity, momentum, energy, and state equations implementing a standard k-ε turbulence model. Individual compressor stage characteristics were generated by performing simulations varying back pressures. Firstly, near choke condition, flow points are run to reduce static back pressure values and later, solutions are restarted with incrementally increasing the static back pressure to compute intermediate points on constant speed line running toward stall.
Inlet boundary conditions: At compressor inlet, a constant total pressure and total temperature conditions are imposed with a turbulence intensity of 1% and flow direction is marked normal to the inlet plane.
Outlet boundary conditions: At the outlet, an average static pressure boundary condition is implemented. Also, a circumferential symmetry condition is imposed on corresponding periodic surfaces with air-fluid medium setup as an ideal gas. A counter-rotating wall boundary is given at the impeller shroud.
5. Results and discussions
5.1. Flow analysis: Velocity triangles
Entry and exit velocity triangles for impeller blades along the radial section with backward swept blades are plotted as shown in Figure 2(a) and (b). Inlet Mach number is within the recommended range of 0.4 to 0.6. A maximum flow angle of 64° at the impeller exit and a Mach number of 1.4 are observed.
5.2. Blade to blade contours
The absolute Mach number plots show no supersonic regions of significant volume. The hub section plot show supersonic flow near suction surface of the blade and the tip Section span plots show a supersonic region on the pressure surface of the blade. The presence of a mild bow shock in front of the leading edge of the blade can also be observed. The Mach number is around 0.6 at the leading edge of the impeller. The following plots shown in Figure 4(a–c) are the relative Mach number contours for the hub, mean and tip.
The Mach number at the leading edge in the relative frame is near 1.2 for the leading edge of the impeller blade. The regions of low Mach number can also be observed, although it is in the relative frame.
5.3. Circumferential contour
5.4. Performance characteristics
Individual component as well as stage performance as a whole can be gauged by different methods.
5.4.1. Impeller performance characteristics
Impeller performance was evaluated by Vista-CCM and their corresponding plots are shown in below Figure 6 represents the variation of pressure ratio with varying mass flow ingested through compressor.
Figure 7 represents the variation of polytropic efficiency with varying mass flow ingested to quantify differential pressure changes occurring through compressor stage.
5.4.2. Stage performance characteristics
Compressor stage performance can be evaluated by obtaining individual characteristic curves as plotted in Figures 8–10. The compressor performance is presented only for a functional design speed of 38,000 rpm. Variables like total pressure ratio, adiabatic efficiency, and power requirement of compressor are plotted against varying mass flow rate as shown in Figures 8–10. The maximum achieved pressure ratio is 5.4 while allowing a fixed mass flow rate of 5.73 kg/s.
A peak efficiency of 75.8% is obtained at a pressure ratio of 5.73. The design point efficiency is predicted to be around 74.5%. Possible reasons for loss in efficiency may be owing to primary and secondary losses.
The impeller works with power a power requirement in the range of 1350–1360 kW on its 100% speed line over the range of total pressure ratio 4.5–5.4.
5.5. Validation studies
The design methodology adapted in current framework uses a backswept angle of 0° and does not follow the design trend established for validation in above figures. Above deviations in the performance trends can be attributed to the structural size and rotational specifications of various stages.
The aero-thermodynamic design of centrifugal compressor turbo-machinery for a small gas turbine engine has been successfully carried out and sequentially investigated for its performance at design speed. The 3D computational analysis infers that the designed centrifugal compressor generates a pressure ratio of 5.4 with an ingested mass flow rate of 5.73 kg/s at a fixed rotational speed of 38,000 rpm. The targeted flow parameters are in well agreement with centrifugal turbo-machinery within the diameter specification of 280 and 340 mm for impeller and diffuser, respectively. The total-to-total efficiency of designed centrifugal compressor is 81% against a target of 82%. The swirl at impeller is reduced from 77 to 27.5° by the diffuser, thus reducing the primary losses viz., flow separation and thereby has the capability of experiencing lesser stall angles.