The smile is there exactly because the model is wrong.
The reason it's used though (despite being wrong) is that it provides a convenient space to look at the underlying - the vol*
The (undiscounted) value of an option is given by:
$$ \int_0^\infty \mathrm{PDF}(s) (s-K)^+ \mathrm{d}s $$
where $\mathrm{PDF}(x)$ is the real probability distribution of the underlying. This is model independent.
Under BS, the value is the following:
$$ \int_0^\infty \mathrm{PDF_{LN}}(F,\sigma)(s) (s-K)^+ \mathrm{d}s $$
where $\mathrm{PDF_{LN}}(F,\sigma)$ is the pdf of a lognormal variable with an expected value of $F$ and vol of $\sigma$.
So now, we just need to make them match - and we have one parameter to solve for - $\sigma$. This gives us the BS vol. If the distribution truly were lognormal, you'd obtain the same vol everywhere. Since it's not, you get a changing vol.
*BS vol.