Source near a wall: finding optimal placement

⚡ Key Takeaways

One-line statement: To model fluid flow near an impermeable wall, place a mirror-image source (the "image") symmetrically on the opposite side of the wall, so the combined flow satisfies the no-penetration boundary condition automatically.

Key result in words: A source of strength $m$ at distance $d$ from a flat wall is exactly equivalent (on the fluid side) to that source plus an identical image source at distance $d$ on the other side; the wall becomes a streamline. "Optimal placement" refers to choosing $d$ — and in richer problems, balancing competing constraints — to achieve a desired flow feature such as a stagnation point, minimal wall velocity, or a target pressure distribution.

Where it applies: Two-dimensional ideal (inviscid, incompressible, irrotational) potential flow; groundwater hydrology (wells near impermeable boundaries); electrostatics and heat conduction by mathematical analogy.

Who/when: The method of images traces to William Thomson (Lord Kelvin) in electrostatics (1848) and was carried into hydrodynamics by Helmholtz, Kirchhoff, and the classical potential-theory tradition of the late 19th century.

💡 The Core Insight

A flat impermeable wall is a plane of symmetry for the flow. If you reflect every singularity across the wall, the normal velocities induced by each pair cancel on the wall while the tangential velocities add. Hence the wall is a streamline of the doubled, unbounded flow. You never solve a PDE with boundary conditions; you solve free-space superposition and let symmetry enforce the wall.

DefinitionTwo-Dimensional Source

A line source of strength $m$ (units: $\text{m}^2/\text{s}$ , area per time per unit depth) located at the origin produces a purely radial velocity field. In polar coordinates $(r,\theta)$ :

1

\quad u_r = \frac{m}{2\pi r}, \qquad u_\theta = 0.

Derivation hint: Conservation of volume — flux through any circle of radius $r$ is $\oint u_r\, r\, d\theta = m$ , independent of $r$ ; solve for $u_r$ . Its complex potential (with $z = x + iy$ ) is

2

\quad w(z) = \phi + i\psi = \frac{m}{2\pi}\ln z,

where $\phi$ is the velocity potential and $\psi$ the stream function. The velocity is recovered from $u - iv = dw/dz = \dfrac{m}{2\pi z}.$ Derivation hint: For a holomorphic $w$ , the Cauchy–Riemann equations make $\phi,\psi$ harmonic conjugates; differentiate $w=\frac{m}{2\pi}\ln z = \frac{m}{2\pi}(\ln r + i\theta)$ and match real/imaginary parts to $u_r,u_\theta$ .

DefinitionNo-Penetration Boundary Condition

On a solid impermeable wall with unit normal $\hat{\mathbf n}$ , the fluid cannot cross:

3

\quad \mathbf u \cdot \hat{\mathbf n} = 0 \quad \text{on the wall}.

Derivation hint: Mass cannot accumulate inside a rigid solid; the normal flux at the surface must be zero (for a stationary wall). Equivalently, the wall is a streamline: $\psi = \text{const}$ along it, since flow runs parallel to streamlines.

TheoremImage of a Source in a Plane Wall

Let a source of strength $m$ be located at $(0, d)$ with $d>0$ , and let the wall be the line $y = 0$ . The flow in the upper half-plane $y>0$ satisfying $\mathbf u\cdot\hat{\mathbf n}=0$ on $y=0$ is given by the superposition of the real source at $(0,d)$ and an identical image source of strength $m$ at $(0,-d)$ :

4

\quad w(z) = \frac{m}{2\pi}\Big[\ln(z - i d) + \ln(z + i d)\Big].

Derivation hint: Place the image by reflection across $y=0$ (i.e. $id \to -id$ ). The combined stream function is symmetric under $y\to -y$ , so $\partial\psi/\partial x = v = 0$ requires checking; instead verify directly that the normal velocity vanishes on $y=0$ . The wall is a streamline, and the upper-half-plane flow is unique.

ProofThe Image Cancels Normal Velocity on the Wall

Step 1 — Write the velocity field. Differentiate [4]:

5

\quad u - iv = \frac{dw}{dz} = \frac{m}{2\pi}\left(\frac{1}{z - i d} + \frac{1}{z + i d}\right).

Derivation hint: $\dfrac{d}{dz}\ln(z\mp id) = \dfrac{1}{z\mp id}$ ; sum the two contributions.

Step 2 — Restrict to the wall $y=0$ , so $z = x$ (real). Then

6

\quad (u - iv)\big|_{y=0} = \frac{m}{2\pi}\left(\frac{1}{x - i d} + \frac{1}{x + i d}\right) = \frac{m}{2\pi}\cdot\frac{2x}{x^2 + d^2}.

Derivation hint: Combine over a common denominator; $(x-id)(x+id)=x^2+d^2$ , numerator $(x+id)+(x-id)=2x$ is real.

Step 3 — Read off the components. The right-hand side of [6] is purely real, so

7

\quad u\big|_{y=0} = \frac{m x}{\pi (x^2 + d^2)}, \qquad v\big|_{y=0} = 0.

Derivation hint: For $u-iv$ real, the imaginary part $-v$ is zero; equate real parts for $u$ .

Step 4 — Conclude. The wall normal is $\hat{\mathbf n} = \hat{\mathbf y}$ , and $v = \mathbf u\cdot\hat{\mathbf n} = 0$ on $y=0$ . Condition [3] is satisfied identically. By uniqueness of solutions to Laplace's equation with given Neumann data on the boundary (and decay at infinity), [4] is the flow. $\blacksquare$

∎

Figure

Source at height

d

above a wall and its mirror image; the wall (dashed) is a streamline.

WT

William Thomson, Lord Kelvin

1824–1907

Kelvin introduced the method of images in 1848 to solve electrostatic problems involving a point charge near a grounded conducting plane, replacing the conductor with a mirror charge of opposite sign. The same geometric trick, transplanted to potential flow, gives the source-near-a-wall result — though with an image of the same sign, because the relevant boundary condition (no normal flow) differs from the electrostatic one (constant potential on the conductor). Kelvin's broad program of treating heat, electricity, and fluid motion through the common language of Laplace's equation made such cross-disciplinary transfers natural and routine.

🧠 Intuition

Imagine the wall as a perfect mirror. Every drop of fluid that the real source pushes toward the wall is met head-on by an identical drop that the image source pushes away from the wall on the other side. Along the mirror line these two impulses are exactly anti-parallel in their normal components, so they cancel — fluid is deflected to slide along the wall rather than through it. The tangential components, by contrast, point the same way and reinforce, which is why the speed along the wall is nonzero and peaks where the geometry is most favorable. The whole construction works because a source's field is perfectly radial and symmetric: reflection maps it onto another genuine source, not some distorted object.

HV

Hermann von Helmholtz

1821–1894

Helmholtz's vortex theorems and his systematic use of potential theory in fluid dynamics (1858 onward) established the framework in which sources, sinks, vortices, and their images became the elementary building blocks of two-dimensional inviscid flow. The superposition principle that makes image methods possible rests on the linearity of Laplace's equation $\nabla^2\phi = 0$ , the governing equation of irrotational incompressible flow that Helmholtz did so much to popularize among physicists.

ExampleWall Velocity and Its Maximum

Given: A source of strength $m = 4\pi\ \text{m}^2/\text{s}$ at height $d = 2\ \text{m}$ above the wall $y=0$ .

Find: (a) the velocity on the wall as a function of $x$ ; (b) the location and value of the maximum wall speed.

Step 1. From [7], the wall velocity is purely tangential:

u(x) = \frac{m x}{\pi(x^2 + d^2)} = \frac{4\pi\, x}{\pi(x^2 + 4)} = \frac{4x}{x^2+4}.

Step 2. Maximize: $\dfrac{du}{dx} = \dfrac{4(x^2+4) - 4x(2x)}{(x^2+4)^2} = \dfrac{4(4 - x^2)}{(x^2+4)^2} = 0 \implies x = \pm 2 = \pm d.$

Step 3. Evaluate at $x = d = 2$ :

u_{\max} = \frac{4(2)}{4+4} = \frac{8}{8} = 1\ \text{m/s}.

> Answer: $u(x) = \dfrac{4x}{x^2+4}\ \text{m/s}$ ; maximum speed $\boxed{u_{\max} = 1\ \text{m/s at } x = \pm d = \pm 2\ \text{m}}$ . The fastest flow along the wall occurs directly beneath a 45° line from the source.

ExampleSource in a Uniform Stream Near a Wall: Stagnation Point

Given: A uniform stream $U = 3\ \text{m/s}$ flows in the $+x$ direction parallel to the wall $y=0$ . A source of strength $m = 12\pi\ \text{m}^2/\text{s}$ sits at $(0,d)$ with $d = 1\ \text{m}$ . Use the image to enforce the wall.

Setup. The complex potential is the stream plus source plus image:

w(z) = U z + \frac{m}{2\pi}\Big[\ln(z - i d) + \ln(z + i d)\Big].

Step 1 — Velocity on the wall ( $y=0$ , $z=x$ ). Add $U$ to the tangential result [7]:

u(x) = U + \frac{m x}{\pi(x^2+d^2)} = 3 + \frac{12\pi\, x}{\pi(x^2+1)} = 3 + \frac{12x}{x^2+1}, \qquad v=0.

Step 2 — Find stagnation points on the wall ( $u=0$ ):

3 + \frac{12x}{x^2+1} = 0 \implies 3(x^2+1) + 12x = 0 \implies x^2 + 4x + 1 = 0.

Step 3 — Solve the quadratic:

x = \frac{-4 \pm \sqrt{16 - 4}}{2} = -2 \pm \sqrt{3} \implies x \approx -0.27\ \text{m} \ \text{ or } \ x \approx -3.73\ \text{m}.

> Answer: Two wall stagnation points at $\boxed{x = -2\pm\sqrt 3 \approx -0.27\ \text{m},\ -3.73\ \text{m}}$ . They lie upstream of the source (negative $x$ ), where the source's leftward-projected outflow opposes the oncoming stream. This is the elementary model of a half-body / Rankine-type fairing sitting on a wall.

⚠️ Common mistake

The image of a source in a wall is a source, not a sink. This trips up students who remember the electrostatic image of a charge near a grounded conductor, which is an opposite charge. The boundary conditions differ: for a conductor the potential is constant (Dirichlet), forcing opposite sign; for an impermeable wall the normal velocity (a Neumann condition on $\phi$ ) is zero, which is satisfied by an image of the same sign and strength. If you mistakenly use a sink, the wall becomes a source/sink line, not a streamline.

⚠️ Common mistake

"Optimal placement" is not a single fixed rule — it depends on the objective. Maximizing wall speed, locating a stagnation point, minimizing peak wall pressure, or keeping a dividing streamline below a given height are different optimization problems with different answers. Always state the objective and constraint before differentiating.

Counterexample

A curved or finite wall is generally NOT solved by a single image. Consider a source near a circular boundary instead of a flat one. The correct construction (the Milne–Thomson circle theorem) requires an image at the inverse point $a^2/\bar{z_0}$ plus a compensating source at the center — a single mirror reflection fails because reflecting a radial field across a circle does not produce another single radial source. Treating a circular wall like a flat mirror gives a flow that does not satisfy $\mathbf u\cdot\hat{\mathbf n}=0$ on the circle.

connection

The construction is the fluid-dynamical sibling of the electrostatic Method of Images and shares its mathematical backbone with Laplace's Equation and harmonic functions. Adding a uniform stream connects directly to Rankine Half-Body and Rankine Oval flows. The complex-potential formalism ([2], [4]) is the entry point to Conformal Mapping, where the wall can be straightened or curved by analytic transformations. When the singularity is a vortex rather than a source, the image picks up a sign flip in circulation, linking to Vortex Dynamics and the Milne-Thomson Circle Theorem.

🔬 Where this appears

⏳ History

- 1828 — George Green's Essay develops the potential function and Green's functions, the analytic foundation for image methods.

- 1848 — William Thomson (Kelvin) introduces the method of images for a point charge near a grounded conducting plane.

- 1858 — Helmholtz publishes his work on vortex motion, embedding singularities and superposition firmly in fluid dynamics.

- 1868 — Kirchhoff and contemporaries extend free-streamline and potential-flow theory; complex-variable methods mature.

- 1880s — W. J. M. Rankine's earlier half-body and oval constructions (1860s–70s) become standard textbook examples of source+stream superposition near and along boundaries.

- 1906 — Kutta and Joukowski formalize complex-potential aerodynamics, generalizing image and conformal techniques.

- 1940 — L. M. Milne-Thomson states the circle theorem, giving the systematic image for curved (circular) boundaries.

- 1960s — The rise of panel methods (Hess & Smith, 1962) turns image-source superposition into a practical numerical tool for arbitrary geometries.

open-problem

The clean single-image result is fragile: it holds exactly only for an infinite, flat, perfectly impermeable wall in inviscid, irrotational, incompressible two-dimensional flow. Real subtleties accumulate when these assumptions break:

Finite or cornered walls demand infinite sequences of images (e.g., a source between two parallel walls needs an infinite ladder of reflections, summable via the digamma/cotangent functions), and the series' convergence must be handled carefully.
Viscosity turns the slip wall velocity [7] into a no-slip boundary layer; the inviscid result then describes only the outer flow and feeds boundary-layer matching.
Three dimensions replace the logarithmic potential with a $1/r$ point source, and "optimal placement" becomes a genuinely 3-D variational problem.
Multi-objective optimality (minimize peak pressure and keep the half-body below a height and fix the discharge) leads to constrained optimization with no closed form — an active design-engineering problem.
A clean open-flavored question: for $N$ sources of fixed total strength near a wall, what placement minimizes the maximum wall shear in the boundary-layer-corrected model? This couples potential flow to viscous theory and resists closed-form solution.

Source near a wall: finding optimal placement

Formal Statement

Why It Holds / Derivation

Building Intuition

Worked Examples

Common Misconceptions

Connections

Real-World Applications

Historical Development

Subtleties & Limits

See Also