Dimensional Analysis and Introduction to Boundary Layers

TL;DR

You'll learn to nondimensionalize the Navier-Stokes equations and see the Reynolds number fall out naturally as the ratio of inertial to viscous forces. You'll then use a scaling argument — not a full solve — to show why viscous effects concentrate in a thin layer near a wall when Re is large. By the end you can estimate boundary layer thickness for any flow just from U, L, and viscosity.

1. The Mental Model

Far from a wall, fluid barely "feels" viscosity — inertia dominates and the flow behaves almost like an ideal fluid. Right at the wall, the no-slip condition forces velocity to zero, and viscosity must matter no matter how small it is. Somewhere in between, at high Reynolds number, viscous effects get squeezed into a razor-thin region. A boundary layer exists because viscosity always wins right at the wall, even when it loses everywhere else.

2. The Core Material

2.1 Nondimensionalizing Navier-Stokes: where Re comes from

Start with the incompressible Navier-Stokes momentum equation, dimensional form:

$$\rho\left(\frac{\partial \mathbf{u}}{\partial t} + \mathbf{u}\cdot abla \mathbf{u}\right) = - abla p + \mu abla^2 \mathbf{u}$$

Pick a characteristic velocity $U$ and a characteristic length $L$ — say, the freestream speed and the plate length in a flow-over-a-plate problem. Define dimensionless variables using these scales:

$$\mathbf{u}^* = \frac{\mathbf{u}}{U}, \quad \mathbf{x}^* = \frac{\mathbf{x}}{L}, \quad t^* = \frac{tU}{L}, \quad p^* = \frac{p}{\rho U^2}$$

Every term in the equation now needs rescaling. The time derivative becomes $\frac{U^2}{L}\frac{\partial \mathbf{u}^*}{\partial t^*}$, the convective term becomes $\frac{U^2}{L}\mathbf{u}^\cdot
abla^\mathbf{u}^$, the pressure gradient becomes $\frac{U^2}{L}
abla^ p^$, and the viscous term becomes $\frac{\mu U}{L^2}
abla^{2}\mathbf{u}^*$. Dividing the whole equation by $\rho U^2/L$ gives:

$$\frac{\partial \mathbf{u}^*}{\partial t^*} + \mathbf{u}^*\cdot abla^* \mathbf{u}^* = - abla^* p^* + \frac{\mu}{\rho U L} abla^{*2} \mathbf{u}^*$$

That coefficient $\frac{\mu}{\rho U L} = \frac{
u}{UL} = \frac{1}{Re}$ is the only free parameter left. Every incompressible flow with the same $Re = UL/
u$ and the same geometry behaves identically once you strip away units — this is why wind tunnel models work.

$$Re = \frac{UL}{ u} = \frac{\text{inertial forces}}{\text{viscous forces}}$$

When $Re$ is large, the equation reads $\frac{\partial \mathbf{u}^}{\partial t^} + \mathbf{u}^\cdot
abla^ \mathbf{u}^ \approx -
abla^ p^*$ — the Euler equation, no viscosity at all. That's a great approximation almost everywhere in the flow. Almost.

2.2 The catch at the wall: Prandtl's scaling argument

Drop viscosity entirely and you get slip at the wall — the inviscid solution allows fluid to slide right past a solid surface. Real fluids don't do that; the no-slip condition demands $u=0$ at $y=0$. So however small $
u$ is, there must be some region near the wall where viscous forces climb back up to match inertia. Prandtl's 1904 insight was to figure out how thick that region has to be, using nothing but order-of-magnitude balancing — no need to solve anything yet.

Set up coordinates for flow over a flat plate: $x$ along the plate (length scale $L$), $y$ normal to it (length scale $\delta$, the unknown boundary layer thickness, with $\delta \ll L$). The streamwise velocity $u$ scales with $U$. From continuity,

$$\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} = 0 \implies \frac{U}{L} \sim \frac{v}{\delta} \implies v \sim \frac{U\delta}{L}$$

So the wall-normal velocity is small — consistent with a thin layer. Now look at the x-momentum balance and compare the inertial term to the viscous term:

$$u\frac{\partial u}{\partial x} \sim \frac{U^2}{L}, \qquad u\frac{\partial^2 u}{\partial y^2} \sim \frac{ u U}{\delta^2}$$

Inside the boundary layer, these two must be the same order of magnitude — that's the whole definition of the boundary layer, the place where viscosity has caught back up to inertia. Setting them equal:

$$\frac{U^2}{L} \sim \frac{ u U}{\delta^2} \implies \delta^2 \sim \frac{ u L}{U} \implies \boxed{\delta \sim \frac{L}{\sqrt{Re_L}}}$$

This single result is the reason boundary layer theory exists: thickness shrinks as $1/\sqrt{Re}$, so at high Reynolds number the viscous region really is thin, and the two-region picture (nearly inviscid outer flow + thin viscous inner layer) is justified.

2.3 The Prandtl boundary layer equations

Given $\delta/L \sim 1/\sqrt{Re} \ll 1$, you can now go back to the full Navier-Stokes equations and throw away every term that's small compared to the dominant ones inside the layer. The y-momentum equation reduces to $\partial p/\partial y \approx 0$ — pressure barely varies across the thin layer, so it's set by the outer inviscid flow and just "imposed" on the boundary layer. The x-momentum equation keeps only the terms that balance at leading order:

$$u\frac{\partial u}{\partial x} + v\frac{\partial u}{\partial y} = -\frac{1}{\rho}\frac{dp}{dx} + u \frac{\partial^2 u}{\partial y^2}$$

$$\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} = 0$$

Notice what's gone: $\partial^2 u/\partial x^2$ (streamwise diffusion) is negligible compared to $\partial^2 u/\partial y^2$ because gradients across the thin layer are far steeper than gradients along it. This is a huge simplification — you've gone from an elliptic system to a parabolic one, which can be marched forward in $x$ like a time-stepping problem.

```mermaid
flowchart TD
A["Full Navier-Stokes equations"] --> B["Nondimensionalize: introduce Re = UL/nu"]
B --> C["High Re limit: outer flow ~ inviscid (Euler)"]
B --> D["No-slip at wall forces viscous region"]
D --> E["Scaling argument: delta ~ L / sqrt