Viscous Pipe Flow

1. Introduction

Viscous pipe flow describes the motion of fluids inside conduits when the influence of viscosity cannot be ignored. Whether the fluid is water in a municipal distribution network, oil in a refinery trunk line, or coolant in a powerplant loop, the pressure drop, flow rate, and required pipe size are all dictated by the interaction between viscous forces and inertia.

Unlike ideal (inviscid) fluids, a viscous fluid experiences shear stresses that develop a velocity profile across the pipe crosssection. The nature of that profile, together with the pipes roughness and geometry, determines the relationship between pressure gradient and volumetric flow rate.

2. Fundamentals of Viscous Flow

2.1 Shear stress and viscosity

Shear stress, \(\tau\), is proportional to the rate of deformation (velocity gradient) through the dynamic viscosity \(\mu\):

\(\tau = \mu \frac{du}{dy}\)

In a circular pipe the gradient is taken radially. The larger the viscosity, the greater the resistance to motion, and the larger the pressure drop needed to sustain a given flow rate.

2.2 Reynolds number

The dimensionless Reynolds number identifies whether the flow is laminar or turbulent:

\(Re = \frac{\rho V D}{\mu} = \frac{V D}{\nu}\)

where \(\rho\) is fluid density, \(V\) the average velocity, \(D\) the pipe diameter, and \(\nu\) the kinematic viscosity. In smooth circular pipes, flow remains laminar for \(Re \lesssim 2100\) and becomes turbulent for \(Re \gtrsim 4000\). The transitional region between these limits depends on pipe roughness, entrance conditions, and disturbances.

3. Laminar Pipe Flow HagenPoiseuille Law

When the flow is fully developed and laminar, the velocity distribution is parabolic, with zero velocity at the wall (noslip condition) and a maximum at the centreline.

\(u(r) = \frac{\Delta P}{4 \mu L}\,(R^{2} - r^{2})\)

Integrating this profile yields the classical HagenPoiseuille equation for the volumetric flow rate \(Q\):

\(Q = \frac{\pi R^{4}}{8 \mu}\,\frac{\Delta P}{L}\)

or, expressed with diameter \(D = 2R\):

\(Q = \frac{\pi D^{4}}{128 \mu}\,\frac{\Delta P}{L}\)

This result shows the strong dependence of flow on pipe diameter (fourthpower law). Small changes in diameter produce large changes in capacity, a fact exploited in pipesizing calculations.

3.1 Pressure drop in laminar flow

The pressure drop per unit length, \(\Delta P/L\), can be written in terms of the Darcy friction factor \(f\):

\(f = \frac{64}{Re}\)

Substituting into the DarcyWeisbach equation,

\(\frac{\Delta P}{L} = f \frac{\rho V^{2}}{2D}\)

gives a consistency check; for laminar flow the friction factor is inversely proportional to Reynolds number.

4. Turbulent Pipe Flow

In turbulent conditions the velocity profile is flatter, with a thin viscous sublayer close to the wall and a fullmixing core. Because turbulence enhances momentum transfer, the friction factor no longer follows the simple \(64/Re\) relationship.

4.1 Empirical frictionfactor correlations

Several correlations are used to estimate \(f\) for turbulent flow. The most widely adopted is the ColebrookWhite equation, which relates \(f\), Reynolds number, and relative roughness \(\varepsilon/D\):

\(\frac{1}{\sqrt{f}} = -2.0\log_{10}\!\Bigg[\frac{\varepsilon/D}{3.7} + \frac{2.51}{Re\sqrt{f}}\Bigg]\)

Because the equation is implicit in \(f\), iterative methods or explicit approximations (e.g., the SwameeJain formula) are normally employed.

4.2 Moody diagram

The Moody chart graphically presents the ColebrookWhite relationship. By locating a point representing the pipes Reynolds number and roughness ratio, the corresponding friction factor can be read directly.

4.3 Pressure drop for turbulent flow

Once \(f\) is known, the pressure drop follows the DarcyWeisbach form:

\(\Delta P = f\,\frac{L}{D}\,\frac{\rho V^{2}}{2}\)

Unlike laminar flow, the dependence on velocity is quadratic, and the sensitivity to pipe roughness can dominate the loss calculation when the Reynolds number is high.

5. Design and Practical Considerations

5.1 Choosing pipe diameter

Design usually begins with a required flow rate \(Q\) and a permissible pressure drop \(\Delta P_{max}\). For a chosen fluid (known \(\mu,\rho\)) and pipe length \(L\), an initial diameter guess is made using the laminar HagenPoiseuille formula. The resulting Reynolds number determines whether the guess leads to laminar or turbulent flow. If turbulent, the friction factor is updated with an appropriate correlation, and the calculation is repeated until convergence.

5.2 Effects of pipe roughness

Roughness \(\varepsilon\) varies with material (e.g., smooth steel, cast iron, PVC). In turbulent flow, a higher \(\varepsilon/D\) raises the friction factor, increasing pressure loss. Correctly accounting for roughness is essential when sizing long runs or when using aged, corroded pipelines.

5.3 Minor losses

In addition to frictional loss, bends, valves, fittings and expansions introduce localized losses. These are expressed as a head loss coefficient \(K\) and added to the major DarcyWeisbach loss:

\(\Delta P_{total}= \Delta P_{friction}+ \sum K\,\frac{\rho V^{2}}{2}\)

5.4 Temperature effects

Viscosity is temperaturedependent; for liquids, \(\mu\) often decreases markedly as temperature rises. Designers therefore must evaluate the fluids temperature profile along the pipe and adjust the viscosity accordingly, especially for heated processes or cooling loops.

5.5 Safety and surge considerations

Rapid changes in flow (pump startup, valve closure) generate pressure surges (water hammer). The wave speed \(c\) is influenced by pipe elasticity and fluid compressibility. Adding surge tanks, air chambers, or slowclosing valves mitigates the transient pressure spikes that could otherwise exceed the pipes design limits.

6. Summary

• Viscous pipe flow is governed by the balance between shear stresses (viscosity) and inertial forces.
• The Reynolds number distinguishes laminar (<2100) from turbulent flow (>4000).
• Laminar flow obeys the HagenPoiseuille law, with a friction factor \(f = 64/Re\) and a pressure drop proportional to the flow rate.
• Turbulent flow requires empirical correlations (ColebrookWhite, SwameeJain) or the Moody diagram to obtain the friction factor; the pressure drop varies with the square of velocity.
• Design involves iterative sizing: select a diameter, compute Reynolds number, determine the appropriate friction factor, and verify that the pressure drop meets system constraints.
• Pipe roughness, minor losses, temperaturedependent viscosity, and transient events such as water hammer must be incorporated for reliable, safe designs.

Understanding the theory and applying the correct correlations allows engineers to predict losses accurately, optimise pipe networks, and ensure longterm operational reliability.