In last weeks blog post "A convection toy problem or can you find the diffusion" I showed the picture of a function propagating through the time. Although the analytical solution of the problem suggests that the shape of the peak remains stable as a function of time, the peak obviously broadens indicating that (numerical) diffusion occurs.
If we go back one step and use a central difference scheme for the discretization of the first order spatial derivative we obtain the following result
We clearly see that during propagation the peak does not only change its shape but that scheme becomes even unstable. If we can analyse this behaviour perhaps we get an answer to our first problem, namely how numerical diffusion is introduced into the stable propagation when using an upwind scheme.
We use a Taylor expansion for the following equation
Using the last but one equation with the last one we obtain