I am learning a bit of fourier analysis, with an interest in physics as well. I originally posted this question on the math stack exchange, but perhaps you physicists have more experience in these kinds of things. The standard equation to solve is the steady state heat equation (Laplace equation) in the plane is
$$ \frac{\partial^2 f}{\partial x^2} + \frac{\partial^2 f}{\partial y^2} = 0 $$
Now I understand that, on functions with a fixed boundary, the solutions to this equation give the steady heat distribution, assuming that the heat at the boundary is a constant temperature. Thus the harmonic solutions make sense.
However, I am not convinced about the physical intuition of a solution on an open subset of the plane. Since heat should disperse equally, shouldn't the only steady state solutions be constant, if there is no boundary at a fixed temperature. Louiville's theorem for harmonic functions (a bounded entire function is constant) tell me that on the whole plane, heat must eventually disperse equally, but I can't figure out intuition for why this isn't the case on bounded open subsets of the plane. Is there a different equation I should be looking at? I can't seem to find any part of the derivation of the heat equation that uses the fact that a boundary must be fixed.