Potential flow

File:Streamlines around a NACA 0012.svg

Potential-flow streamlines around a NACA 0012 airfoil at 11° angle of attack, with upper and lower streamtubes identified. The flow is two-dimensional and the airfoil has infinite span.

In fluid dynamics, potential flow or irrotational flow refers to a description of a fluid flow with no vorticity in it. Such a description typically arises in the limit of vanishing viscosity, i.e., for an inviscid fluid and with no vorticity present in the flow.

Potential flow describes the velocity field as the gradient of a scalar function: the velocity potential. As a result, a potential flow is characterized by an irrotational velocity field, which is a valid approximation for several applications. The irrotationality of a potential flow is due to the curl of the gradient of a scalar always being equal to zero.

In the case of an incompressible flow the velocity potential satisfies Laplace's equation, and potential theory is applicable. However, potential flows also have been used to describe compressible flows and Hele-Shaw flows. The potential flow approach occurs in the modeling of both stationary as well as nonstationary flows.

Applications of potential flow include: the outer flow field for aerofoils, water waves, electroosmotic flow, and groundwater flow. For flows (or parts thereof) with strong vorticity effects, the potential flow approximation is not applicable. In flow regions where vorticity is known to be important, such as wakes and boundary layers, potential flow theory is not able to provide reasonable predictions of the flow.<ref name=B_378_380>Batchelor (1973) pp. 378–380.</ref> However, there are often large regions of a flow in which the assumption of irrotationality is valid, allowing the use of potential flow for various applications; these include flow around aircraft, groundwater flow, acoustics, water waves, and electroosmotic flow.<ref name=Kirby>Template:Citation</ref>

Description and characteristics

File:Construction of a potential flow.svg

A potential flow is constructed by adding simple elementary flows and observing the result.

File:Potential cylinder.svg

Streamlines for the incompressible potential flow around a circular cylinder in a uniform onflow.

In potential or irrotational flow, the vorticity vector field is zero, i.e.,

<math display="block">\boldsymbol\omega \equiv \nabla\times\mathbf v=0,</math>

where <math>\mathbf v(\mathbf x,t)</math> is the velocity field and <math>\boldsymbol\omega(\mathbf x,t)</math> is the vorticity field. Like any vector field having zero curl, the velocity field can be expressed as the gradient of certain scalar, say <math>\varphi(\mathbf x,t)</math> which is called the velocity potential, since the curl of the gradient is always zero. We therefore have<ref name=B_99_101>Batchelor (1973) pp. 99–101.</ref>

<math display="block"> \mathbf{v} = \nabla \varphi.</math>

The velocity potential is not uniquely defined since one can add to it an arbitrary function of time, say <math>f(t)</math>, without affecting the relevant physical quantity which is <math>\mathbf v</math>. The non-uniqueness is usually removed by suitably selecting appropriate initial or boundary conditions satisfied by <math>\varphi</math> and as such the procedure may vary from one problem to another.

In potential flow, the circulation <math>\Gamma</math> around any simply-connected contour <math>C</math> is zero. This can be shown using the Stokes theorem,

<math display="block">\Gamma \equiv \oint_C \mathbf v\cdot d\mathbf l = \int \boldsymbol\omega\cdot d\mathbf f=0</math>

where <math>d\mathbf l</math> is the line element on the contour and <math>d\mathbf f</math> is the area element of any surface bounded by the contour. In multiply-connected space (say, around a contour enclosing solid body in two dimensions or around a contour enclosing a torus in three-dimensions) or in the presence of concentrated vortices, (say, in the so-called irrotational vortices or point vortices, or in smoke rings), the circulation <math>\Gamma</math> need not be zero. In the former case, Stokes theorem cannot be applied and in the later case, <math>\boldsymbol\omega</math> is non-zero within the region bounded by the contour. Around a contour encircling an infinitely long solid cylinder with which the contour loops <math>N</math> times, we have

<math display="block">\Gamma = N \kappa</math>

where <math>\kappa</math> is a cyclic constant. This example belongs to a doubly-connected space. In an <math>n</math>-tuply connected space, there are <math>n-1</math> such cyclic constants, namely, <math>\kappa_1,\kappa_2,\dots,\kappa_{n-1}.</math>

Incompressible flow

In case of an incompressible flow — for instance of a liquid, or a gas at low Mach numbers; but not for sound waves — the velocity Template:Math has zero divergence:<ref name=B_99_101/>

<math display="block">\nabla \cdot \mathbf{v} =0 \,,</math>

Substituting here <math>\mathbf v = \nabla\varphi</math> shows that <math>\varphi</math> satisfies the Laplace equation<ref name=B_99_101/>

<math display="block">\nabla^2 \varphi = 0 \,,</math>

where Template:Math is the Laplace operator (sometimes also written Template:Math). Since solutions of the Laplace equation are harmonic functions, every harmonic function represents a potential flow solution. As evident, in the incompressible case, the velocity field is determined completely from its kinematics: the assumptions of irrotationality and zero divergence of flow. Dynamics in connection with the momentum equations, only have to be applied afterwards, if one is interested in computing pressure field: for instance for flow around airfoils through the use of Bernoulli's principle.

In incompressible flows, contrary to common misconception, the potential flow indeed satisfies the full Navier–Stokes equations, not just the Euler equations, because the viscous term

<math display="block">\mu\nabla^2\mathbf v = \mu\nabla(\nabla\cdot\mathbf v)-\mu\nabla\times\boldsymbol\omega=0</math>

is identically zero. It is the inability of the potential flow to satisfy the required boundary conditions, especially near solid boundaries, makes it invalid in representing the required flow field. If the potential flow satisfies the necessary conditions, then it is the required solution of the incompressible Navier–Stokes equations.

In two dimensions, with the help of the harmonic function <math>\varphi</math> and its conjugate harmonic function <math>\psi</math> (stream function), incompressible potential flow reduces to a very simple system that is analyzed using complex analysis (see below).

Compressible flow

Steady flow

Potential flow theory can also be used to model irrotational compressible flow. The derivation of the governing equation for <math>\varphi</math> from Eulers equation is quite straightforward. The continuity and the (potential flow) momentum equations for steady flows are given by

<math display="block">\rho \nabla\cdot\mathbf v + \mathbf v\cdot\nabla \rho = 0, \quad (\mathbf v \cdot\nabla)\mathbf v= -\frac{1}{\rho}\nabla p = -\frac{c^2}{\rho}\nabla \rho</math>

where the last equation follows from the fact that entropy is constant for a fluid particle and that square of the sound speed is <math>c^2=(\partial p/\partial\rho)_s</math>. Eliminating <math>\nabla\rho</math> from the two governing equations results in

<math display="block">c^2\nabla\cdot\mathbf v - \mathbf v\cdot (\mathbf v \cdot \nabla)\mathbf v=0.</math>

The incompressible version emerges in the limit <math>c\to\infty</math>. Substituting here <math>\mathbf v=\nabla\varphi</math> results in<ref name="landau">Template:Cite book</ref><ref name=Anderson>Template:Cite book</ref>

<math display="block">(c^2-\varphi_x^2)\varphi_{xx}+(c^2-\varphi_y^2)\varphi_{yy}+(c^2-\varphi_z^2)\varphi_{zz}-2(\varphi_x\varphi_y\varphi_{xy}+\varphi_y\varphi_z\varphi_{yz}+\varphi_z\varphi_x\phi_{zx})=0</math>

where <math>c=c(v)</math> is expressed as a function of the velocity magnitude <math>v^2=(\nabla\phi)^2</math>. For a polytropic gas, <math>c^2 = (\gamma-1)(h_0-v^2/2)</math>, where <math>\gamma</math> is the specific heat ratio and <math>h_0</math> is the stagnation enthalpy. In two dimensions, the equation simplifies to

<math display="block">(c^2-\varphi_x^2)\varphi_{xx}+(c^2-\varphi_y^2)\varphi_{yy}-2\varphi_x\varphi_y\varphi_{xy}=0.</math>

Validity: As it stands, the equation is valid for any inviscid potential flows, irrespective of whether the flow is subsonic or supersonic (e.g. Prandtl–Meyer flow). However in supersonic and also in transonic flows, shock waves can occur which can introduce entropy and vorticity into the flow making the flow rotational. Nevertheless, there are two cases for which potential flow prevails even in the presence of shock waves, which are explained from the (not necessarily potential) momentum equation written in the following form

<math display="block">\nabla (h+v^2/2) - \mathbf v\times\boldsymbol\omega = T \nabla s</math>

where <math>h</math> is the specific enthalpy, <math>\boldsymbol\omega</math> is the vorticity field, <math>T</math> is the temperature and <math>s</math> is the specific entropy. Since in front of the leading shock wave, we have a potential flow, Bernoulli's equation shows that <math>h+v^2/2</math> is constant, which is also constant across the shock wave (Rankine–Hugoniot conditions) and therefore we can writeTemplate:R

<math display="block">\mathbf v\times\boldsymbol\omega = -T \nabla s</math>

1) When the shock wave is of constant intensity, the entropy discontinuity across the shock wave is also constant i.e., <math>\nabla s=0</math> and therefore vorticity production is zero. Shock waves at the pointed leading edge of two-dimensional wedge or three-dimensional cone (Taylor–Maccoll flow) has constant intensity. 2) For weak shock waves, the entropy jump across the shock wave is a third-order quantity in terms of shock wave strength and therefore <math>\nabla s</math> can be neglected. Shock waves in slender bodies lies nearly parallel to the body and they are weak.

Nearly parallel flows: When the flow is predominantly unidirectional with small deviations such as in flow past slender bodies, the full equation can be further simplified. Let <math>U\mathbf{e}_x</math> be the mainstream and consider small deviations from this velocity field. The corresponding velocity potential can be written as <math>\varphi = x U + \phi</math> where <math>\phi</math> characterizes the small departure from the uniform flow and satisfies the linearized version of the full equation. This is given by

<math display="block">(1-M^2) \frac{\partial^2\phi}{\partial x^2} + \frac{\partial^2\phi}{\partial y^2} + \frac{\partial^2\phi}{\partial z^2} =0</math>

where <math>M=U/c_\infty</math> is the constant Mach number corresponding to the uniform flow. This equation is valid provided <math>M</math> is not close to unity. When <math>|M-1|</math> is small (transonic flow), we have the following nonlinear equationTemplate:R

<math display="block">2\alpha_*\frac{\partial\phi}{\partial x} \frac{\partial^2\phi}{\partial x^2} = \frac{\partial^2\phi}{\partial y^2} + \frac{\partial^2\phi}{\partial z^2}</math>

where <math>\alpha_*</math> is the critical value of Landau derivative <math>\alpha = (c^4/2\upsilon^3)(\partial^2 \upsilon/\partial p^2)_s</math><ref>1942, Landau, L.D. "On shock waves" J. Phys. USSR 6 229-230</ref><ref>Thompson, P. A. (1971). A fundamental derivative in gasdynamics. The Physics of Fluids, 14(9), 1843-1849.</ref> and <math>\upsilon=1/\rho</math> is the specific volume. The transonic flow is completely characterized by the single parameter <math>\alpha_*</math>, which for polytropic gas takes the value <math>\alpha_*=\alpha=(\gamma+1)/2</math>. Under hodograph transformation, the transonic equation in two-dimensions becomes the Euler–Tricomi equation.

Unsteady flow

The continuity and the (potential flow) momentum equations for unsteady flows are given by

<math display="block">\frac{\partial\rho}{\partial t} + \rho \nabla\cdot\mathbf v + \mathbf v\cdot\nabla \rho = 0, \quad \frac{\partial\mathbf v}{\partial t}+ (\mathbf v \cdot\nabla)\mathbf v= -\frac{1}{\rho}\nabla p =-\frac{c^2}{\rho}\nabla \rho=-\nabla h.</math>

The first integral of the (potential flow) momentum equation is given by

<math display="block">\frac{\partial\varphi}{\partial t} + \frac{v^2}{2} + h = f(t), \quad \Rightarrow \quad \frac{\partial h}{\partial t} = -\frac{\partial^2\varphi}{\partial t^2} - \frac{1}{2}\frac{\partial v^2}{\partial t} + \frac{df}{dt}</math>

where <math>f(t)</math> is an arbitrary function. Without loss of generality, we can set <math>f(t)=0</math> since <math>\varphi</math> is not uniquely defined. Combining these equations, we obtain

<math display="block">\frac{\partial^2\varphi}{\partial t^2} + \frac{\partial v^2}{\partial t}=c^2\nabla\cdot\mathbf v - \mathbf v\cdot (\mathbf v \cdot \nabla)\mathbf v.</math>

Substituting here <math>\mathbf v=\nabla\varphi</math> results in

<math display="block">\varphi_{tt} + (\varphi_x^2+ \varphi_y^2+ \varphi_z^2)_t= (c^2-\varphi_x^2)\varphi_{xx}+(c^2-\varphi_y^2)\varphi_{yy}+(c^2-\varphi_z^2)\varphi_{zz}-2(\varphi_x\varphi_y\varphi_{xy}+\varphi_y\varphi_z\varphi_{yz}+\varphi_z\varphi_x\phi_{zx}).</math>

Nearly parallel flows: As in before, for nearly parallel flows, we can write (after introudcing a recaled time <math>\tau=c_\infty t</math>)

<math display="block">\frac{\partial^2\phi}{\partial \tau^2} + 2M \frac{\partial^2\phi}{\partial x\partial\tau}= (1-M^2) \frac{\partial^2\phi}{\partial x^2} + \frac{\partial^2\phi}{\partial y^2} + \frac{\partial^2\phi}{\partial z^2}</math>

provided the constant Mach number <math>M</math> is not close to unity. When <math>|M-1|</math> is small (transonic flow), we have the following nonlinear equationTemplate:R

<math display="block">\frac{\partial^2\phi}{\partial \tau^2} + 2\frac{\partial^2\phi}{\partial x\partial\tau} = -2\alpha_*\frac{\partial\phi}{\partial x} \frac{\partial^2\phi}{\partial x^2} + \frac{\partial^2\phi}{\partial y^2} + \frac{\partial^2\phi}{\partial z^2}.</math>

Sound waves: In sound waves, the velocity magntiude <math>v</math> (or the Mach number) is very small, although the unsteady term is now comparable to the other leading terms in the equation. Thus neglecting all quadratic and higher-order terms and noting that in the same approximation, <math>c</math> is a constant (for example, in polytropic gas <math>c^2=(\gamma-1)h_0</math>), we have<ref>Lamb (1994) §287, pp. 492–495.</ref>Template:R

<math display="block">\frac{\partial^2 \varphi}{\partial t^2} = c^2 \nabla^2 \varphi,</math>

which is a linear wave equation for the velocity potential Template:Mvar. Again the oscillatory part of the velocity vector Template:Math is related to the velocity potential by Template:Math, while as before Template:Math is the Laplace operator, and Template:Mvar is the average speed of sound in the homogeneous medium. Note that also the oscillatory parts of the pressure Template:Mvar and density Template:Mvar each individually satisfy the wave equation, in this approximation.

Applicability and limitations

Potential flow does not include all the characteristics of flows that are encountered in the real world. Potential flow theory cannot be applied for viscous internal flows,<ref name=B_378_380/> except for flows between closely spaced plates. Richard Feynman considered potential flow to be so unphysical that the only fluid to obey the assumptions was "dry water" (quoting John von Neumann).<ref>Template:Citation, p. 40-3. Chapter 40 has the title: The flow of dry water.</ref> Incompressible 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.<ref name=B_404_405>Batchelor (1973) pp. 404–405.</ref> More precisely, potential flow cannot account for the behaviour of flows that include a boundary layer.<ref name=B_378_380/> Nevertheless, understanding potential flow is important in many branches of fluid mechanics. In particular, simple potential flows (called elementary 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.{{ safesubst:#invoke:Unsubst||date=__DATE__ |$B= Template:Fix }}

Analysis for two-dimensional incompressible flow

{{#invoke:Labelled list hatnote|labelledList|Main article|Main articles|Main page|Main pages}} Potential flow in two dimensions is simple to analyze using conformal mapping, by the use of transformations of the complex plane. However, use of complex numbers is not required, as for example in the classical analysis of fluid flow past a cylinder. It is not possible to solve a potential flow using complex numbers in three dimensions.<ref name=B_106_108>Batchelor (1973) pp. 106–108.</ref>

The basic idea is to use a holomorphic (also called analytic) or meromorphic function Template:Mvar, which maps the physical domain Template:Math to the transformed domain Template:Math. While Template:Mvar, Template:Mvar, Template:Mvar and Template:Mvar are all real valued, it is convenient to define the complex quantities

<math display="block">\begin{align}

 z &= x + iy \,, \text{ and } &
 w &= \varphi + i\psi \,.

\end{align}</math>

Now, if we write the mapping Template:Mvar as<ref name=B_106_108/>

<math display="block">\begin{align}

 f(x + iy) &= \varphi + i\psi \,, \text{ or } &
      f(z) &= w \,.

\end{align}</math>

Then, because Template:Mvar is a holomorphic or meromorphic function, it has to satisfy the Cauchy–Riemann equations<ref name=B_106_108/>

<math display="block">\begin{align}

 \frac{\partial\varphi}{\partial x} &=  \frac{\partial\psi}{\partial y} \,, &
 \frac{\partial\varphi}{\partial y} &= -\frac{\partial\psi}{\partial x} \,.

\end{align}</math>

The velocity components Template:Math, in the Template:Math directions respectively, can be obtained directly from Template:Mvar by differentiating with respect to Template:Mvar. That is<ref name=B_106_108/>

<math display="block">\frac{df}{dz} = u - iv</math>

So the velocity field Template:Math is specified by<ref name=B_106_108/>

<math display="block">\begin{align}

 u &= \frac{\partial\varphi}{\partial x} = \frac{\partial\psi}{\partial y}, &
 v &= \frac{\partial\varphi}{\partial y} = -\frac{\partial\psi}{\partial x} \,.

\end{align}</math>

Both Template:Mvar and Template:Mvar then satisfy Laplace's equation:<ref name=B_106_108/>

<math display="block">\begin{align}

 \Delta\varphi &= \frac{\partial^2\varphi}{\partial x^2} + \frac{\partial^2\varphi}{\partial y^2} = 0 \,,\text{ and } &
    \Delta\psi &= \frac{\partial^2\psi}{\partial x^2} + \frac{\partial^2\psi}{\partial y^2} = 0 \,.

\end{align}</math>

So Template:Mvar can be identified as the velocity potential and Template:Mvar is called the stream function.<ref name=B_106_108/> Lines of constant Template:Mvar are known as streamlines and lines of constant Template:Mvar are known as equipotential lines (see equipotential surface).

Streamlines and equipotential lines are orthogonal to each other, since<ref name=B_106_108/>

<math display="block">

 \nabla \varphi \cdot \nabla \psi =
 \frac{\partial\varphi}{\partial x} \frac{\partial\psi}{\partial x} + \frac{\partial\varphi}{\partial y} \frac{\partial\psi}{\partial y} =
 \frac{\partial \psi}{\partial y} \frac{\partial \psi}{\partial x} - \frac{\partial \psi}{\partial x} \frac{\partial \psi}{\partial y} =
 0 \,.

</math>

Thus the flow occurs along the lines of constant Template:Mvar and at right angles to the lines of constant Template:Mvar.<ref name=B_106_108/>

Template:Math is also satisfied, this relation being equivalent to Template:Math. So the flow is irrotational. The automatic condition Template:Math then gives the incompressibility constraint Template:Math.

Examples of two-dimensional incompressible flows

{{#invoke:Labelled list hatnote|labelledList|Main article|Main articles|Main page|Main pages}} Any differentiable function may be used for Template:Mvar. 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.

Power laws

In case the following power-law conformal map is applied, from Template:Math to Template:Math:<ref name=B_409_413>Batchelor (1973) pp. 409–413.</ref>

<math display="block">w=Az^n \,,</math>

then, writing Template:Mvar in polar coordinates as Template:Math, we have<ref name=B_409_413/>

<math display="block">\varphi=Ar^n\cos n\theta \qquad \text{and} \qquad \psi=Ar^n\sin n\theta \,.</math>

In the figures to the right examples are given for several values of Template:Mvar. The black line is the boundary of the flow, while the darker blue lines are streamlines, and the lighter blue lines are equi-potential lines. Some interesting powers Template:Mvar are:<ref name=B_409_413/>

• Template:Math: this corresponds with flow around a semi-infinite plate,
• Template:Math: flow around a right corner,
• Template:Math: a trivial case of uniform flow,
• Template:Math: flow through a corner, or near a stagnation point, and
• Template:Math: flow due to a source doublet

The constant Template:Mvar is a scaling parameter: its absolute value Template:Math determines the scale, while its argument Template:Math introduces a rotation (if non-zero).

Power laws with Template:Math: uniform flow

If Template:Math, that is, a power law with Template:Math, the streamlines (i.e. lines of constant Template:Mvar) are a system of straight lines parallel to the Template:Mvar-axis. This is easiest to see by writing in terms of real and imaginary components:

<math display="block">f(x+iy) = A\, (x+iy) = Ax + i Ay </math>

thus giving Template:Math and Template:Math. This flow may be interpreted as uniform flow parallel to the Template:Mvar-axis.

Power laws with Template:Math

If Template:Math, then Template:Math and the streamline corresponding to a particular value of Template:Mvar are those points satisfying

<math display="block">\psi=Ar^2\sin 2\theta \,,</math>

which is a system of rectangular hyperbolae. This may be seen by again rewriting in terms of real and imaginary components. Noting that Template:Math and rewriting Template:Math and Template:Math it is seen (on simplifying) that the streamlines are given by

<math display="block">\psi=2Axy \,.</math>

The velocity field is given by Template:Math, or

<math display="block">\begin{pmatrix} u \\ v \end{pmatrix} = \begin{pmatrix} \frac{\partial \varphi}{\partial x} \\[2px] \frac{\partial \varphi}{\partial y} \end{pmatrix} = \begin{pmatrix} + {\partial \psi \over \partial y} \\[2px] - {\partial \psi \over \partial x} \end{pmatrix} = \begin{pmatrix} +2Ax \\[2px] -2Ay \end{pmatrix} \,.</math>

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 Template:Math at Template:Math). The Template:Math streamline is particularly interesting: it has two (or four) branches, following the coordinate axes, i.e. Template:Math and Template:Math. As no fluid flows across the Template:Mvar-axis, it (the Template:Mvar-axis) may be treated as a solid boundary. It is thus possible to ignore the flow in the lower half-plane where Template:Math and to focus on the flow in the upper halfplane. 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) Template:Math are ignored.

Power laws with Template:Math

If Template:Math, the resulting flow is a sort of hexagonal version of the Template:Math case considered above. Streamlines are given by, Template:Math and the flow in this case may be interpreted as flow into a 60° corner.

Power laws with Template:Math: doublet

If Template:Math, the streamlines are given by

<math display="block">\psi = -\frac{A}{r}\sin\theta.</math>

This is more easily interpreted in terms of real and imaginary components: <math display="block">\begin{align}

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

\end{align}</math>

Thus the streamlines are circles that are tangent to the x-axis at the origin. The circles in the upper half-plane thus flow clockwise, those in the lower half-plane flow anticlockwise. Note that the velocity components are proportional to Template:Math; and their values at the origin is infinite. This flow pattern is usually referred to as a doublet, or dipole, and can be interpreted as the combination of a source-sink pair of infinite strength kept an infinitesimally small distance apart. The velocity field is given by

<math display="block">(u,v)=\left( \frac{\partial \psi}{\partial y}, - \frac{\partial \psi}{\partial x} \right) = \left(A\frac{y^2-x^2}{\left(x^2+y^2\right)^2},-A\frac{2xy}{\left(x^2+y^2\right)^2}\right) \,.</math>

or in polar coordinates:

<math display="block">(u_r, u_\theta)=\left( \frac{1}{r} \frac{\partial \psi}{\partial \theta}, - \frac{\partial \psi}{\partial r} \right) = \left(-\frac{A}{r^2}\cos\theta, -\frac{A}{r^2}\sin\theta\right) \,.</math>

Power laws with Template:Math: quadrupole

If Template:Math, the streamlines are given by

<math display="block">\psi=-\frac{A}{r^2}\sin 2 \theta \,.</math>

This is the flow field associated with a quadrupole.<ref>Template:Cite book</ref>

Line source and sink

A line source or sink of strength <math>Q</math> (<math>Q>0</math> for source and <math>Q<0</math> for sink) is given by the potential

<math display="block">w = \frac{Q}{2\pi} \ln z</math>

where <math>Q</math> in fact is the volume flux per unit length across a surface enclosing the source or sink. The velocity field in polar coordinates are

<math display="block">u_r = \frac{Q}{2\pi r},\quad u_\theta=0</math>

i.e., a purely radial flow.

Line vortex

A line vortex of strength <math>\Gamma</math> is given by

<math display="block">w=\frac{\Gamma}{2\pi i}\ln z</math>

where <math>\Gamma</math> is the circulation around any simple closed contour enclosing the vortex. The velocity field in polar coordinates are

<math display="block">u_r = 0,\quad u_\theta=\frac{\Gamma}{2\pi r}</math>

i.e., a purely azimuthal flow.

Analysis for three-dimensional incompressible flows

For three-dimensional flows, complex potential cannot be obtained.

Point source and sink

The velocity potential of a point source or sink of strength <math>Q</math> (<math>Q>0</math> for source and <math>Q<0</math> for sink) in spherical polar coordinates is given by

<math display="block">\phi = -\frac{Q}{4\pi r}</math>

where <math>Q</math> in fact is the volume flux across a closed surface enclosing the source or sink. The velocity field in spherical polar coordinates are

<math display="block">u_r = \frac{Q}{4\pi r^2}, \quad u_\theta=0, \quad u_\phi = 0.</math>

See also

• Potential flow around a circular cylinder
• Aerodynamic potential-flow code
• Conformal mapping
• Darwin drift
• Flownet
• Laplacian field
• Laplace equation for irrotational flow
• Potential theory
• Stream function
• Velocity potential
• Helmholtz decomposition

References

• Template:Citation
• Template:Citation
• Template:Citation
• Template:Citation

Further reading

• Template:Citation
• Template:Citation

External links

• {{#invoke:citation/CS1|citation

|CitationClass=web }}

• {{#invoke:citation/CS1|citation

|CitationClass=web }} — Java applets for exploring conformal maps

• Potential Flow Visualizations - Interactive WebApps