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.
[1]
T. Blesgen,et al.
A generalization of the Navier-Stokes equations to two-phase flows
,
1999
.
[2]
Roland Span,et al.
Equations of State for Technical Applications. I. Simultaneously Optimized Functional Forms for Nonpolar and Polar Fluids
,
2003
.
[3]
Gerhart Eigenberger,et al.
Dynamic Numerical Simulation of Gas-Liquid Two Phase Flows, Euler-Euler versus Euler-Lagrange
,
1997
.
[4]
Masoud Ziabasharhagh,et al.
A numerical investigation of two-phase steam flow around a 2-D turbine's rotor tip☆
,
2009
.
[5]
E. Feher.
SUPERCRITICAL THERMODYNAMIC POWER CYCLE.
,
1967
.
[6]
J. Yoo,et al.
Direct numerical simulation of turbulent supercritical flows with heat transfer
,
2005
.
[7]
Stefano Rebay,et al.
Real-gas effects in Organic Rankine Cycle turbine nozzles
,
2008
.
[8]
J. Quirk.
A Contribution to the Great Riemann Solver Debate
,
1994
.
[9]
V. Venkatakrishnan.
Convergence to steady state solutions of the Euler equations on unstructured grids with limiters
,
1995
.
[10]
P. Roe.
Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes
,
1997
.
[11]
M. Liou.
A Sequel to AUSM
,
1996
.
[12]
F. Menter.
Two-equation eddy-viscosity turbulence models for engineering applications
,
1994
.
[13]
G. Angelino.
Carbon Dioxide Condensation Cycles For Power Production
,
1968
.
[14]
Jonathan Richard Shewchuk,et al.
Triangle: Engineering a 2D Quality Mesh Generator and Delaunay Triangulator
,
1996,
WACG.
[15]
Hiroshige Kikura,et al.
Research on Flow Characteristics of Supercritical CO2 Axial Compressor Blades by CFD Analysis
,
2010
.
[16]
Heinz Pitsch,et al.
Reynolds-Averaged Navier-Stokes Simulations of the HyShot II Scramjet
,
2012
.
[17]
Roland Span,et al.
Equations of State for Technical Applications. II. Results for Nonpolar Fluids
,
2003
.