Continuity Equation: Fluid Mechanics Explained

Conservation of mass for a control volume — derived from the Reynolds Transport Theorem by setting one symbol to 1.

Conservation of MassMechanical EngineeringFree preview
⏱️ About 16 min

In the final lesson of Module 3 we built a single generic tool — the Reynolds Transport Theorem — and made a promise: set its intensive property <em>b</em> to 1 and out drops the continuity equation, the most-used equation in all of fluid mechanics. This lesson collects on that promise. We will not re-derive RTT; we will specialize it to mass, read off the result, and then simplify it down to the compact <em>A₁V₁ = A₂V₂</em> that engineers actually carry in their heads. Along the way you will see exactly why a nozzle speeds a fluid up (and why that fact has <em>nothing</em> to do with pressure).

💡
The big idea: Setting <em>b = 1</em> (so <em>B = mass</em>) in the Reynolds Transport Theorem, and using that a system's mass never changes (<em>dm/dt = 0</em>), yields the <strong>integral continuity equation</strong>: the mass inside a control volume changes only by the net mass flowing across its surface — <em>∂/∂t ∫_CV ρ dV + ∫_CS ρ(V·n̂) dA = 0</em>. For steady flow the accumulation term vanishes, leaving <em>Σṁ_out = Σṁ_in</em>. For steady <em>incompressible</em> flow (constant ρ), the density cancels and mass conservation becomes volume conservation, <em>ΣQ_out = ΣQ_in</em>; for one inlet and one outlet this is <em>A₁V₁ = A₂V₂</em>, with <em>ṁ = ρAV</em> the mass flow rate through any port.
🎯 By the end, you'll be able to
  • Derive the integral continuity equation as the b = 1 (mass) case of the Reynolds Transport Theorem
  • State the port forms ṁ = ρAV and Q = AV and relate them by ṁ = ρQ
  • Reduce the general equation for steady flow and for steady incompressible flow
  • Apply A₁V₁ = A₂V₂ to a single-inlet–single-outlet duct and explain why a nozzle accelerates the fluid
📎 Helpful to know first
  • The Reynolds Transport Theorem

Continuity is the Reynolds Transport Theorem with b = 1

In the final lesson of Module 3 we derived a single, generic equation — the Reynolds Transport Theorem (RTT) — and then made a promise: set its intensive property b to 1, and out drops the continuity equation. This lesson collects on that promise. We will not re-derive RTT; we will simply specialize it to mass and read off the most-used equation in fluid mechanics.

Recall the move. The extensive property B = mass m has intensive property b = dm/dm = 1. Plugging b = 1 into RTT, the left side becomes dm/dt of a system — and a system, by definition, always contains the same matter, so dm/dt = 0. Setting the right-hand side to zero gives the integral continuity equation: the mass inside any control volume can change only because mass flows across its boundary. Nothing else. Mass is neither created nor destroyed; it is only transported.

That one line — the accumulation of mass inside the control volume balances the net mass flowing across its surface — is the whole of mass conservation written for a region of space. Every special case we use, from a contracting nozzle to a branching pipe network, is this same equation with some terms dropped. The rest of this lesson is about which terms to drop, and when.

\[ \frac{\partial}{\partial t}\int_{CV}\rho\, dV\;+\;\int_{CS}\rho\,(\mathbf{V}\cdot\hat{\mathbf{n}})\, dA\;=\;0 \]
The integral continuity equation — mass conservation for a control volume, obtained from RTT by setting b = 1 (so B = m) and using dm/dt of a system = 0. The first term is the rate at which mass accumulates inside the CV; the second is the net mass flux out across the control surface (n̂ is the outward normal, so V·n̂ is the outflow speed). Mass inside the CV changes only by crossing the boundary.
🔑 The payoff Module 3 promised

This is the moment the Reynolds Transport Theorem was built for. Last module we carried the generic theorem all the way to the statement 'with b = 1, the accumulation of mass plus the net mass outflow is zero' — and stopped there, one line short. This lesson simply names that line the continuity equation and starts using it. The same template, with b = V or b = e, will hand us the momentum equation (M6) and the energy equation (M5) later — continuity is just the first of three siblings.

A fixed control volume (dashed box) with one inlet port on the left and one outlet port on the right. Mass flows in at the inlet as m-dot-1 = rho A1 V1 and out at the outlet as m-dot-2 = rho A2 V2. The accumulation term (rate of change of the mass integral inside the box) plus the net flux across the boundary sum to zero. inlet (1) ṁ₁ = ρ A₁ V₁ (flows IN) outlet (2) ṁ₂ = ρ A₂ V₂ (flows OUT) Control Volume accumulation: ∂/∂t ∫ ρ dV (mass stored inside) Control Surface (CS) — the dashed boundary, crossed only by the flux terms Steady: accumulation = 0, so ṁ₁ (in) = ṁ₂ (out).

A fixed control volume shown as a dashed box with one inlet port on the left and one outlet port on the right. Mass flows in at the inlet as m-dot-1 = rho A1 V1 and out at the outlet as m-dot-2 = rho A2 V2. Inside the box is the accumulation term (rate of change of the mass integral). For steady flow the accumulation term is zero and the inflow equals the outflow.

A control volume and its control surface. The two parts of the continuity equation are visible: mass accumulating inside the box, and mass crossing the boundary at each port as ṁ = ρAV. For steady flow the inside stops changing and what flows in must equal what flows out.

The port form: ṁ = ρAV and Q = AV

Real hardware — pumps, valves, lengths of pipe — connects to the rest of a system through a small number of openings, or ports. When the flow through each port is steady and roughly uniform across the opening, the surface integral collapses into a sum over ports, and each port contributes a single number: the mass flow rate

where ρ is the density, A the port's cross-sectional area, and V the average speed normal to that area. This port form was introduced in the RTT lesson; here it becomes our working currency. Pair it with the volume flow rate Q = A V (cubic metres per second), and note that ṁ = ρ Q — mass flow is just density times volume flow.

A useful bookkeeping habit concerns sign. At an outlet the velocity points outward (V·n̂ > 0) so the port's ρAV is positive; at an inlet the velocity points inward (V·n̂ < 0) so the same product enters with a minus sign. Keeping those signs straight is what turns a jumble of connected pipes into a solvable mass balance.

\[ \dot m=\rho\, A\, V\quad(\text{mass flow rate}),\qquad Q=A\, V\quad(\text{volume flow rate}),\qquad \dot m=\rho\, Q \]
The port forms used in almost every engineering calculation. ṁ = ρAV is the mass crossing a port per unit time; Q = AV is the volume crossing per unit time; the two are related by ṁ = ρQ. For a fixed control volume with several ports, the integral continuity equation becomes Σ_out(ρAV) − Σ_in(ρAV) plus any accumulation.

Two special cases that do most of the work

The integral form is completely general but is rarely used as-is, because two simplifications apply to almost every engineering flow.

The first is steady flow. 'Steady' means nothing changes with time at any fixed point, so the accumulation term ∂/∂t ∫_CV ρ dV vanishes. What remains says: the net mass flowing out across the control surface is zero — equivalently, the total mass flow rate in equals the total mass flow rate out. In port form, Σṁ_out = Σṁ_in. This steady balance is the form used for the great majority of design calculations.

The second is steady plus incompressible. A liquid such as water is, for practical purposes, incompressible — its density ρ is constant and cancels out of the balance. The mass balance then becomes a volume balance: the total volume flow in equals the total volume flow out, ΣQ_out = ΣQ_in. For the commonest geometry of all — a single inlet and a single outlet — this collapses to the beautifully compact

which says that the product of area and velocity is the same at every cross-section. That single equation is the continuity equation as most engineers actually use it. Squeeze the area (a nozzle) and the velocity rises to keep A·V constant; widen it (a diffuser) and the velocity falls. The duct simulator below lets you feel that trade-off directly.

\[ \text{steady: }\sum_{\text{out}}\dot m=\sum_{\text{in}}\dot m\qquad\big|\qquad \text{steady, incompressible: }\sum_{\text{out}}Q=\sum_{\text{in}}Q\qquad\big|\qquad A_{1}V_{1}=A_{2}V_{2}\ \ (\text{one in, one out}) \]
The three working forms, each more specialized than the last. Steady flow drops the accumulation term (Σṁ_out = Σṁ_in). Adding incompressibility cancels ρ and turns it into a volume balance (ΣQ_out = ΣQ_in). For a single inlet and single outlet this reduces to A₁V₁ = A₂V₂ — the velocity at any section is inversely proportional to its area.
🎮 Interactive: a variable-area duct — feel continuity live LIVE
Predict first: Drag the outlet area A₂ down below the inlet area A₁ and watch the outlet velocity V₂ climb to keep Q = A·V constant. Then toggle 'show mass flow' to confirm ṁ₁ = ρA₁V₁ equals ṁ₂ = ρA₂V₂ at both stations — mass is conserved even as the velocity swings.

An interactive variable-area duct. Three sliders set the inlet area, outlet area, and inlet velocity; animated flow markers move faster through the narrow section, and a live readout shows the outlet velocity V2 = A1 V1 / A2 and the volume flow rate Q, with an optional mass-flow label rho A V shown equal at both stations.

The duct_continuity simulator draws a converging or diverging duct whose shape you control with three sliders (inlet area A₁, outlet area A₂, inlet velocity V₁). Animated flow markers speed up where the duct narrows and slow where it widens, exactly as A₁V₁ = A₂V₂ demands, and the readout reports V₂ and the volume flow Q live. The mass-flow toggle labels ṁ = ρAV at both stations to show they match.
⚠️ Continuity tells you nothing about pressure

Here is the trap that catches almost everyone first. You squeeze a nozzle, the fluid speeds up, and the tempting conclusion is that pressure is pushing it harder. It is not — and continuity cannot even discuss pressure, because pressure never appears in the equation. Continuity is purely kinematic: given the areas and one velocity, it fixes the other velocity, full stop. The fluid speeds up in a nozzle because the area falls and A·V must stay constant — no more, no less. Why the fluid obligingly accelerates, and what happens to the pressure when it does, is the energy equation's job (Module 5). Until then, do not let continuity tempt you into a pressure claim.

📝 Worked example: Water (ρ = 1000 kg/m³) flows steadily through a contracting pipe with inlet area A₁ = 0.008 m², outlet area A₂ = 0.002 m², and inlet velocity V₁ = 3 m/s. Find (a) the outlet velocity V₂, (b) the volume flow rate Q, and (c) the mass flow rate ṁ.
  1. Steady, incompressible, one inlet and one outlet: A₁V₁ = A₂V₂.
  2. Solve for V₂: V₂ = A₁V₁/A₂ = (0.008)(3)/(0.002) = 0.024/0.002 = 12 m/s. The area fell 4:1, so the velocity rose 4:1.
  3. Volume flow rate (same at every section): Q = A₁V₁ = (0.008)(3) = 0.024 m³/s (= 24 L/s).
  4. Check with the outlet: Q = A₂V₂ = (0.002)(12) = 0.024 m³/s. ✓ Same value.
  5. Mass flow rate: ṁ = ρQ = (1000)(0.024) = 24 kg/s. Equivalently ṁ = ρA₁V₁ = ρA₂V₂ = 24 kg/s at both ports.
✓ (a) V₂ = 12 m/s; (b) Q = 0.024 m³/s (24 L/s); (c) ṁ = 24 kg/s.
✏️ Practice: Water flows steadily through a nozzle with inlet area A₁ = 0.005 m², outlet area A₂ = 0.001 m² (a 5:1 contraction), and inlet speed V₁ = 2 m/s. Find the outlet speed V₂, in m/s.
m/s
Solution
  1. One inlet, one outlet, incompressible: A₁V₁ = A₂V₂.
  2. V₂ = A₁V₁/A₂ = (0.005)(2)/(0.001) = 0.010/0.001 = 10 m/s.
  3. The 5:1 area contraction produces a 5:1 velocity increase.
✏️ Practice: Air (ρ = 1.2 kg/m³) flows steadily through a duct of area A = 0.02 m² at an average speed V = 5 m/s. Find the mass flow rate ṁ = ρAV through the duct, in kg/s.
kg/s
Solution
  1. Mass flow rate through a port: ṁ = ρ A V.
  2. ṁ = (1.2)(0.02)(5) = 0.12 kg/s.
  3. This is the quantity that must balance (in = out) at any steady junction — the subject of the next lesson.

Check your understanding

1. The integral continuity equation is obtained from the Reynolds Transport Theorem by setting the intensive property b to:
Mass per unit mass is 1, so b = 1 and B = mass. With dm/dt of a system equal to 0, RTT gives accumulation + net mass outflow = 0 — the continuity equation. The same template with b = V gives momentum; b = e gives energy.
2. For a steady, incompressible flow through a single inlet and a single outlet, continuity reduces to:
Steady drops the accumulation term; incompressibility cancels the density, leaving a volume balance Q = A₁V₁ = A₂V₂. Velocity is inversely proportional to area: halve the area and the velocity doubles.
3. A fluid speeds up as it passes through a contracting nozzle. The correct physical explanation is that:
Continuity fixes the velocity from the area change alone; pressure never enters the equation. Why the fluid accelerates, and what happens to pressure, is the energy equation's job (Module 5). Do not let continuity tempt you into a pressure claim.
✅ Key takeaways
  • Continuity is the Reynolds Transport Theorem with b = 1 (so B = mass): ∂/∂t ∫_CV ρ dV + ∫_CS ρ(V·n̂) dA = 0.
  • The port forms are ṁ = ρAV (mass flow rate) and Q = AV (volume flow rate), related by ṁ = ρQ.
  • Steady flow drops the accumulation term (Σṁ_out = Σṁ_in); steady incompressible flow cancels ρ and becomes a volume balance (ΣQ_out = ΣQ_in).
  • For one inlet and one outlet, A₁V₁ = A₂V₂. A nozzle accelerates the fluid because its area falls — continuity is kinematic and silent on pressure.
➡️ You now hold the continuity equation in all its useful forms. The next lesson applies them to real hardware — nozzles and diffusers (where area alone sets the velocity), branching pipe junctions (where several flows must balance), and filling or draining tanks (where the unsteady accumulation term finally does some work). Every one of them is the same equation with different terms kept.
Want to test yourself on this? Try the Mechanical Aptitude test →
🎓 Go deeper: external courses & trusted references