Add pure diffusion case.

 * Both solver and a convergence test.
 * Write a separate function to plot the flow field.
 * Triple the mesh size automatically.
Check convergence of linearly interpolated field.