Open access peer-reviewed chapter

Numerical Simulations of a High-Resolution RANS-FVDM Scheme for the Design of a Gas Turbine Centrifugal Compressor

By Chellapilla V K N S N Moorthy and Vadapalli Srinivas

Submitted: September 15th 2017Reviewed: October 31st 2017Published: February 6th 2018

DOI: 10.5772/intechopen.72098

Downloaded: 231

Abstract

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.

Keywords

  • numerical simulations
  • Reynolds-averaged Navier-Stokes equations
  • finite volume discretization method
  • compressor design

1. Introduction

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.

Figure 1.

Velocity triangles of impeller for a centrifugal compressor.

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

  1. At the inlet of the impeller, very high relative velocity of the shroud Wsh,1should be avoided.

  2. Separation of the flow in the passages of the impeller should be minimized so that the losses are minimized.

  3. To meet the above limiting condition, the de Haller impeller ratio should be W2WSh,1>0.75.

  4. The stability of the diffuser strongly depends on αc,2and for both vaneless and vaned diffusers, αc,270°. Further, it should be taken from 60to 65°for better design.

  5. For the backswept impellers, the ratio of the relative velocities should be increased for betterment in the design.

  6. 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.

  7. Work coefficient for radial impellers with zero back sweep is around 0.9 and 0.5 for backswept impellers.

  8. 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.

  9. The above two conditions lead to decrease in the maximum possible work per stage for a backswept impellers.

  10. Hence, if the maximum work per stage is not required, blade swept angles up to 50°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

u2Cu,2u1Cu,1=gcCpgT0,1p0,2p0,1RCp1ηp,c1E1

2.3. Slip factor

In general for the design of the compressor, the inlet swirl is considered as zero and hence, Cu,1=0. Furthermore, the Wiesner’s correlation for slip factor is given by

σw=slipfactor=Cu,2,acCu,2,tl=1cosβ2Z0.7E2

where number of blades is denoted by Z and β2represents the angle between the radial direction and tangent to the rotor blade at the periphery.

αc,2should be less than 60°such that the downstream diffuser does not get prone to stall when C2and Cr,2are reduced by keeping very large αc,2.

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

ψ=Cu,2,scu2=tanβ2tanαc,2+1σw1E3

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 45°blade angle at outlet are taken.

2.5. Blade peripheral speed

The peripheral speed of the impeller is calculated using the following correlations

h0=ψu2E4

where increase in the enthalpy is given by

h0=CpT0,2T0,1E5

The isentropic law can be applied, which is given by

T0,2T0,1=p0,2p0,1RCp1ηp,cE6

The polytropic efficiency ηp,ccan be calculated by choosing isentropic efficiency ηs,cas 0.83 and optimum specific speed Nsas 0.8 for the better design and using the following correlation

ηs,c=rRCp1rRCpηp,c1E7

Other blade parameters can be calculated using the regular trigonometric correlations from the velocity triangles of the blade like the following

Cu,2,ac=ψu2E8
Cr,2=Cu,2,actanαc,2E9
C2=Cu,2,acsinαc,2E10

2.6. Density at the inlet of the blade

The Mach number at the inlet can be calculated using the following correlation

CgcRT0=2CpR[11+M22CpR11]E11

Static density can be calculated from the following correlations

ρ0ρst=1+M22CpR1CpR1E12

2.7. Rotational speed

The rotational speed can be calculated using the following correlation

N=60Ns2πgch034V̇E13

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

ṁ=Cr,2ρst,2πd2b2E14

where,d2, the rotor diameter can be calculated using the correlation

u2=πd2N60E15

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 [4].

ρuit+ρuiujxj=pxi+xjμuixj+ujxixi23μukxk+SiE16

Energy equation is given by

ρϵt+ρujexj=pujxj+xjkTxj++SiE17

where

=μuixjuixj+ujxi23μukxk2E18

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 [4]

ρui¯t+ρui¯uj¯xj=ρfi¯+xjp¯δij+2μSij¯ρuiuj¯E19

where Sij¯=ui¯xj+uj¯xiis 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 [4].

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 [4].

3.2.1.1. 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)

ρkt+ρkuixi=xj[μ+μtσkkxj+Pk+PbρεYM+SkE20

The length scale is represented by ε, which is the rate of dissipation and can be calculated from the Eq. (4)

ρϵt+ρϵuixi=xj[μ+μtσϵϵxj+C1ϵϵkPk+C3ϵPbC2ϵρε2k+SϵE21

where Pkrepresents the turbulence kinetic energy due to the mean velocity gradients, Pbrepresents the turbulence kinetic energy due to buoyancy, and YMrepresents 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:

  1. The physical problem is discretized into control volumes.

  2. Each control volume is formulated with balance equations involving integrals.

  3. Numerical integration is applied to approximate the integral equations.

  4. Nodal value interpolations are used to approximate the derivatives and function values.

  5. 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 [5] 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. NoName of the parameterValue
1Mass flow rate5.73 kg/s
2Pressure ratio5.4:1
3Efficiency82%
4Maximum impeller diameter280 mm
5Maximum diffuser diameter340 mm
6Rotational speed38,000 rpm
7Inlet pressure101,325 Pa
8Inlet temperature288.15 K

Table 1.

Geometrical and functional specifications of analysed turbo-machinery.

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 [6]. 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 [6].

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 [7]. 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 β2=0°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.

Figure 2.

(a) Inlet Mach triangles for impeller blades only in radial section with backward swept blades, β2=0° with zero swirl at entry; (b) exit Mach triangles for impeller blades only in radial section with backward swept blades, β2=0° with zero swirl at entry.

5.2. Blade to blade contours

Blade to blade contours for absolute Mach number at hub, mean and tip sections are presented in below Figure 3 (a) and (c).

Figure 3.

(a) Hub section absolute Mach number contour, (b) mean section absolute Mach number contour, and (c) tip section absolute Mach number contour.

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.

Figure 4.

(a) Hub section relative Mach number contour, (b) mean section relative Mach number contour, and (c) tip section relative Mach number contour.

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

Circumferential relative Mach number contours at the impeller inlet and outlet are presented in Figure 5(a) and (b). The higher color temperature indicates higher values of relative Mach number.

Figure 5.

(a) Inlet relative Mach contour and (b) outlet relative Mach contour.

From above Figure 5(a), corresponding inlet relative Mach number contours shows that passage area is lying in the supersonic regime. The Figure 5(b) shows the effect of lean angle of impeller blades.

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 6.

Mass flow versus pressure ratio by Vista CCM.

Figure 7 represents the variation of polytropic efficiency with varying mass flow ingested to quantify differential pressure changes occurring through compressor stage.

Figure 7.

Mass flow versus polytropic efficiency by Vista CCM.

5.4.2. Stage performance characteristics

Compressor stage performance can be evaluated by obtaining individual characteristic curves as plotted in Figures 810. 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 810. The maximum achieved pressure ratio is 5.4 while allowing a fixed mass flow rate of 5.73 kg/s.

Figure 8.

Pressure ratio versus mass flow.

Figure 9.

Efficiency versus mass flow.

Figure 10.

Power versus mass flow.

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 obtained results presented in the Figures 11 and 12 shows state-of-the-art validity of current design with design performance work stated in the literature.

Figure 11.

Comparison with compressor in [8].

Figure 12.

Comparison with trend established in paper [9].

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.

6. Conclusion

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.

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Chellapilla V K N S N Moorthy and Vadapalli Srinivas (February 6th 2018). Numerical Simulations of a High-Resolution RANS-FVDM Scheme for the Design of a Gas Turbine Centrifugal Compressor, Numerical Simulations in Engineering and Science, Srinivas P. Rao, IntechOpen, DOI: 10.5772/intechopen.72098. Available from:

chapter statistics

231total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Susceptibility of Group-IV and III-V Semiconductor-Based Electronics to Atmospheric Neutrons Explored by Geant4 Numerical Simulations

By Daniela Munteanu and Jean-Luc Autran

Related Book

First chapter

Monte Carlo Simulations in NDT

By Frank Sukowski and Norman Uhlmann

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us