No history yet

Reynolds Equation Derivation

From Oceans to Bearings

The Navier-Stokes equations are the master equations of fluid dynamics, describing everything from ocean currents to airflow over a wing. But using them in their full form to analyze the lubricant in a tiny gap between machine parts is like using a sledgehammer to crack a nut. The equations are notoriously complex. To understand lubrication, we need a simplified model tailored to the unique environment of a thin fluid film.

This is where the Reynolds equation comes in. It's a special version of the Navier-Stokes equations, adapted for situations where a fluid is squeezed into a very thin layer. To get there, we start with the full equations and make a few key assumptions based on the geometry of a typical bearing.

First, because the film thickness, hh, is tiny compared to its length and width, we assume fluid properties don't change across this thin gap. This means pressure is constant in the direction perpendicular to the surface. Any pressure variation happens along the direction of flow, not through the film's thickness. Second, we assume the fluid flow is slow and orderly (laminar), so inertial forces—the tendency of the fluid to resist changes in velocity—are negligible compared to viscous forces. Essentially, the fluid's stickiness dominates its behavior.

Simplifying the Math

With these assumptions, the mighty Navier-Stokes equations for an incompressible fluid simplify dramatically. We start with the x-momentum equation:

ρ(uux+vuy+wuz)=px+μ(2ux2+2uy2+2uz2)\rho \left( u \frac{\partial u}{\partial x} + v \frac{\partial u}{\partial y} + w \frac{\partial u}{\partial z} \right) = -\frac{\partial p}{\partial x} + \mu \left( \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} + \frac{\partial^2 u}{\partial z^2} \right)

Our 'thin film' assumptions allow us to discard most of these terms. The inertial terms on the left are zero because we assume they are negligible. On the right, the velocity gradients along the flow directions (xx and zz) are much smaller than the gradient across the thin film (yy). This leaves a much simpler relationship:

px=μ2uy2\frac{\partial p}{\partial x} = \mu \frac{\partial^2 u}{\partial y^2}

This is a differential equation we can actually solve. Integrating twice with respect to yy (the direction across the film) gives us the velocity profile of the fluid. We use boundary conditions to find the constants of integration: the fluid velocity is zero at the stationary surface and matches the velocity UU of the moving surface. The result is a velocity profile that has two parts: a linear part from the shearing motion (Couette flow) and a parabolic part driven by the pressure gradient (Poiseuille flow).

Generating Pressure

Knowing the velocity profile is key, but it doesn't tell us how pressure is generated. For that, we need to invoke the principle of mass conservation, expressed by the continuity equation. This equation states that for an incompressible fluid, the amount of fluid flowing into a small volume must equal the amount flowing out.

By substituting our new velocity profile into the continuity equation and integrating across the film thickness, we arrive at the one-dimensional Reynolds equation.

x(h312μpx)=U2hx\frac{\partial}{\partial x} \left( \frac{h^3}{12\mu} \frac{\partial p}{\partial x} \right) = \frac{U}{2} \frac{\partial h}{\partial x}

This equation is the core of hydrodynamic lubrication. It shows that for a pressure gradient to exist, the film thickness hh must vary with xx. In other words, you need a converging wedge. If the surfaces are parallel ( rac{\partial h}{\partial x} = 0), the right side is zero, and no pressure is generated. The fluid is simply dragged along. But with a wedge, the fluid is forced into a narrower space, creating the pressure that lifts the surfaces apart.

For realistic bearings that are finite in width, we must also consider flow in the z-direction (side leakage). This expands the equation into its two-dimensional form:

x(h312μpx)+z(h312μpz)=U2hx\frac{\partial}{\partial x} \left( \frac{h^3}{12\mu} \frac{\partial p}{\partial x} \right) + \frac{\partial}{\partial z} \left( \frac{h^3}{12\mu} \frac{\partial p}{\partial z} \right) = \frac{U}{2} \frac{\partial h}{\partial x}

To solve this equation for a specific bearing, we apply boundary conditions. Typically, the pressure at the edges of the bearing is assumed to be ambient pressure (gauge pressure of zero). Solving the equation with these conditions gives the complete pressure distribution across the bearing surface, which can then be integrated to find the total load the bearing can support.

Let's check your understanding of these foundational principles.

Quiz Questions 1/5

What is the primary purpose of the Reynolds equation in fluid dynamics?

Quiz Questions 2/5

Which of the following is NOT a key assumption used to simplify the Navier-Stokes equations into the Reynolds equation?