Skip to main content

Control-Volume Mass and Energy Analysis

Most engineering devices are open systems: turbines, compressors, nozzles, diffusers, pumps, boilers, condensers, throttling valves, mixing chambers, and heat exchangers all have mass crossing their boundaries. Control-volume analysis keeps the boundary fixed around the device and accounts for the energy carried with flowing matter.

The key new property is enthalpy. Flow work PvPv is automatically included in h=u+Pvh=u+Pv, so the steady-flow energy equation can be written compactly. Cengel's method is device-oriented: after writing the general balance, simplify it using what is known about a particular device. A turbine emphasizes shaft work; a nozzle emphasizes kinetic energy; a heat exchanger emphasizes heat transfer between streams; a throttle is nearly isenthalpic.

Definitions

  • A control volume is a selected region in space through which mass may flow. Its boundary is the control surface.
  • Mass flow rate is m˙=ρVnA\dot m=\rho V_n A, where VnV_n is the velocity component normal to the area.
  • Volume flow rate is V˙=VnA\dot V=V_nA for uniform flow, and m˙=ρV˙=V˙/v\dot m=\rho\dot V=\dot V/v.
  • Steady flow means properties at any fixed point in the control volume do not change with time. Mass and energy stored in the control volume are then constant.
  • Flow work is the work needed to push mass into or out of a control volume. It appears as PvPv per unit mass and combines with uu to form enthalpy.
  • A nozzle increases velocity at the expense of enthalpy. A diffuser increases pressure by slowing a stream.
  • A turbine produces shaft work from a pressure and enthalpy drop. A compressor or pump consumes shaft work to raise pressure.
  • A throttling valve is a restriction with negligible heat transfer, work, and kinetic-potential changes; it is commonly modeled as h1h2h_1\approx h_2.
  • A mixing chamber combines streams, usually with no shaft work and negligible kinetic and potential energy changes.
  • A heat exchanger transfers heat between streams without shaft work; the outer boundary may be nearly adiabatic while heat flows internally between fluids.

The reliable sequence is mass balance first, energy balance second, device assumptions third. The mass balance may determine an unknown flow rate or mixing fraction before the energy equation can be evaluated. In transient filling and emptying problems, storage terms remain and the steady-flow simplification is not valid. For this topic, a complete engineering model should state the boundary, the time basis, the property model, and the sign convention before any numbers are substituted. In control-volume mass and energy analysis, that habit is especially important because several formulas look similar while answering different physical questions. A closed-system expression, a steady-flow expression, an ideal-gas relation, and a property-table interpolation may all contain pressure, temperature, or enthalpy, but they do not have the same assumptions. The safest workflow is to write the general balance or defining relation first, cancel terms with a written reason, and only then insert table values or constants.

The second modeling habit is to keep the basis visible. Some calculations are per unit mass, some per mole, some per kg dry air, and some per unit time. A correct formula on the wrong basis is a common source of errors that look numerically plausible. When a table gives kJ/kg\mathrm{kJ/kg}, multiply by m˙\dot m to get kW\mathrm{kW}; when a reaction is balanced in kmol, convert to mass only after the element balance is complete; when a mixture property uses mole fraction, do not substitute mass fraction without conversion.

Key results

The conservation of mass for a control volume is

dmCVdt=m˙inm˙out.\frac{dm_{\mathrm{CV}}}{dt} =\sum \dot m_{\mathrm{in}}-\sum \dot m_{\mathrm{out}}.

At steady state,

m˙in=m˙out.\sum \dot m_{\mathrm{in}}=\sum \dot m_{\mathrm{out}}.

The general steady-flow energy equation is

Q˙W˙+inm˙(h+V22+gz)=outm˙(h+V22+gz).\dot Q-\dot W +\sum_{\mathrm{in}}\dot m\left(h+\frac{V^2}{2}+gz\right) = \sum_{\mathrm{out}}\dot m\left(h+\frac{V^2}{2}+gz\right).

For one-inlet, one-exit steady devices, divide by m˙\dot m:

qw=(h2h1)+V22V122+g(z2z1).q-w=(h_2-h_1)+\frac{V_2^2-V_1^2}{2}+g(z_2-z_1).

Common reductions are:

turbine: wouth1h2,compressor: winh2h1,nozzle: h1h2V22V122,throttle: h1h2.\begin{aligned} \text{turbine: } & w_{\mathrm{out}}\approx h_1-h_2, \\ \text{compressor: } & w_{\mathrm{in}}\approx h_2-h_1, \\ \text{nozzle: } & h_1-h_2\approx \frac{V_2^2-V_1^2}{2}, \\ \text{throttle: } & h_1\approx h_2. \end{aligned}

Remember that velocity terms in SI are in m2/s2\mathrm{m^2/s^2}; divide by 10001000 to convert to kJ/kg\mathrm{kJ/kg}. These results should be read as a hierarchy rather than a list of isolated equations. Conservation of mass and energy set the allowed accounting; property relations supply the missing state data; the second law or equilibrium criterion decides direction, limits, and losses. A numerical answer is not finished until it passes three checks: the units reduce to the requested quantity, the sign matches the stated energy or entropy transfer direction, and the magnitude is reasonable compared with a limiting case. Useful limiting cases include zero heat transfer, reversible operation, incompressible behavior, ideal-gas behavior, saturated-liquid or saturated-vapor endpoints, and equal reservoir temperatures.

Because the textbook often moves between exact laws and engineering approximations, the approximation should be named in the solution. Examples include constant specific heats, negligible kinetic energy, negligible pump work, adiabatic devices, isentropic turbomachinery, ideal-gas mixtures, dry-air approximations, and linear interpolation. Naming the approximation makes later refinement straightforward: replace the approximate property model or restore the neglected term without rebuilding the whole analysis.

Visual

DeviceDominant conversionUsual simplificationsCore result
Nozzleenthalpy to kinetic energyq0q\approx0, w=0w=0, Δz0\Delta z\approx0V2V_2 rises as hh falls
Diffuserkinetic energy to pressure/enthalpyq0q\approx0, w=0w=0velocity falls
Turbineenthalpy to shaft workadiabatic, small ΔKE\Delta KE, ΔPE\Delta PEW˙m˙(h1h2)\dot W\approx\dot m(h_1-h_2)
Compressorshaft work to enthalpyadiabatic, small ΔKE\Delta KEW˙inm˙(h2h1)\dot W_{in}\approx\dot m(h_2-h_1)
Throttlepressure drop at nearly constant enthalpyq0q\approx0, w=0w=0h1h2h_1\approx h_2
Heat exchangerheat between streamsouter boundary adiabaticenthalpy lost by hot stream equals gained by cold stream

Worked example 1: adiabatic turbine power

Problem. Steam enters an adiabatic turbine at h1=3200 kJ/kgh_1=3200\ \mathrm{kJ/kg} and exits at h2=2500 kJ/kgh_2=2500\ \mathrm{kJ/kg}. The mass flow rate is 2.0 kg/s2.0\ \mathrm{kg/s}. Kinetic and potential energy changes are negligible. Find the turbine power output.

Method.

  1. Steady mass balance gives the same mass flow rate at inlet and exit.
  2. For an adiabatic turbine with negligible kinetic and potential changes:
W˙out=m˙(h1h2).\dot W_{\mathrm{out}}=\dot m(h_1-h_2).
  1. Substitute:
W˙out=(2.0)(32002500)=1400 kJ/s.\dot W_{\mathrm{out}} =(2.0)(3200-2500) =1400\ \mathrm{kJ/s}.
  1. Convert kJ/s\mathrm{kJ/s} to kW\mathrm{kW}:
W˙out=1400 kW=1.4 MW.\dot W_{\mathrm{out}}=1400\ \mathrm{kW}=1.4\ \mathrm{MW}.

Checked answer. The positive output is expected because the steam enthalpy drops through the turbine.

Worked example 2: exit velocity from an adiabatic nozzle

Problem. Air enters a nozzle with h1=430 kJ/kgh_1=430\ \mathrm{kJ/kg} and V1=40 m/sV_1=40\ \mathrm{m/s}. It exits with h2=390 kJ/kgh_2=390\ \mathrm{kJ/kg}. Heat transfer, work, and elevation change are negligible. Find the exit velocity.

Method.

  1. The nozzle energy balance is
h1+V122=h2+V222h_1+\frac{V_1^2}{2}=h_2+\frac{V_2^2}{2}

when both kinetic terms are in the same units as enthalpy.

  1. Convert the inlet kinetic energy:
V122=4022=800 m2/s2=0.800 kJ/kg.\frac{V_1^2}{2}=\frac{40^2}{2}=800\ \mathrm{m^2/s^2}=0.800\ \mathrm{kJ/kg}.
  1. Solve for exit kinetic energy:
V222=(h1h2)+0.800=40.8 kJ/kg.\frac{V_2^2}{2} =(h_1-h_2)+0.800 =40.8\ \mathrm{kJ/kg}.
  1. Convert back to m2/s2\mathrm{m^2/s^2} and solve:
V2=2(40.8)(1000)=286 m/s.V_2=\sqrt{2(40.8)(1000)}=286\ \mathrm{m/s}.

Checked answer. The exit velocity is about 286 m/s286\ \mathrm{m/s}. The velocity increase is large because a 40 kJ/kg40\ \mathrm{kJ/kg} enthalpy drop is a large kinetic-energy change.

Code

import math

def turbine_power(mdot, h1, h2):
return mdot * (h1 - h2) # kW

def nozzle_exit_velocity(h1, h2, V1):
ke2_kJkg = (h1 - h2) + V1**2 / 2000.0
return math.sqrt(2.0 * ke2_kJkg * 1000.0)

print(turbine_power(2.0, 3200.0, 2500.0))
print(nozzle_exit_velocity(430.0, 390.0, 40.0))

Common pitfalls

  • Using internal energy instead of enthalpy for steady-flow devices with flow work.
  • Forgetting the factor of 10001000 when converting V2/2V^2/2 to kJ/kg\mathrm{kJ/kg}.
  • Applying steady-flow equations to startup, shutdown, filling, or emptying without storage terms.
  • Assuming every valve or nozzle is adiabatic without checking heat transfer or residence time.
  • Dropping mass balance in heat exchangers and mixing chambers before determining flow rates.
  • Starting from a special-case equation before checking that its assumptions actually hold. Write the general balance or definition first, then reduce it.
  • Leaving property-table values unlabeled. Record the substance, phase region, pressure or temperature row, interpolation fraction, and units so the result can be audited.
  • Rounding intermediate states too aggressively. Keep extra digits through property lookup, quality calculation, and efficiency ratios, then round the final answer to justified precision.
  • Skipping a limiting-case check. Test the result against reversible operation, zero pressure drop, saturated endpoints, ideal-gas behavior, or equal-temperature reservoirs when those limits are meaningful.
  • Treating a numerical solver or chart as a substitute for physical reasoning. Software can return a precise-looking number even when the selected phase, reference state, or boundary model is wrong.
  • Forgetting to state whether the reported answer is specific, total, or rate based.

Connections