3.1 Thermodynamics Fundamentals and Fluid Flow Principles

3.1 Thermodynamics Fundamentals and Fluid Flow Principles

Quick Answer: The First Law of Thermodynamics governs energy conservation across closed systems ($Q - W = \Delta U$) and steady open control volumes ($\dot{Q} - \dot{W}_s = \dot{m}[\Delta h + \Delta(v^2/2) + g\Delta z]$). The Second Law sets the upper bound on thermal efficiency via the Carnot limit ($\eta = 1 - T_L / T_H$, requiring absolute temperatures in Kelvin or Rankine). Fluid flow couples mass continuity ($Q = A_1 v_1 = A_2 v_2$) with the extended Bernoulli energy equation, where viscous pipe friction head loss is evaluated via the Darcy-Weisbach formula ($h_f = f \frac{L}{D} \frac{v^2}{2g}$), with $f = 64/Re$ in laminar flow ($Re < 2,100$).

Industrial engineers regularly analyze thermodynamic cycles in power generation, compressed air distribution, industrial drying, facility heating, ventilation, and air conditioning (HVAC), as well as fluid transport in chemical manufacturing and assembly plant piping. Mastery of foundational energy balances, gas behavior, heat transfer, and pipe hydraulics is directly assessed on the NCEES FE exam.

---

1. First and Second Laws of Thermodynamics

Thermodynamic systems are classified by their boundaries:

• Closed System (Control Mass): No mass crosses the boundary; energy enters or leaves as heat ($Q$) or work ($W$).
• Open System (Control Volume): Both mass and energy cross the control surface (e.g., pumps, turbines, compressors, nozzles).
• Isolated System: Neither mass nor energy crosses the system boundary.

The First Law of Thermodynamics (Conservation of Energy)

For a closed system undergoing a finite process between state 1 and state 2, the First Law is expressed as:

$$Q - W = \Delta U + \Delta KE + \Delta PE$$

Where:

• $Q$ = net heat transfer added to the system (positive if heat enters; negative if rejected)
• $W$ = net work done by the system on its surroundings (positive for expansion work; negative for work input)
• $\Delta U = m c_v (T_2 - T_1)$ = change in internal energy for an ideal gas with constant specific heats
• $\Delta KE = \frac{1}{2} m (v_2^2 - v_1^2)$ = change in kinetic energy
• $\Delta PE = m g (z_2 - z_1)$ = change in gravitational potential energy

In stationary industrial systems, kinetic and potential energy changes are typically negligible ($\Delta KE \approx 0$, $\Delta PE \approx 0$), simplifying the First Law to $Q - W = \Delta U$.

For an open system operating at steady state, steady flow (SSSF), mass conservation requires $\sum \dot{m}_{in} = \sum \dot{m}_{out} = \dot{m}$. The steady-flow energy equation on a rate basis is:

$$\dot{Q} - \dot{W}_s = \dot{m} \left[ \left(h_2 - h_1\right) + \frac{v_2^2 - v_1^2}{2} + g\left(z_2 - z_1\right) \right]$$

Where:

• $\dot{W}_s$ = shaft power delivered by the control volume (e.g., turbine output is positive; compressor/pump input is negative)
• $h = u + P v$ = specific enthalpy (kJ/kg or Btu/lbm)

The Second Law of Thermodynamics and Carnot Efficiency

The Second Law dictates the natural direction of physical processes:

• Kelvin-Planck Statement: It is impossible for any heat engine operating in a thermodynamic cycle to receive heat from a single thermal reservoir and deliver an equivalent amount of work without rejecting heat to a lower-temperature reservoir.
• Clausius Statement: It is impossible to construct a cyclic refrigerator or heat pump that transfers heat from a colder body to a hotter body without an external input of work.

The theoretical maximum thermal efficiency of any heat engine operating between a high-temperature thermal reservoir at $T_H$ and a low-temperature sink at $T_L$ is given by the Carnot efficiency:

$$\eta_{th,max} = \eta_{Carnot} = 1 - \frac{T_L}{T_H} = \frac{T_H - T_L}{T_H}$$

CRITICAL EXAM RULE: Reservoir temperatures $T_H$ and $T_L$ in the Carnot efficiency formula MUST ALWAYS be expressed in absolute temperature units: Kelvin ($T_K = T_{^\circ C} + 273.15$) for SI units, or Rankine ($T_R = T_{^\circ F} + 459.67$) for US Customary units. Substituting Celsius or Fahrenheit directly is the single most common distractor on the FE exam.

For refrigerators and heat pumps operating on reversed Carnot cycles, performance is rated by the Coefficient of Performance (COP):

$$COP_{refrig, Carnot} = \frac{Q_L}{W_{in}} = \frac{T_L}{T_H - T_L}$$ $$COP_{HP, Carnot} = \frac{Q_H}{W_{in}} = \frac{T_H}{T_H - T_L} = COP_{refrig, Carnot} + 1$$

---

2. State Variables and Pure Substance Properties

Thermodynamic state is defined by intensive properties (independent of system size: pressure $P$, temperature $T$, specific volume $v$, specific internal energy $u$, specific enthalpy $h$, specific entropy $s$) and extensive properties (dependent on total mass: volume $V$, internal energy $U$).

The Vapor Dome and Phase Quality

For pure substances such as water or refrigerants, phase change occurs under the saturation vapor dome:

• Subcooled (Compressed) Liquid: Liquid at a temperature below the saturation temperature for the given pressure ($T < T_{sat}(P)$).
• Saturated Liquid (subscript $f$): Liquid on the verge of boiling at $T = T_{sat}(P)$; quality $x = 0$.
• Saturated Vapor (subscript $g$): Vapor on the verge of condensation at $T = T_{sat}(P)$; quality $x = 1$.
• Two-Phase Liquid-Vapor Mixture: Liquid and vapor coexisting in equilibrium at $T_{sat}(P)$ and $P_{sat}(T)$.

The dryness fraction or vapor quality ($x$) is defined strictly within the two-phase mixture region as the ratio of vapor mass to total mixture mass:

$$x = \frac{m_{vapor}}{m_{total}} = \frac{m_g}{m_f + m_g}, \quad 0 \le x \le 1$$

Any intensive thermodynamic property $y$ (where $y$ can represent $v, u, h,$ or $s$) within the saturated liquid-vapor mixture region is computed using the saturation table values:

$$y = y_f + x y_{fg} = y_f + x (y_g - y_f)$$

Where $y_{fg} = y_g - y_f$ is the property difference between saturated vapor and saturated liquid (e.g., $h_{fg}$ is the latent heat of vaporization).

---

3. Ideal Gas Laws and Polytropic Processes

At temperatures sufficiently above the critical temperature and pressures well below the critical pressure, real gases obey the Ideal Gas Law:

$$P v = R T \quad \iff \quad P V = m R T = n \bar{R} T$$

Where:

• $P$ = absolute pressure (Pa or $\text{N/m}^2$, $\text{lbf/ft}^2$)
• $v = V/m$ = specific volume ($\text{m}^3\text{/kg}$, $\text{ft}^3\text{/lbm}$)
• $T$ = absolute temperature (K or $^\circ\text{R}$)
• $\bar{R} = 8.314\text{ kJ/(kmol}\cdot\text{K)} = 1,545\text{ ft}\cdot\text{lbf/(lbmol}\cdot^\circ\text{R)}$ = universal gas constant
• $R = \frac{\bar{R}}{M}$ = specific gas constant, where $M$ is molar mass (for air: $M \approx 28.97\text{ kg/kmol}$, yielding $R_{air} \approx 0.2870\text{ kJ/(kg}\cdot\text{K)} = 53.34\text{ ft}\cdot\text{lbf/(lbm}\cdot^\circ\text{R)}$)

Specific Heat Relationships

For an ideal gas, internal energy and enthalpy depend strictly on temperature:

$$du = c_v dT, \quad dh = c_p dT$$ $$c_p - c_v = R, \quad k = \frac{c_p}{c_v}$$ $$c_v = \frac{R}{k - 1}, \quad c_p = \frac{k R}{k - 1}$$

For standard atmospheric air at room temperature ($300\text{ K}$): $k \approx 1.4$, $c_p \approx 1.005\text{ kJ/(kg}\cdot\text{K)}$, and $c_v \approx 0.718\text{ kJ/(kg}\cdot\text{K)}$.

Polytropic Processes ($P v^n = \text{constant}$)

Many closed-system compression and expansion processes follow a polytropic path $P v^n = C$:

---

4. Heat Transfer Fundamentals

Thermal energy transfers across system boundaries via three distinct mechanisms:

Conduction (Fourier's Law)

Heat transfer through a stationary solid or fluid by microscopic molecular vibrations and electron motion:

$$\dot{Q}_{cond} = -k A \frac{dT}{dx} = \frac{k A (T_1 - T_2)}{L} = \frac{T_1 - T_2}{R_{cond}}$$

Where $k$ is thermal conductivity ($\text{W/(m}\cdot\text{K)}$ or $\text{Btu/(hr}\cdot\text{ft}\cdot^\circ\text{F)}$), $A$ is heat transfer surface area, $L$ is slab thickness, and $R_{cond} = \frac{L}{k A}$ is conductive thermal resistance ($\text{K/W}$). For composite industrial walls in series:

$$R_{total} = \sum R_i = \frac{L_1}{k_1 A} + \frac{L_2}{k_2 A} + \cdots + \frac{L_n}{k_n A}$$

Convection (Newton's Law of Cooling)

Heat transfer between a solid surface and an adjacent moving fluid:

$$\dot{Q}_{conv} = h A (T_s - T_\infty) = \frac{T_s - T_\infty}{R_{conv}}$$

Where $h$ is the convective heat transfer coefficient ($\text{W/(m}^2\cdot\text{K)}$ or $\text{Btu/(hr}\cdot\text{ft}^2\cdot^\circ\text{F)}$), $T_s$ is surface temperature, $T_\infty$ is bulk fluid temperature, and $R_{conv} = \frac{1}{h A}$.

Radiation (Stefan-Boltzmann Law)

Electromagnetic radiation emitted by all matter at non-zero absolute temperature:

$$\dot{Q}_{rad} = \epsilon \sigma A \left(T_s^4 - T_{surr}^4\right)$$

Where:

• $\epsilon$ = surface emissivity ($0 \le \epsilon \le 1$; $\epsilon = 1$ for an ideal blackbody)
• $\sigma = 5.670 \times 10^{-8}\text{ W/(m}^2\cdot\text{K}^4) = 0.1714 \times 10^{-8}\text{ Btu/(hr}\cdot\text{ft}^2\cdot^\circ\text{R}^4)$ = Stefan-Boltzmann constant
• Temperatures $T_s$ and $T_{surr}$ must be in Kelvin or Rankine.

---

5. Fluid Properties and Hydrostatics

Fundamental Fluid Properties

• Density ($\rho$): Mass per unit volume ($\text{kg/m}^3$ or $\text{slug/ft}^3$, $\text{lbm/ft}^3$). Water at $4^\circ\text{C}$: $\rho_{H_2O} = 1,000\text{ kg/m}^3 = 1.94\text{ slugs/ft}^3 = 62.4\text{ lbm/ft}^3$.
• Specific Weight ($\gamma$): Weight per unit volume: $\gamma = \rho g$. For water: $\gamma_{H_2O} = 9,810\text{ N/m}^3 = 62.4\text{ lbf/ft}^3$.
• Specific Gravity ($SG$): Dimensionless density ratio relative to water at $4^\circ\text{C}$: $SG = \frac{\rho}{\rho_{H_2O}} = \frac{\gamma}{\gamma_{H_2O}}$.
• Dynamic (Absolute) Viscosity ($\mu$): Fluid resistance to shear deformation ($\text{Pa}\cdot\text{s} = \text{N}\cdot\text{s/m}^2$ or $\text{lbf}\cdot\text{s/ft}^2$). Newtonian shear stress is $\tau = \mu \frac{du}{dy}$.
• Kinematic Viscosity ($\nu$): Ratio of dynamic viscosity to density: $\nu = \frac{\mu}{\rho}$ ($\text{m}^2\text{/s}$ or $\text{ft}^2\text{/s}$). 1 stoke = $1\text{ cm}^2\text{/s} = 10^{-4}\text{ m}^2\text{/s}$.

Hydrostatic Pressure Distribution

In a static fluid of constant density under gravity, pressure depends solely on vertical depth $h$:

$$P = P_0 + \rho g h = P_0 + \gamma h$$ $$\Delta P = \gamma \Delta h$$

• Absolute vs. Gage Pressure: $$P_{abs} = P_{gage} + P_{atm}$$ Where standard atmospheric pressure at sea level is $P_{atm} = 101.325\text{ kPa} = 1.01325\text{ bar} = 14.696\text{ psia}$. When $P_{abs} < P_{atm}$, vacuum pressure is $P_{vac} = P_{atm} - P_{abs}$.

Hydrostatic Force on Submerged Plane Surfaces

The total resultant hydrostatic force acting on a flat submerged surface of area $A$ is:

$$F_R = \bar{P} A = \left(P_0 + \gamma \bar{h}\right) A$$

Where $\bar{h}$ is the vertical depth from the free surface to the centroid of the area. The center of pressure ($y_{cp}$), where the resultant force acts along the inclined surface, lies strictly below the area centroid:

$$y_{cp} = \bar{y} + \frac{I_{xc}}{\bar{y} A}$$

Where $I_{xc}$ is the area moment of inertia about the horizontal centroidal axis (for a rectangle of width $b$ and height $h$: $I_{xc} = \frac{b h^3}{12}$).

---

6. Continuity and Bernoulli's Principles

The Continuity Equation (Conservation of Mass)

For one-dimensional steady flow through a conduit with cross-sectional areas $A_1$ and $A_2$:

$$\dot{m} = \rho_1 A_1 v_1 = \rho_2 A_2 v_2 = \text{constant}$$

For an incompressible fluid (constant density $\rho_1 = \rho_2$):

$$Q = A_1 v_1 = A_2 v_2 = \text{constant}$$

Where $Q$ is volumetric flow rate ($\text{m}^3\text{/s}$ or $\text{ft}^3\text{/s}$, gpm), and $v$ is average flow velocity. Because circular pipe area is $A = \frac{\pi D^2}{4}$, velocity scales inversely with the square of diameter:

$$v_2 = v_1 \left(\frac{D_1}{D_2}\right)^2$$

Bernoulli's Equation

Along a streamline in steady, incompressible, frictionless (inviscid) flow with no shaft work and no heat transfer:

$$\frac{P_1}{\gamma} + \frac{v_1^2}{2g} + z_1 = \frac{P_2}{\gamma} + \frac{v_2^2}{2g} + z_2 = \text{constant}$$

Each term possesses units of length (meters or feet), representing an energy head:

• Pressure Head: $\frac{P}{\gamma}$
• Velocity (Dynamic) Head: $\frac{v^2}{2g}$
• Elevation (Potential) Head: $z$

Graphically, the Hydraulic Grade Line (HGL) represents the sum of elevation and pressure heads ($HGL = z + \frac{P}{\gamma}$), while the Energy Grade Line (EGL) includes the velocity head ($EGL = HGL + \frac{v^2}{2g}$). In frictionless flow with no pumps or turbines, EGL is strictly horizontal.

---

7. Viscous Pipe Flow: Head Loss and Industrial Pumping

Real industrial fluids experience viscous shear friction against pipe walls, causing mechanical energy dissipation into thermal internal energy.

Reynolds Number and Flow Regimes

The dimensionless Reynolds Number ($Re$) characterizes the ratio of inertial forces to viscous forces:

$$Re = \frac{\rho v D}{\mu} = \frac{v D}{\nu}$$

Where $D$ is the internal pipe diameter. In internal circular pipe flow:

• Laminar Flow: $Re < 2,100$ (fluid moves in parallel, smooth concentric laminas without macroscopic mixing)
• Transitional Flow: $2,100 \le Re \le 4,000$
• Turbulent Flow: $Re > 4,000$ (chaotic, fluctuating eddy currents; enhanced wall shear)

Darcy-Weisbach Equation for Major Head Loss

The major friction head loss ($h_f$) in a straight pipe of length $L$ and diameter $D$ is calculated via:

$$h_f = f \frac{L}{D} \frac{v^2}{2g}$$

Where $f$ is the dimensionless Darcy friction factor:

• For Laminar Flow ($Re < 2,100$): Friction depends strictly on Reynolds number and is completely independent of pipe surface roughness: $$f = \frac{64}{Re}$$
• For Turbulent Flow ($Re > 4,000$): The friction factor depends on both the Reynolds number and relative roughness ($\epsilon / D$), obtained from the Moody diagram or the Swamee-Jain explicit approximation: $$f = \frac{0.25}{\left[\log_{10}\left(\frac{\epsilon/D}{3.7} + \frac{5.74}{Re^{0.9}}\right)\right]^2}$$

Minor Losses

Turbulent dissipation through pipe fittings, valves, bends, elbows, contractions, and expansions is computed using loss coefficients ($K_L$):

$$h_m = \sum K_L \frac{v^2}{2g}$$

The Extended Energy Equation with Pumps and Turbines

Incorporating shaft work and total head losses ($h_L = h_f + h_m$) into the energy equation between upstream station 1 and downstream station 2:

$$\frac{P_1}{\gamma} + \alpha_1 \frac{v_1^2}{2g} + z_1 + h_{pump} = \frac{P_2}{\gamma} + \alpha_2 \frac{v_2^2}{2g} + z_2 + h_{turbine} + h_L$$

Where $h_{pump}$ is the head added to the fluid by a pump, $h_{turbine}$ is head extracted by a turbine, and kinetic energy correction factor $\alpha \approx 1.0$ for turbulent flow ($\alpha = 2.0$ for fully developed laminar flow).

Pump Hydraulic and Electrical Power

The mechanical power delivered directly to the liquid by the pump impeller is:

$$\dot{W}_{fluid} = \gamma Q h_{pump} = \rho g Q h_{pump} = \Delta P_{pump} Q$$

The electric motor input power required to drive a pump with efficiency $\eta_{pump}$ is:

$$\dot{W}_{brake} = \frac{\dot{W}_{fluid}}{\eta_{pump}} = \frac{\gamma Q h_{pump}}{\eta_{pump}}$$

---

8. Step-by-Step Worked Engineering Examples

Example 1: Steady-Flow Air Compressor Power

Problem: An industrial rotary screw air compressor in an assembly facility compresses air steadily from $P_1 = 100\text{ kPa}$, $T_1 = 300\text{ K}$ to $P_2 = 800\text{ kPa}$, $T_2 = 580\text{ K}$. The mass flow rate of air is $\dot{m} = 0.50\text{ kg/s}$. During compression, heat is rejected to ambient cooling air at a rate of $\dot{Q}_{out} = 25\text{ kW}$. Assuming air behaves as an ideal gas with constant $c_p = 1.005\text{ kJ/(kg}\cdot\text{K)}$, and negligible kinetic and potential energy changes, calculate the required shaft power input ($\dot{W}_{in}$) to the compressor.

Step-by-Step Solution:

1. Identify the system and energy equation: The compressor is an open steady-flow control volume. With $\Delta KE = 0$ and $\Delta PE = 0$: $$\dot{Q} - \dot{W}_s = \dot{m} (h_2 - h_1)$$
2. Determine enthalpy change for an ideal gas: $$\Delta h = h_2 - h_1 = c_p (T_2 - T_1)$$ $$\Delta h = 1.005\text{ kJ/(kg}\cdot\text{K)} \times (580\text{ K} - 300\text{ K}) = 1.005 \times 280 = 281.4\text{ kJ/kg}$$
3. Substitute sign conventions: Heat is rejected, so $\dot{Q} = -25\text{ kW} = -25\text{ kJ/s}$. $$-25\text{ kJ/s} - \dot{W}_s = 0.50\text{ kg/s} \times 281.4\text{ kJ/kg} = 140.7\text{ kW}$$
4. Solve for shaft work: $$-\dot{W}_s = 140.7 - (-25) = 165.7\text{ kW} \implies \dot{W}_s = -165.7\text{ kW}$$ The negative sign denotes work entering the system. The power input required is $\dot{W}_{in} = 165.7\text{ kW}$.

---

Example 2: Pipe Friction Head Loss and Pumping Power

Problem: Water ($\rho = 1,000\text{ kg/m}^3$, dynamic viscosity $\mu = 1.00 \times 10^{-3}\text{ Pa}\cdot\text{s}$) is pumped through a horizontal commercial steel pipe ($D = 0.15\text{ m}$, length $L = 300\text{ m}$) at a flow rate of $Q = 0.03534\text{ m}^3\text{/s}$. The Darcy friction factor is determined to be $f = 0.020$. Minor losses total $\sum K_L = 4.5$. The pump has an efficiency of $\eta = 75\%$. Find the required pump input brake power.

Step-by-Step Solution:

1. Compute flow velocity: $$A = \frac{\pi D^2}{4} = \frac{\pi (0.15\text{ m})^2}{4} = 0.01767\text{ m}^2$$ $$v = \frac{Q}{A} = \frac{0.03534\text{ m}^3\text{/s}}{0.01767\text{ m}^2} = 2.00\text{ m/s}$$
2. Compute velocity head: $$\frac{v^2}{2g} = \frac{(2.00\text{ m/s})^2}{2 \times 9.81\text{ m/s}^2} = \frac{4.00}{19.62} = 0.2039\text{ m}$$
3. Calculate major head loss ($h_f$) and minor head loss ($h_m$): $$h_f = f \frac{L}{D} \frac{v^2}{2g} = 0.020 \times \left(\frac{300}{0.15}\right) \times 0.2039 = 0.020 \times 2,000 \times 0.2039 = 8.155\text{ m}$$ $$h_m = \sum K_L \frac{v^2}{2g} = 4.5 \times 0.2039 = 0.918\text{ m}$$ $$h_L = h_f + h_m = 8.155 + 0.918 = 9.073\text{ m}$$
4. Apply extended energy equation: Because the pipe is horizontal ($z_1 = z_2$) and constant diameter ($v_1 = v_2$), with open discharge to atmospheric pressure ($P_1 = P_2 = P_{atm}$): $$h_{pump} = h_L = 9.073\text{ m}$$
5. Calculate pump input power: $$\dot{W}_{brake} = \frac{\rho g Q h_{pump}}{\eta} = \frac{1,000 \times 9.81 \times 0.03534 \times 9.073}{0.75} = \frac{3,145.5\text{ W}}{0.75} = 4,194\text{ W} = 4.19\text{ kW}$$

---

9. Common FE Exam Traps in Thermodynamics and Fluid Flow