Stabilized FIF/FEM formulation for multidimensional advection-diffusion-reaction problems

A stabilized finite method (FEM) for the multidimensional steady state advection-diffusion-absorption equations is presented. The stabilized formulations is based on the modified governing differential equations derived via the Finite Calculus (FIC) method. For 1D problems the stabilization terms act as a nonlinear additional difffusion governed by a single stabilization parameter. It is shown that for multidimensional problems an orthotropic stabilizing diffusion must be added along the principal directions of curvature of the solution. A simple iterative algorithm yielding a stable and accurate solution for all the range of physical parameters and boundary conditions is described. Numerical results for 1D and 2D problems with sharp gradients are presented showing the effectiveness and accuracy of the new stabilized formulation