Editing Stream function

{{Short description|Function for incompressible divergence-free flows in two dimensions}}
{{for|three-dimensional flows with [[axisymmetry]]|Stokes stream function}}
[[File:Stream function.png|thumb|For an incompressible-flow velocity vector field (red, top), its [[Streamlines, streaklines, and pathlines|streamlines]] (dashed) can be computed as the contours of the stream function (bottom).]]
In [[fluid dynamics]], two types of '''stream function''' (or '''streamfunction''') are defined:
* The two-dimensional (or Lagrange) stream function, introduced by [[Joseph Louis Lagrange]] in 1781,<ref>{{citation |first=J.-L. |last=Lagrange |author-link=Joseph-Louis Lagrange |url=http://gallica.bnf.fr/ark:/12148/bpt6k229223s/f697.image |contribution=Mémoire sur la théorie du mouvement des fluides (in: Nouveaux Mémoires de l'Académie Royale des Sciences et Belles-Lettres de Berlin, année 1781) |title=Oevres de Lagrange |year=1868 |volume=Tome IV |pages=695–748 }}</ref> is defined for [[incompressible flow|incompressible]] ([[divergence-free]]), two-dimensional [[fluid flow|flows]].
* The [[Stokes stream function]], named after [[George Gabriel Stokes]],<ref>{{citation | first=G.G. | last=Stokes | authorlink=George Gabriel Stokes | year= 1842 | title= On the steady motion of incompressible fluids | journal= Transactions of the Cambridge Philosophical Society | volume= 7 | pages= 439–453 | bibcode=1848TCaPS...7..439S }}<br>Reprinted in: {{citation | first= G.G. | last=Stokes | year= 1880 | title= Mathematical and Physical Papers, Volume I | publisher= Cambridge University Press | pages= 1–16 | url= https://archive.org/details/mathphyspapers01stokrich }} </ref> is defined for incompressible, three-dimensional flows with [[axisymmetry]].
The properties of stream functions make them useful for analyzing and graphically illustrating flows.

The remainder of this article describes the two-dimensional stream function.

== Two-dimensional stream function ==

=== Assumptions ===
The two-dimensional stream function is based on the following assumptions:
* The flow field can be described as two-dimensional plane flow, with velocity vector
:<math>
\quad \mathbf{u} = \begin{bmatrix}
u (x,y,t) \\
v (x,y,t) \\
0
\end{bmatrix}.
</math>
* The velocity satisfies the [[continuity equation]] for incompressible flow:
:<math>
\quad \nabla \cdot \mathbf{u} = 0.
</math>
* The domain has no holes, or only has holes that have no net flux inwards or outwards.

Although in principle the stream function doesn't require the use of a particular coordinate system, for convenience the description presented here uses a right-handed [[Cartesian coordinate system]] with coordinates <math>(x, y, z)</math>.

=== Derivation ===
==== The test surface ====
Consider two points <math>A</math> and <math>P</math> in the <math>xy</math> plane, and a continuous curve <math>AP</math>, also in the <math>xy</math> plane, that connects them. Then every point on the curve <math>AP</math> has <math>z</math> coordinate <math>z = 0</math>. Let the total length of the curve <math>AP</math> be <math>L</math>.

Suppose a ribbon-shaped surface is created by extending the curve <math>AP</math> upward to the horizontal plane <math>z = b</math> <math>(b>0)</math>, where <math>b</math> is the thickness of the flow. Then the surface has length <math>L</math>, width <math>b</math>, and area <math>b\, L</math>. Call this the ''test surface''.

==== Flux through the test surface ====
[[File:Stream function definition.svg|thumb|right|The volume [[flux]] through the test surface connecting the points <math>A</math> and <math>P.</math>]]
The total [[flux|volumetric flux]] through the test surface is
:<math>
Q (x, y, t) = \int_0^b \int_0^L \mathbf{u} \cdot \hat\mathbf{n}\, \mathrm{d}s\, \mathrm{d}z
</math>
where <math>s</math> is an arc-length parameter defined on the curve <math>AP</math>, with <math>s = 0</math> at the point <math>A</math> and <math>s = L</math> at the point <math>P</math>.
Here <math>\hat\mathbf{n}</math> is the unit vector perpendicular to the test surface, i.e.,
:<math>
\hat\mathbf{n}\, \mathrm{d}s = -R\, \mathrm{d} \mathbf{r} = \begin{bmatrix}
\mathrm{d}y \\
- \mathrm{d}x \\
0
\end{bmatrix}
</math>
where <math>R</math> is the <math>3 \times 3</math> [[rotation matrix]] corresponding to a <math>90^\circ</math> anticlockwise rotation about the positive <math>z</math> axis:
:<math>
R = R_z(90^\circ) = \begin{bmatrix}
0 & -1 & 0 \\  
1 & 0 & 0 \\
0 & 0 & 1 
\end{bmatrix}.
</math>
The integrand in the expression for <math>Q</math> is independent of <math>z</math>, so the outer integral can be evaluated to yield
:<math>
Q (x, y, t) = b\, \int_A^P \left( u\, \mathrm{d} y - v\, \mathrm{d} x \right)
</math>

==== Classical definition ====
[[Horace Lamb|Lamb]] and [[George Batchelor|Batchelor]] define the stream function <math>\psi</math> as follows.<ref>{{harvtxt|Lamb|1932|pages=62–63}} and {{harvtxt|Batchelor|1967|pages=75–79}}</ref>   

:<math>\psi(x,y,t) = \int_A^P \left( u\, \mathrm{d} y - v\, \mathrm{d} x \right)</math>

Using the expression derived above for the total volumetric flux, <math>Q</math>, this can be written as

:<math>\psi(x,y,t) = \frac{Q(x, y, t)}{b}</math>.

In words, the stream function <math>\psi</math> is the volumetric flux through the test surface per unit thickness, where thickness is measured perpendicular to the plane of flow.

The point <math>A</math> is a reference point that defines where the stream function is identically zero. Its position is chosen more or less arbitrarily and, once chosen, typically remains fixed.

An [[infinitesimal]] shift <math>\mathrm{d} P=(\mathrm{d} x,\mathrm{d} y)</math> in the position of point <math>P</math> results in the following change of the stream function:

:<math>\mathrm{d} \psi = u\, \mathrm{d} y - v\, \mathrm{d} x</math>.

From the [[exact differential]] 

:<math>\mathrm{d} \psi = \frac{\partial \psi}{\partial x}\, \mathrm{d} x + \frac{\partial \psi}{\partial y}\, \mathrm{d} y,</math>

so the flow velocity components in relation to the stream function <math>\psi</math> must be

:<math>
u= \frac{\partial \psi}{\partial y},
\qquad
v = -\frac{\partial \psi}{\partial x}.
</math>

Notice that the stream function is [[linear]] in the velocity. Consequently if two incompressible flow fields are superimposed, then the stream function of the resultant flow field is the algebraic sum of the stream functions of the two original fields.

=== Effect of shift in position of reference point ===
Consider a shift in the position of the reference point, say from <math>A</math> to <math>A'</math>. Let <math>\psi '</math> denote the stream function relative to the shifted reference point <math>A'</math>:
:<math>
\psi '(x,y,t) = \int_{A '}^P \left( u\, \mathrm{d} y - v\, \mathrm{d} x \right).
</math>
Then the stream function is shifted by
:<math>\begin{align}
\Delta \psi ( t ) &= \psi '(x,y,t) - \psi (x,y,t) \\
&= \int_{A '}^A \left( u\, \mathrm{d} y - v\, \mathrm{d} x \right),
\end{align}
</math>
which implies the following:
* A shift in the position of the reference point effectively adds a constant (for steady flow) or a function solely of time (for nonsteady flow) to the stream function <math>\psi</math> at every point <math>P</math>.
* The shift in the stream function, <math>\Delta \psi</math>, is equal to the total volumetric flux, per unit thickness, through the continuous surface that extends from point <math>A'</math> to point <math>A</math>. Consequently <math>\Delta \psi = 0</math> if and only if <math>A</math> and <math>A'</math> lie on the same streamline.

=== In terms of vector rotation ===
The velocity <math>\mathbf{u}</math> can be expressed in terms of the stream function <math>\psi</math> as
:<math>
\mathbf{u} = - R \, \nabla \psi
</math>
where <math>R</math> is the <math>3 \times 3</math> rotation matrix corresponding to a <math>90^\circ</math> anticlockwise rotation about the positive <math>z</math> axis. Solving the above equation for <math>\nabla \psi</math> produces the equivalent form
:<math>
\nabla \psi = R\, \mathbf{u}.
</math>
From these forms it is immediately evident that the vectors <math>\mathbf{u}</math> and <math>\nabla \psi</math> are
* perpendicular:  <math> \mathbf{u} \cdot \nabla \psi = 0 </math>
* of the same length:  <math> |\mathbf{u}| = |\nabla \psi| </math>.
Additionally, the compactness of the rotation form facilitates manipulations (e.g., see [[Stream function#Condition of existence|Condition of existence]]).

=== In terms of vector potential and stream surfaces ===
In general, a divergence-free field like <math>\mathbf{u}</math>, also known as a [[solenoidal vector field]], can always be represented as the curl of some [[vector potential]] <math>\boldsymbol{A}</math>:
:<math>
\mathbf{u}= \nabla \times  \boldsymbol{A}.
</math>

The stream function <math>\psi</math> can be understood as providing the strength of a vector potential that is directed perpendicular to the plane:<ref>{{cite book|doi=10.1016/B978-0-12-815489-2.00005-8 |chapter=Viscous Fluid Flow |title=Free-Surface Flow |date=2019 |last1=Katopodes |first1=Nikolaos D. |pages=324–426 |isbn=978-0-12-815489-2 }}</ref>
:<math>
\boldsymbol{A}(x,y,t) = \begin{bmatrix}
0 \\
0 \\
\psi(x,y,t)
\end{bmatrix},
</math>
in other words <math>\boldsymbol{A} = \psi \hat\mathbf{z}</math>, where <math>\hat\mathbf{z}</math> is the unit vector pointing in the positive <math>z</math> direction.

This can also be written as the vector cross product
:<math>
\mathbf{u} = \nabla \psi \times \hat\mathbf{z}
</math>
where we've used the [[vector calculus identity]]
:<math>
\nabla \times \left( \psi \hat\mathbf{z} \right) = \psi \nabla \times \hat\mathbf{z} + \nabla \psi \times \hat\mathbf{z}.
</math>
Noting that <math>\hat\mathbf{z} = \nabla z</math>, and defining <math>\phi = z</math>, one can express the velocity field as
:<math>
\mathbf{u} = \nabla \psi \times \nabla \phi .
</math>
This form shows that the level surfaces of <math>\psi</math> and the level surfaces of <math>z</math> (i.e., horizontal planes) form a system of orthogonal [[Streamsurface|stream surfaces]].

=== Alternative (opposite sign) definition ===
An alternative definition, sometimes used in [[meteorology]] and [[oceanography]], is
:<math>
\psi' = - \psi.
</math>

=== Relation to vorticity ===
{{see also|Vorticity}}
In two-dimensional plane flow, the [[vorticity]] vector, defined as <math> \boldsymbol{\omega} = \nabla \times \mathbf{u}</math>, reduces to <math>\omega \, \hat\mathbf{z}</math>, where
:<math>\omega = - \nabla^2 \psi </math>
or
:<math>\omega = + \nabla^2 \psi' </math>
These are forms of [[Poisson's equation]].

=== Relation to streamlines ===

Consider two-dimensional plane flow with two infinitesimally close points <math>P = (x,y,z)</math> and <math>P ' = (x+dx,y+dy,z)</math> lying in the same horizontal plane. From calculus, the corresponding infinitesimal difference between the values of the stream function at the two points is

:<math>\begin{align}
\mathrm{d} \psi (x, y, t) &= \psi (x + \mathrm{d} x, y + \mathrm{d} y, t) - \psi(x, y, t) \\
&= {\partial \psi \over \partial x} \mathrm{d} x + {\partial \psi \over \partial y} \mathrm{d} y \\
&= \nabla \psi \cdot \mathrm{d} \mathbf{r}
\end{align}
</math>

Suppose <math>\psi</math> takes the same value, say <math>C</math>, at the two points <math>P</math> and <math>P '</math>. Then this gives

:<math>
0 = \nabla \psi \cdot \mathrm{d} \mathbf{r} ,
</math>

implying that the vector <math>\nabla \psi</math> is normal to the surface <math>\psi = C</math>. Because <math>\mathbf{u} \cdot \nabla \psi = 0</math> everywhere (e.g., see [[Stream function#Condition of existence|In terms of vector rotation]]), each streamline corresponds to the intersection of a particular stream surface and a particular horizontal plane. Consequently, in three dimensions, unambiguous identification of any particular streamline requires that one specify corresponding values of both the stream function and the elevation (<math>z</math> coordinate).

The development here assumes the space domain is three-dimensional. The concept of stream function can also be developed in the context of a two-dimensional space domain. In that case level sets of the stream function are curves rather than surfaces, and streamlines are level curves of the stream function. Consequently, in two dimensions, unambiguous identification of any particular streamline requires that one specify the corresponding value of the stream function only.

=== Condition of existence ===

It's straightforward to show that for two-dimensional plane flow <math>\mathbf{u}</math> satisfies the curl-divergence equation
:<math>
(\nabla \cdot \mathbf{u})\, \hat\mathbf{z} = - \nabla \times (R \, \mathbf{u})
</math>
where <math>R</math> is the <math>3 \times 3</math> [[rotation matrix]] corresponding to a <math>90^\circ</math> anticlockwise rotation about the positive <math>z</math> axis. This equation holds regardless of whether or not the flow is incompressible.

If the flow is incompressible (i.e., <math>\nabla \cdot \mathbf{u} = 0</math>), then the curl-divergence equation gives
:<math>\mathbf{0} = \nabla \times (R\, \mathbf{u})</math>.
Then by [[Stokes' theorem]] the line integral of <math>R\, \mathbf{u}</math> over every closed loop vanishes
:<math>
\oint_{\partial\Sigma} (R\, \mathbf{u}) \cdot \mathrm{d}\mathbf{\Gamma}= 0.
</math>
Hence, the line integral of <math>R\, \mathbf{u}</math> is path-independent. Finally, by the converse of the [[gradient theorem]], a scalar function <math>\psi (x,y,t)</math> exists such that
:<math>R\, \mathbf{u} = \nabla \psi</math>.
Here <math>\psi</math> represents the stream function.

Conversely, if the stream function exists, then <math>R\, \mathbf{u} = \nabla \psi</math>. Substituting this result into the curl-divergence equation yields <math>\nabla \cdot \mathbf{u} = 0</math> (i.e., the flow is incompressible).

In summary, the stream function for two-dimensional plane flow exists if and only if the flow is incompressible.

=== Potential flow ===
For two-dimensional [[potential flow]], streamlines are perpendicular to [[equipotential]] lines. Taken together with the [[velocity potential]], the stream function may be used to derive a [[complex number|complex]] potential. In other words, the stream function accounts for the [[solenoidal]] part of a two-dimensional [[Helmholtz decomposition]], while the velocity potential accounts for the [[irrotational]] part.

=== Summary of properties ===
The basic properties of two-dimensional stream functions can be summarized as follows:
# The ''x''- and ''y''-components of the flow velocity at a given point are given by the [[partial derivative]]s of the stream function at that point.
# The value of the stream function is constant along every [[Streamlines, streaklines, and pathlines|streamline]] (streamlines represent the trajectories of particles in steady flow). That is, in two dimensions each streamline is a [[level curve]] of the stream function.
# The difference between the stream function values at any two points gives the volumetric flux through the vertical surface that connects the two points.

== Two-dimensional stream function for flows with time-invariant density==
If the fluid density is time-invariant at all points within the flow, i.e.,

:<math> \frac{\partial \rho}{\partial t} = 0 </math>,

then the [[continuity equation]] (e.g., see [[Continuity equation#Fluid dynamics]]) for two-dimensional plane flow becomes

:<math> \nabla \cdot ( \rho\, \mathbf{u} ) = 0. </math>

In this case the stream function <math>\psi</math> is defined such that

:<math>\rho\, u = \frac{\partial \psi}{\partial y}, \quad \rho\, v = - \frac{\partial \psi}{\partial x}</math>

and represents the mass flux (rather than volumetric flux) per unit thickness through the test surface.

==See also==
*[[Elementary flow]]

== References ==

===Citations===
{{reflist}}

===Sources===
{{refbegin}}
* {{citation |first=G. K. |last=Batchelor |author-link=George Batchelor |title=An Introduction to Fluid Dynamics |year=1967 |publisher=Cambridge University Press |isbn=0-521-09817-3 }}
* {{citation |first=H. |last=Lamb |author-link=Horace Lamb |title=Hydrodynamics |publisher=Cambridge University Press, republished by Dover Publications |year=1932 |edition=6th |isbn=0-486-60256-7 }} 
* {{citation |first1=B. S. |last1=Massey| first2=J. |last2=Ward-Smith |title=Mechanics of Fluids |edition=7th |publisher=Nelson Thornes |location=UK |year=1998 }}
* {{citation |first=F. M. |last=White |title=Fluid Mechanics |edition=5th |publisher=McGraw-Hill |location=New York |year=2003 }} 
* {{citation |first=T. W. |last=Gamelin |title=Complex Analysis |publisher=Springer |location=New York |year=2001 |isbn=0-387-95093-1 }}
* {{citation |publisher=[[American Meteorological Society]] |work=AMS Glossary of Meteorology |url=http://glossary.ametsoc.org/wiki/Streamfunction |title=Streamfunction |accessdate=2014-01-30 }}
{{refend}}

[[Category:Continuum mechanics]]
[[Category:Fluid dynamics]]

==External links==

* [http://prj.dimanov.com/ Joukowsky Transform Interactive WebApp]