Numerical Investigation of Convection Diffusion Reaction Systems This work is concerned with the numerical solution of selected convection-diffusion-reaction (CDR) type mathematical models with dominating convective and reactive terms, coupled with some algebraic equations. Five established CDR-type models are analyzed namely, the gas-solid reaction, chemotaxis, liquid chromatography, radiation hydrodynamical, and hy- perbolic heat condition models. These models are encountered in various scientific and engi- neering fields, such as chemical engineering, biological systems, astrophysics, heat transfer, and fluid dynamics. The Laplace transformation is applied as a basic tool to find the ana- lytical solutions of linear CDR models for different types of boundary conditions. However, for the nonlinear models, numerical techniques are the only tools to get physical solutions. The nonlinear transport and stiff source (reaction) terms of the governing differential equa- tions produce discontinuities and narrow peaks in the solution. It is difficult to capture steep variations in the solution through a less accurate numerical scheme. Therefore, ef- ficient and accurate numerical methods are needed to obtain physically reliable solutions in acceptable computational time. The objective of this thesis project is to develop and implement simpler, robust, and accurate numerical frameworks for the solution of one and two-dimensional CDR type systems. The space-time CE/SE-method, the discontinuous Galerkin (DG) finite element method , and different high resolution finite volume schemes (FVSs) are proposed to numerically approximate the solution of these models. Several case studies are carried out. The validity and performance of the suggested numerical techniques are revealed through test problems and by comparing their results with each other, analytical solutions, and the results of some available finite volume schemes in the literature.