Abstract
We consider the Mizukami-Hughes method for the numerical solution of scalar two-dimensional steady convection-diffusion equations using conforming triangular piecewise linear finite elements. We propose several modifications of this method to eliminate its shortcomings.
The improved method still satisfies the discrete maximum principle and gives very accurate discrete solutions in convection-dominated regime, which is illustrated by several numerical experiments.