Stream function

The stream function is defined for two-dimensional flows of various kinds. The stream function can be used to plot streamlines, which represent the trajectories of particles in a steady flow. Streamlines are perpendicular to equipotential lines. In most cases, the stream function is the imaginary part of the complex potential, while the potential function is the real part.

Considering the particular case of fluid dynamics, the difference between the stream function values at any two points gives the volumetric flow rate (or flux) through a line connecting the two points.

Since streamlines are tangent to the velocity vector of the flow, the value of the stream function must be constant along a streamline. If there were a flux across a line, it would necessarily not be tangent to the flow, hence would not be a streamline.

The usefulness of the stream function lies in the fact that the velocity components in the "x"- and "y"- directions at a given point are given by the partial derivatives of the stream function at that point. A stream function may be defined for any flow of dimensions greater than two, however the two dimensional case is generally the easiest to visualize and derive.

Taken together with the velocity potential, the stream function may be used to derive a complex potential for a potential flow.

Two dimensional stream function

The stream function $psi$ for a two dimensional flow is defined such that the flow velocity can be expressed as::

Where if the velocity vector $mathbf\{u\}\; =\; (u,v,0)$.

In Cartesian coordinate system this is equivalent to:$u=\; frac\{partialpsi\}\{partial\; y\},qquadv=\; -frac\{partialpsi\}\{partial\; x\}$ Where $u$ and $v$ are the velocities in the $x$ and $y$ coordinate directions, respectively.

This formulation of the stream function satisfies the two dimensional continuity equation::$frac\{partial\; u\}\{partial\; x\}\; +\; frac\{partial\; v\}\{partial\; y\}\; =\; 0$

Derivation of the two dimensional stream function

Consider two points A and B in two dimensional plane flow. If the distance between these two points is very small: δn, and a stream of flow passes between these points with an average velocity, q perpendicular to the line AB, the volume flow rate per unit thickness, δΨ is given by::$delta\; psi\; =\; q\; delta\; n,$

As δn → 0, rearranging this expression, we get::$q\; =\; frac\{partial\; psi\}\{partial\; n\},$

Now consider two dimensional plane flow with reference to a coordinate system. Suppose an observer looks along an arbitrary axis in the direction of increase and sees flow crossing the axis from "left to right". A sign convention is adopted such that the velocity of the flow is "positive".

Flow in Cartesian coordinates

By observing the flow into an elemental square in an x-y Cartesian coordinate system, we have::$delta\; psi\; =\; u\; delta\; y,$:$delta\; psi\; =\; -v\; delta\; x,$

where u is the velocity parallel to and in the direction of the x-axis, and v is the velocity parallel to and in the direction of the y-axis. Thus, as δn → 0 and by rearranging, we have::$u\; =\; frac\{partial\; psi\}\{partial\; y\},$:$v\; =\; -\; frac\{partial\; psi\}\{partial\; x\},$

Flow in Polar coordinates

By observing the flow into an elemental square in an r-θ Polar coordinate system, we have::$delta\; psi\; =\; v\_r\; (\; r\; delta\; heta\; ),$:$delta\; psi\; =\; -v\_\; heta\; delta\; r,$

where v_r is the radial velocity component (parallel to the r-axis), and v_θ is the tangential velocity component (parallel to the θ-axis). Thus, as δn → 0 and by rearranging, we have::$v\_r\; =\; frac\{1\}\{r\}\; frac\{partial\; psi\}\{partial\; heta\},$:$v\_\; heta\; =\; -frac\{partial\; psi\}\{partial\; r\},$

Continuity: The Derivation

Consider two dimensional plane flow within a Cartesian coordinate system. Continuity states that if we consider incompressible flow into an elemental square, the flow into that small element must equal the flow out of that element.

The total flow into the element is given by::$delta\; psi\_\{in\}\; =\; u\; delta\; y\; +\; v\; delta\; x.,$

The total flow out of the element is given by::$delta\; psi\_\{out\}\; =\; left(\; u\; +\; frac\{partial\; u\}\{partial\; x\}delta\; x\; ight)\; delta\; y\; +\; left(\; v\; +\; frac\{partial\; v\}\{partial\; y\}delta\; y\; ight)\; delta\; x.,$

Thus we have::$delta\; psi\_\{in\}\; =\; delta\; psi\_\{out\},$:$u\; delta\; y\; +\; v\; delta\; x\; =\; left(\; u\; +\; frac\{partial\; u\}\{partial\; x\}delta\; x\; ight)\; delta\; y\; +\; left(\; v\; +\; frac\{partial\; v\}\{partial\; y\}delta\; y\; ight)\; delta\; x,$

and simplifying to::$frac\{partial\; u\}\{partial\; x\}\; +\; frac\{partial\; v\}\{partial\; y\}\; =\; 0.$

Substituting the expressions of the stream function into this equation, we have:

:$frac\{partial^2\; psi\}\{partial\; x\; partial\; y\}\; -\; frac\{partial^2\; psi\}\{partial\; y\; partial\; x\}\; =\; 0.$

Vorticity

In Cartesian coordinates, the stream function can be found from vorticity using the following Poisson's equation::$abla\; ^2\; psi\; =\; -omega$

where $vec\; omega\; =\; (\; 0,\; 0,\; omega\; )$ and $vec\; omega\; =\; vec\; abla\; imes\; vec\; v\; .$

See also

*Potential flow

References

* B. S. Massey and J. Ward-Smith, "Mechanics of Fluids", 7^th ed., Nelson Thornes, UK (1998).
* F. M. White, "Fluid Mechanics", 5^th ed., McGraw-Hill, New York (2003).
* T. W. Gamelin, "Complex Analysis", Springer, New York (2001). ISBN 0-387-95093-1.