Velocity Potential and Stream Function

In fluid mechanics, describing complex flow patterns can be daunting. However, for a special class of ideal flows, we can replace complicated vector velocity fields with much simpler scalar functions. Mastering the velocity potential and the stream function provides you with powerful analytical tools to model everything from airflow over a wing to water flow around a bridge pier, transforming a vector problem into a more manageable scalar one.

From Vector Fields to Scalars: Core Definitions

To begin, we must understand the specific flow conditions that allow these scalar functions to exist. The velocity potential, denoted by $\varphi$, is defined for irrotational flows. A flow is irrotational when the vorticity, or local rotation of fluid elements, is zero everywhere. Mathematically, this means $\nabla \times \vec{V}=0$, where $\vec{V}$ is the velocity vector. This condition is crucial because a fundamental vector identity states that the curl of the gradient of any scalar function is always zero. Therefore, we can define the velocity potential such that the velocity field is its gradient: $$\vec{V}=\nabla \varphi$$ This single scalar function $\varphi (x,y,z,t)$ contains all the information about the three-dimensional velocity field. A major consequence is that for an incompressible fluid (where density is constant), substituting $\vec{V}=\nabla \varphi$ into the conservation of mass equation, $\nabla \cdot \vec{V}=0$, leads to Laplace's equation: $${\nabla }^2\varphi =0$$ Thus, finding the velocity potential for an incompressible, irrotational flow reduces to solving Laplace's equation, a well-studied linear partial differential equation.

Conversely, the stream function, traditionally denoted by $\psi$, is primarily defined for two-dimensional, incompressible flows. The requirement for a two-dimensional plane (like the xy-plane) is key. For such a flow, the conservation of mass equation is $\partial u/\partial x+\partial v/\partial y=0$. This condition allows us to define a function $\psi (x,y,t)$ such that the velocity components are given by its partial derivatives: $$u=\frac{\partial \psi }{\partial y},v=-\frac{\partial \psi }{\partial x}$$ You can verify that this definition automatically satisfies the two-dimensional incompressible continuity equation. The most powerful feature of the stream function is that lines of constant $\psi$ are streamlines. Along a streamline, the value of the stream function does not change, providing a direct visual map of the flow field.

Orthogonal Families and the Power of Superposition

When a flow is both two-dimensional, incompressible, and irrotational, both the velocity potential $\varphi$ and the stream function $\psi$ exist. This scenario is the heart of potential flow theory. In this case, the lines of constant $\varphi$ (called equipotential lines or potential lines) and the lines of constant $\psi$ (streamlines) form two families of curves that are everywhere orthogonal to each other. This orthogonality arises because the conditions for $\varphi$ and $\psi$ become the Cauchy-Riemann equations from complex analysis: $$\frac{\partial \varphi }{\partial x}=\frac{\partial \psi }{\partial y}=u,\frac{\partial \varphi }{\partial y}=-\frac{\partial \psi }{\partial x}=v$$ These equations confirm that $\varphi$ and $\psi$ are harmonic conjugate functions. This relationship is not just mathematically elegant; it allows us to construct complex flow patterns by simply adding together simpler, known solutions—a principle called superposition. Because Laplace's equation is linear, the sum of two or more solutions is also a solution.

For example, you can combine a uniform flow potential ($\varphi =U_{\infty }x$) with a doublet potential (which represents flow around a circular cylinder). The resulting flow field, described by both a new $\varphi$ and a new $\psi$, accurately models the ideal, frictionless flow around a cylinder. This ability to enable analytical solutions for ideal flow around bodies is the primary engineering application of these functions. By strategically combining basic potentials (uniform flow, source, sink, vortex, doublet), you can model flow over airfoils, through contractions, and around buildings.

Worked Example: Combining a Source and Uniform Flow

Let's see superposition in action. Consider a two-dimensional, incompressible, irrotational flow. We will superimpose a uniform flow from left to right and a line source placed at the origin.

• Uniform Flow: ${\psi }_{uf}=U_{\infty }y$, ${\varphi }_{uf}=U_{\infty }x$.
• Source Flow (strength $m$): ${\psi }_{src}=\frac{m}{2\pi }\theta$, ${\varphi }_{src}=\frac{m}{2\pi }\ln r$, where $r$ and $\theta$ are polar coordinates.

The combined stream function is simply the sum: $${\psi }_{total}=U_{\infty }y+\frac{m}{2\pi }\theta =U_{\infty }r\sin \theta +\frac{m}{2\pi }\theta$$

To find the shape of the body in the flow (a streamline that represents a solid boundary), we look for a stagnation point where velocity is zero and trace the streamline that passes through it. In this combined flow, a stagnation point forms on the x-axis upstream of the source. The streamline ($\psi =\text{constant}$) that passes through this point and encloses the source represents the shape of a semi-infinite body, often called a Rankine half-body. By plotting lines of constant $\psi$ and constant $\varphi$, you would see the orthogonal network of streamlines and potential lines defining this flow geometry.

Common Pitfalls

1. Applying the stream function to 3D flows: A common error is trying to use the classic stream function $\psi$ for general three-dimensional flows. Remember, the $\psi$ defined here is strictly for two-dimensional or axisymmetric flows. For fully 3D flows, other approaches, like a vector potential, are required.
2. Assuming all flows have a velocity potential: The velocity potential $\varphi$ exists only if the flow is irrotational ($\nabla \times \vec{V}=0$). This is a strong constraint. Real flows with significant shear, boundary layers, or wakes are rotational and cannot be described by a velocity potential. Potential flow theory is an idealization that neglects viscosity.
3. Confusing streamlines and pathlines: While constant $\psi$ lines are streamlines (instantaneous snapshots of flow direction), they are not necessarily pathlines (the actual trajectory of a fluid particle over time), except in steady flow. Always ensure the flow is steady if you are using streamlines to trace particle paths.
4. Misinterpreting the physical meaning of $\varphi$: The velocity potential $\varphi$ itself does not have a direct, intuitive physical meaning like pressure or height. Its primary utility is as a mathematical construct whose gradient gives velocity. However, for steady, incompressible, irrotational flow, Bernoulli's equation can be applied between different $\varphi$ contours to find pressure distribution.

Summary

• The velocity potential ($\varphi$) is a scalar function used for irrotational flows, where $\vec{V}=\nabla \varphi$. For incompressible flow, it satisfies Laplace's equation (${\nabla }^2\varphi =0$).
• The stream function ($\psi$) is a scalar function for two-dimensional, incompressible flows, defined such that $u=\partial \psi /\partial y$ and $v=-\partial \psi /\partial x$. Constant values of $\psi$ lie along streamlines.
• When a flow is both 2D, incompressible, and irrotational, both functions exist. Their contours—potential lines and streamlines—form orthogonal families of curves, related by the Cauchy-Riemann equations.
• The linearity of Laplace's equation allows solutions to be built via superposition. This property enables analytical solutions for ideal flow around bodies by adding together basic flow elements like uniform flow, sources, sinks, and vortices.
• These concepts form the foundation of potential flow theory, a powerful analytical tool for modeling frictionless flows, though care must be taken to understand its limitations regarding viscosity and flow rotation.