Stream function
From Free net encyclopedia
In fluid dynamics, the stream function is defined for two-dimensional flows. The difference between the values of the stream function at any two points gives the volumetric flow rate (or flux) through a line connecting the two points.
Note that since streamlines are tangent to the flow, the value of the stream function must be the same 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 fluid flow.
Contents |
Two Dimensional Stream Function: Overview
The stream function <math>\psi</math> for a two dimensional flow is defined such that the flow velocity can be expressed as:
- <math>
\mathbf{u}= \nabla \times \boldsymbol{\psi} </math>
Where <math>\boldsymbol{\psi} = (0,0,\psi)</math> if the velocity vector <math>\mathbf{u} = (u,v,0)</math> or <math>\boldsymbol{\psi} = (0,\psi,0)</math> if <math>\mathbf{u} = (u,0,w)</math>.
In Cartesian coordinate system this is equivalant to
- <math>
u= - \frac{\partial\psi}{\partial y},\qquad v=\frac{\partial\psi}{\partial x} </math> Where <math>u</math> and <math>v</math> are the velocities in the <math>x</math> and <math>y</math> directions, respectively:
- <math>
u = \frac{\partial x}{\partial t},\qquad v = \frac{\partial y}{\partial t} </math>
This formulation of the stream function satisfies the two dimensional continuity equation:
- <math>
\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} = 0 </math>
Two Dimensional Stream Function: The Derivation
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:
- <math>\delta \psi = q \delta n\,</math> (1.1)
As δn → 0, rearranging this expression, we get:
- <math>q = \frac{\partial \psi}{\partial n}\,</math> (1.2)
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 right to left. A sign convention is adopted such that the velocity of the flow is positive. However, this sign convention is not universal and thus ought be taken with caution. A couple of examples ought clarify this point:
Flow in Cartesian coordinates (x,y)
By observing the flow into an elemental square in an x-y coordinate system, we have:
- <math>\delta \psi = -u \delta y\,</math> (1.3)
- <math>\delta \psi = v \delta x\,</math> (1.4)
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. In this case, we have a negative velocity of -u since the flow crosses the y-axis from left to right. Thus as δn → 0 and by rearranging, we have:
- <math>u = - \frac{\partial \psi}{\partial y}\,</math> (1.5)
- <math>v = \frac{\partial \psi}{\partial x}\,</math> (1.6)
Flow in Polar coordinates (r,θ)
By observing the flow into an elemental square in an r-θ coordinate system (where r expands outward from the origin and θ is positive anticlockwise), we have:
- <math>\delta \psi = -v_r ( r \delta \theta )\,</math> (1.7)
- <math>\delta \psi = v_\theta \delta r\,</math> (1.8)
where vr is the velocity parallel to and in the direction of the r-axis, and vθ is the velocity parallel to and in the direction of the θ-axis. In this case, we have a negative velocity of -vr since the flow crosses the θ-axis from left to right. Thus as δn → 0 and by rearranging, we have:
- <math>v_r = - \frac{1}{r} \frac{\partial \psi}{\partial \theta}\,</math> (1.9)
- <math>v_\theta = \frac{\partial \psi}{\partial r}\,</math> (1.10)
Continuity: The Derivation
Consider two dimensional plane flow within a Cartesian coordinate system. Continuity states that if we consider 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:
- <math>\delta \psi_{in} = u \delta y + v \delta x\,</math> (2.1)
The total flow out of the element is given by:
- <math>\delta \psi_{out} = \left( u + \frac{\partial u}{\partial x}\delta x \right) \delta y + \left( v + \frac{\partial v}{\partial y}\delta y \right) \delta x\,</math> (2.2)
Thus we have:
- <math>\delta \psi_{in} = \delta \psi_{out}\,</math> (2.3)
- <math> u \delta y + v \delta x\ = \left( u + \frac{\partial u}{\partial x}\delta x \right) \delta y + \left( v + \frac{\partial v}{\partial y}\delta y \right) \delta x\,</math> (2.4)
and simplifying to:
- <math>\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} = 0</math> (2.5)
(NB: Whilst we have ignored the earlier sign convention for Ψ, the mathematics ought work out fine but this proof is perhaps slightly tidier and easier to follow). Substituting the expressions of the stream function into this equation, we have:
- <math>\frac{\partial^2 \psi}{\partial x \partial y} - \frac{\partial^2 \psi}{\partial y \partial x} = 0</math> (2.6)
Vorticity
It can be proven that the stream function equation for vorticity satisifies Laplace's equation:
- <math>\frac{\partial^2 \psi}{\partial x^2} + \frac{\partial^2 \psi}{\partial y^2} = 0</math> (3.1)
See also
References
- B. Massey, "Mechanics of Fluids (Seventh Edition)"
- M. Oldfield, "Potential Flow", Lecture Notes (2005), University of Oxford