Fundamentals of Stream Function and Velocity Potential Function

When analyzing complex fluid dynamics, dealing directly with vector velocity fields can quickly become mathematically cumbersome. To streamline these calculations, fluid dynamicists often turn to scalar formulations. For two-dimensional, incompressible flows, the Stream Function ($\psi$) and the Velocity Potential Function ($\phi$) serve as two of the most powerful mathematical tools in classical mechanics.

The stream function is specifically engineered to satisfy the continuity equation for incompressible flows without requiring simultaneous tracking of both velocity components.

1. Mathematical Formulation

Consider a standard two-dimensional incompressible flow, where the continuity equation is expressed as:
$$\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} = 0$$

To identically satisfy this conservation law, we introduce a scalar function $\psi(x, y)$ such that the velocity components $u$ and $v$ are defined as:
$$u = \frac{\partial \psi}{\partial y}, \quad v = -\frac{\partial \psi}{\partial x}$$

By substituting these definitions into the continuity equation, the expression holds true regardless of the specific form of $\psi$.

2. Physical Significance

• Streamlines: Lines of constant stream function ($\psi = \text{constant}$) correspond directly to streamlines. Fluid particles travel along these contours, and the velocity vector $\mathbf{u}$ is always tangent to them.
• Volumetric Flow Rate: A major practical advantage of the stream function is how it quantifies flow rate. The volume flow rate $Q$ per unit thickness passing between any two streamlines, $\psi_1$ and $\psi_2$, is simply given by their numerical difference:
  $$Q = \psi_2 - \psi_1$$

The Velocity Potential Function ($\phi$)

While the stream function focuses on mass conservation, the velocity potential function is rooted in the kinematic property of rotationality.

1. Mathematical Formulation

If a flow field is entirely irrotational—meaning its vorticity or curl is zero ($\nabla \times \mathbf{u} = 0$)—vector calculus dictates that the velocity field can be represented as the gradient of a single scalar function:
$$\mathbf{u} = \nabla \phi \implies u = \frac{\partial \phi}{\partial x}, \quad v = \frac{\partial \phi}{\partial y}$$

Here, $\phi$ is designated as the velocity potential.

2. Conditions and Applicability

The velocity potential does not exist for every fluid motion; it requires the flow to be irrotational. In practical engineering scenarios, high-Reynolds-number flows outside of boundary layers can often be reliably approximated as irrotational regions.

3. Physical Significance

Contours where $\phi = \text{constant}$ are known as equipotential lines. Because the velocity vector is defined as the gradient of $\phi$, the flow velocity is always perpendicular to these equipotential lines at any given point.

The Intersection of Both Worlds: Potential Flow

When a flow field is simultaneously incompressible and irrotational, it is classified as a Potential Flow. Under these conditions, both $\psi$ and $\phi$ coexist and exhibit a profound mathematical coupling.

1. The Cauchy-Riemann Equations

Comparing the velocity definitions for both functions yields:
$$u = \frac{\partial \phi}{\partial x} = \frac{\partial \psi}{\partial y}$$
$$v = \frac{\partial \phi}{\partial y} = -\frac{\partial \psi}{\partial x}$$

These identities match the famous Cauchy-Riemann equations from complex analysis. Consequently, $\phi$ and $\psi$ act as conjugate harmonic functions.

2. Governing Laplace Equations

Substituting these relations back into the flow governing conditions shows that both functions independently satisfy the Laplace equation:
$$\nabla^2 \phi = 0, \quad \nabla^2 \psi = 0$$

This elegant reduction transforms complex fluid mechanics problems into linear boundary-value problems, allowing engineers to use superposition principles (such as combining sources, sinks, and vortices) to model intricate flow patterns.

3. Orthogonal Grids

In a potential flow, streamlines ($\psi = \text{constant}$) and equipotential lines ($\phi = \text{constant}$) intersect at right angles everywhere. This orthogonal grid provides an intuitive, visually powerful method for mapping fluid behavior.

Illustrative Example: Uniform Flow

To see these concepts in action, consider a steady, two-dimensional uniform flow moving parallel to the $x$-axis, where $u = U$ (a constant) and $v = 0$.

1. Deriving the Velocity Potential ($\phi$):
   $$\frac{\partial \phi}{\partial x} = U \implies \phi = Ux + f(y)$$
   $$\frac{\partial \phi}{\partial y} = 0 \implies f'(y) = 0 \implies \phi = Ux + C$$

2. Deriving the Stream Function ($\psi$):
   $$\frac{\partial \psi}{\partial y} = U \implies \psi = Uy + g(x)$$
   $$-\frac{\partial \psi}{\partial x} = 0 \implies g'(x) = 0 \implies \psi = Uy + C$$

3. Key Takeaways:

   • The equipotential lines ($\phi = Ux$) form straight lines perpendicular to the $x$-axis.
   • The streamlines ($\psi = Uy$) form straight lines parallel to the $x$-axis.
   • Both functions satisfy $\nabla^2 = 0$ and maintain a mutually orthogonal geometry.

Summary and Comparison