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

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:

$$\text{Specific Work } w = U_1 V_{u1} - U_2 V_{u2}$$ $$\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 = \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}\cdot \text{K)} = 53.35\text{ ft}\cdot \text{lbf/(lbm}\cdot \text{°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$:

$$\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:

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

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.

$$\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}$$ $$\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}$):

$$N_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}$):

$${}^{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:

$$P_{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:

$$C_{p,max} = \frac{16}{27} \approx 0.5926$$ $$\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:

$$\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.