Analyze thermodynamic cycles and heat transfer in pump systems
This skill provides a systematic approach to thermodynamic analysis of pump systems, including energy balance, entropy analysis, heat transfer, and cycle calculations. The methodology follows classical thermodynamics principles applied specifically to pumping systems and fluid machinery.
Define the Control Volume:
The first step in any thermodynamic analysis is to clearly define the system boundaries.
For Pump Systems:
System Classification:
Assumptions to Document:
Coordinate System:
Thermodynamic State Definition:
For a pure substance, two independent intensive properties define the state.
State Point 1 (Inlet/Suction):
Identify and measure:
Properties at inlet:
- Pressure: P₁ (Pa, bar, psi)
- Temperature: T₁ (K, °C)
- Velocity: c₁ (m/s)
- Elevation: z₁ (m)
- Phase: Liquid, vapor, two-phase
State Point 2 (Outlet/Discharge):
Identify and measure:
Properties at outlet:
- Pressure: P₂ (Pa, bar, psi)
- Temperature: T₂ (K, °C)
- Velocity: c₂ (m/s)
- Elevation: z₂ (m)
- Phase: Liquid, vapor, two-phase
Additional State Points:
For complex systems:
Property Diagrams:
Plot state points on:
Fluid Property Sources:
For Pure Substances:
For Incompressible Liquids:
Specific volume: v ≈ constant
Specific heat: cp ≈ constant
Internal energy: du = cp·dT
Enthalpy: dh = cp·dT + v·dP
For Ideal Gases:
Equation of state: P·v = R·T
Specific heats: cp, cv = f(T)
Internal energy: du = cv·dT
Enthalpy: dh = cp·dT
Entropy: ds = cp·dT/T - R·dP/P
Property Evaluation Methods:
Method 1: Using CoolProp (Recommended)
from CoolProp.CoolProp import PropsSI
# Water at 25°C, 2 bar
h = PropsSI('H', 'T', 298.15, 'P', 200000, 'Water') # J/kg
s = PropsSI('S', 'T', 298.15, 'P', 200000, 'Water') # J/kg-K
rho = PropsSI('D', 'T', 298.15, 'P', 200000, 'Water') # kg/m³
Method 2: Property Tables
Method 3: Empirical Correlations
Critical Properties to Evaluate:
At each state point:
General Energy Equation for Open Systems:
For a control volume in steady flow:
Q̇ - Ẇ = ṁ[(h₂ - h₁) + (c₂² - c₁²)/2 + g(z₂ - z₁)]
Where:
For Pump Systems:
Work is done ON the fluid (Ẇ < 0 by convention), so:
Ẇshaft = ṁ[(h₂ - h₁) + (c₂² - c₁²)/2 + g(z₂ - z₁)] + Q̇loss
Simplified Forms:
Incompressible Liquid (Constant Density):
Ẇshaft/ṁ = (P₂ - P₁)/ρ + (c₂² - c₁²)/2 + g(z₂ - z₁) + q_loss
This can be written in terms of head:
H = (P₂ - P₁)/(ρg) + (c₂² - c₁²)/(2g) + (z₂ - z₁) + h_loss
Adiabatic Pump (No Heat Loss):
Ẇshaft = ṁ[(h₂ - h₁) + (c₂² - c₁²)/2 + g(z₂ - z₁)]
Negligible Kinetic and Potential Energy Changes:
Ẇshaft = ṁ(h₂ - h₁) + Q̇loss
Energy Balance Components:
Enthalpy Change:
Kinetic Energy Change:
Potential Energy Change:
Heat Transfer:
Temperature Rise in Pumped Fluid:
The temperature rise due to inefficiency:
ΔT = (1 - η)·Ẇshaft / (ṁ·cp)
Where η is the pump efficiency.
For incompressible liquids:
ΔT = (1 - η)·g·H / cp
Example: Water pump, H = 100 m, η = 0.80, cp = 4186 J/kg-K
ΔT = (1 - 0.80)·9.81·100 / 4186 = 0.47°C
Entropy Balance:
For a control volume in steady flow:
ṁ(s₂ - s₁) = Q̇/T_boundary + Ṡgen
Where:
Entropy Generation:
Measures irreversibility in the process:
Ṡgen = ṁ(s₂ - s₁) - Q̇/T_boundary
For adiabatic pump:
Ṡgen = ṁ(s₂ - s₁)
Isentropic Efficiency:
Compares actual process to ideal isentropic process:
η_isentropic = (h₂s - h₁) / (h₂ - h₁)
Where:
For Incompressible Liquids:
The isentropic work is:
w_isentropic = ∫v·dP = v·(P₂ - P₁) = (P₂ - P₁)/ρ
Exergy Analysis:
Exergy (availability) represents maximum useful work:
ψ = (h - h₀) - T₀(s - s₀) + c²/2 + gz
Where subscript 0 denotes dead state (ambient conditions).
Exergy destruction:
Ẋdestroyed = T₀·Ṡgen
Second Law Efficiency:
η_II = Ẋuseful / Ẋinput = 1 - Ẋdestroyed / Ẋinput
This represents how effectively the pump uses available energy.
Loss Breakdown by Entropy Generation:
Different loss mechanisms contribute to entropy generation:
Thermodynamic Cycles Relevant to Pumps:
While pumps are typically single components, they operate within larger thermodynamic cycles.
Components:
Pump Analysis in Rankine Cycle:
State points:
Pump work (isentropic):
w_pump,s = v₃·(P₄ - P₃) = h₄s - h₃
Actual pump work:
w_pump = w_pump,s / η_pump
Impact on Cycle Efficiency:
η_cycle = (w_turbine - w_pump) / q_in
Pump work is typically 1-3% of turbine work in Rankine cycles.
Not directly applicable to liquid pumps, but relevant for:
Key Differences from Rankine:
Relevant for:
Components:
Coefficient of Performance (COP):
COP = Q_useful / W_input
Where pump work contributes to W_input.
Physical Mechanism:
Inefficiencies in pumps convert mechanical energy to thermal energy, causing temperature rise.
Calculation Method:
From Energy Balance:
ΔT = (h₂ - h₁) / cp - v·(P₂ - P₁) / cp
Simplified for Liquids:
ΔT = (1 - η)·v·(P₂ - P₁) / cp
Or in terms of head:
ΔT = (1 - η)·g·H / cp
Factors Affecting Temperature Rise:
Practical Implications:
Example Calculations:
Water pump:
Hydraulic oil:
High viscosity fluid:
Overall Efficiency:
η = η_h · η_vol · η_mech
Component Efficiencies:
Hydraulic Efficiency (η_h):
Volumetric Efficiency (η_vol):
Mechanical Efficiency (η_mech):
Thermodynamic Loss Analysis:
Each loss mechanism increases entropy:
Friction Loss:
Ṡ_friction = ṁ·g·h_friction / T_avg
Leakage Loss:
Ṡ_leakage = ṁ_leak·(s_discharge - s_suction)
Heat Transfer Loss:
Ṡ_heat = Q̇_loss / T_boundary - Q̇_loss / T_fluid
Total Entropy Generation:
Ṡ_total = Ṡ_friction + Ṡ_leakage + Ṡ_heat + Ṡ_mechanical
Exergy Efficiency:
η_exergy = 1 - T₀·Ṡ_total / Ẇshaft
This gives a more complete picture of thermodynamic performance.
Heat Transfer Modes in Pumps:
Conduction:
Convection:
Radiation:
Heat Loss Estimation:
From Casing Walls:
Q̇_loss = h·A·(T_fluid - T_ambient)
Where:
From Energy Balance:
Q̇_loss = Ẇshaft - ṁ·(h₂ - h₁) - ṁ·(c₂² - c₁²)/2 - ṁ·g·(z₂ - z₁)
Impact on Performance:
Thermal Management:
Cooling Requirements:
Insulation:
Staging Rationale:
High head requirements divided among multiple impellers:
H_total = H_stage1 + H_stage2 + ... + H_stage_n
Thermodynamic Analysis:
Between Stages:
State 1 → Impeller 1 → State 2 → Diffuser 1 → State 3 →
Impeller 2 → State 4 → ... → Final discharge
Temperature Rise per Stage:
ΔT_stage = (1 - η_stage)·g·H_stage / cp
Total Temperature Rise:
ΔT_total = Σ ΔT_stage = (1 - η_overall)·g·H_total / cp
Interstage Cooling:
For high-head or high-temperature applications:
Q̇_cooling = ṁ·cp·(T_beforecooling - T_aftercooling)
Benefits:
Pressure Compounding Effects:
As pressure increases through stages:
Energy Distribution:
Not all stages contribute equally:
Steam Power Plant Configuration:
Boiler → Turbine → Condenser → Feedwater Pump → Boiler
↓ ↓ ↓ ↓
High P Power Low P High P
High T Output Low T Moderate T
State Point Analysis:
State 1: Turbine exhaust entering condenser
State 2: Condenser outlet / Pump inlet
State 3: Pump outlet / Boiler inlet
State 4: Boiler outlet / Turbine inlet
Pump Work Calculation:
Isentropic work:
w_pump,s = h₃s - h₂ = v₂·(P₃ - P₂)
For water at 40°C (v ≈ 0.001008 m³/kg):
w_pump,s = 0.001008·(250e5 - 0.08e5) ≈ 2.52 kJ/kg
Actual work:
w_pump = w_pump,s / η_pump = 2.52 / 0.85 ≈ 2.96 kJ/kg
Temperature rise:
ΔT = (h₃ - h₂) / cp ≈ 0.44 / 4.18 ≈ 0.11°C (negligible)
Cycle Efficiency:
η_Rankine = (w_turbine - w_pump) / q_in
Typical values:
Pump work is small (~1-2% of turbine work) but critical for operation.
Vapor Compression Cycle with Liquid Pump:
Some industrial refrigeration systems use liquid pumps instead of or in addition to compressors.
State Points:
Pump Application:
COP Calculation:
COP = Q_evaporator / (W_compressor + W_pump)
Pump work usually small compared to compressor work.
Applications:
Thermodynamic Analysis:
Heating Mode:
Q̇_delivered = ṁ·cp·(T_supply - T_return)
W_pump = ṁ·g·H / η_pump
COP_heating = Q̇_delivered / W_pump
Cooling Mode:
Q̇_removed = ṁ·cp·(T_return - T_supply)
COP_cooling = Q̇_removed / W_pump
Optimization:
Thermodynamic analysis of pump systems requires systematic application of:
Key outputs:
See equations-reference.md for detailed equations and examples.md for worked examples.