Here is my understanding from a pure math point of view. But I am not sure if this is what you want for the ''intuition''.
In analysis and PDE, that $\Delta u=f$ says that $u$ ''gain 2 derivatives than'' $f$.
Assuming $u$ and $f$ are both defined on $\mathbb R^n$ and has at most polynomial growth, if we take Fourier transform on both sides then we have $-4\pi^2|\xi|^2\hat u(\xi)=\hat f(\xi)$. Basically it says that if $\hat f(\xi)=O(|\xi|^{-m}$ then $\hat u(\xi)=O(|\xi|^{-m-2})$ decay two more degree faster than $\hat f$ as $\xi\to\infty$.
The decay of Fourier coefficients is equivalent as saying the smoothness of the original functions. By taking Fourier transform back we see that this is saying $u$ is smoother than $f$ with order 2.
To get a better idea, we can consider $L^2$ Sobolev spaces. By the so-called Plancherel identity if $f\in W^{k,2}(\mathbb R^n)$ if and only if $\int (1+|\xi|^2)^{k/2}|\hat f(\xi)|^2d\xi<\infty$. If we igonore the "+1" in the coefficient, i.e. the low frequency, this is saying that $u\in W^{k,2}(\mathbb R^n)$.
In practice we can also make the same description locally, namely if $f$ ''looks nice near $x_0$'' then $u$ ''looks 2-order nicer near $x_0$''. This is what the Schauder's estimate comes into play.
Certainly thinking that, if $f$ is order $k$ differentiable (i.e. $C^k$) near $x_0$ then $u$ is order $k+2$ differentiable (i.e. $C^{k+2}$) near $x_0$, is a good starting point. Although that is not true in general. If you replace $C^k$ by Sobolev regularity $W^{k,p}$ for $1<p<\infty$ or H"older regularity $C^{k,\alpha}$ for $0<\alpha<1$, it will be correct.
I think checking the word elliptic differential operator can be a further path to get the intuition.