Potential flow

From Wikipedia, the free encyclopedia

A potential flow is characterized by an irrotational velocity field. This property allows the description of the velocity field as the gradient of a scalar function (because taking the curl of the gradient is equivalent to make the cross product of two parallel vectors which is always zero). The flow can be compressible or incompressible, stationnary or unstationnary.

The free-vortex flow is solution of the vorticity transport equation (obtained by taking the curl of the Euler momentum equation) only if the sources terms are zero, that is that the flow is barotropic (the pressure varies only with the density) and that the external forces are negligible or derives from a potential (the force fields are described by taking the gradient of scalar functions).

In incompressible, steady or unsteady, fluid dynamics, potential flow obeys the following equations

\nabla \times \mathbf{v} = 0 (zero rotation)
\nabla \cdot \mathbf{v} = 0 (zero divergence = volume conservation)

Equivalently,

\mathbf{v} = \nabla \Phi \;
        

where:

  • \mathbf{v} is the fluid velocity vector
  • Φ is the fluid flow potential, scalar
  • "\nabla \times" is curl
  • "\nabla \cdot" is divergence.
  • "\nabla" is gradient

Note that "\nabla \cdot \mathbf{v}" is something different from "\nabla \mathbf{v}".

The equations above imply \nabla^2 \Phi=0, or Laplace's equation, holds.

Together with the Navier-Stokes equations or the Euler equations, these equations can be used to calculate solutions to many practical flow situations. In two dimensions, potential flow reduces to a very simple system that is analysed using complex numbers (see below).

Potential flow does not include all the characteristics of flows that are encountered in the real world. For example, potential flow excludes turbulence, which is commonly encountered in nature. Richard Feynman considered potential flow to be so unphysical that the only fluid to obey the assumptions was "dry water".

Potential flow also makes a number of invalid predictions, such as d'Alembert's paradox, which states that the drag on any object moving through an infinite fluid otherwise at rest is zero.

More precisely, potential flow cannot account for the behaviour of flows that include a boundary layer.

Nevertheless, understanding potential flow is important in many branches of fluid mechanics. In particular, simple potential flows (called elemental flows) such as the free vortex and the point source possess ready analytical solutions. These solutions can be superposed to create more complex flows satisfying a variety of boundary conditions. These flows correspond closely to real-life flows over the whole of fluid mechanics; in addition, many valuable insights arise when considering the deviation (often slight) between an observed flow and the corresponding potential flow.

Potential flow finds many applications in fields such as aircraft design. For instance, in computational fluid dynamics, one technique is to couple a potential flow solution outside the boundary layer to a solution of the boundary layer equations inside the boundary layer.

The absence of boundary layer effects means that any streamline can be replaced by a solid boundary with no change in the flow field, a technique used in many aerodynamic design approaches. Another technique would be the use of Riabouchinsky solids.

Contents

[edit] Analysis

Potential flow in two dimensions is simple to analyse using complex numbers, viewed for convenience on the Argand diagram.

The basic idea is to define a holomorphic function f. If we write

f(x + iy) = φ + iψ

then the Cauchy-Riemann equations show that

\frac{\partial\phi}{\partial x}=\frac{\partial\psi}{\partial y}, \qquad \frac{\partial\phi}{\partial y}=-\frac{\partial\psi}{\partial x}.

(it is conventional to regard all symbols as real numbers; and to write z = x + iy and w = φ + iψ).

Also notice that \frac{df}{dz}=u-iv (Some steps are missing, try calculating yourself)

The velocity field \underline{u}=(u,v), specified by

u=\frac{\partial\phi}{\partial x},\qquad v=\frac{\partial\phi}{\partial y}

then satisfies the requirements for potential flow:

\nabla\cdot\underline{u}= \nabla^2\phi= \frac{\partial^2\phi}{\partial x^2}+\frac{\partial^2\phi}{\partial y^2}= {\partial \over \partial x} {\partial \psi \over \partial y} + {\partial \over \partial y} \left( - {\partial \psi \over \partial x} \right) = 0

and

\left|\nabla\times\underline{u}\right|= \frac{\partial v}{\partial x}-\frac{\partial u}{\partial y}= \frac{\partial^2\phi}{\partial x\partial y}- \frac{\partial^2\phi}{\partial y\partial x}=0.

ψ is defined as the stream function. Lines of constant ψ are known as streamlines and lines of constant φ are known as equipotential lines (see equipotential surface).

The two sets of curves intersect at right angles, for

\nabla \phi \cdot \nabla \psi = \frac{\partial\phi}{\partial x}\frac{\partial\psi}{\partial x}+ \frac{\partial\phi}{\partial y}\frac{\partial\psi}{\partial y}= {\partial \psi \over \partial y} {\partial \psi \over \partial x} - {\partial \psi \over \partial x} {\partial \psi \over \partial y} = 0.

[edit] Examples: general considerations

Any differentiable function may be used for f. The examples that follow use a variety of elementary functions; special functions may also be used.

Note that multi-valued functions such as the natural logarithm may be used, but attention must be confined to a single Riemann surface.

[edit] Examples: Power laws

If

w = Azn

then, writing x + iy = reiθ, we have

φ = Arncosnθ

and

ψ = Arnsinnθ

[edit] Power law with n = 1

If w = Az1, that is, a power law with n = 1, the streamlines (ie lines of constant ψ) are a system of straight lines parallel to the x-axis. This is easiest to see by writing in terms of real and imaginary components:

f(x+iy)=A\times(x+iy)=Ax+i\cdot Ay

thus giving φ = Ax and ψ = Ay.

[edit] Power law with n = 2

If n = 2, then w = Az2 and the streamline corresponding to a particular value of ψ are those points satisfying

ψ = Ar2sin2θ

which is a system of rectangular hyperbolae. This may be seen by again rewriting in terms of real and imaginary components. Noting that \sin 2\theta=2\sin\theta\,\cos\theta and rewriting sinθ = y / r and cosθ = x / r it is seen (on simplifying) that the streamlines are given by

ψ = 2Axy.

The velocity field is given by \nabla\phi, or

(u,v)= \left( {\partial \phi \over \partial x}, {\partial \phi \over \partial y} \right) =  \left( {\partial \psi \over \partial y}, - {\partial \psi \over \partial x} \right) = \left(2Ax,-2Ay\right)

In fluid dynamics, the flowfield near the origin corresponds to a stagnation point. Note that the fluid at the origin is at rest (this follows on differentiation of f(z) = z2 at z = 0).

The ψ = 0 streamline is particularly interesting: it has two (or four) branches, following the coordinate axes, ie x = 0 and y = 0.

As no fluid flows across the x-axis, it (the x-axis) may be treated as a solid boundary. It is thus possible to ignore the flow in the lower half-plane where y < 0 and to focus on the flow in the upper half-plane.

With this interpretation, the flow is that of a vertically directed jet impinging on a horizontal flat plate.

The flow may also be interpreted as flow into a 90 degree corner if the regions specified by (say) x < 0 and y < 0 are ignored.

[edit] Power law with n = 3

If n = 3 the resulting flow is a sort of hexagonal version of the n = 2 case considered above. Streamlines are given by 3x2yy3 = ψ.

[edit] Power law with n = − 1

if n=\,-1, the streamlines are given by

\psi=-\frac{A}{r}\sin\theta.

This is more easily interpreted in terms of real and imaginary components:

\psi = {-A y \over r^2} = {-A y \over x^2 + y^2},
x^2 + y^2 + {A y \over \psi} = 0,
x^2+\left(y+\frac{A}{2\psi}\right)^2=\left(\frac{A}{2\psi}\right)^2

Thus the streamlines are circles that are tangent to the x-axis at the origin. The velocity field is given by

(u,v)=\left( {\partial \psi \over \partial y}, - {\partial \psi \over \partial x} \right) =  \left(A\frac{y^2-x^2}{(x^2+y^2)^2},-A\frac{2xy}{(x^2+y^2)^2}\right)

The circles in the upper half-plane thus flow clockwise, those in the lower half-plane flow anticlockwise. Note that speeds go as r − 2; and the speed at the origin is infinite.

[edit] Power law with n equals minus 2

{this section is to be completed}

[edit] See also

[edit] External links

In other languages