5.5 Turbomachinery, Compressible Flow & Nozzles, Renewables & Energy Storage

Key Takeaways

  • The Euler Turbine Equation $\dot{W} = \dot{m}(U_1 V_{u1} - U_2 V_{u2})$ governs angular momentum exchange in all dynamic turbomachines.
  • A converging nozzle chokes when back pressure falls to or below the critical ratio $P^{\ast}/P_0 = (2/(k+1))^{k/(k-1)} = 0.528$ for air ($k = 1.40$); beyond choking, mass flow responds only to reservoir pressure, temperature, or throat area — and supersonic exit flow requires a converging–diverging (de Laval) nozzle with sonic throat.
  • Hydroelectric turbine selection is governed by specific speed ($N_s$): Pelton wheels for high head ($>300\text{ m}$), Francis turbines for medium head ($30-300\text{ m}$), and Kaplan propeller turbines for low head ($<30\text{ m}$).
  • The Betz Limit establishes that an ideal open-rotor wind turbine can extract a maximum of $16/27 \approx 59.3\%$ of kinetic energy from wind: $P_{max} = \frac{1}{2} \rho A V^3 C_{p,max}$.
  • Plant Capacity Factor ($CF$) measures annual energy generated against maximum continuous rated nameplate output over 8,760 hours per year.
Last updated: August 2026

Turbomachinery Principles, Renewable Energy & Energy Storage

Modern power engineering integrates conventional turbomachinery with renewable energy systems (hydroelectric, wind, solar) and grid-scale energy storage. The NCEES PE Mechanical exam tests velocity triangle analysis, turbine selection by specific speed, nuclear core fundamentals, aerodynamic wind power extraction, and energy storage capacity factor economics.


1. Turbomachinery Principles & Velocity Triangles

Turbomachinery transfers energy between a continuous fluid stream and a rotating blade rotor through dynamic angular momentum exchange.

+-----------------------------------------------------------------------------------------+
|                              TURBOMACHINERY VELOCITY VECTORS                            |
|                                                                                         |
|   Blade (Tangential) Speed:     U = omega * r = 2 * pi * N * r                          |
|   Absolute Fluid Velocity:      V  (Velocity relative to stationary casing)             |
|   Relative Fluid Velocity:      W  (Velocity relative to moving blade rotor)           |
|   Vector Sum:                   V_vector = U_vector + W_vector                          |
|   Whirl (Tangential) Component: V_u = V * cos(alpha)                                    |
|   Axial / Meridional Component: V_a = V * sin(alpha)                                    |
+-----------------------------------------------------------------------------------------+
          Rotor Inlet Velocity Triangle              Rotor Exit Velocity Triangle
                V_u1                                       V_u2
            +----------+                               +----------+
            |         /|                               |         /|
            |        / |                               |        / |
        V_a1|    V1 /  | W1                        V_a2|    V2 /  | W2
            |      /   |                               |      /   |
            | alpha1\  | beta1                         | alpha2\  | beta2
            +--------+-+                               +--------+-+
               U1                                         U2

The Euler Turbine Equation

Applying the principle of conservation of angular momentum across a rotating runner yields the Euler Turbomachinery Equation:

Specific Work w=U1Vu1U2Vu2\text{Specific Work } w = U_1 V_{u1} - U_2 V_{u2}

Shaft Power W˙=m˙(U1Vu1U2Vu2)\text{Shaft Power } \dot{W} = \dot{m} \left(U_1 V_{u1} - U_2 V_{u2}\right)

Impulse vs. Reaction Blading

+-----------------------------------------------------------------------------------------+
|                        IMPULSE VS REACTION STAGE COMPARISON                             |
|                                                                                         |
|   IMPULSE STAGE (Degree of Reaction R = 0):                                             |
|   - Complete fluid pressure drop occurs entirely across stationary nozzles/vanes.       |
|   - Pressure across moving rotor blades remains CONSTANT (P_inlet = P_exit).            |
|   - Relative velocity magnitude is constant across rotor: W_2 = W_1 (ignoring friction).|
|   - Optimum blade speed ratio: (U / V_1)_opt = cos(alpha_1) / 2  (~0.45 - 0.50).        |
|                                                                                         |
|   50% REACTION STAGE (Parsons Stage, Degree of Reaction R = 0.50):                      |
|   - Pressure drops equally across stationary stators and rotating rotor blades.         |
|   - Rotor passages act as moving converging nozzles, accelerating relative flow W_2 > W_1|
|   - Symmetric velocity triangles: V_1 = W_2,  W_1 = V_2,  alpha_1 = beta_2.             |
|   - Optimum blade speed ratio: (U / V_1)_opt = cos(alpha_1)  (~0.85 - 0.90).            |
|   - Requires twice as many stages as impulse, but achieves higher peak efficiency.      |
+-----------------------------------------------------------------------------------------+

Compressor Operational Limits: Surge, Stall & Choke

In dynamic compressors (axial and centrifugal), flow stability is bounded by three critical aerodynamic limits:

+-----------------------------------------------------------------------------------------+
|                         COMPRESSOR INSTABILITY MECHANISMS                               |
|                                                                                         |
|   SURGE:           Violent, low-frequency system-wide flow reversal and pressure        |
|                    oscillations across the entire compressor and duct network.          |
|                    Occurs at low mass flow rates; causes catastrophic mechanical damage.|
|                                                                                         |
|   ROTATING STALL:  Localized aerodynamic boundary layer separation on blade suction     |
|                    surfaces that rotates around the annulus at 30-70% rotor speed.      |
|                    Causes severe blade fatigue vibration without total flow reversal.   |
|                                                                                         |
|   CHOKE (STONEWALL):Mass flow rate reaches maximum limit when local velocity reaches     |
|                    Mach 1.0 (sonic velocity) at the blade throat passages.              |
+-----------------------------------------------------------------------------------------+

2. Compressible Gas Flow: Mach Number, Nozzles & Diffusers

When gas velocity approaches the local speed of sound, density can no longer be treated as constant, and the incompressible Bernoulli/Darcy machinery of Chapter 4 fails. The Thermal and Fluid Systems specification explicitly lists Mach number, shock properties, and isentropic flow through compressors, nozzles, and diffusers.

The Mach Number & Speed of Sound

c=kRTM=Vcc = \sqrt{k R T} \qquad M = \frac{V}{c}

where $k = c_p / c_v$ (1.40 for air), $R$ is the gas constant ($287\text{ J/(kg·K)} = 53.35\text{ ft·lbf/(lbm·°R)}$ for air), and $T$ is absolute temperature. Regimes: $M < 1$ subsonic, $M = 1$ sonic, $M > 1$ supersonic.

Isentropic Stagnation (Total) Property Relations

For isentropic flow relative to stagnation temperature $T_0$ and stagnation pressure $P_0$:

T0T=1+k12M2P0P=(1+k12M2)kk1ρ0ρ=(1+k12M2)1k1\frac{T_0}{T} = 1 + \frac{k-1}{2} M^2 \qquad \frac{P_0}{P} = \left(1 + \frac{k-1}{2} M^2\right)^{\frac{k}{k-1}} \qquad \frac{\rho_0}{\rho} = \left(1 + \frac{k-1}{2} M^2\right)^{\frac{1}{k-1}}

Stagnation pressure is the pressure reached if the fluid is brought to rest isentropically. Stagnation temperature is constant through adiabatic machinery — even across a shock wave — but stagnation pressure drops irreversibly through a shock.

Critical Pressure Ratio & Choked Flow

As a converging nozzle exhausts from a stagnation reservoir, mass flux rises as the back pressure $P_b$ falls — until the throat reaches sonic velocity ($M = 1$). Beyond that point the flow is choked: further reducing $P_b$ cannot increase mass flow. The choking threshold is the critical pressure ratio:

PP0=(2k+1)kk1=0.528for air (k=1.40)\frac{P^{\ast}}{P_0} = \left(\frac{2}{k+1}\right)^{\frac{k}{k-1}} = 0.528 \quad \text{for air } (k = 1.40)

ConditionConsequence
$P_b / P_0 > 0.528$Throat flow subsonic; mass flow depends on back pressure
$P_b / P_0 \le 0.528$Choked; $\dot{m} = \dot{m}_{max}$ fixed by $P_0$, $T_0$, and throat area only

[!WARNING] Classic trap: once a nozzle is choked, lowering the downstream pressure does not increase mass flow. Only raising the reservoir pressure $P_0$, increasing the throat area, or lowering $T_0$ increases $\dot{m}_{max}$.

Converging-Diverging (de Laval) Nozzles — the Area Rule Inverts

For subsonic flow, converging passages accelerate flow and diverging passages decelerate it (a diffuser). Once the throat goes sonic, the area–velocity coupling inverts: a supersonic gas accelerates in a diverging passage. A converging–diverging nozzle therefore reaches supersonic exit velocity only if the throat is at $M = 1$ — and between the design and first-critical back pressures, a normal shock stands inside the divergent section instead.

Normal Shock Waves (Qualitative Facts the Exam Tests)

  • Upstream Mach $M_1 > 1$; downstream Mach $M_2 < 1$ — flow through a normal shock is always supersonic-to-subsonic.
  • Static pressure, temperature, and density rise discontinuously across the shock.
  • The shock is adiabatic but irreversible: stagnation temperature is unchanged while stagnation pressure falls and entropy rises ($\Delta s > 0$).

Nozzles vs. Diffusers from the Steady-Flow Energy Equation

With $\dot{Q} \approx 0$ and $\dot{W} = 0$, SFEE reduces to $h_1 + V_1^2/2 = h_2 + V_2^2/2$:

  • Nozzle: pressure drops, velocity rises, enthalpy falls ($h_2 < h_1$). Isentropic nozzle efficiency $\eta_N = \frac{h_1 - h_{2a}}{h_1 - h_{2s}}$.
  • Diffuser: pressure rises, velocity falls, enthalpy rises. Diffuser effectiveness is assessed against the ideal pressure-recovery path.

3. Hydroelectric Power & Hydraulic Turbine Selection

Hydroelectric stations convert gravitational potential energy of water into shaft work.

Total Dynamic Head Hnet=(z1z2)+P1P2ρg+V12V222ghloss\text{Total Dynamic Head } H_{net} = \left(z_1 - z_2\right) + \frac{P_1 - P_2}{\rho g} + \frac{V_1^2 - V_2^2}{2g} - h_{loss}

Hydroelectric Power Output W˙hydro=ηoverallρgQHnet\text{Hydroelectric Power Output } \dot{W}_{hydro} = \eta_{overall} \rho g Q H_{net}

Where $\rho \approx 1000\text{ kg/m}^3$ ($62.4\text{ lbm/ft}^3$), $g = 9.81\text{ m/s}^2$ ($32.174\text{ ft/s}^2$), $Q$ is volumetric flow rate ($\text{m}^3/\text{s}$ or $\text{ft}^3/\text{s}$), and $\eta_{overall} = \eta_{turbine} \cdot \eta_{generator}$.

+-----------------------------------------------------------------------------------------+
|                         HYDRAULIC TURBINE SELECTION MATRIX                              |
|                                                                                         |
|   TURBINE TYPE       HEAD REGIME (H)         FLOW RATE (Q)     SPECIFIC SPEED (N_s)     |
|   ----------------   ---------------------   ---------------   --------------------     |
|   PELTON (Impulse)   High (>300 m / 1000 ft) Low Flow          Low (2 - 20)             |
|   FRANCIS (Reaction) Medium (30 - 300 m)     Medium Flow       Medium (20 - 120)        |
|   KAPLAN (Axial)     Low (<30 m / 100 ft)    High Flow         High (100 - 250+)        |
+-----------------------------------------------------------------------------------------+

Specific Speed Formula

In US Customary units ($N$ in $\text{rpm}$, $P$ in $\text{bhp}$, $H$ in $\text{ft}$):

Ns=NPH5/4N_s = \frac{N \sqrt{P}}{H^{5/4}}


4. Nuclear Power Engineering: Core Fundamentals & Plant Types

Nuclear power generates thermal heat via the controlled neutron-induced fission of uranium-235 ($^{235}\text{U}$) or plutonium-239 ($^{239}\text{Pu}$):

92235U+01nthermalFission Fragments+2.43  01nfast+200 MeV{}^{235}_{92}\text{U} + {}^1_0\text{n}_{thermal} \longrightarrow \text{Fission Fragments} + 2.43 \; {}^1_0\text{n}_{fast} + 200\text{ MeV}

+-----------------------------------------------------------------------------------------+
|                             NUCLEAR REACTOR CORE COMPONENTS                             |
|                                                                                         |
|   FUEL RODS:     Ceramic UO2 pellets (enriched to 3-5% U-235) sealed in Zircaloy tubes. |
|   MODERATOR:     Light water (H2O), heavy water (D2O), or graphite that slows high-speed|
|                  fast neutrons (2 MeV) to thermal energy levels (0.025 eV) via elastic  |
|                  collisions to sustain fission chain reactions.                         |
|   CONTROL RODS:  Neutron absorbers (Boron carbide B4C, Cadmium, Silver-Indium-Cadmium) |
|                  inserted to regulate reactor power or SCRAM (emergency shutdown).      |
|   COOLANT:       Fluid (pressurized water) circulating through the core to extract heat.|
+-----------------------------------------------------------------------------------------+
       Pressurized Water Reactor (PWR)               Boiling Water Reactor (BWR)
   +---------------------------------------+     +---------------------------------------+
   | TWO-LOOP SYSTEM:                      |     | SINGLE-LOOP DIRECT SYSTEM:            |
   | 1. Primary Loop at high pressure      |     | 1. Water boils directly inside the    |
   |    (~15.5 MPa / 2250 psia) prevents   |     |    reactor core at ~7.0 MPa           |
   |    bulk boiling in reactor core.      |     |    (1020 psia).                       |
   | 2. Primary coolant transfers heat to  |     | 2. Steam separators send steam        |
   |    secondary loop in Steam Generator  |     |    directly to the turbine.           |
   |    U-tube heat exchangers.            |     | 3. High cycle efficiency, but turbine |
   | 3. Turbine steam is non-radioactive.  |     |    requires radiation shielding.      |
   +---------------------------------------+     +---------------------------------------+

5. Renewable Energy: Wind Power & Solar Systems

Wind Power Mechanics & Betz Limit

The kinetic power available in an upstream wind column of cross-sectional area $A = \frac{\pi}{4} D^2$ is:

Pwind=12m˙V2=12(ρAV)V2=12ρAV3P_{wind} = \frac{1}{2} \dot{m} V^2 = \frac{1}{2} (\rho A V) V^2 = \frac{1}{2} \rho A V^3

By momentum conservation across an actuator disk, the theoretical upper limit of aerodynamic power extraction is established by the Betz Limit:

Cp,max=16270.5926C_{p,max} = \frac{16}{27} \approx 0.5926

Actual Electrical Power Output Pelec=12ρAV3Cpηmechηgen\text{Actual Electrical Power Output } P_{elec} = \frac{1}{2} \rho A V^3 C_p \eta_{mech} \eta_{gen}

Where $C_p \approx 0.40 - 0.48$ for modern high-efficiency three-bladed horizontal axis wind turbines (HAWT).

+-----------------------------------------------------------------------------------------+
|                            WIND TURBINE OPERATING REGIMES                               |
|                                                                                         |
|   Cut-In Wind Speed (V_cut-in ~ 3-4 m/s):   Rotor starts generating power               |
|   Rated Wind Speed  (V_rated ~ 11-14 m/s):  Turbine reaches maximum nameplate capacity  |
|   Cut-Out Wind Speed (V_cut-out ~ 25 m/s):  Blades pitch to feather (shut down) for     |
|                                             structural protection during storms         |
+-----------------------------------------------------------------------------------------+

Solar Photovoltaic (PV) & Concentrating Solar Thermal (CSP)

  1. Photovoltaic (PV): Converts photon energy directly into DC electricity via the photoelectric effect in semiconductor $p-n$ junctions. System performance is rated at Standard Test Conditions (STC: $1000\text{ W/m}^2$ irradiance, $25^\circ\text{C}$ cell temperature, AM 1.5 spectrum). Inverter efficiency and temperature derating ($-0.35% \text{ to } -0.45% / ^\circ\text{C}$ above $25^\circ\text{C}$) are critical design factors.
  2. Concentrated Solar Power (CSP): Parabolic troughs or central power towers focus solar radiation onto receivers heating thermal oils or molten salts ($60% \text{ NaNO}_3 / 40% \text{ KNO}_3$ to $565^\circ\text{C}$), generating steam to drive a Rankine power block.

6. Grid-Scale Energy Storage & Plant Capacity Factor

Grid energy storage decouples energy generation from real-time customer demand:

+-----------------------------------------------------------------------------------------+
|                           GRID-SCALE ENERGY STORAGE COMPARISON                          |
|                                                                                         |
|   TECHNOLOGY             ROUND-TRIP EFFICIENCY   DISCHARGE DURATION   RESPONSE TIME     |
|   --------------------   ---------------------   ------------------   -------------     |
|   Pumped Hydro (PHS)     75% - 82%               6 - 24+ hours        Minutes           |
|   Compressed Air (CAES)  50% - 70%               4 - 12 hours         Minutes           |
|   Battery (BESS Li-ion)  85% - 92%               1 - 4 hours          Milliseconds      |
|   Thermal (Molten Salt)  85% - 90% (thermal)     4 - 16 hours         Minutes           |
+-----------------------------------------------------------------------------------------+

Capacity Factor Calculation

The Capacity Factor ($CF$) measures the actual electrical energy produced by a power facility over a given period (typically 1 calendar year = 8,760 hours) relative to its theoretical maximum continuous nameplate output:

Capacity Factor CF=Actual Net Energy Generated (MWh)Nameplate Rated Capacity (MW)×8760 hours/year×100\text{Capacity Factor } CF = \frac{\text{Actual Net Energy Generated (MWh)}}{\text{Nameplate Rated Capacity (MW)} \times 8760\text{ hours/year}} \times 100

  • Nuclear Power: $CF \approx 90-95%$ (continuous baseload)
  • Combined Cycle Natural Gas: $CF \approx 50-80%$ (intermediate/baseload)
  • Wind Power: $CF \approx 30-45%$ (intermittent resource)
  • Solar PV: $CF \approx 18-28%$ (diurnal resource)

7. Step-by-Step Worked Problem: Hydroelectric & Wind Power Calculations

Problem:

  1. Hydro Power: A hydroelectric facility operates with an effective net head of $H = 80\text{ m}$ ($262.5\text{ ft}$) and a volumetric flow rate of $Q = 25\text{ m}^3/\text{s}$. The combined turbine-generator efficiency is $\eta = 88%$. Calculate the electrical power generation in $\text{MW}$.
  2. Wind Power: A utility wind turbine with a rotor diameter of $D = 90\text{ m}$ operates in atmospheric air with density $\rho = 1.225\text{ kg/m}^3$. At a wind speed of $V = 12\text{ m/s}$, the aerodynamic power coefficient is $C_p = 0.44$ and the drivetrain/generator efficiency is $\eta_{mech-gen} = 92%$. Calculate the electrical power output in $\text{MW}$.
+-----------------------------------------------------------------------------------------+
|                        RENEWABLE POWER CALCULATION STEPS                                |
|                                                                                         |
|   PART 1: HYDROELECTRIC POWER GENERATION                                                |
|           Formula: P_hydro = eta * rho * g * Q * H                                      |
|           Given: eta = 0.88, rho = 1000 kg/m^3, g = 9.81 m/s^2, Q = 25 m^3/s, H = 80 m   |
|           P_hydro = 0.88 * (1000 kg/m^3) * (9.81 m/s^2) * (25 m^3/s) * (80 m)           |
|           P_hydro = 0.88 * 19,620,000 W = 17,265,600 W = 17.27 MW                       |
|                                                                                         |
|   PART 2: WIND TURBINE ELECTRICAL OUTPUT                                                |
|           Rotor Swept Area:                                                             |
|           A = (pi / 4) * D^2 = (pi / 4) * (90 m)^2 = 6,361.73 m^2                       |
|                                                                                         |
|           Available Wind Kinetic Power:                                                 |
|           P_wind = 0.5 * rho * A * V^3                                                  |
|           P_wind = 0.5 * (1.225 kg/m^3) * (6,361.73 m^2) * (12 m/s)^3                   |
|           P_wind = 0.5 * 1.225 * 6361.73 * 1728 = 6,733,255 W = 6.733 MW                |
|                                                                                         |
|           Actual Electrical Output:                                                     |
|           P_elec = P_wind * C_p * eta_mech-gen                                          |
|           P_elec = 6.733 MW * 0.44 * 0.92 = 2.726 MW                                    |
+-----------------------------------------------------------------------------------------+

8. Common Exam Traps & PE Pro-Tips

  • Trap 1 — Velocity Triangles Whirl Direction: In Euler's equation $w = U_1 V_{u1} - U_2 V_{u2}$, the whirl velocity $V_u$ is the tangential component of absolute velocity in the direction of blade motion $U$. If exit swirl opposes rotor rotation, $V_{u2}$ is negative, making $- U_2 V_{u2}$ additive.
  • Trap 2 — Hydraulic Head Units: In hydroelectric formulas $\dot{W} = \rho g Q H$, ensure consistent SI units ($Q$ in $\text{m}^3/\text{s}$, $H$ in $\text{m}$, $\rho$ in $\text{kg/m}^3$, $g = 9.81\text{ m/s}^2$ yields Watts) or US Customary units ($P = \gamma Q H / 550$ where $\gamma = 62.4\text{ lbf/ft}^3$, $Q$ in $\text{cfs}$, $H$ in $\text{ft}$ yields horsepower).
  • Trap 3 — Capacity Factor Calendar Hours: Annual capacity factor calculations must divide by $8,760\text{ hours}$ ($24 \times 365$), not $8,000$ or operating hours. The denominator represents theoretical continuous 24/7/365 production.
Test Your Knowledge

Compressed air ($k = 1.40$) at stagnation pressure $P_0 = 600\text{ kPa}$ feeds a converging nozzle discharging into a vessel. Below approximately what vessel back pressure will the nozzle be choked?

A
B
C
D
Test Your Knowledge

A proposed hydroelectric power project has an effective head of 450 meters and a design flow rate of 4.0 m³/s. Which type of hydraulic turbine is most appropriate for this installation?

A
B
C
D
Test Your Knowledge

In a 50% reaction steam turbine stage (Parsons stage), which of the following statements accurately characterizes the thermodynamic state and flow behavior across the moving rotor blades?

A
B
C
D
Test Your Knowledge

A 100 MW rated nameplate wind farm generates a total of 315,360 MWh of net electrical energy during a full non-leap calendar year (8,760 hours). What is the annual capacity factor of this wind farm?

A
B
C
D
Test Your Knowledge

What fundamental physical limitation is defined by the Betz Limit for open-rotor wind turbines?

A
B
C
D