Main Content

Compressor (2P)

Two-phase fluid compressor in a thermodynamic cycle

Since R2022a

Libraries:
Simscape / Fluids / Two-Phase Fluid / Fluid Machines

Description

The Compressor (2P) block represents a dynamic compressor, such as a centrifugal or axial compressor, in a two-phase fluid network. You can use the Parameterization parameter to parameterize the block analytically based on the design point or by using a tabulated compressor map. A positive rotation of port R relative to port C causes fluid to flow from port A to port B. Port R and port C are mechanical rotational conserving ports associated with the compressor shaft and casing, respectively.

The surge margin is the ratio between the surge pressure ratio at a given mass flow rate and the operating point pressure ratio minus 1. When you set Parameterization to Tabulated, the block outputs the surge margin at port SM.

The design point is the design operational pressure ratio across and mass flow rate through the compressor during simulation. The compressor operating point and the point of maximum efficiency do not need to coincide.

The Compressor (2P) block assumes that superheated fluid enters the inlet. You can use the Report when fluid is not fully vapor parameter to choose what the block does when the fluid does not meet superheated conditions.

Compressor Map

The compressor map plots the isentropic efficiency of the compressor and the lines of constant corrected shaft speed between the two extremes of choked flow and surge flow. To visualize the block map, right-click the block and select Fluids > Plot Compressor Map. Each time you modify the block settings, click Reload Data on the figure window.

Compressor map

Each corrected speed line shows how the pressure ratio varies with the corrected mass flow rate when the shaft spins at the corresponding corrected speed. The variable β indicates the relative position along the corrected speed lines between the two extremes. Choked flow corresponds to β = 0 and surge flow corresponds to β = 1. The map also plots contours of isentropic efficiency as a function of pressure ratio and corrected mass flow rate, which provides a relative indication of how much power the compressor needs to operate at various combinations of pressure ratio and corrected mass flow rate.

Corrected Mass Flow Rate

Due to the large changes in pressure and temperature inside a compressor, the compressor map plots performance in terms of the corrected mass flow rate. The map adjusts the corrected mass flow rate from the inlet mass flow rate by using the reference pressure and reference temperature,

m˙corr=m˙A(TATref/pApref),

where:

  • A is the mass flow rate at port A.

  • TA is the temperature at port A.

  • Tref is the value of the Reference temperature for corrected flow parameter. When you set Parameterization to Analytical, this is the inlet temperature at the design operating condition.

  • m˙corr is the corrected mass flow rate. When you set Parameterization to Analytical, the block uses the Corrected mass flow rate at design point parameter. When you set Parameterization to Tabulated, the block uses the Corrected mass flow rate table, mdot(N,beta) parameter.

  • pA is the pressure at port A.

  • pref is the Reference pressure for corrected flow parameter. When using the analytical parameterization, this is the inlet pressure at the design operating condition.

The block derives TA from the specific internal energy, uA, and specific pressure, pA.

Corrected Speed

The block also adjusts the shaft speed, ω, according to the reference temperature, such that the corrected shaft speed is

ωcorr=ωTATref.

Shaft Torque

The block calculates the shaft torque, τ, as

τ=m˙AΔhtotalηmω,

where:

  • Δhtotal is the change in specific total enthalpy.

  • ηm is the compressor Mechanical efficiency.

  • ω is the relative shaft angular velocity, ωR - ωC.

The block relates the efficiency in the compressor map as

Δhtotal=Δhisenηisen,

where:

  • Δhisen is the isentropic change in specific total enthalpy.

  • ηisen is the isentropic efficiency.

A threshold region when flow approaches zero ensures that the compressor generates no torque when the flow rate is near zero or reversed.

Analytical Parameterization

You can generate the compressor map analytically by setting Parameterization to Analytical. The block fits a model of the compressor map based on [1] to the specified values for the Corrected speed at design point, Pressure ratio at design point, and Corrected mass flow rate at design point parameters. This method does not use β lines and the block does not report a surge margin.

Pressure Ratio

The block finds the pressure ratio at a given shaft speed and mass flow rate as

π=1+(πD1)[N˜ab+2N˜kln(1m˜N˜bk)],

where:

  • πD is the Pressure ratio at design point parameter.

  • N˜ is the normalized corrected shaft speed,

    NND,

    where ND is the Corrected speed at design point parameter.

  • m˜ is the normalized corrected mass flow rate,

    m˙corrm˙D,

    where m˙D is the Corrected mass flow rate at design point parameter.

  • a is the Spine shape, a parameter.

  • b is the Speed line spread, b parameter.

  • k is the Speed line roundness, k parameter.

The spine refers to the black line where the isentropic efficiency contours start to bend. The map speed lines are the shaft constant-speed lines that intersect the spine perpendicularly.

Compressor analytical map

Isentropic Efficiency Parameterization

When you set Efficiency specification to Analytical, the block models variable compressor efficiency as

η=η0(1C|p˜m˜a+Δa1m˜|cD|m˜m˜01|d),

where:

  • η0 is the value of the Maximum isentropic efficiency parameter.

  • C is the value of the Efficiency contour gradient orthogonal to spine, C parameter.

  • D is the value of the Efficiency contour gradient along spine, D parameter.

  • c is the value of the Efficiency peak flatness orthogonal to spine, c parameter.

  • d is the value of the Efficiency peak flatness along spine, d parameter.

  • p˜ is the normalized corrected pressure ratio,

    π1πD1,

    where πD is the Corrected pressure ratio at design point parameter.

  • m˜0 is the normalized corrected mass flow rate at which the compressor reaches the value of the Maximum isentropic efficiency parameter.

You can adjust the efficiency variables for different performance characteristics. Alternatively, you can choose a constant efficiency value by using the Constant isentropic efficiency parameter.

Tabulated Data Parameterization

When you set Parameterization to Tabulated, the isentropic efficiency, pressure ratio, and corrected mass flow rate of the compressor are functions of the corrected speed, N, and the map index, β. The block uses linear interpolation between data points for the efficiency, pressure ratio, and corrected mass flow rate.

If β exceeds 1, compressor surge occurs and the block assumes the pressure ratio remains at β = 1, while the mass flow rate continues to change. If the simulation conditions fall below β = 0, the block models the effects of choked flow, where the mass flow rate remains at its value at β = 0 while the pressure ratio continues to change. To constrain the compressor performance in the map boundaries, the block extrapolates isentropic efficiency to the nearest point.

You can choose to be notified when the operating point pressure ratio exceeds the surge pressure ratio. Set Report when surge margin is negative to Warning to receive a warning or to Error to stop the simulation when this occurs.

Continuity Equations

The block conserves mass such that

m˙A+m˙B=0,

where m˙B is the mass flow rate at port B.

The block computes the energy balance equation as

ΦA+ΦB+Pfluid=0,

where:

  • ΦA is the energy flow rate at port A.

  • ΦB is the energy flow rate at port B.

  • Pfluid is the hydraulic power delivered to the fluid,which the block determines from the change in specific enthalpy Pfluid=m˙AΔhtotal.

Assumptions and Limitations

  • The block assumes that superheated fluid enters at A.

  • The block only defines compressor map flow from port A to port B. Reverse flow results may not be accurate.

  • The block only represents dynamic compressors.

Ports

Output

expand all

Physical signal output port associated with the surge margin at a given mass flow rate. The block calculates the surge margin as

SM(m˙corr)=pr,surge(m˙corr)pr(m˙corr)1.

Dependencies

To enable this port, set Parameterization to Tabulated.

Conserving

expand all

Two-phase fluid conserving port associated with the compressor inlet.

Two-phase fluid conserving port associated with the compressor outlet.

Mechanical rotational conserving port associated with the compressor case.

Mechanical rotational conserving port associated with the compressor shaft.

Parameters

expand all

Compressor Map

Compressor performance parameterization. Choose the analytical or tabulated parameterization:

  • Analytical: A pressure ratio–corrected mass flow rate curve defines the peak compressor performance. You can choose to model isentropic efficiency as a constant or analytically.

  • Tabulated: A user-supplied compressor map defines the compressor performance. The block determines the compressor operating points by using linear interpolation between the corrected mass flow rate, pressure ratio, and isentropic efficiency tables at given points in the user-provided corrected shaft speed and β vectors. The default map comes from data reported in [3].

Shaft speed at the intended compressor pressure ratio and corrected mass flow rate, corrected for temperature.

Dependencies

To enable this parameter, set Parameterization to Analytical.

Outlet-to-inlet pressure ratio at the intended compressor corrected mass flow rate.

Dependencies

To enable this parameter, set Parameterization to Analytical.

Mass flow rate at the intended compressor pressure ratio, corrected for temperature and pressure.

Dependencies

To enable this parameter, set Parameterization to Analytical.

Isentropic efficiency model type. Choose a constant or analytical model. Use Analytical to specify variable values.

Dependencies

To enable this parameter, set Parameterization to Analytical.

Maximum compressor isentropic efficiency. Isentropic efficiency is the ratio of isentropic change to actual change.

Dependencies

To enable this parameter, set Parameterization to Analytical and Efficiency specification to Analytical.

Minimum compressor isentropic efficiency. Isentropic efficiency is the ratio of isentropic change to actual change.

Dependencies

To enable this parameter, set Parameterization to Analytical and Efficiency specification to Analytical.

Mass flow rate at maximum efficiency, corrected for temperature and pressure. The point of maximum efficiency does not necessarily coincide with the compressor design point.

Dependencies

To enable this parameter, set Parameterization to Analytical and Efficiency specification to Analytical.

Pressure ratio at maximum efficiency. The point of maximum efficiency does not necessarily coincide with the compressor design point.

Dependencies

To enable this parameter, set Parameterization to Analytical and Efficiency specification to Analytical.

Value of the constant isentropic efficiency.

Dependencies

To enable this parameter, set Parameterization to Analytical and Efficiency specification to Constant.

Vector of corrected shaft speeds.

Dependencies

To enable this parameter, set Parameterization to Tabulated.

Vector of relative positions along the corrected speed lines. Choked flow occurs when β = 0 and surge flow occurs at β = 1. β lines are perpendicular to the compressor shaft constant speed lines, N.

Dependencies

To enable this parameter, set Parameterization to Tabulated.

M-by-N matrix of compressor outlet-to-inlet pressure ratios at the specified corrected shaft speed and β value. The block employs linear interpolation between table elements. M and N are the sizes of the corresponding vectors:

  • M is the number of vector elements in the Corrected speed index vector, N parameter.

  • N is the number of vector elements in the Beta index vector, beta parameter.

The default value is:

[1.1814, 1.2385, 1.2792, 1.3017, 1.3057, 1.3057; 1.3648, 1.5157, 1.6298, 1.691, 1.7011, 1.7032; 1.587, 1.8357, 2.011, 2.0925, 2.1006, 2.1006; 1.8255, 2.28, 2.4777, 2.5225, 2.5633, 2.5389]

Dependencies

To enable this parameter, set Parameterization to Tabulated.

M-by-N matrix of corrected mass flow rates at the specified corrected shaft speed and β value. The block employs linear interpolation between table elements. M and N are the sizes of the corresponding vectors:

  • M is the number of vector elements in the Corrected speed index vector, N parameter.

  • N is the number of vector elements in the Beta index vector, beta parameter.

The default value is:

[.1503, .1313, .11, .0887, .0668, .0502; .2265, .1999, .1649, .1322, .0992, .0639; .2869, .2545, .216, .1798, .1424, .108; .3275, .2846, .2466, .2285, .209, .1661] kg/s

Dependencies

To enable this parameter, set Parameterization to Tabulated.

M-by-N matrix of compressor isentropic efficiencies at the specified corrected shaft speed and β value. The block employs linear interpolation between table elements. M and N are the sizes of the corresponding vectors:

  • M is the number of vector elements in the Corrected speed index vector, N parameter.

  • N is the number of vector elements in the Beta index vector, beta parameter.

The default value is:

[.56, .66, .71, .695, .659, .635; .558, .682, .755, .743, .697, .638; .57, .705, .765, .752, .712, .652; .552, .718, .755, .752, .736, .67]

Dependencies

To enable this parameter, set Parameterization to Tabulated.

Whether the block does nothing, generates a warning, or generates an error when it detects a negative surge margin.

Dependencies

To enable this parameter, set Parameterization to Tabulated.

Map Coefficients

To enable the Map Coefficients parameters, set Parameterization to Analytical.

Exponent in the analytical parameterization of the compressor map that characterizes the spine shape.

Dependencies

To enable this parameter, set Parameterization to Analytical.

Exponent in the analytical parameterization of the compressor map that characterizes the constant speed line spacing.

Dependencies

To enable this parameter, set Parameterization to Analytical.

Coefficient in the analytical parameterization of the compressor map that characterizes the constant speed line shape.

Dependencies

To enable this parameter, set Parameterization to Analytical.

Exponent associated with the compressor map peak flatness orthogonal to the spine.

Dependencies

To enable this parameter, set Parameterization to Analytical and Efficiency specification to Analytical.

Exponent associated with the efficiency map peak flatness along the spine.

Dependencies

To enable this parameter, set Parameterization to Analytical and Efficiency specification to Analytical.

Coefficient associated with the efficiency contour gradient orthogonal to the spine.

Dependencies

To enable this parameter, set Parameterization to Analytical and Efficiency specification to Analytical.

Coefficient associated with the efficiency contour gradient along the spine.

Dependencies

To enable this parameter, set Parameterization to Analytical and Efficiency specification to Analytical.

Reference data

Reference inlet pressure for the compressor map. When you set Parameterization to Tabulated, the data supplier specifies this value with respect to the tabulated compressor data. When you set Parameterization to Analytical, this value is the inlet pressure at the design operating condition.

Reference inlet temperature for the compressor map. When you set Parameterization to Tabulated, the data supplier specifies this value with respect to the tabulated compressor data. When you set Parameterization to Analytical, this value is the inlet temperature at the design operating condition.

Ratio of the power delivered to the fluid flow to the power driving the mechanical shaft.

Compressor inlet cross-sectional area.

Compressor outlet cross-sectional area.

Whether the block does nothing, generates a warning, or generates an error when it detects fluid that is not fully vapor.

References

[1] Greitzer, E. M. et al. “N+3 Aircraft Concept Designs and Trade Studies. Volume 2: Appendices – Design Methodologies for Aerodynamics, Structures, Weight, and Thermodynamic Cycles.” NASA Technical Report, 2010.

[2] Kurzke, Joachim. "How to Get Component Maps for Aircraft Gas Turbine Performance Calculations." Volume 5: Manufacturing Materials and Metallurgy; Ceramics; Structures and Dynamics; Controls, Diagnostics and Instrumentation; Education; General, American Society of Mechanical Engineers, 1996, p. V005T16A001.

[3] Plencner, Robert M. “Plotting component maps in the Navy/NASA Engine Program (NNEP): A method and its usage.” NASA Technical Memorandum, 1989.

Extended Capabilities

C/C++ Code Generation
Generate C and C++ code using Simulink® Coder™.

Version History

Introduced in R2022a