11.3 Viscous Flow in Pipes, Boundary Layer & Minor Losses

Viscous Flow in Pipes, Boundary Layer & Minor Losses

In Coal India mining complexes, fluid transport encompasses massive pipeline systems: deep sump mine dewatering mains, sand-stowing slurry lines, hydraulic transport of coal fines, and pressurized mine fire-fighting networks. Analyzing viscous dissipation, pressure drops, pipe network equivalencies, and boundary layer mechanics is vital for sizing pumps, pipes, and ventilation conduits.

---

1. Reynolds Number and Flow Regimes

The Reynolds Number ($Re$) represents the ratio of dynamic inertia forces to viscous shear forces:

$$Re = \frac{\text{Inertia Force}}{\text{Viscous Force}} = \frac{\rho V D}{\mu} = \frac{V D}{\nu}$$

Where:

• $V$ is mean flow velocity ($\text{m/s}$)
• $D$ is internal pipe diameter ($\text{m}$) or characteristic hydraulic diameter $D_h = \frac{4A}{P_w}$
• $\mu$ is dynamic viscosity ($\text{Pa}\cdot\text{s}$), $\nu$ is kinematic viscosity ($\text{m}^2/\text{s}$)

Critical Reynolds Numbers in Engineering Geometry:

---

2. Hagen-Poiseuille Laminar Flow in Circular Pipes

For steady, fully developed, incompressible laminar flow through a straight horizontal circular pipe of radius $R$ and diameter $D = 2R$:

               Pipe Centerline (r = 0, u = u_max, tau = 0)
      ---------------------------------------------------------
       r ^        Parabolic Velocity Profile u(r)
         |           /-------------------\
         |          /                     \
         |         /                       \
         +--------+-------------------------+-----------------> x
         |         \                       /
         |          \                     /
         v           \-------------------/
      ---------------------------------------------------------
               Pipe Wall (r = R, u = 0, tau = tau_w = max)

2.1 Shear Stress Distribution ($\tau(r)$)

Considering a coaxial cylindrical fluid element of radius $r$ and length $dx$, force equilibrium yields:

$$\tau(r) = -\frac{r}{2} \left(\frac{dP}{dx}\right) = \tau_w \left(\frac{r}{R}\right)$$

• Linear Distribution: $\tau = 0$ at pipe centerline ($r = 0$), and reaches maximum value $\tau_w$ at the wall ($r = R$): $$\tau_w = -\frac{R}{2} \left(\frac{dP}{dx}\right) = \frac{\Delta P \cdot D}{4L}$$
• Note: The linear variation of shear stress across the pipe radius holds true for both laminar and turbulent flows because it is derived purely from static force equilibrium.

2.2 Parabolic Velocity Profile ($u(r)$)

Using $\tau = -\mu \frac{du}{dr}$ and integrating with no-slip boundary condition $u(R) = 0$:

$$\nu(r) = -\frac{1}{4\mu} \left(\frac{dP}{dx}\right) (R^2 - r^2) = u_{\max} \left(1 - \frac{r^2}{R^2}\right)$$

• Maximum Velocity ($u_{\max}$ at $r = 0$): $$u_{\max} = -\frac{R^2}{4\mu} \left(\frac{dP}{dx}\right) = \frac{\Delta P \cdot R^2}{4\mu L} = \frac{\Delta P \cdot D^2}{16\mu L}$$

2.3 Mean / Average Velocity ($V_{\text{avg}}$ or $\bar{u}$)

Integrating discharge $Q = \int_0^R u(r) \cdot 2\pi r \, dr$ and dividing by area $A = \pi R^2$:

$$V_{\text{avg}} = \frac{u_{\max}}{2} = -\frac{R^2}{8\mu}\left(\frac{dP}{dx}\right) = \frac{\Delta P \cdot D^2}{32\mu L}$$

Core GATE / PSU Relation: In fully developed laminar pipe flow, $u_{\max} = 2 \times V_{\text{avg}}$.

• Radial Position where Local Velocity equals Average Velocity: $$u(r) = V_{\text{avg}} = \frac{u_{\max}}{2} \implies u_{\max}\left(1 - \frac{r^2}{R^2}\right) = \frac{u_{\max}}{2}$$ $$1 - \frac{r^2}{R^2} = 0.5 \implies \frac{r^2}{R^2} = 0.5 \implies r = \frac{R}{\sqrt{2}} \approx 0.707 R$$

2.4 Hagen-Poiseuille Pressure Drop and Discharge Equations

$$\Delta P = P_1 - P_2 = \frac{32 \mu V_{\text{avg}} L}{D^2} = \frac{128 \mu Q L}{\pi D^4}$$ $$Q = \frac{\pi \Delta P D^4}{128 \mu L}$$

Kinetic Energy and Momentum Correction Factors:

• Kinetic Energy Correction Factor ($\alpha = \frac{1}{A}\int (u/V)^3 dA$):
  • For Laminar circular pipe flow: $\alpha = 2.0$
  • For Turbulent circular pipe flow: $\alpha \approx 1.03 - 1.06$
• Momentum Correction Factor ($\beta = \frac{1}{A}\int (u/V)^2 dA$):
  • For Laminar circular pipe flow: $\beta = \frac{4}{3} \approx 1.33$
  • For Turbulent circular pipe flow: $\beta \approx 1.01 - 1.03$

---

3. Darcy-Weisbach Equation & Pipe Friction Factors

The fundamental equation for major frictional head loss $h_f$ in any conduit is the Darcy-Weisbach Formula:

$$h_f = \frac{f L V^2}{2g D} = \frac{4 f' L V^2}{2g D} = \frac{f L Q^2}{12.1 D^5}$$

Where:

• $f$ is the dimensionless Darcy Friction Factor (or Darcy-Weisbach friction coefficient)
• $f'$ is the Fanning Friction Factor ($f = 4f'$)

3.1 Friction Factor in Laminar Flow

Equating Darcy-Weisbach loss $h_f = \frac{\Delta P}{\rho g} = \frac{32 \mu V L}{\rho g D^2}$ to $\frac{f L V^2}{2g D}$:

$$\frac{f L V^2}{2g D} = \frac{32 \mu V L}{\rho g D^2} \implies f = \frac{64 \mu}{\rho V D} = \frac{64}{Re}$$ $$f_{\text{laminar}} = \frac{64}{Re}, \quad f'_{\text{Fanning}} = \frac{16}{Re}$$

Critical Note: In laminar flow, the friction factor depends only on Reynolds number and is completely independent of pipe relative surface roughness $\epsilon/D$.

3.2 Friction Factor in Turbulent Flow (Moody Chart)

1. Smooth Pipes ($4000 < Re < 10^5$ - Blasius Formula): $$f = \frac{0.3164}{Re^{1/4}} = \frac{0.3164}{Re^{0.25}}$$
2. Rough Pipes in Fully Rough / Wholly Turbulent Regime: Viscous sublayer thickness $\delta_v \approx \frac{11.6\nu}{u_{\tau}}$ becomes much smaller than roughness asperities $k_s$. Friction factor becomes independent of $Re$ and depends purely on relative roughness (Nikuradse equation): $$\frac{1}{\sqrt{f}} = 2.0 \log_{10}\left(\frac{D}{2\epsilon}\right) + 1.74$$

---

4. Minor Losses in Pipe Systems

Minor losses occur due to local flow disturbances, geometry transitions, valves, and fittings:

$$h_L = k_L \frac{V^2}{2g}$$

     1. Sudden Enlargement:             2. Sudden Contraction:
     V_1                       V_2      V_1                    V_2
     =========\                 ======  ======                /=======
               \               /              \   Vena       /
      Recirc.   \-------------/                \-Contracta--/
      Eddies                                     Recirc. Eddies

4.1 Sudden Expansion (Borda-Carnot Formula)

Due to momentum deceleration and eddy dissipation in the corners:

$$h_{\text{exp}} = \frac{(V_1 - V_2)^2}{2g} = \frac{V_1^2}{2g}\left(1 - \frac{A_1}{A_2}\right)^2$$

• Pipe Exit Loss (Discharge into large sump/reservoir where $A_2 \to \infty, V_2 = 0$): $$h_{\text{exit}} = \frac{V_1^2}{2g}$$

4.2 Sudden Contraction

Flow converges to a vena contracta ($A_c = C_c A_2$) before expanding to fill the smaller pipe:

$$h_c = \frac{(V_c - V_2)^2}{2g} = \left(\frac{1}{C_c} - 1\right)^2 \frac{V_2^2}{2g} = k_c \frac{V_2^2}{2g}$$

For typical water flow with $C_c \approx 0.62$:

$$k_c = \left(\frac{1}{0.62} - 1\right)^2 \approx (0.613)^2 \approx 0.375 - 0.50 \implies h_c \approx 0.5 \frac{V_2^2}{2g}$$

4.3 Summary of Minor Loss Coefficients:

---

5. Pipe Networks: Series, Parallel, and Siphons

5.1 Pipes in Series & Dupuit's Equivalent Pipe Equation

When pipes of different lengths ($L_1, L_2, \dots$) and diameters ($D_1, D_2, \dots$) are connected in series:

1. Discharge is constant: $Q = Q_1 = Q_2 = Q_3 = \dots$
2. Total friction head loss is additive: $H_L = h_{f1} + h_{f2} + h_{f3} + \dots$

Using Darcy-Weisbach $h_f = \frac{f L Q^2}{12.1 D^5}$ and assuming a uniform friction factor $f$:

$$\frac{f L_{\text{eq}} Q^2}{12.1 D_{\text{eq}}^5} = \frac{f L_1 Q^2}{12.1 D_1^5} + \frac{f L_2 Q^2}{12.1 D_2^5} + \frac{f L_3 Q^2}{12.1 D_3^5} + \dots$$ $$\frac{L_{\text{eq}}}{D_{\text{eq}}^5} = \frac{L_1}{D_1^5} + \frac{L_2}{D_2^5} + \frac{L_3}{D_3^5} + \dots = \sum_{i=1}^n \frac{L_i}{D_i^5}$$

This is Dupuit's Equation for equivalent pipelines.

5.2 Pipes in Parallel

When a main pipe branches into two or more parallel lines that rejoin downstream:

1. Head loss across each parallel branch is identical: $$h_{f1} = h_{f2} = h_{f3} = \dots = h_{f,\text{total}}$$ $$\frac{f_1 L_1 Q_1^2}{12.1 D_1^5} = \frac{f_2 L_2 Q_2^2}{12.1 D_2^5}$$
2. Total discharge is the sum of branch flows: $$Q_{\text{total}} = Q_1 + Q_2 + Q_3 + \dots$$

---

5.3 Hydraulic Gradient Line (HGL) and Total Energy Line (TEL)

    TEL (Total Energy Line = P/gamma + V^2/2g + z)  -- Slopes downward in flow direction
     \------------------------------------------------------\ 
      \  V^2/2g (Dynamic Velocity Head)                      \
       \------------------------------------------------------\ HGL (P/gamma + z)
        \                                                      \
         \ Pipe Conduit                                         \
    =================================================================
          Datum (z = 0)

• Total Energy Line (TEL): Graphical representation of the total head ($H = \frac{P}{\gamma} + \frac{V^2}{2g} + z$) along the pipeline. In real fluid flow without energy input (pumps), TEL must always slope downward in the direction of flow due to continuous frictional dissipation.
• Hydraulic Gradient Line (HGL): Graphical representation of piezometric head ($h_{\text{piezo}} = \frac{P}{\gamma} + z$).
• Vertical Separation: The vertical distance between TEL and HGL at any point is identically equal to the velocity head $\frac{V^2}{2g}$.
• Siphon and Vacuum Risk: If the pipeline profile rises above the HGL, $\frac{P}{\gamma}$ becomes negative (gauge vacuum). In water pipelines, if the absolute pressure drops below the saturation vapor pressure ($P_v \approx 2.34\text{ kPa}$ abs at $20^{\circ}\text{C} \implies -7.5\text{ to } -7.7\text{ m of water gauge}$), cavitation and vapor lock occur, breaking the siphon action.

---

6. Prandtl's Boundary Layer Theory

Introduced by Ludwig Prandtl (1904), boundary layer theory divides high Reynolds number flow around a body into two domains:

1. Thin Boundary Layer Region: A very thin layer adjacent to solid boundaries where velocity gradients $\frac{\partial u}{\partial y}$ are large, and viscous shear stresses cannot be neglected.
2. Outer Potential Flow Region: Outside the boundary layer, velocity gradients are near zero; fluid can be treated as frictionless/inviscid and irrotational (governed by Euler/Bernoulli equations).

                  U_inf (Free Stream Velocity)
          --------------------------------------------->
          
          y ^         delta(x) Boundary Layer Thickness
            |         .-'''''''''''''''''''''-.   Outer Inviscid Flow
            |       .'                         `.
            |      /                             \ (u = 0.99 U_inf)
            |     /  Boundary Layer (Viscous)     \
            +----+---------------------------------+-------------> x
            Leading Edge (x = 0)                   Plate Surface (y = 0, u = 0)

6.1 Characteristic Boundary Layer Thickness Parameters

1. Nominal Boundary Layer Thickness ($\delta$ or $\delta_{99}$): The normal distance from the surface where local velocity $u(y)$ reaches $99\%$ of the free stream velocity ($u = 0.99 U_{\infty}$).
2. Displacement Thickness ($\delta^*$): The distance by which the solid boundary must be shifted outward to compensate for the reduction in mass flow rate caused by boundary layer friction: $$\delta^* = \int_0^\delta \left(1 - \frac{u}{U_{\infty}}\right) dy$$
3. Momentum Thickness ($\theta$): The distance by which the boundary must be shifted to compensate for the loss of momentum flux in the boundary layer: $$\theta = \int_0^\delta \frac{u}{U_{\infty}}\left(1 - \frac{u}{U_{\infty}}\right) dy$$
4. Energy Thickness ($\delta^{**}$): The distance compensating for the kinetic energy deficit: $$\delta^{**} = \int_0^\delta \frac{u}{U_{\infty}}\left(1 - \frac{u^2}{U_{\infty}^2}\right) dy$$
5. Shape Factor ($H$): Dimensionless ratio of displacement thickness to momentum thickness: $$H = \frac{\delta^*}{\theta}$$
   • For Blasius laminar boundary layer: $H = \frac{1.72}{0.664} \approx 2.59$
   • For turbulent boundary layer ($1/7\text{th}$ power law): $H \approx 1.286$

---

6.2 Blasius Exact Solution for Laminar Boundary Layer on Flat Plate

For a smooth flat plate of length $L$ parallel to uniform flow $U_{\infty}$ ($Re_x = \frac{\rho U_{\infty} x}{\mu} < 5\times 10^5$):

$$\frac{\delta}{x} = \frac{5.0}{\sqrt{Re_x}}$$ $$\frac{\delta^*}{x} = \frac{1.72}{\sqrt{Re_x}}, \quad \frac{\theta}{x} = \frac{0.664}{\sqrt{Re_x}}$$

• Local Skin Friction Coefficient ($C_{fx}$): $$C_{fx} = \frac{\tau_w}{\frac{1}{2}\rho U_{\infty}^2} = \frac{0.664}{\sqrt{Re_x}}$$
• Average Drag Coefficient ($C_D$ or $\bar{C}_f$ over total plate length $L$): $$C_D = \frac{1}{L}\int_0^L C_{fx} \, dx = \frac{1.328}{\sqrt{Re_L}} = 2 \times C_{f(x=L)}$$
• Total Skin Friction Drag Force on One Side of Plate ($F_D$): $$F_D = C_D \cdot \left(\frac{1}{2}\rho U_{\infty}^2\right) \cdot (b \cdot L) = \frac{1.328}{\sqrt{Re_L}} \left(\frac{1}{2}\rho U_{\infty}^2\right) (b L)$$

---

7. Boundary Layer Separation & Drag Forces

     Accelerating Flow (dP/dx < 0):      Decelerating Flow (dP/dx > 0):
     Favorable Gradient                  Adverse Gradient -> SEPARATION!
     Boundary Layer Attached             (du/dy)_wall = 0 at separation point S

     ----->                              ----->              Recirculation
     ------->  (Attached)                ------>       S    / Backflow
     --------->                          ----->       *    (~~~~~~~~~~~)
     =========================           ============================

7.1 Boundary Layer Separation Mechanics

Separation occurs when fluid layers adjacent to the wall lose kinetic energy due to wall shear friction and are forced to advance against an adverse pressure gradient ($\frac{dP}{dx} > 0$, decelerating flow).

Wall Velocity Gradient Criteria at $y = 0$:

1. $\left(\frac{\partial u}{\partial y}\right)_{y=0} > 0$: Flow remains firmly attached; shear stress $\tau_w > 0$.
2. $\left(\frac{\partial u}{\partial y}\right)_{y=0} = 0$: Point of Imminent Separation; wall shear stress drops to zero ($\tau_w = 0$).
3. $\left(\frac{\partial u}{\partial y}\right)_{y=0} < 0$: Flow has separated; backflow, vortices, and wake zone develop.

7.2 Drag and Lift on Submerged Bodies

• Total Drag ($F_D$): Component of total fluid dynamic force parallel to free-stream velocity: $$F_D = \underbrace{F_{D,\text{friction}}}_{\text{Skin Friction Drag}} + \underbrace{F_{D,\text{pressure}}}_{\text{Form / Pressure Drag}} = \frac{1}{2} C_D \rho A U_{\infty}^2$$
• Streamlined Body (e.g., Airfoil, Hydrofoil): Separation is delayed until the extreme trailing edge; Form Drag $\approx 10\%$, Skin Friction Drag $\approx 90\%$.
• Bluff Body (e.g., Sphere, Cylinder, Mine Haul Truck): Early boundary layer separation creates a massive low-pressure wake; Form Drag $\approx 90\%$, Skin Friction Drag $\approx 10\%$.

---

8. Worked Engineering Examples

Problem 1: Hagen-Poiseuille Viscous Mine Pipeline Flow

A lubricating oil of dynamic viscosity $\mu = 0.29\text{ Pa}\cdot\text{s} = 0.29\text{ N}\cdot\text{s/m}^2$ and relative density $S = 0.90$ ($\rho = 900\text{ kg/m}^3$) is pumped through a horizontal steel pipe of internal diameter $D = 100\text{ mm} = 0.10\text{ m}$ and length $L = 500\text{ m}$ at a volumetric flow rate of $Q = 0.015\text{ m}^3/\text{s}$. Determine:

1. The Reynolds number and verify the flow regime.
2. The pressure drop ($\Delta P$) across the pipeline.
3. The maximum flow velocity ($u_{\max}$) and the wall shear stress ($\tau_w$).
4. The theoretical pumping power required.

Solution:

1. Mean Flow Velocity ($V_{\text{avg}}$) and Reynolds Number: $$A = \frac{\pi D^2}{4} = \frac{\pi \times (0.10)^2}{4} = 0.007854\text{ m}^2$$ $$V_{\text{avg}} = \frac{Q}{A} = \frac{0.015}{0.007854} = 1.910\text{ m/s}$$ $$Re = \frac{\rho V D}{\mu} = \frac{900 \times 1.910 \times 0.10}{0.29} = \frac{171.9}{0.29} = 592.76$$ Since $Re = 592.76 < 2000$, the flow is strictly in the laminar flow regime.
2. Pressure Drop ($\Delta P$): $$\Delta P = \frac{128 \mu Q L}{\pi D^4} = \frac{128 \times 0.29 \times 0.015 \times 500}{\pi \times (0.10)^4} = \frac{278.4}{\pi \times 10^{-4}} = 886,175\text{ Pa} = 886.18\text{ kPa}$$
3. Maximum Velocity ($u_{\max}$) and Wall Shear Stress ($\tau_w$): $$u_{\max} = 2 \times V_{\text{avg}} = 2 \times 1.910 = 3.82\text{ m/s}$$ $$\tau_w = \frac{\Delta P \cdot D}{4L} = \frac{886175 \times 0.10}{4 \times 500} = \frac{88617.5}{2000} = 44.31\text{ N/m}^2 = 44.31\text{ Pa}$$
4. Pumping Power Required ($P_{\text{pump}}$): $$P_{\text{pump}} = Q \cdot \Delta P = 0.015\text{ m}^3/\text{s} \times 886,175\text{ N/m}^2 = 13,292.6\text{ W} = 13.29\text{ kW}$$

---

Problem 2: Dupuit's Equivalent Pipe Sizing

Three pipes in series connect two mine drainage reservoirs: Pipe 1 ($L_1 = 400\text{ m}, D_1 = 200\text{ mm}$), Pipe 2 ($L_2 = 300\text{ m}, D_2 = 150\text{ mm}$), and Pipe 3 ($L_3 = 200\text{ m}, D_3 = 100\text{ mm}$). Determine the equivalent length of a uniform pipe of diameter $D_{\text{eq}} = 200\text{ mm}$.

Solution:

Applying Dupuit's Equation with $D_{\text{eq}} = D_1 = 0.20\text{ m}$:

$$L_{\text{eq}} = D_{\text{eq}}^5 \left(\frac{L_1}{D_1^5} + \frac{L_2}{D_2^5} + \frac{L_3}{D_3^5}\right)$$ $$L_{\text{eq}} = L_1 + L_2\left(\frac{D_{\text{eq}}}{D_2}\right)^5 + L_3\left(\frac{D_{\text{eq}}}{D_3}\right)^5$$ $$\frac{D_{\text{eq}}}{D_2} = \frac{200}{150} = \frac{4}{3} \implies \left(\frac{4}{3}\right)^5 = \frac{1024}{243} \approx 4.214$$ $$\frac{D_{\text{eq}}}{D_3} = \frac{200}{100} = 2.0 \implies 2^5 = 32.0$$ $$L_{\text{eq}} = 400 + 300(4.214) + 200(32.0) = 400 + 1264.2 + 6400 = 8064.2\text{ m}$$

The compound pipeline is equivalent to a single $200\text{ mm}$ diameter pipe of length $8064.2\text{ m}$.