Enriched Finite Elements Method for Convevtion-Diffusion Problems

Abstract

In this thesis, we consider stabilization techniques for linear convection-diffusionreaction (CDR) problems. The survey begins with two stabilization techniques: streamline upwind Petrov-Galerkin method (SUPG) and Residual-free bubbles method (RFB). We briefly recall the general ideas behind them, trying to underline their potentials and limitations. Next, we propose a stabilization technique for one-dimensional CDR problems based on the RFB method and particularly designed to treat the most interesting case of small diffusion. We replace the RFB functions by their cheap, yet efficient approximations which retain the same qualitative behavior. The approximate bubbles are computed on a suitable sub-grid, the choice of whose nodes are critical and determined by minimizing the residual of a local problem. The resulting numerical method has similar stability features with the RFB method for the whole range of problem parameters. We also note that the location of the sub-grid nodes suggested by the strategy herein coincides with the one described by Brezzi and his coworkers. Next, the approach in one-dimensional case is extended to two-dimensional CDR problems. Based on the numerical experiences gained with this work, the pseudo RFBs retain the stability features of RFBs for the whole range of problem parameters. Finally, a numerical scheme for one-dimensional time-dependent CDR problem is studied. A numerical approximation with the Crank-Nicolson operator for time and a recent method suggested by Neslitürk and his coworkers for the space discretization is constructed. Numerical results confirm the good performance of the method.