The method of images models transient diffusion near an impermeable boundary by mirroring the real source across the wall. For a point solute source at (z=a), an equal fictitious source is placed at (z=-a), where the impermeable boundary is at (z=0). The two transient diffusion fields are superimposed, producing a concentration distribution whose normal gradient is zero at the wall:

[ \left.\frac{\partial C}{\partial z}\right|_{z=0}=0 ]

An impermeable boundary does not remove diffusing mass; it reflects the diffusion field. An equal-strength image source mathematically creates this reflection and enforces the required zero-flux condition.

Why the Boundary Condition Matters

Impermeable walls permit no normal mass flux

For one-dimensional diffusion normal to a solid boundary, Fick's first law gives

[ N_z=-D\frac{\partial C}{\partial z} ]

At an impermeable wall, (N_z=0). Assuming a constant diffusivity, this requires

[ \left.\frac{\partial C}{\partial z}\right|_{z=0}=0 ]

The concentration itself does not need to be zero at the wall. That distinction separates an impermeable boundary from an absorbing or fixed-concentration boundary.

The physical wall is replaced mathematically

The method of images removes the boundary from the mathematical domain by extending the concentration field into a fictitious region. A mirror source is then introduced so that the combined solution has the correct symmetry at the actual wall.

For a no-flux boundary, the image source has the same sign and magnitude as the real source. This produces an even concentration field about the boundary plane.

How the Image Source Satisfies No Flux

Superposition follows from diffusion-equation linearity

For constant diffusivity, transient diffusion is governed by Fick's second law:

[ \frac{\partial C}{\partial t}

D\frac{\partial^2 C}{\partial z^2} ]

Because this equation is linear, solutions can be added. The concentration generated by the real source and the concentration generated by the image source therefore form another valid solution.

For a point source, the free-space Green's function in one dimension has the form

[ G(z,t;a)

\frac{1}{\sqrt{4\pi Dt}} \exp\left[-\frac{(z-a)^2}{4Dt}\right] ]

The wall-bounded solution is

[ C(z,t)

M\left[G(z,t;a)+G(z,t;-a)\right] ]

where (M) represents the source strength, with the appropriate concentration or mass units for the problem.

The concentration gradient cancels at the wall

At (z=0), the real and image sources are equally distant from the boundary. Their concentration gradients have equal magnitude but opposite direction, so they cancel:

[ \left. \frac{\partial C}{\partial z} \right|_{z=0}

0 ]

Consequently,

[ N_z(0,t)

-D\left.\frac{\partial C}{\partial z}\right|_{z=0}

0 ]

The image source does not represent additional physical solute. It is a mathematical construction that reproduces the reflection caused by the impermeable boundary.

What the Model Predicts in a Pilot Unit

Concentrations increase near the boundary

In an unbounded medium, solute can diffuse away from the real source in both directions. Near an impermeable wall, the image source supplies the reflected portion of that diffusion field.

As a result, the concentration near the wall is generally higher than it would be in an equivalent unbounded system. The wall changes the spatial distribution without allowing mass to cross the wall.

The disturbance spreads over a diffusion length

Transient diffusion affects a characteristic distance of order

[ \ell_D\sim (Dt)^{1/2} ]

For step-change or semi-infinite diffusion problems, the concentration profile is commonly expressed with an error function or complementary error function. A typical similarity variable is

[ \eta=\frac{x}{2\sqrt{Dt}} ]

The method of images addresses the boundary condition, while the error-function solution describes how the transient profile evolves with time.

Flux remains zero at the wall but not nearby

The normal flux is exactly zero at the impermeable surface. However, concentration gradients and mass fluxes can be substantial at locations just inside the fluid or porous medium.

This distinction is important in pilot measurements: a wall probe may record elevated concentration even though the wall itself receives no mass flux.

Application to Environmental Unit Operations

Soil contaminant transport rigs

In a soil column or packed-bed pilot unit, an impermeable layer can represent a compacted soil horizon, rock surface, liner, or sealed base. A contaminant source near that layer can be modeled using a reflected image source.

The resulting solution estimates transient concentration profiles near the interface without requiring a numerical mesh to extend through the solid layer.

Diffusion and mass-transfer columns

In a diffusion column with non-reactive walls, the method can approximate localized solute spreading when radial or transverse transport is important. It is particularly useful for interpreting concentration buildup near walls and checking whether wall effects influence measured breakthrough curves.

For absorption, extraction, or membrane-related experiments, the same transient diffusion principles help identify whether the observed response is controlled by early penetration or by longer-term transport through a finite layer.

Pilot-scale interpretation

The model can be used to compare measurements with an analytical benchmark. If concentration data near a wall deviate strongly from the image-source prediction, possible causes include wall adsorption, reaction, leakage, convection, nonuniform porosity, or an incorrect diffusivity.

This makes the method useful not only for prediction but also for diagnosing whether the assumed boundary condition is physically valid.

Understanding the Trade-offs

The method assumes a simple geometry

The basic image construction is exact for an infinite planar boundary with constant properties. It becomes less direct for curved walls, irregular vessel shapes, finite-width columns, or several nearby boundaries.

Multiple parallel impermeable walls can sometimes be handled with repeated image sources, but the resulting series may be cumbersome. Numerical methods are usually more practical for complex geometries.

Diffusion must dominate the transport

The standard solution assumes that molecular diffusion or effective porous-medium diffusion controls the transient. Significant advection, dispersion, density-driven flow, or turbulence can invalidate a purely diffusive image-source model.

In environmental systems, the effective diffusivity may also depend on porosity, tortuosity, saturation, and sorption. These effects must be incorporated separately rather than assumed to be represented by the image source.

Reactive or absorbing walls require a different image

An equal-strength image source enforces a zero-gradient, no-flux condition. It does not correctly describe a wall that consumes solute, maintains a fixed concentration, or permits partial transfer.

For a fixed-concentration boundary, an image source of opposite sign is commonly used in the simplest planar problems. Mixed or reactive boundary conditions generally require another analytical solution or a numerical treatment.

Finite domains eventually depart from the semi-infinite model

At short times, the diffusion zone may be small compared with the pilot unit's dimensions, so the semi-infinite image solution is often appropriate. At longer times, solute can reach other walls, outlets, interfaces, or the far side of a membrane.

The relevant transition can be estimated by comparing ((Dt)^{1/2}) with the distance to the next boundary. Once those distances are comparable, a finite-domain model is needed.

Making the Right Choice for Your Goal

The method is most useful when the wall is planar, impermeable, and close enough to affect the transient concentration field.

If your primary focus is wall-adjacent contaminant dispersion: Use an equal-strength mirrored source to enforce (\partial C/\partial z=0) and evaluate the resulting transient concentration profile.
If your primary focus is pilot-unit design: Compare the diffusion length ((Dt)^{1/2}) with the unit dimensions to determine whether a semi-infinite image solution is adequate.
If your primary focus is membrane or interfacial mass transfer: Combine the boundary treatment with an error-function transient solution to distinguish early penetration from later finite-film behavior.
If your primary focus is realistic environmental prediction: Include advection, dispersion, reaction, sorption, variable properties, and additional boundaries when they are significant.
If your primary focus is experimental validation: Use the analytical image-source solution as a benchmark for testing measured concentration profiles and identifying non-ideal wall effects.

By mirroring the source, the method of images turns an impermeable boundary into a mathematically exact zero-flux condition for transient diffusion in simple geometries.

Summary Table:

Key aspect	Practical insight
Boundary condition	An impermeable wall requires zero normal flux: (N_z=0) and (\partial C/\partial z=0).
Image source	An equal-strength mirror source creates the reflected diffusion field.
Governing model	Superposition of transient Green's functions solves the planar no-flux problem.
Wall effect	Concentrations may increase near the wall, although no mass crosses it.
Diffusion scale	The affected distance grows approximately as (\ell_D\sim(Dt)^{1/2}).
Best applications	Soil columns, packed beds, diffusion columns, and environmental pilot studies.
Key limitations	Complex geometry, advection, reaction, sorption, variable properties, and finite domains may require numerical modeling.

Planning a pilot study or upgrading a mass-transfer training system? Contact LABPARK to discuss Educational and Vocational Unit Operations Pilot Plants tailored to your application. LABPARK supports universities, research institutes, and enterprises with practical systems for chemical engineering, bioprocess and biotech, and environmental and water treatment. Our team can help you evaluate diffusion behavior, model impermeable or reactive boundaries, select suitable instrumentation, and build a reliable platform for teaching, research, process validation, and scale-up.