Thermodynamics and Heat Transfer Explained: Energy Conversion, Transport Phenomena & Thermal Systems
- Mastering the strict governing distinction between equilibrium states (Thermodynamics) and non-equilibrium transport rates (Heat Transfer).
- Formulating closed/open control volume exergy destruction balances alongside 3D multi-mode transport (Fourier, Navier-Stokes convection, Stefan-Boltzmann).
- Bridging classical logarithmic mean temperature difference (LMTD) & ε-NTU methods with Physics-Informed Neural Networks (PINNs) and multi-scale thermal digital twins.
Thermodynamics and heat transfer form the foundational energy language of mechanical and thermal systems engineering. Classical thermodynamics defines the strict boundaries of what nature permits: how internal energy is transformed, how open and closed systems exchange work and heat, and how the inexorable generation of entropy limits maximum thermal conversion efficiency. Heat transfer provides the complementary temporal and spatial mechanisms: establishing the quantitative rates at which thermal energy moves across finite temperature gradients via conduction, convection, and electromagnetic radiation. Together, these disciplines enable engineers to model, simulate, and optimize power generation plants, aerospace thermal protection shields, microelectronics cooling architectures, and advanced electric vehicle battery management packs.
IDEF0 System Functional Decomposition
The IDEF0 (Integration Definition for Function Modeling) diagram formalized above abstracts thermal engineering pedagogy into a rigorous transformational control volume:
- Inputs (I): Unprocessed thermodynamic principles (energy conservation, state postulates, entropy) and raw heat transport physical properties (thermal conductivity k, fluid dynamic viscosity μ, surface emissivity ε, and specific heat capacities cp, cv).
- Controls (C): Curricular accreditation benchmarks (Washington Accord, ABET Criteria 3, EUR-ACE framework standards), dimensional homogeneity constraints, second-law irreversibility bounds, and environmental safety regulations.
- Mechanisms (M): High-performance computing clusters (CFD / OpenFOAM), calibrated subsonic wind tunnels, laser flash thermal diffusivity analyzers, and Physics-Informed Neural Network (PINN) inference engines.
- Outputs (O): Certified mechanical engineers equipped with advanced diagnostic judgment, verified thermal management solutions, and thermodynamic cycles optimized for minimum exergy destruction.
Mechanical & Materials Engineering Curriculum Cluster
Select any foundational pillar below to explore related modules across the PREP4UNI.ONLINE Mechanical Engineering series:
Mechanical Engineering
Sets the First and Second Law energy boundaries that dictate work conversion limits in mechanical power machinery.
Automotive Engineering
Sizes bottom cold plate heat exchangers to keep EV battery cells within 20–35°C and suppresses thermal propagation.
Biomechanical Engineering
Solves Pennes' bioheat transfer equation to model targeted freezing in cryosurgical probes and whole-body metabolic heat dissipation.
Control Systems in Mechanical Engineering
Tunes PID feedback controllers for chiller expansion valves, furnace fuel injection, and air handling mixing boxes.
Fluid Mechanics and Hydraulics
Establishes boundary layer velocity profiles that determine convective heat transfer coefficients (h) and pressure drop penalties.
Manufacturing and Production Engineering
Tracks latent heat dissipation during metal casting solidification and thermal cooling curves during metal heat treatment.
Mechanical Design and CAD
Transforms calculated heat transfer areas (A) into compact plate-fin cores, shell-and-tube baffles, and heat sink fins in CAD.
Nanotechnology and Advanced Materials
Applies low-diffusivity ceramic thermal barrier coatings (YSZ) to turbine blades and engineers high-conductivity carbon spreaders.
Robotics and Automation in Mechanical Engineering
Removes Joule heat generated in high-duty robot joint motors and safeguards electronics operating in hot foundry environments.
Solid Mechanics: Stress, Strain & Structural Analysis
Evaluates thermal expansion stresses (σ = EαΔT) in constrained structures and calculates high-temperature material creep rupture.
Thermodynamics and Heat Transfer
State postulates, exergy destruction, multi-mode Fourier/Newton/Stefan transport, LMTD/NTU methods, and PINN thermal models.
Vibrations and Acoustics
Analyzes heat-driven acoustic pressure pulsations (Rayleigh criterion) in combustion chambers and thermoacoustic refrigerators.
1. Deep Technical Foundations: Thermodynamic Systems & Governing Laws
Thermodynamic state analysis requires establishing a well-defined control surface dividing the universe into system and surroundings. The state of a simple compressible system is completely defined by two independent intensive properties according to the state postulate.
1.1 The First Law of Thermodynamics: Conservation of Energy
For a general open control volume with mass inflows and outflows, kinetic, potential, and flow work terms, the transient conservation of energy is formulated as:
$$ \frac{dE_{\text{cv}}}{dt} = \dot{Q} - \dot{W}_{\text{shaft}} + \sum_{\text{in}} \dot{m}_i \left( h_i + \frac{V_i^2}{2} + g z_i \right) - \sum_{\text{out}} \dot{m}_e \left( h_e + \frac{V_e^2}{2} + g z_e \right) $$where Ecv is the total stored internal, kinetic, and potential energy within the control volume, Q̇ is the net heat transfer rate across the boundary, Ẇshaft is the shaft power delivered to the environment, ṁ is mass flow rate, and h represents specific enthalpy defined by:
$$ h = u + P v $$For an ideal gas with constant specific heats, enthalpy and internal energy depend strictly on absolute temperature: dh = cpdT and du = cvdT, related via the specific gas constant R = cp − cv.
1.2 The Second Law of Thermodynamics, Entropy & Exergy Destruction
The Second Law dictates that the total entropy of an isolated system must always increase over time, reaching a maximum at thermodynamic equilibrium. The differential entropy generation rate Ṡgen is defined by the Clausius inequality extended to control volumes:
$$ \frac{dS_{\text{cv}}}{dt} = \sum_j \frac{\dot{Q}_j}{T_j} + \sum_{\text{in}} \dot{m}_i s_i - \sum_{\text{out}} \dot{m}_e s_e + \dot{S}_{\text{gen}} $$where Ṡgen ≥ 0. The equality holds strictly for reversible processes; all real physical transformations incur irreversibilities (fluid friction, unrestrained expansion, chemical reaction, finite-gradient heat transfer), enforcing Ṡgen > 0.
To evaluate the true thermodynamic utility of energy, engineers perform exergy analysis. Exergy (available work potential relative to an environmental dead state at T0, P0) is destroyed directly proportionally to entropy generation through the Gouy-Stodola theorem:
$$ \dot{X}_{\text{destroyed}} = T_0 \dot{S}_{\text{gen}} $$1.3 Ideal Power & Refrigeration Cycles
The maximum thermal efficiency of any heat engine operating between a high-temperature thermal reservoir at TH and a low-temperature sink at TL is bounded by the reversible Carnot efficiency:
$$ \eta_{\text{th, Carnot}} = 1 - \frac{T_L}{T_H} $$In industrial implementations, vapor and gas power cycles deviate from Carnot due to pump irreversibilities, finite transfer coefficients, and moisture content constraints:
- Rankine Cycle (Steam Turbines): Employs isentropic compression in a liquid pump, constant-pressure heat addition in a boiler, isentropic expansion in a steam turbine, and constant-pressure condensation. Superheating and reheating stages are deployed to prevent turbine blade moisture erosion.
- Brayton Cycle (Gas Turbines / Jet Propulsion): Operates via adiabatic compression, isobaric combustion, adiabatic gas expansion, and exhaust. Thermal efficiency is driven by the pressure ratio rp = P2 / P1: $$ \eta_{\text{Brayton}} = 1 - \frac{1}{r_p^{(\gamma - 1)/\gamma}} $$ where γ = cp / cv is the specific heat ratio.
- Vapor-Compression Refrigeration Cycle: Utilizes an evaporator, compressor, condenser, and isenthalpic expansion valve. The Coefficient of Performance (COP) is governed by: $$ \text{COP}_{\text{R}} = \frac{q_L}{w_{\text{in}}} = \frac{h_1 - h_4}{h_2 - h_1} $$
2. Heat Transfer Phenomena: Conduction, Convection & Radiation
2.1 Conduction Heat Transfer: Fourier's Law & The Heat Diffusion Equation
Conduction occurs via lattice vibrational waves (phonons) and free electron migration. The local vector heat flux q̋'' is proportional to the negative spatial temperature gradient:
$$ \mathbf{q}'' = -k \nabla T = -k \left( \frac{\partial T}{\partial x}\hat{\mathbf{i}} + \frac{\partial T}{\partial y}\hat{\mathbf{j}} + \frac{\partial T}{\partial z}\hat{\mathbf{k}} \right) $$Applying differential energy conservation across an infinitesimal control volume yields the 3D transient heat diffusion equation with volumetric thermal generation q̇:
$$ \rho c_p \frac{\partial T}{\partial t} = \frac{\partial}{\partial x}\left(k \frac{\partial T}{\partial x}\right) + \frac{\partial}{\partial y}\left(k \frac{\partial T}{\partial y}\right) + \frac{\partial}{\partial z}\left(k \frac{\partial T}{\partial z}\right) + \dot{q} $$Assuming homogeneous and isotropic thermal conductivity, this simplifies to:
$$ \frac{\partial T}{\partial t} = \alpha \nabla^2 T + \frac{\dot{q}}{\rho c_p} $$where α = k / (ρcp) is the thermal diffusivity (m2/s), measuring the rate at which thermal disturbances propagate relative to thermal energy storage capacity.
2.2 Convective Heat Transfer & Boundary Layer Dynamics
Convection combines molecular conduction at the solid boundary with bulk fluid advection. Local heat flux is governed by Newton's law of cooling:
$$ q'' = h (T_s - T_\infty) $$The heat transfer coefficient h is not an intrinsic fluid property; it is an empirical and analytical function of velocity profile, boundary layer development, turbulence, surface morphology, and fluid thermo-physical properties.
The boundary layer equations couple the continuity, momentum (Navier-Stokes), and thermal energy equations. Across a flat plate, the hydrodynamic boundary layer thickness δ(x) and thermal boundary layer thickness δt(x) scale according to the dimensionless Prandtl number (Pr):
$$ \text{Pr} = \frac{\nu}{\alpha} = \frac{\mu c_p}{k} $$$$ \frac{\delta(x)}{\delta_t(x)} \approx \text{Pr}^{1/3} $$Convective behavior is governed by key dimensionless similarity groups:
- Reynolds Number: ReL = (ρ U∞ L) / μ = (U∞ L) / ν (ratio of inertial to viscous forces).
- Nusselt Number: NuL = (h L) / kf (ratio of convective to pure molecular conductive heat transfer across boundary thickness L).
- Grashof Number: GrL = [g β (Ts − T∞) L3] / ν2 (ratio of buoyancy forces to viscous forces in natural convection, where β is the volumetric thermal expansion coefficient).
- Rayleigh Number: RaL = GrL · Pr (primary governing stability parameter for natural convection plume development).
2.3 Thermal Radiation: Electromagnetic Energy Transport
Thermal radiation is electromagnetic emission spanning wavelengths from λ ≈ 0.1 μm to 100 μm. An ideal blackbody emitter radiates energy according to the Stefan-Boltzmann law:
$$ E_b = \sigma T^4 $$where σ = 5.670374 × 10−8 W/(m2·K4). For a diffuse, gray real surface with emissivity ε (0 ≤ ε ≤ 1), the emissive power is:
$$ E = \epsilon \sigma T^4 $$Between two gray, diffuse, opaque surfaces forming an enclosure, the net radiation heat exchange is evaluated using the surface and space resistance network formulation:
$$ \dot{Q}_{12} = \frac{\sigma (T_1^4 - T_2^4)}{\frac{1 - \epsilon_1}{\epsilon_1 A_1} + \frac{1}{A_1 F_{12}} + \frac{1 - \epsilon_2}{\epsilon_2 A_2}} $$where F12 is the geometric view factor (fraction of radiation leaving surface 1 intercepted directly by surface 2), satisfying the reciprocity theorem A1F12 = A2F21.
3. Contemporary Computational Paradigms in Thermal Engineering
While classical 1D analytical models provide essential intuition, contemporary aerospace, nuclear, and automotive designs rely on advanced computational methods to solve complex multi-phase, conjugate heat transfer (CHT) problems.
Physics-Informed Neural Networks (PINNs) in Thermal Field Inversion
Traditional Finite Volume Methods (FVM) require complex computational grids and high processing budgets when resolving multi-scale boundary layers. Physics-Informed Neural Networks (PINNs) circumvent boundary meshing by embedding the non-linear heat diffusion and Navier-Stokes equations directly into the loss function of deep neural architectures:
$$ \mathcal{L}_{\text{total}}(\theta) = \mathcal{L}_{\text{data}}(\theta) + \lambda_{\text{PDE}} \mathcal{L}_{\text{PDE}}(\theta) + \lambda_{\text{BC}} \mathcal{L}_{\text{BC}}(\theta) $$$$ \mathcal{L}_{\text{PDE}}(\theta) = \frac{1}{N_f} \sum_{i=1}^{N_f} \left\| \frac{\partial \hat{T}}{\partial t} - \alpha \nabla^2 \hat{T} - \frac{\dot{q}}{\rho c_p} \right\|^2_{(x_i, y_i, z_i, t_i)} $$By training the network weights θ via automatic differentiation, PINNs enable real-time temperature field reconstruction from sparse point-sensor arrays (such as fiber Bragg grating sensors in battery modules), solving inverse heat transfer problems where boundary heat fluxes are unknown.
In high-throughput thermal design optimization, Fourier Neural Operators (FNOs) learn mappings between function spaces (such as arbitrary heat generation topologies and continuous steady-state temperature profiles) up to three orders of magnitude faster than conventional numerical solvers, making them ideal for generative heat exchanger topological optimization.
4. Systems Engineering Trade-Off Matrix: Thermal Management Architectures
Thermal systems design requires continuous optimization between thermal performance, pumping parasitics, structural mass, operational envelope, and reliability failure modes:
| Cooling Strategy | Primary Advantage | Bottlenecks & Failure Modes | Mitigation Strategy |
|---|---|---|---|
| Forced Air & Extruded Heat Sinks | Zero liquid leakage risk; lowest system integration capital cost; simple mechanical controls. | Limited to heat fluxes < 50 W/cm2; acoustic fan noise; dust/particulate fouling choking fin passages. | Integrate vapor chamber heat spreaders; aerodynamic fin profiling (serrated/pin fins); electrostatic filtration. |
| Single-Phase Liquid Cold Plates | High convective transfer coefficients (h > 5,000 W/m2·K); accommodates fluxes up to 300 W/cm2. | Pumping power penalties (ΔP); galvanic corrosion at fittings; seal degradation causing fluid leaks. | Ethylene glycol/water (EGW) with corrosion inhibitors; micro-channel manifold routing; burst-resistant brazing. |
| Two-Phase Evaporative / Heat Pipes | Nearly isothermal heat transport utilizing latent heat of vaporization (Δhvap); passive capillary operation. | Capillary pumping limits; boiling dry-out limits; non-condensable gas (NCG) generation degrading condenser zone. | Sintered copper powder wicks; strict chemical passivation procedures; working fluid degassing during vacuum bakeout. |
| Direct Immersion Liquid Cooling | Eliminates thermal interface materials (TIM); uniform temperature profile across dense electronic modules. | Expensive dielectric fluids (fluorochemicals/synthetic hydrocarbons); material chemical compatibility leaching. | Material compatibility screening (polymers, elastomers); sealed modular tank cassettes; active secondary heat rejection loops. |
5. Industrial Systems Case Applications
Aerospace: Gas Turbine Film Cooling
Modern turbine inlet temperatures (Tinlet > 1,600°C) significantly exceed the melting threshold of nickel-based superalloys (Tmelt ≈ 1,350°C). High-pressure compressor bleed air is routed through internal serpentine passages with rib turbulators, then ejected through angled shaped surface holes to create a protective, insulating thin film of cooler air across the blade aerodynamic profile.
Automotive: EV Battery Pack Thermal Runaway
Lithium-ion pouch and cylindrical cells generate internal Joulean (I2Rinternal) and entropic/reversible electrochemical reaction heat. Thermal management requires keeping the pack within 20–35°C to prevent lithium plating and premature capacity fade, while incorporating bottom cold plates and phase-change aerogels to arrest cell-to-cell thermal runaway propagation during catastrophic internal shorting.
Semiconductors: Microchannel Cold Plates
High-performance compute clusters and AI GPUs generate localized heat fluxes exceeding 150 W/cm2. Microchannel heat sinks with hydraulic diameters Dh < 200 μm exploit thin boundary layers where h ∝ k / Dh. Maintaining laminar flow prevents excessive pumping pressure drop penalties while achieving thermal resistances below 0.05 K/W.
6. Interactive Simulation Engine: Counter-Flow Heat Exchanger Modeling
Counter-Flow Heat Exchanger Numerical Simulator
7. Frequently Asked Questions
What fundamental boundary separates thermodynamics from heat transfer?
Why is exergy destruction preferable to energy balance alone in optimization?
What is the physical significance of the Biot number (Bi) in transient conduction?
How does the Prandtl number govern the relationship between velocity and thermal boundary layers?
8. Comprehensive Student Practice Modules
Section 1: Interactive Knowledge Self-Assessment
1. Why is enthalpy (h) rather than internal energy (u) the primary state variable used in open steady-flow devices?
2. What physical mechanism prevents a real thermal power plant from operating at the Carnot cycle efficiency?
3. What is the fundamental difference between thermal conductivity and thermal diffusivity?
4. Why do counter-flow heat exchangers achieve higher thermal effectiveness than parallel-flow configurations?
5. Under what condition does adding external insulation to a circular pipe increase heat loss?
6. Why is radiation often neglected in low-temperature electronics cooling, yet dominant in spacecraft thermal design?
7. What is the boundary layer separation mechanism and how does it influence convection?
8. How does nucleate boiling achieve heat transfer coefficients an order of magnitude higher than single-phase liquid convection?
Section 2: Analytical & Scenario-Based Case Studies
Case 1: Spacecraft Re-entry Thermal Protection System (TPS) Degradation
Scenario: During atmospheric re-entry at Mach 22, an orbital spacecraft incurs severe hypersonic shock-layer stagnation temperatures exceeding 8,000 K. The leading edge relies on reinforced carbon-carbon (RCC) and sintered silica ceramic tiles. Post-flight inspection reveals surface micro-spallation and localized sintering.
Engineering Analysis: At hypersonic stagnation, air dissociates into a reactive plasma. Heat transfer transitions from pure molecular conduction to catalytic surface recombination of nascent oxygen and nitrogen atoms. The failure mode stems from non-catalytic coating erosion, which spiked exothermic catalytic surface recombination heat flux. Mitigation involves applying high-emissivity silicon-hexaboride anti-oxidation coatings that simultaneously minimize surface catalyticity and maximize radiative surface reradiation: q̇rad = εσTs4.
Case 2: EV Fast-Charging Station Liquid Cold Plate Thermal Choking
Scenario: A 350 kW DC fast charger utilizes a water-glycol cooled copper cold plate to regulate internal silicon carbide (SiC) MOSFET inverter modules. During peak summer ambient operations (42°C), junction temperatures trip the emergency 175°C threshold despite maximum coolant flow rate.
Engineering Analysis: Pumping coolant at excessive velocities can induce laminar-turbulent transition delays in narrow microchannels, while driving hydraulic friction pressure losses beyond the pump's head capacity. Crucially, thermal interface material (TIM) pump-out occurred due to cyclic CTE mismatch expansion. The solution requires transitioning from a grease TIM to a phase-change metallic pad (indium alloy) and designing multi-pass manifold cold plates to cut channel length in half, reducing ΔP while elevating local Nusselt numbers.
Case 3: Data Center Closed-Loop Water-Cooled Rack Exergy Audit
Scenario: An enterprise compute facility consumes 10 MW of electrical power. The cooling loop supplies chilled water at 7°C and returns it at 13°C, requiring 2.5 MW of chiller compressor electrical power.
Engineering Analysis: Generating 7°C chilled water when ambient air is 20°C creates massive exergy destruction across the vapor-compression chillers. Raising the rack cooling supply water temperature to 35°C (allowable with direct-to-chip liquid cold plates) permits 100% dry-cooler "free cooling" without compressor engagement for 80% of the annual operating hours, slashing cooling overhead from 25% of facility energy to under 4%.
Case 4: Combined Cycle Gas Turbine (CCGT) Heat Recovery Steam Generator (HRSG) Pinch Point Analysis
Scenario: An engineer attempts to maximize steam production by decreasing the HRSG pinch-point temperature difference ΔTpinch between exhaust gas and evaporator saturation temperature from 15 K to 2 K.
Engineering Analysis: While reducing ΔTpinch lowers thermodynamic exergy destruction during heat transfer, the required heat transfer area scales asymptotically: A ∝ 1 / ΔTLMTD → ∞. Sizing for 2 K demands an economically non-viable capital expenditure increase in finned tubing, alongside severe exhaust-side gas pressure drops that back-pressure the upstream gas turbine, dropping its cycle work output.
Case 5: Cryogenic Liquid Hydrogen Fuel Tank Boil-off Mitigation
Scenario: A commercial aviation hydrogen storage vessel operating at 20 K (liquid phase at 2 bar) exhibits unacceptably high passive venting boil-off rates during tarmac turnaround.
Engineering Analysis: Heat leak is dominated by two vectors: residual gas conduction across vacuum jacket annular spaces due to vacuum degradation (loss of cryopumping), and thermal bridging across structural titanium load-bearing struts. Redesign mandates incorporating multilayer insulation (MLI) with aluminized Mylar radiation shields separated by low-conductivity silk netting, supplemented by vapor-cooled shielding (VCS) that routes escaping cold boil-off gas through intermediate structural supports to intercept parasitic heat before it reaches the core liquid.
Case 6: Nuclear Reactor Containment Core Catcher Passive Thermal Management
Scenario: In Generation III+ pressurized water reactors, a hypothetical core-melt accident results in liquid corium relocation into a sub-vessel core catcher. Heat must be continuously extracted without active AC electrical pumps.
Engineering Analysis: The containment relies on passive natural circulation driven by thermal siphon loops. The density differential between cold water downcomers and boiling chimney headers (Δρ = ρcold − ρhot) establishes a driving buoyancy head ΔPbuoyant = Δρ g H that balances two-phase frictional pressure drops, sustaining steady-state passive cooling heat flux margins safely below the Critical Heat Flux (CHF) limit.
Case 7: Supercritical CO2 (sCO2) Power Cycle Recuperator Transients
Scenario: In next-generation concentrated solar power sCO2 Brayton cycles, rapid solar cloud transients cause severe localized thermal stresses within printed circuit heat exchangers (PCHEs).
Engineering Analysis: Near the critical point (Tc = 31.0°C, Pc = 7.39 MPa), carbon dioxide exhibits dramatic property non-linearities: specific heat capacity cp spikes toward infinity. This local property fluctuation distorts the thermal capacitance ratio across the heat exchanger, inducing severe localized temperature gradients. Mitigating thermal fatigue requires utilizing bypass control valves and training PINN real-time predictive controllers to anticipate cloud-induced fluid property shifts and modulate mass flow rates proactively.
Section 3: Step-by-Step Quantitative Worked Examples
Worked Problem 1: Steady 1D Composite Wall Conduction & Interface Temperature Determination
An industrial furnace wall consists of an inner fireclay refractory layer (L1 = 0.20 m, k1 = 1.2 W/(m·K)) and an outer rockwool insulation layer (L2 = 0.10 m, k2 = 0.05 W/(m·K)). The inner gas temperature is T∞,1 = 1,100°C with internal convection coefficient h1 = 60 W/(m2·K). The ambient environment is at T∞,2 = 25°C with external convection coefficient h2 = 12 W/(m2·K). Determine the heat flux per unit area q'' and the critical interface temperature Tint between the refractory brick and insulation.
Internal convection: R''conv,1 = 1 / h1 = 1 / 60 = 0.01667 m2·K/W
Refractory layer: R''cond,1 = L1 / k1 = 0.20 / 1.2 = 0.16667 m2·K/W
Insulation layer: R''cond,2 = L2 / k2 = 0.10 / 0.05 = 2.00000 m2·K/W
External convection: R''conv,2 = 1 / h2 = 1 / 12 = 0.08333 m2·K/W
Step 2: Calculate total thermal resistance:
R''total = 0.01667 + 0.16667 + 2.00000 + 0.08333 = 2.26667 m2·K/W
Step 3: Determine heat flux:
ΔToverall = 1,100 − 25 = 1,075 K
q'' = ΔToverall / R''total = 1,075 / 2.26667 = 474.26 W/m2
Step 4: Compute interface temperature Tint:
Tint = T∞,1 − q'' · (R''conv,1 + R''cond,1)
R''upstream = 0.01667 + 0.16667 = 0.18334 m2·K/W
ΔTdrop = 474.26 × 0.18334 = 86.95 K
Tint = 1,100 − 86.95 = 1,013.05°C
Worked Problem 2: Open Brayton Gas Turbine Compressor Exergy Analysis
Air enters the adiabatic axial compressor of a gas turbine at T1 = 300 K, P1 = 100 kPa with mass flow rate ṁ = 25 kg/s. The compressor pressure ratio is rp = 16, and its isentropic efficiency is ηc = 0.86. Taking air properties cp = 1.005 kJ/(kg·K) and γ = 1.40, evaluate: (a) actual exit temperature T2, (b) required mechanical shaft power Ẇcomp, and (c) rate of exergy destruction Ẋdest assuming dead state T0 = 300 K, P0 = 100 kPa.
T2s = T1 · (P2 / P1)(γ − 1)/γ = 300 × (16)(0.4 / 1.4) = 300 × (16)0.2857 = 300 × 2.2085 = 662.56 K
Step 2: Determine actual discharge temperature T2:
ηc = (T2s − T1) / (T2 − T1) ⇒ (T2 − 300) = (662.56 − 300) / 0.86 = 362.56 / 0.86 = 421.58 K
T2 = 300 + 421.58 = 721.58 K (448.43°C)
Step 3: Calculate shaft power input:
Ẇcomp = ṁ · cp · (T2 − T1) = 25 × 1.005 × 421.58 = 10,592.20 kW ≈ 10.59 MW
Step 4: Compute entropy generation and exergy destruction:
s2 − s1 = cp ln(T2 / T1) − R ln(P2 / P1)
R = cp · (γ − 1) / γ = 1.005 × 0.4 / 1.4 = 0.2871 kJ/(kg·K)
cp ln(721.58 / 300) = 1.005 × ln(2.4053) = 1.005 × 0.87767 = 0.88206 kJ/(kg·K)
R ln(16) = 0.2871 × 2.77259 = 0.79601 kJ/(kg·K)
Δs = 0.88206 − 0.79601 = 0.08605 kJ/(kg·K)
Ṡgen = ṁ Δs = 25 × 0.08605 = 2.1513 kW/K
Ẋdest = T0 · Ṡgen = 300 × 2.1513 = 645.38 kW
Worked Problem 3: Forced Convection & Pressure Drop in a Rectangular Duct
Transformer oil at bulk mean temperature Tb = 40°C flows through a flat cooling channel of length L = 2.0 m, width W = 0.20 m, and height H = 0.02 m. The surface of the channel is maintained at uniform temperature Ts = 80°C. Fluid properties: ρ = 876 kg/m3, kf = 0.144 W/(m·K), μ = 0.035 Pa·s, cp = 1,970 J/(kg·K), Pr = 480. If mean velocity is Um = 0.75 m/s, evaluate the convective heat transfer coefficient h and the total thermal dissipation rate Q̇.
Ac = W · H = 0.20 × 0.02 = 0.004 m2
P = 2(W + H) = 2(0.20 + 0.02) = 0.44 m
Dh = 4Ac / P = (4 × 0.004) / 0.44 = 0.03636 m
ReDh = (ρ · Um · Dh) / μ = (876 × 0.75 × 0.03636) / 0.035 = 23.89 / 0.035 = 682.5
Step 2: Characterize flow regime & select Nusselt correlation:
Because ReDh = 682.5 < 2,300, the flow is strictly laminar.
Thermal entry length: Lt ≈ 0.05 · Re · Pr · Dh = 0.05 × 682.5 × 480 × 0.03636 = 595.6 m.
Since actual length L = 2.0 m ≪ Lt, the flow is thermally developing (Graetz problem). Sieder-Tate correlation for laminar developing entry:
NuL = 1.86 · [ (Re · Pr · Dh) / L ]1/3 · (μb / μs)0.14
Assuming property viscosity ratio term ≈ 1.15:
(Re · Pr · Dh) / L = (682.5 × 480 × 0.03636) / 2.0 = 11,911.6 / 2.0 = 5,955.8
(5,955.8)1/3 = 18.126
NuL = 1.86 × 18.126 × 1.15 = 38.77
Step 3: Calculate heat transfer coefficient:
h = (NuL · kf) / Dh = (38.77 × 0.144) / 0.03636 = 5.583 / 0.03636 = 153.55 W/(m2·K)
Step 4: Calculate total surface area & heat transfer rate:
As = P · L = 0.44 × 2.0 = 0.88 m2
ΔT = Ts − Tb = 80 − 40 = 40 K
Q̇ = h · As · ΔT = 153.55 × 0.88 × 40 = 5,404.96 W ≈ 5.40 kW
Worked Problem 4: Radiation Enclosure Exchange with a Re-radiating Surface
A vacuum furnace enclosure consists of two parallel circular coaxial discs of diameter D = 1.0 m separated by distance L = 0.50 m. Base 1 is heated to T1 = 1,000 K with emissivity ε1 = 0.80. Base 2 is maintained at T2 = 500 K with emissivity ε2 = 0.50. The cylindrical side refractory wall (Surface 3) is well-insulated, behaving as an ideal adiabatic reradiating surface. The direct view factor between the disks is F12 = 0.38. Calculate the net radiative heat exchange rate Q̇12 between the two plates.
A1 = A2 = (π / 4) · (1.0)2 = 0.7854 m2
Step 2: Determine view factors to reradiating wall:
By summation rule: F11 + F12 + F13 = 1.0. For flat disk, F11 = 0 ⇒ F13 = 1 − 0.38 = 0.62.
By symmetry: F23 = 0.62.
Step 3: Construct radiation resistance network:
Surface resistance 1: R1 = (1 − ε1) / (ε1 A1) = (1 − 0.8) / (0.8 × 0.7854) = 0.2 / 0.6283 = 0.3183 m−2
Surface resistance 2: R2 = (1 − ε2) / (ε2 A2) = (1 − 0.5) / (0.5 × 0.7854) = 0.5 / 0.3927 = 1.2732 m−2
Direct space resistance: R12 = 1 / (A1 F12) = 1 / (0.7854 × 0.38) = 1 / 0.29845 = 3.3506 m−2
Reradiating path resistance: Rrerad = [1 / (A1 F13)] + [1 / (A2 F23)] = 2 / (0.7854 × 0.62) = 2 / 0.4869 = 4.1072 m−2
Step 4: Equivalent space resistance:
Rspace,eq = (R12 · Rrerad) / (R12 + Rrerad) = (3.3506 × 4.1072) / (3.3506 + 4.1072) = 13.7615 / 7.4578 = 1.8453 m−2
Step 5: Total equivalent resistance & net heat rate:
Rtotal = R1 + Rspace,eq + R2 = 0.3183 + 1.8453 + 1.2732 = 3.4368 m−2
Blackbody emissive powers (σ = 5.67037 × 10−8 W/(m2·K4)):
Eb1 = 5.67037 × 10−8 × (1,000)4 = 56,703.7 W/m2
Eb2 = 5.67037 × 10−8 × (500)4 = 5.67037 × 10−8 × (6.25 × 1010) = 3,543.98 W/m2
ΔEb = 56,703.7 − 3,543.98 = 53,159.72 W/m2
Q̇12 = ΔEb / Rtotal = 53,159.72 / 3.4368 = 15,467.79 W ≈ 15.47 kW
Worked Problem 5: Transient Quenching & Lumped Capacitance Validity
A solid spherical ball bearing made of stainless steel (ρ = 7,850 kg/m3, cp = 480 J/(kg·K), k = 15.5 W/(m·K)) has diameter D = 12 mm. It is removed from an annealing furnace at T0 = 850°C and quenched into an agitated oil bath maintained at T∞ = 40°C with convection coefficient h = 450 W/(m2·K). Calculate: (a) the Biot number to confirm lumped capacitance validity, and (b) the time t required for the sphere center to reach 100°C.
For a sphere: Lc = V / As = [(π/6)D3] / [πD2] = D / 6
Lc = 0.012 / 6 = 0.002 m (2.0 mm)
Step 2: Evaluate the Biot Number (Bi):
Bi = (h · Lc) / k = (450 × 0.002) / 15.5 = 0.90 / 15.5 = 0.05806
Since Bi = 0.0581 < 0.10, internal temperature gradients are negligible (< 3%), confirming lumped capacitance validity.
Step 3: Calculate thermal time constant τt:
τt = (ρ · V · cp) / (h · As) = (ρ · cp · Lc) / h
τt = (7,850 × 480 × 0.002) / 450 = 7,536 / 450 = 16.747 s
Step 4: Solve for elapsed quench time t:
(T(t) − T∞) / (T0 − T∞) = exp(−t / τt)
(100 − 40) / (850 − 40) = 60 / 810 = 0.07407
ln(0.07407) = −2.6027
−t / 16.747 = −2.6027 ⇒ t = 2.6027 × 16.747 = 43.59 seconds
Worked Problem 6: Shell-and-Tube Heat Exchanger Sizing via LMTD and ε-NTU
A counter-flow shell-and-tube heat exchanger is designed to cool engine lube oil (cp = 2.10 kJ/(kg·K)) from Th,in = 100°C to Th,out = 60°C at a rate of ṁh = 2.0 kg/s. Cooling water (cp = 4.18 kJ/(kg·K)) enters at Tc,in = 20°C with mass flow rate ṁc = 1.5 kg/s. The overall heat transfer coefficient is estimated at U = 450 W/(m2·K). Calculate: (a) total thermal duty Q̇, (b) cooling water exit temperature Tc,out, (c) log-mean temperature difference ΔTLMTD, and (d) required surface area A.
Ch = ṁh · cp,h = 2.0 × 2.10 = 4.20 kW/K = 4,200 W/K
Cc = ṁc · cp,c = 1.5 × 4.18 = 6.27 kW/K = 6,270 W/K
Q̇ = Ch · (Th,in − Th,out) = 4,200 × (100 − 60) = 168,000 W = 168.0 kW
Step 2: Determine cooling water discharge temperature Tc,out:
Q̇ = Cc · (Tc,out − Tc,in) ⇒ 168,000 = 6,270 × (Tc,out − 20)
Tc,out − 20 = 168,000 / 6,270 = 26.80 K
Tc,out = 20 + 26.80 = 46.80°C
Step 3: Calculate Log-Mean Temperature Difference (ΔTLMTD):
For counter-flow:
ΔT1 = Th,in − Tc,out = 100 − 46.80 = 53.20 K
ΔT2 = Th,out − Tc,in = 60 − 20 = 40.00 K
ΔTLMTD = (ΔT1 − ΔT2) / ln(ΔT1 / ΔT2) = (53.20 − 40.00) / ln(53.20 / 40.00)
ΔTLMTD = 13.20 / ln(1.330) = 13.20 / 0.28518 = 46.29 K
Step 4: Compute required heat exchanger surface area A:
Q̇ = U · A · ΔTLMTD ⇒ A = Q̇ / (U · ΔTLMTD)
A = 168,000 / (450 × 46.29) = 168,000 / 20,830.5 = 8.065 m2