---
title: "Lubrication theory"
sort_title: "Lubrication theory"
date: 2021-06-03
categories:
- Physics
- Fluid mechanics
- Fluid dynamics
layout: "concept"
---

**Lubricants** are widely used to reduce friction between two moving surfaces.
In fluid mechanics, **lubrication theory**
is the study of fluids that are tightly constrained in one dimension,
especially those in small gaps between moving surfaces.

For simplicity, we limit ourselves to 2D
by assuming that everything is constant along the $$z$$-axis.
Consider a gap of width $$d$$ (along $$y$$) and length $$L$$ (along $$x$$),
where $$d \ll L$$, containing the fluid.
Outside the gap, the lubricant has a
[Reynolds number](/know/concept/reynolds-number/) $$\mathrm{Re} \approx U L / \nu$$.

Inside the gap, the Reynolds number $$\mathrm{Re}_\mathrm{gap}$$ is different.
This is because advection will dominate along the $$x$$-axis (gap length),
and viscosity along the $$y$$-axis (gap width).
Therefore:

$$\begin{aligned}
    \mathrm{Re}_\mathrm{gap}
    \approx \frac{|(\va{v} \cdot \nabla) \va{v}|}{|\nu \nabla^2 \va{v}|}
    \approx \frac{U^2 / L}{\nu U / d^2}
    \approx \frac{d^2}{L^2} \mathrm{Re}
\end{aligned}$$

If $$d$$ is small enough compared to $$L$$, then $$\mathrm{Re}_\mathrm{gap} \ll 1$$.
More formally, we need $$d \ll L / \sqrt{\mathrm{Re}}$$,
so we are inside the boundary layer,
in the realm of the [Prandtl equations](/know/concept/prandtl-equations/).

Let $$\mathrm{Re}_\mathrm{gap} \ll 1$$.
We are then dealing with *Stokes flow*, in which case
the [Navier-Stokes equations](/know/concept/navier-stokes/equations/)
can be reduced to the following *Stokes equations*:

$$\begin{aligned}
    \pdv{p}{x}
    = \eta \: \Big( \pdvn{2}{v_x}{x} + \pdvn{2}{v_x}{y} \Big)
    \qquad \qquad
    \pdv{p}{y}
    = \eta \: \Big( \pdvn{2}{v_y}{x} + \pdvn{2}{v_y}{y} \Big)
\end{aligned}$$

Let the $$y = 0$$ plane be an infinite flat surface
(a good approximation because $$d \ll L$$),
sliding in the positive $$x$$-direction at a constant velocity $$U$$.
On the other side of the gap, an arbitrary surface
is described by a height function $$h(x)$$.

Since the gap is so narrow,
and the surfaces' movements cause large shear stresses inside it,
$$v_y$$ is negligible compared to $$v_x$$.
Furthermore, because the gap is so long,
we assume that $$\ipdv{v_x}{x}$$ is negligible compared to $$\ipdv{v_x}{y}$$.
This reduces the Stokes equations to:

$$\begin{aligned}
    \pdv{p}{x}
    = \eta \pdvn{2}{v_x}{y}
    \qquad \qquad
    \pdv{p}{y}
    = 0
\end{aligned}$$

This result could also be derived from the Prandtl equations.
In any case, it tells us that $$p$$ only depends on $$x$$,
allowing us to integrate the former equation:

$$\begin{aligned}
    v_x
    = \frac{p'}{2 \eta} y^2 + C_1 y + C_2
\end{aligned}$$

Where $$C_1$$ and $$C_2$$ are integration constants.
At $$y = 0$$, the viscous *no-slip* condition demands that $$v_x = U$$, so $$C_2 = U$$.
Likewise, at $$y = h(x)$$, we need $$v_x = 0$$, leading us to:

$$\begin{aligned}
    v_x
    = \frac{p'}{2 \eta} y^2 - \Big( \frac{p'}{2 \eta} h + \frac{U}{h} \Big) y + U
\end{aligned}$$

The moving bottom surface drags fluid in the $$x$$-direction
at a volumetric rate $$Q(x)$$, given by:

$$\begin{aligned}
    Q(x)
    &= \int_0^{h(x)} v_x(x, y) \dd{y}
    \\
    &= \bigg[ \frac{p'}{6 \eta} y^3 - \frac{p'}{4 \eta} h y^2 - \frac{U}{2 h} y^2 + U y \bigg]_0^{h}
    \\
    &= - \frac{p'}{12 \eta} h^3(x) + \frac{U}{2} h(x)
\end{aligned}$$

Let us assume that the lubricant is incompressible,
meaning that the same volume of fluid must be both leaving and entering the gap.
In that case, $$Q$$ is independent of $$x$$,
which allows us to write $$p'(x)$$ in terms of
measurable constants and the known function $$h(x)$$:

$$\begin{aligned}
    \boxed{
        p'
        = 6 \eta \: \Big( \frac{U}{h^2} - \frac{2 Q}{h^3} \Big)
    }
\end{aligned}$$

Then we insert this into our earlier expression for $$v_x$$, yielding:

$$\begin{aligned}
    v_x
    &= 3 y (y - h) \Big( \frac{U}{h^2} - \frac{2 Q}{h^3} \Big) - \frac{U h}{h^2} y + \frac{U h^2}{h^2}
\end{aligned}$$

Which, after some rearranging, can be written in the following form:

$$\begin{aligned}
    \boxed{
        v_x
        = U \frac{(3 y - h) (y - h)}{h^2} - Q \frac{6 y (y - h)}{h^3}
    }
\end{aligned}$$

With this, we can find $$v_y$$ by exploiting incompressibility,
i.e. the continuity equation states:

$$\begin{aligned}
    \pdv{v_y}{y}
    = - \pdv{v_x}{x}
    = - 2 h' \frac{U h - 3 Q}{h^4} \big( 2 h y - 3 y^2 \big)
\end{aligned}$$

Integrating with respect to $$y$$ therefore leads to
the following transverse velocity $$v_y$$:

$$\begin{aligned}
    \boxed{
        v_y
        = - 2 h' \frac{U h - 3 Q}{h^4} y^2 (h - y)
    }
\end{aligned}$$

Usually, the lubricant is not getting pumped through the system.
In that case, although the pressure gradient $$p'$$ need not be zero in all points
(i.e. there may be complex dynamics inside the gap),
we do expect that its integral across the gap vanishes
(because both sides are at the same pressure):

$$\begin{aligned}
    0
    = \int_L p'(x) \dd{x}
    = 6 \eta U \int_L \frac{1}{h(x)^2} \dd{x} - 12 \eta Q \int_L \frac{1}{h(x)^3} \dd{x}
\end{aligned}$$

Isolating this for $$Q$$, and defining $$q$$ as below, yields a simple equation:

$$\begin{aligned}
    Q
    = \frac{1}{2} U q
    \qquad \qquad
    q
    \equiv \frac{\int_L h^{-2} \dd{x}}{\int_L h^{-3} \dd{x}}
\end{aligned}$$

We substitute this into $$v_x$$ and rearrange to get an interesting expression:

$$\begin{aligned}
    v_x
    &= U \frac{3 y^2 - h y - 3 h y + h^2}{h^2} - U q \frac{3 y^2 - 3 h y}{h^3}
    \\
    &= U \Big( 1 - \frac{y}{h} \Big) \Big( 1 - \frac{3 y (h - q)}{h^2} \Big)
\end{aligned}$$

The first factors are always positive, but the last one can be negative,
if any $$y$$-values satisfy:

$$\begin{aligned}
    h^2 < 3 y (h - q)
    \qquad \implies \qquad
    y > \frac{h^2}{3 (h - q)}
\end{aligned}$$

Since $$h \le y$$, such $$y$$-values will only exist
if $$h$$ is larger than some threshold:

$$\begin{aligned}
    3 (h - q) > h
    \qquad \implies \qquad
    h > \frac{3}{2} q
\end{aligned}$$

If this condition is satisfied, there will be some flow reversal:
rather than just getting dragged by the shearing motion,
the lubricant instead "rolls" inside the gap.
This is confirmed by $$v_y$$:

$$\begin{aligned}
    v_y
    = - U h' \frac{2 h - 3 q}{h^4} y^2 (h - y)
\end{aligned}$$



## References
1.  B. Lautrup,
    *Physics of continuous matter: exotic and everyday phenomena in the macroscopic world*, 2nd edition,
    CRC Press.