Modelling of Real Gas Compression Using the Finite Volume Thermofluid Network Methodology and Mean-Line Approach for Centrifugal Compressors

Abstract

A novel compressible real gas 1D finite volume thermofluid network modelling methodology is proposed that can be applied in the analysis of supercritical CO2 (sCO2) power cycles. It is applied to simulate centrifugal compressor performance via a mean line approach that directly solves the mass, energy, and momentum balance equations for the discretized flow channels within the machine along with the associated loss models, component characteristics, and applicable compressible real gas equation of state. This can be integrated seamlessly with the rest of the power cycle components like pipes, valves, heat exchangers, etc., thereby eliminating the need for static performance maps and the associated difficulties to ensure similarity over the whole range of operation. This direct simulation of the internal flows also makes it possible to do integrated cycle design optimization that includes the geometry of the compressors, and to simulate fast transients for dynamic analysis. The methodology was applied to simulate well-documented air and sCO2 centrifugal compressors and compared the generated results with that of conventional thermofluid network incompressible and compressible ideal gas approaches. The results indicate that the proposed compressible real gas approach can capture the stagnation pressure changes and isentropic efficiencies with good accuracy for both cases over the whole range of mass flow rates and shaft rotational speeds. In contrast, the conventional thermofluid network incompressible and compressible ideal gas approaches are only applicable in selected regions of the fluid property regime and suffer from discretization errors and inappropriate treatment of the relationship between stagnation and static pressures. Furthermore, comparative analyses between the newly proposed approach and the traditional mechanical energy balance method are performed, showing that the latter approach needs to be discretized more finely through the compression process to correctly capture the polytropic compression path shape, whereas the former more accurately captures the stagnation pressure changes through the analyzed compressors with only a single increment per compressor component.