Here is a genuinely powerful, and genuinely intuitive, fact about steady-state temperature, steady electric potential, or any other harmonic quantity: it can never have a local maximum or minimum sitting strictly inside its domain. The hottest and coldest spots of a metal plate that has settled into equilibrium always occur somewhere on the plate's edge — never buried in the middle. That sounds almost too clean to be a real mathematical theorem, but it has real teeth, and we can prove it outright.
Think about why this has to be physically sensible first. If some interior point really were strictly hotter than every point around it, heat would immediately start flowing away from it in all directions — down the temperature gradient — and the temperature there would drop. A genuine steady state, one that isn't changing anymore, simply can't sustain a hidden interior hot spot. The maximum principle turns that physical intuition into an ironclad mathematical guarantee.
The proof idea is short and comes directly from
Now suppose, for contradiction, that some interior point
Let
Take a circular metal plate whose rim is welded to a ring of thermostats reading
anywhere from
The maximum principle answers this immediately, without solving a single equation: the interior
can be at most as hot as the hottest edge reading —
The picture below shows a simpler, fully solvable case: rim data
The same guarantee works on any shape, not just circles. Picture a square metal plate whose bottom
edge is heated to a steady
This is also a genuine, practical sanity-check for numerical simulations. If you solve a steady-state heat problem on a computer and the output shows an interior grid cell reading hotter than every boundary cell, that is not a legitimate physical result — it is a red flag that a bug exists somewhere in the code (a sign error, a bad mesh cell, a solver that hasn't converged). The maximum principle is a free, automatic check that costs nothing to run.
The maximum principle is not just a curiosity — it is the cleanest route to uniqueness
for the Dirichlet problem. Suppose two harmonic functions
The same argument gives stability: if two boundary datasets differ by at most
This is exactly why the sanity check above works: two slightly different numerical runs of the same physical setup, with boundary data agreeing to within a tiny tolerance, are guaranteed to agree everywhere inside to within that same tolerance. A simulation that drifts far from a known solution despite matching boundary conditions closely is, again, telling you something is wrong.
It's tempting to assume every "nice" equation from physics obeys the same clean maximum principle. It doesn't. This particular argument leans on the mean-value property, which is special to Laplace's equation (and close relatives, the general class of elliptic equations). It does not carry over in the same simple form to:
The steady-state, no-time-dependence, "everything has already settled down" character of Laplace's equation is exactly what makes the simple version of the maximum principle work.
Physicists and engineers who run numerical solvers for steady-state problems lean on the maximum principle constantly, often without naming it. It's one of the cheapest correctness checks in numerical computing: after solving for a steady-state temperature or potential field, scan the computed interior values and compare them against the range of the boundary values.
If every interior value falls neatly between the smallest and largest boundary values, that's consistent with — though of course not proof of — a correct solution. But if even a single interior grid point pokes above the hottest boundary value or below the coldest one, you don't need to hunt for a subtle physics explanation. Something in the code is broken: a mismatched sign, a misapplied boundary condition, a solver that stopped iterating too early. Mathematics itself is telling you the answer is impossible.