Potential Flows: Velocity Potential, Source, Sink and Vortex, Superposition and Flow past a Cylinder
2. Uniform flow, source, sink and free vortex
Four building blocks solve Laplace's equation in the plane. Let U be the stream speed along +x, m the source strength (volume flow per unit depth, taken positive for a source, negative for a sink) and Γ the circulation of a vortex (counter-clockwise positive).
| Flow | Velocity | φ | ψ |
|---|---|---|---|
| Uniform | u = U, v = 0 | Ux | Uy |
| Source / sink | v_r = m/(2πr), v_θ = 0 | (m/2π) ln r | (m/2π) θ |
| Free vortex | v_θ = Γ/(2πr), v_r = 0 | Γθ/2π | −(Γ/2π) ln r |
The source speed v_r falls as 1/r because the same flow m crosses a circle of circumference 2πr. The free vortex is irrotational everywhere except the singular centre, though the fluid circulates: the circulation round any loop that encloses the centre is Γ, while the circulation round a loop that does not enclose it is zero. Velocity is infinite at the centre of both, so these flows describe only the region away from the origin.
Pressure in a free vortex. Bernoulli applies between any two points of an irrotational flow. With v_θ = Γ/(2πr), p + ½ρΓ²/(4π²r²) is constant, so pressure rises outward: p(r) = p_∞ − ρΓ²/(8π²r²). For water, a vortex with 8 m/s at r = 0.5 m has 2 m/s at r = 2 m, so p₂ − p₁ = ½ × 1000 × (8² − 2²) = 30 kPa.
3. Superposition: source-sink pair and the half-body
Laplace's equation is linear, so the sum of any solutions is a solution: φ = Σφᵢ and ψ = Σψᵢ, and velocity vectors add. Pressure does not add, because Bernoulli is quadratic in velocity; compute the total velocity first and then find p. Superposition lets simple flows build flows past shapes: the streamline ψ = 0 (or another chosen value) that separates fluid from the source is taken as a solid surface, because no fluid crosses a streamline.
Source in a uniform stream (Rankine half-body). ψ = Uy + (m/2π)θ = U r sinθ + (m/2π)θ. The source sends fluid upstream and the stream pushes it back, so they balance at a stagnation point on the axis upstream of the source at distance a = m/(2πU). The streamline through it is a body shape that is open downstream. Far downstream its total width is m/U (half-width m/2U), since the source's whole flow m must pass through a strip of width w moving at U.
Worked example. U = 10 m/s and m = 12 m²/s. The stagnation point lies a = 12/(2π × 10) = 0.191 m upstream of the source, and the body is 12/10 = 1.2 m wide far downstream. Changing m to 24 doubles both the distance and the width.
A source-sink pair of equal strength, separated by 2b, has streamlines that leave the source and enter the sink (the flow is that of a dipole as b → 0). Placed in a uniform stream, the source-sink pair gives a closed oval, the Rankine oval, which is a closed body; the half-body above is its open-ended cousin.
4. Doublet and flow past a circular cylinder
A doublet is the limit of a source-sink pair brought together (b → 0) while the product m·2b = κ is held fixed. Its stream function is ψ = −(κ/2π) sinθ/r. Add it to a uniform stream: ψ = U sinθ (r − R²/r), with R² = κ/(2πU). On r = R, ψ = 0, so the circle r = R is a streamline and can be a solid cylinder of radius R. The velocity at the cylinder surface is purely tangential:
- v_r = U cosθ (1 − R²/r²) = 0 on r = R, and v_θ = −U sinθ (1 + R²/r²) = −2U sinθ there.
- The two stagnation points are at θ = 0 and θ = 180° (front and rear), where the speed is zero.
- The surface speed is largest at θ = 90° and 270°, where it is 2U, twice the free stream.
The pressure coefficient is C_p = (p − p_∞)/(½ρU²) = 1 − (V/U)², by Bernoulli between the surface point and far upstream. With V = 2U sinθ on the surface, C_p = 1 − 4 sin²θ.
| θ | V/U | C_p | Meaning |
|---|---|---|---|
| 0°, 180° | 0 | +1 | stagnation, maximum pressure |
| 30° | 1 | 0 | pressure equals p_∞ |
| 60° | 1.73 | −2 | suction |
| 90° | 2 | −3 | maximum speed, minimum pressure |
Worked example. Air (ρ = 1.2 kg/m³) at U = 40 m/s meets a cylinder. At the shoulder (θ = 90°), p − p_∞ = −3 × ½ × 1.2 × 40² = −2880 Pa, and at the front stagnation point p − p_∞ = +960 Pa. The pressure difference between the two is 4 × 960 = 3840 Pa.
5. Pressure distribution and the d'Alembert paradox
The potential-flow pressure distribution C_p = 1 − 4 sin²θ is symmetric about the flow axis (front to back) and about the stream direction (top to bottom). Integrating pressure over the surface, the forces along and across the stream both vanish. A body in steady, inviscid, irrotational, unseparated flow of unbounded extent therefore feels no drag. This is d'Alembert's paradox, since real bodies certainly experience drag.
The paradox is resolved by the viscous boundary layer (Section XE1.9). A real boundary layer cannot follow the pressure recovery on the rear half, where the pressure rises, so the flow separates and leaves a low-pressure wake. The pressure on the rear is then not restored, front and rear differ, and pressure (form) drag arises, along with a small skin-friction drag. Potential flow is therefore a good model of the outer flow and of the front part of the body, and a poor model of the wake.
6. Rotating cylinder, circulation and the Kutta-Joukowski lift
Add a free vortex of circulation Γ at the centre of the cylinder: ψ = U sinθ (r − R²/r) − (Γ/2π) ln(r/R). Still ψ = 0 on r = R, so the cylinder remains a streamline. On the surface the speed is q = |2U sinθ − Γ/(2πR)|. The stagnation points satisfy sinθ = Γ/(4πUR) (with Γ taken counter-clockwise): they move toward each other along the surface as Γ grows, meet at the top (or bottom) of the cylinder when Γ = 4πUR, and leave the surface for a larger Γ.
With circulation the pressure distribution is no longer symmetric top to bottom: the surface speed is higher on the side where the rotation adds to the stream, so the pressure is lower there, and a net transverse force results. The Kutta-Joukowski theorem gives the lift per unit span on any cylinder of arbitrary cross-section in steady, ideal, irrotational flow: L′ = ρUΓ. The force is perpendicular to the free stream, directed toward the side where the surface moves with the stream, and the drag is still zero.
Worked example. A spinning cylinder of R = 0.10 m has surface speed ωR = 15 m/s. Assuming the circulation equals that of the surface, Γ = 2πR(ωR) = 9.42 m²/s. In air (ρ = 1.2 kg/m³) at U = 20 m/s the lift per metre is L′ = 1.2 × 20 × 9.42 = 226 N/m. This is the Magnus effect; in real flow, viscosity and separation make the lift smaller and some drag remains.
Key takeaways
- A velocity potential exists for irrotational flow; a stream function exists for any 2D incompressible flow. Both satisfy Laplace's equation only for 2D incompressible irrotational flow, where the two families of lines cross at right angles.
- Source: v_r = m/(2πr); free vortex: v_θ = Γ/(2πr); both fall as 1/r. Pressure in a free vortex rises outward, p = p_∞ − ρΓ²/(8π²r²).
- Velocities add under superposition; pressures do not. A source in a uniform stream gives a half-body with its stagnation point at m/(2πU) upstream and far-downstream width m/U.
- Doublet plus uniform stream is flow past a cylinder: surface speed 2U sinθ, C_p = 1 − 4 sin²θ, maximum speed 2U and minimum C_p = −3 at θ = 90°.
- The symmetric potential-flow pressure gives zero drag (d'Alembert's paradox); real drag comes from separation and the wake, which need viscosity.
- Adding circulation Γ gives lift L′ = ρUΓ per unit span (Kutta-Joukowski), perpendicular to the stream, with zero drag in ideal flow.
Practice questions (17)
Attempt each one before opening the answer. Every explanation names the tempting wrong option as well as the right one, because that is where marks are lost.
A velocity potential function φ with V = ∇φ can be defined for a flow only when the flow is
Show answer
Answer: C — irrotational, since then ∇ × V = 0
A curl-free field is the gradient of a scalar, so irrotationality is the condition. Incompressibility is the condition for a stream function (2D) and decides whether φ satisfies Laplace's equation; steadiness and dimension do not decide the existence of φ.In a two-dimensional, incompressible, irrotational flow, which statement about φ and ψ is correct?
Show answer
Answer: A — Both satisfy Laplace's equation, and lines of constant φ cross lines of constant ψ at right angles
Irrotationality gives ∇²ψ = 0 and incompressibility gives ∇²φ = 0, and ∇φ·∇ψ = 0 follows from u = ∂φ/∂x = ∂ψ/∂y and v = ∂φ/∂y = −∂ψ/∂x. The statement that ψ never satisfies Laplace's equation is wrong for exactly this flow, and the claim that ψ is only for irrotational and φ only for rotational flow reverses the truth: φ needs irrotationality.Which of the following stream functions, with u = ∂ψ/∂y and v = −∂ψ/∂x, represents an irrotational flow?
Show answer
Answer: D — ψ = 2xy
Vorticity is ω_z = −∇²ψ. For 2xy, ∇²ψ = 0, so the flow is irrotational. For x² + y² it is 4, for x²y it is 2y and for y³ it is 6y, none of which vanishes everywhere. x² + y² is the tempting one because of its symmetry; it is a solid-body rotation, not a potential flow.A two-dimensional source of strength m = 4π m²/s (volume flow per unit depth) is placed in still water. The radial speed at r = 2.0 m is
Show answer
Answer: D — 1.0 m/s
v_r = m/(2πr) = 4π/(2π × 2.0) = 1.0 m/s. Using m/r (2π m/s ≈ 6.3) or m/(4πr) = 0.5 m/s mistakes the geometry: the flow crosses a circle of circumference 2πr, not a sphere of area 4πr².In a free vortex in an incompressible fluid, with gravity effects ignored, the pressure
Show answer
Answer: C — rises as r increases, since speed falls as 1/r
Bernoulli holds across the whole irrotational flow: p + ½ρv² is constant with v = Γ/(2πr), so p rises outward and is lowest near the core. Uniform pressure would mean no centripetal force. The options with speed rising with r or the greatest pressure at the centre reverse the Bernoulli relation: higher speed goes with lower pressure.A source placed in a uniform stream gives the flow past a semi-infinite body that is open downstream. This body is called
Show answer
Answer: B — a Rankine half-body
A source plus a uniform stream gives the half-body. The oval needs a source and a sink of equal strength, the circular cylinder comes from a doublet, and an aerofoil needs a transformation with circulation.For potential flow of speed U past a circular cylinder, the largest speed on the cylinder surface is
Show answer
Answer: A — 2U, at θ = 90° from the front stagnation point
The surface speed is 2U|sinθ|, which is largest at θ = 90°. The number 3 is the magnitude of the minimum pressure coefficient, not a speed ratio, and at the stagnation points the speed is zero.The potential-flow solution for a steady flow past a non-rotating circular cylinder predicts zero drag. The real drag arises mainly because
Show answer
Answer: B — viscosity makes the boundary layer separate, leaving a low-pressure wake
The potential-flow pressure is symmetric fore and aft, so the force cancels. Real boundary layers separate in the adverse pressure gradient on the rear half, so the rear pressure is not recovered. Compressibility is not needed (the paradox holds at low speed), and the stagnation pressure is correctly p_∞ + ½ρU².For flow past a rotating cylinder of radius R, the surface speed is q = |2U sinθ − Γ/(2πR)|. When Γ = 4πUR, the stagnation points
Show answer
Answer: A — merge into one point on the surface, at 90° to the stream direction
q = 0 needs sinθ = Γ/(4πUR) = 1, so θ = 90°, a single point. For Γ below 4πUR there are two stagnation points on the surface and for Γ above it they move off into the fluid. The lift ρUΓ is not zero; it is proportional to Γ.Which of the following statements about the velocity potential and stream function are correct?
Show answer
Answer: A — A velocity potential exists whenever the flow is irrotational; B — A stream function exists for any two-dimensional incompressible flow, rotational or not
Irrotationality gives φ, and incompressibility in 2D gives ψ. But ψ satisfies Laplace's equation only if the flow is also irrotational, since ω_z = −∇²ψ. In uniform flow along x, φ = Ux, so constant-φ lines are vertical and cross the horizontal streamlines at right angles.Which of the following are correct for the superposition of potential flows?
Show answer
Answer: A — The velocity potentials of the component flows add, because Laplace's equation is linear; C — The velocity vectors of the component flows add at each point
Linearity means φ and ψ add, and so do the velocities. Pressure follows Bernoulli, p = p₀ − ½ρ|V|², which is quadratic in the total velocity, so it must be computed from the summed velocity. Stream functions add, they do not multiply.For potential flow of speed U past a circular cylinder, which of the following statements are correct?
Show answer
Answer: A — The pressure coefficient at a stagnation point is +1; B — The minimum pressure coefficient on the surface is −3; D — The surface pressure at θ = 30° equals the free-stream pressure
C_p = 1 − 4 sin²θ gives +1 at θ = 0 and −3 at θ = 90°. At θ = 30°, sinθ = 0.5, so the speed is 2U × 0.5 = U (not 2U), C_p = 0 and the pressure equals p_∞.A two-dimensional source of strength m = 12 m²/s (volume flow per unit depth) is placed in a uniform stream of speed U = 10 m/s. The stagnation point lies on the axis upstream of the source. The distance of this stagnation point from the source, in metres to three decimal places, is ____.
Numerical answer — type the value.
Show answer
Answer: 0.191
At the stagnation point the source speed m/(2πa) points upstream and cancels U, so U − m/(2πa) = 0. So a = m/(2πU) = 12/(2π × 10) = 0.19099 m, which is 0.191 m. Forgetting the 2π gives 1.2 m, which is the far-downstream width m/U.Air of density 1.2 kg/m³ flows steadily past a long circular cylinder with an undisturbed speed of 30 m/s. Using ideal (potential) flow, find the difference in pressure between the front stagnation point and the point on the surface at 90° from it, in pascals, as a whole number.
Numerical answer — type the value.
Show answer
Answer: 2160
With ½ρU² = 0.5 × 1.2 × 900 = 540 Pa, C_p = +1 at the front and −3 at 90°. The difference is (1 − (−3)) × 540 = 4 × 540 = 2160 Pa. Equivalently, at 90° the speed is 60 m/s, so the difference is ½ × 1.2 × 60² = 2160 Pa.A free vortex in water (ρ = 1000 kg/m³) has a tangential speed of 8.0 m/s at r = 0.50 m. Ignoring gravity, the pressure at r = 2.0 m exceeds the pressure at r = 0.50 m by ____ kPa (answer as a whole number).
Numerical answer — type the value.
Show answer
Answer: 30
v ∝ 1/r gives v₂ = 8 × 0.5/2.0 = 2.0 m/s. Bernoulli: p₂ − p₁ = ½ρ(v₁² − v₂²) = 0.5 × 1000 × (64 − 4) = 30,000 Pa = 30 kPa. Taking v₂ = 4 m/s (treating it as a forced vortex) would give 24 kPa.A two-dimensional flow has the stream function ψ = 3xy (in m²/s, with x and y in metres), and u = ∂ψ/∂y, v = −∂ψ/∂x. The speed at the point (x, y) = (2, 1) m, in m/s to two decimal places, is ____.
Numerical answer — type the value.
Show answer
Answer: 6.71
u = ∂ψ/∂y = 3x = 6 m/s and v = −∂ψ/∂x = −3y = −3 m/s at (2, 1). Speed = √(36 + 9) = √45 = 6.708, which is 6.71 m/s. Swapping the formulas for u and v gives the same speed here, but not the direction.A long cylinder of radius 0.10 m rotates in a uniform air stream (ρ = 1.2 kg/m³, U = 20 m/s). Treat the circulation as Γ = 2πR(ωR) for a surface speed ωR = 15 m/s in ideal flow. Using the Kutta-Joukowski result L′ = ρUΓ, the lift per metre of span is approximately
Show answer
Answer: A — 226 N/m
Γ = 2π × 0.10 × 15 = 9.42 m²/s, and L′ = 1.2 × 20 × 9.42 = 226 N/m. Using Γ = ωR (15) gives 360 N/m, using 2ρUΓ gives 452 N/m and using ½ρUΓ gives 113 N/m; the theorem has no such factor.