---
title: Darcy–Weisbach Friction Factor
url: https://www.emergentmind.com/topics/darcy-weisbach-friction-factor
type: topic
---

# Darcy–Weisbach Friction Factor

The Darcy–Weisbach friction factor, denoted \( f \), is a central dimensionless parameter quantifying energy loss due to wall shear in internal and open-channel flows, linking the pressure drop or head loss to the mean velocity, geometry, and, critically, fluid–wall interaction. It forms the basis for the Darcy–Weisbach equation, which is universally adopted in both analytical models and simulation frameworks across hydraulic, environmental, and industrial fluid mechanics.

## 1. Fundamental Definitions and Canonical Formulations

The classical Darcy–Weisbach equation connects the head loss per unit length or total pressure drop in a conduit to the mean flow parameters via
\[
h_f = f\,\frac{L}{D}\,\frac{V^2}{2g}, \quad \text{or equivalently,} \quad \Delta p = f\,\frac{L}{D}\,\frac{\rho V^2}{2}
\]
where \( f \) is the Darcy–Weisbach friction factor, \( L \) is the pipe length, \( D \) diameter, \( V \) mean velocity, and \( g \) gravitational acceleration. For general geometries and flow regimes, modifications appear (wetted perimeter, hydraulic radius, effective velocities), but the core definition remains.

The friction factor \( f \) encapsulates all wall-shear effects and depends on Reynolds number \( Re \), relative roughness \( \epsilon/D \), and, for non-Newtonian fluids, the rheological response. In laminar (fully developed) flow of Newtonian fluids in a circular pipe, \( f \) adopts the exact form \( f = 64/Re \), where \( Re = VD/\nu \), and for turbulent regimes or open channels, it is determined empirically or via semi-theoretical correlations [1612.08540, 2110.14966, 1609.04711].


## 2. Regime-Specific Correlations: Laminar, Transitional, and Turbulent Laws

Explicit expressions and regime boundaries for \( f \) depend on flow state, wall roughness, and geometry. For Newtonian fluids in pipes:
- **Laminar (Re \( < \sim\!2000 \))**:  \( f = 64/Re \) (Hagen–Poiseuille law). This holds for both smooth and rough pipes if \( Re \) remains below the transition, with roughness effects negligible due to the dominance of viscous sublayers [1905.12479, 1912.13388].
- **Turbulent (Re \( > \sim\!3000 \))**, smooth/rough:
    - *Blasius law* for hydraulically smooth pipes: \( f = 0.3164\,Re^{-1/4} \) (\( 3\times10^3 \lesssim Re \lesssim 10^5 \)) [1612.08540].
    - *General turbulent regime*: Implicit relations, e.g., Colebrook–White,
      \[
    \frac{1}{\sqrt{f}} = -2.0\, \log_{10}\left(\frac{\epsilon/D}{3.7} + \frac{2.51}{Re\sqrt{f}}\right)
      \]
      or explicit forms such as the Haaland approximation [2110.14966]. Roughness-limited flows (\( Re\rightarrow\infty \)) asymptote to \( f \sim (\epsilon/D)^{1/3} \) (“Strickler regime”).

- **Transitional regime** (\( 1.5\times 10^3 < Re < 2.7\times10^3 \)): Laboratory measurements reveal that \( f \) does not interpolate linearly between laminar and turbulent values, but follows reproducible, protocol-dependent curves depending on driving conditions (see Table 1 below) [2110.14966].

| \( Re \) | \( f_{\text{CPD}} \) (gravity-driven) | \( f_{\text{CMF}} \) (mass-flux) |
|---------|--------------------|------------------|
| 1600    | 0.0530             | 0.0660           |
| 2000    | 0.0380             | 0.0510           |
| 2400    | 0.0310             | 0.0440           |
| 2700    | 0.0280             | —                |

The appropriate regime and correlation must be selected based on the specific \( Re \), relative roughness, and flow control protocol [2110.14966].

## 3. Dependency on Surface Roughness and Extension to Arbitrary Geometries

Wall roughness effects are incorporated via \( \epsilon/D \) or equivalent sand-grain roughness \( k_s/D \), modifying \( f \) and shifting regime boundaries:
- In turbulent pipe flow, measured friction factors for rough pipes deviate systematically from smooth curves as \( Re \) increases, saturating at a roughness-controlled asymptote [1612.08540, 2602.17082].
- In laminar regimes, the scaling \( f \propto 1/Re \) typically persists even with periodic/bar roughness, but the prefactor increases quadratically with roughness height and decreases linearly with pitch [1905.12479]:
    \[
    f = \frac{64 + 520(k/D) + 3000(k/D)^2 - 35(p/D)}{Re}
    \]
    Valid for pipes with periodic square-bar roughness, 200 \(\leq\) Re \(\leq\) 2000.

For open-channel and shallow flows (e.g., rainfall-driven overland flow), analogs of the Darcy–Weisbach law are derived, with the friction factor adapted to local Reynolds numbers and flow depth [1609.04711].

## 4. Extensions: Non-Newtonian Fluids and Multiphase Systems

The Darcy–Weisbach friction factor is explicitly generalized for non-Newtonian, non-single phase, and viscoelastic flows:
- **Non-Newtonian turbulent pipe flow** is governed by a generalized Reynolds number (\( Re_g \)) and wall-consistency index, yielding an implicit log-law [1009.0968]:
    \[
    \frac{1}{\sqrt{f}} = 4.06\log_{10}(Re_g \sqrt{f}) + 2.10 - \frac{2.81}{n'}
    \]
    which must be solved iteratively for \( f \) for given rheological parameters.
- **Laminar viscoelastic fluids** (e.g., FENE-P model): The friction factor is strictly bounded by the values for pure solvent and total zero-shear viscosity
    \[
    \frac{64\,\beta^2}{Re} \leq f_{VE} \leq \frac{64}{Re}
    \]
    with the reduction below \( 64/Re \) due to strain-rate redistribution (shear-thinning) [1912.13388].
- **Multiphase/aerated flows**: In high-Froude, air–water mixtures, the effective friction factor is decomposed into uniform, spatial, and temporal (evolving) contributions,
    \[
    f_e = f_{\rm uniform} + f_{\rm spatial} + f_{\rm temporal}
    \]
    incorporating Favre (density-weighted) averaging and vertical structure corrections for momentum and pressure terms [2601.05523].

## 5. Modern Computational and Symbolic Approaches to Friction-Factor Correlations

The need for accurate, explicit, and robust friction factor expressions for both routine engineering and simulation code integration has driven recent developments:
- **Continuous crossover formulations** (e.g., Churchill equation) are implemented in modern open-source solvers like openKARST, blending laminar and turbulent regimes seamlessly and removing the need for iterative solution of implicit equations such as Colebrook–White [2503.22320].
- **Data-based symbolic regression**: Physics-informed symbolic regression, constrained by order-of-magnitude momentum and turbulence scaling, yields explicit expressions matching reference data (Nikuradse, Superpipe) to \(\sim\)1% RMS error across \(10^3 \leq Re \leq 10^7\) and all roughness levels. These models encode both viscous and turbulent mechanisms, recover the characteristic inflection (“minimum”) in \( f(Re) \) for rough pipes, and impose physically justified sensitivity exponents for velocity, roughness, viscosity, and density [2602.17082].

Example explicit expression (“Candidate 1”):
\[
f(Re,\epsilon/D) = T_1 + T_2 + T_3 + T_4 + 0.046
\]
where each \( T_i \) is an explicit function of \( Re \) and \( \epsilon/D \), designed to recover the appropriate asymptotic scaling in all regimes [2602.17082].


## 6. Sensitivity, Physical Constraints, and Influence of Protocol

Quantitative sensitivity envelopes—velocity, roughness, viscosity, and density exponents—serve as physical priors for correlation development and validation, ensuring that \( f \) behaves consistently with underlying momentum and energy balances [2602.17082]. In transitional flows, laboratory evidence has established that friction factor curves depend on protocol (constant pressure vs. constant mass flux), initial conditions, and flow history, with transitional \( f \) curves occupying a distinct region in \( f \)–\( Re \) space not interpolated by classical laminar–turbulent laws [2110.14966].

Moreover, in flows with geometric or physical complexity (e.g., helical pipes), rigorous upper bounds on \( f \) can be constructed using variational “background” methods, which account for curvature and torsion effects but may be quantitatively loose relative to direct numerical/experimental estimates at high \( Re \) [2009.01322].


## 7. Practical Applications, Implementation, and Limitations

The Darcy–Weisbach friction factor underpins pressure-loss, head-loss, and flow-capacity computations across hydraulics, pipeline engineering, subsurface fracture propagation, open-channel design, and even transient multi-regime network solvers. Successive generations of correlations—empirical, semi-theoretical, and data-driven—have systematically expanded \( f \)’s applicability, precision, and computational convenience. However, accurate prediction still depends on correct identification of the flow regime, careful characterization of wall properties (roughness, geometry), protocol adaptation (especially in the transitional regime), and suitable consideration of fluid rheology and multiphase effects.

Current best practice is to employ explicit, physics-consistent friction-factor formulas validated on the target flow regime and geometry, and to employ modern continuous or data-driven models where legacy iterative formulas (e.g., Colebrook–White) introduce numerical instability or excess computational cost [2503.22320, 2602.17082].

---
**References:**  
[1609.04711], [1612.08540], [2110.14966], [1905.12479], [1912.13388], [1009.0968], [2601.05523], [2503.22320], [2602.17082], [2009.01322], [1609.01237]

Source: https://www.emergentmind.com/topics/darcy-weisbach-friction-factor