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.
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:
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
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$:
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:
| Condition | Consequence |
|---|---|
| $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.
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}$):
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}$):
+-----------------------------------------------------------------------------------------+
| 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:
By momentum conservation across an actuator disk, the theoretical upper limit of aerodynamic power extraction is established by the Betz Limit:
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)
- 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.
- 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:
- 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:
- 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}$.
- 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.
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 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?
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 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?
What fundamental physical limitation is defined by the Betz Limit for open-rotor wind turbines?