FDM/FEM for nonlinear convection-diffusion-reaction equations with Neumann boundary conditions—Convergence analysis for smooth and nonsmooth solutions

This paper aims to present in a systematic form the stability and convergence analysis of a numerical method defined in nonuniform grids for nonlinear elliptic and parabolic convection-diffusion-reaction equations with Neumann boundary conditions. The method proposed can be seen simultaneously as a finite difference scheme and as a fully discrete piecewise linear finite element method. We establish second convergence order with respect to a discrete H1-norm which shows that the method is simultaneously supraconvergent and superconvergent. Numerical results to illustrate the theoretical results are included.