The output shows the oscillation introduced by large Peclet number. Essentially it is caused by unsymmetry of the equation, induced by advection velocity terms. I’ll try to revisit this phenomenon later.