Research Papers: Gas Turbines: Microturbines and Small Turbomachinery

Computational Fluid Dynamics of a Radial Compressor Operating With Supercritical CO2

[+] Author and Article Information
Rene Pecnik

e-mail: r.pecnik@tudelft.nl

Piero Colonna

Process and Energy Department,
Delft University of Technology,
Leeghwaterstraat 44, 2628 CA Delft,
The Netherlands

1Corresponding author.

Contributed by International Gas Turbine Institute (IGTI) of ASME for publication in the JOURNAL OF ENGINEERING FOR GAS TURBINES AND POWER. Manuscript received June 19, 2012; final manuscript received July 5, 2012; published online October 11, 2012. Editor: Dilip R. Ballal.

J. Eng. Gas Turbines Power 134(12), 122301 (Oct 11, 2012) (8 pages) doi:10.1115/1.4007196 History: Received June 19, 2012; Revised July 05, 2012

The merit of using supercritical CO2(scCO2) as the working fluid of a closed Brayton cycle gas turbine is now widely recognized, and the development of this technology is now actively pursued. scCO2 gas turbine power plants are an attractive option for solar, geothermal, and nuclear energy conversion. Among the challenges that must be overcome in order to successfully bring the technology to the market is that the efficiency of the compressor and turbine operating with the supercritical fluid should be increased as much as possible. High efficiency can be reached by means of sophisticated aerodynamic design, which, compared to other overall efficiency improvements, like cycle maximum pressure and temperature increase, or increase of recuperator effectiveness, does not require an increase in equipment cost, but only an additional effort in research and development. This paper reports a three-dimensional computational fluid dynamics (CFD) study of a high-speed centrifugal compressor operating with CO2 in the thermodynamic region slightly above the vapor–liquid critical point. The investigated geometry is the compressor impeller tested in the Sandia scCO2 compression loop facility. The fluid dynamic simulations are performed with a fully implicit parallel Reynolds-averaged Navier–Stokes code based on a finite volume formulation on arbitrary polyhedral mesh elements. In order to account for the strongly nonlinear variation of the thermophysical properties of supercritical CO2, the CFD code is coupled with an extensive library for the computation of properties of fluids and mixtures. A specialized look-up table approach and a meshing technique suited for turbomachinery geometries are also among the novelties introduced in the developed methodology. A detailed evaluation of the CFD results highlights the challenges of numerical studies aimed at the simulation of technically relevant compressible flows occurring close to the liquid–vapor critical point. The data of the obtained flow field are used for a comparison with experiments performed at the Sandia scCO2 compression-loop facility.

Copyright © 2012 by ASME
Your Session has timed out. Please sign back in to continue.


Angelino, G., 1969, “Carbon Dioxide Condensation Cycles for Power Production,” ASME J. Eng. Power Trans., 10, pp. 272–287.
Feher, E. G., 1968, “The Supercritical Thermodynamic Power Cycle,” Energ. Convers., 8, pp. 85–90. [CrossRef]
Dostal, V., Driscoll, M. J., and Hejzlar, P., 2004, “A Supercritical Carbon Dioxide Cycle for Next Generation Nuclear Reactors,” MIT-ANP-Series, Massachusetts Institute of Technology, Boston, MA, MIT-ANP-TR-100.
Wright, S. A., Radel, R. F., Vernon, M. E., Rochau, G. E., and Pickard, P. S., 2010, “Operation and Analysis of a Supercritical CO2 Brayton Cycle,” Sandia Report 2010–0171.
Colonna, P., Harinck, J., Rebay, S., and Guardone, A., 2008, “Real-Gas Effects in Organic Rankine Cycle Turbine Nozzles,” J. Propul. Power, 24(2), pp. 282–294. [CrossRef]
Harinck, J., Colonna, P., Guardone, A., and Rebay, S., 2010, “Influence of Thermodynamic Models in 2D Flow Simulations of Turboexpanders,” ASME J. Turbomach., 132(1), p. 011001. [CrossRef]
Takagi, K., Muto, Y., Ishizuka, T., Kikura, H., and Aritomi, M., 2010, “Research on Flow Characteristics of Supercritical CO2 Axial Compressor Blades by CFD Analysis,” J. Power Energ. Syst., 4(1), pp. 138–149. [CrossRef]
Pecnik, R., Terrapon, V. E., Ham, F., Iaccarino, G., and Pitsch, H., 2012, “Reynolds-Averaged Navier–Stokes Simulations of the Hyshot II Scramjet,” AIAA J., 50(8), pp. 1717–1732. [CrossRef]
Pečnik, R., Witteveen, J. A. S., and Iaccarino, G., 2011, “Uncertainty Quantification for Laminar-Turbulent Transition Prediction in RANS Turbomachinery Applications,” 49th AIAA Aerospace Sciences Meeting, Orlando, FL, Jan. 4–7, AIAA Paper No. 2011–660.
Spalart, P. R., and Allmaras, S. R., 1992, “A One-Equation Turbulence Model for Aerodynamic Flows,” AIAA Paper No. 92-0439.
Menter, F. R., 1994, “Two-Equation Eddy-Viscosity Turbulence Models for Engineering Applications,” AIAA J., 32(8), pp. 1598–1605. [CrossRef]
Satish, B., Buschelman, K., Eijkhout, V., Gropp, W. D., Kaushik, D., Knepley, M. G., McInnes, L. C., Smith, B. F., and Zhang, H., 2009, PETSc Web page, http://www.mcs.anl.gov/petsc
Berger, M., Aftosmis, M. J., and Murman, S. M., 2005, “Analysis of Slope Limiters on Irregular Grids,” 43rd AIAA Aerospace Sciences Meeting, Reno, NV, Jan. 10-13, AIAA Paper No. 2005-0490.
Barth, T. J., and Jespersen, D. C., 1989, “The Design and Application of Upwind Schemes on Unstructured Meshes,” AIAA Paper No. 89-0366.
Venkatakrishnan, V., 1995, “Convergence to Steady State Solutions of the Euler Equations on Unstructures Grids With Limiters,” J. Comput. Phys., 118(1), pp. 120–130. [CrossRef]
Roe, P. L., 1981, “Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes,” J. Comput. Phys., 43(2), pp. 357–372. [CrossRef]
Candler, G. V., Barnhardt, M. D., Drayna, T. M., Nompelis, I., Peterson, D. M., and Subbareddy, P., 2007, “Unstructured Grid Approaches for Accurate Aeroheating Simulations,” AIAA Paper No. 2007-3959.
Quirk, J. J., 1994, “A Contribution to the Great Riemann Solver Debate,” Int. J. Numer. Meth. Fluids, 18(6), pp. 555–574. [CrossRef]
Van Leer, B., Lee, W. T., and Powell, K. G., 1989, “Sonic-Point Capturing,” AIAA Paper No. 89-1945.
Nompelis, I., Drayna, T. W., and Candler, G. V., 2005, “A Parallel Unstructured Implicit Solver for Hypersonic Reacting Flow Simulation,” AIAA Paper No. 2005–4867.
Liou, M.-S., 1996, “A Sequel to AUSM: AUSM+,” J. Comput. Phys., 129(2), pp. 364–382. [CrossRef]
Ham, F., and Iaccarino, G., 2004, “Energy Conservation in Collocated Discretization Schemes on Unstructured Meshes,” Annual Research Briefs 2004, Center for Turbulence Research, NASA Ames/Stanford University, Stanford, CA.
Colonna, P., and Van der Stelt, T., 2005, “Fluidprop: A Program for the Estimation of Thermophysical Properties of Fluids,” Energy Technology Section, Delft University of Technology, Delft, The Netherlands, available at: http://www.fluidprop.com
Lemmon, E. W., McLinden, M. O., and Huber, M. L., 2002, “NIST Standard Reference Database 23: Reference Fluid Thermodynamic and Transport Properties-REFPROP, Version 8.0,” National Institute of Standards and Technology, Standard Reference Data Program, Gaithersburg, MD.
Span, R., and Wagner, W., 2003, “Equations of State for Technical Applications. I. Simultaneously Optimized Functional Forms for Nonpolar and Polar Fluids,” Int. J. Thermophys., 24, pp. 1–39. [CrossRef]
Span, R., and Wagner, W., 2003, “Equations of State for Technical Applications. II. Results for Nonpolar Fluids,” Int. J. Thermophys., 24, pp. 41–109. [CrossRef]
ANSYS, 2009, “ANSYS BladeGen, Release 13.0 Users Guide,” ANSYS, Inc., Canonsburg, PA.
Shewchuk, J., 1996, “Triangle: Engineering a 2D Quality Mesh Generator and Delaunay Triangulator,” Applied Computational Geometry Towards Geometric Engineering (Lecture Notes in Computer Science), Vol. 1148, M.Lin and D.Manocha, eds., Springer, Berlin/Heidelberg, pp. 203–222.
Fuller, R. L., and Eisemann, K., 2011, “Centrifugal Compressor Off-Design Performance for Super-Critical CO2,” Supercritical CO2 Power Cycle Symposium, Boulder, CO, May.
Bae, J., Yoo, J., and Choi, H., 2005, “Direct Numerical Simulation of Turbulent Supercritical Flows With Heat Transfer,” Phys. Fluids, 17(10), pp. 1–24. [CrossRef]
Sokolichin, A., Eigenberger, G., Lapin, A., and Lübbert, A., 1997, “Dynamic Numerical Simulation of Gas-Liquid Two-Phase Flows: Euler/Euler Versus Euler/Lagrange,” Chem. Eng. Sci., 52(4), pp. 611–626. [CrossRef]
Blesgen, T., 1999, “Generalization of the Navier–Stokes Equations to Two-Phase Flows,” J. Phys. D Appl. Phys., 32(10), pp. 1119–1123. [CrossRef]
Nikkhahi, B., Shams, M., and Ziabasharhagh, M., 2009, “A Numerical Investigation of Two-Phase Steam Flow Around a 2-D Turbine’s Rotor Tip,” Int. Commun. Heat Mass Transf., 36(6), pp. 632–639. [CrossRef]


Grahic Jump Location
Fig. 1

The tabulated region in the temperature-density plane and schematic of the interpolation; (a) the tabulated region in the temperature-density plane, with the two tables which are separated by the saturation lines; (b) schematic of the interpolation applied in the look-up table

Grahic Jump Location
Fig. 2

Look-up table convergence and interpolation error using the density ρ and the internal energy e. The number of intervals is the same for the two independent variables. Dashed line represents a bilinear interpolation scheme, while the solid line shows the convergence rate for the interpolation scheme using the gradients at each surrounding node.

Grahic Jump Location
Fig. 3

Centrifugal compressor geometry and computational mesh

Grahic Jump Location
Fig. 4

Simulation results for P2 = 90 bar and shaft speed = 50 krpm plotted in the T-s diagram. (a) Averaged states of the compression process. The starting and the end state (black squares) of the compression are obtained by averaging the inlet and the outlet plane of the computational domain. (b) The point could represents all control volumes in the computational domain plotted in the T-s plane.

Grahic Jump Location
Fig. 5

The red surfaces show the two phase regions at the blade tip and the blade trailing edges for P2 = 90 bar and shaft speed = 50 krpm

Grahic Jump Location
Fig. 6

Pressure contour plots in planes at different span locations for P2 = 90 bar and shaft speed = 50 krpm

Grahic Jump Location
Fig. 7

Impeller performance for 50 krpm shaft speed



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In