Discontinuous propagating fronts: linear and nonlinear systems The aim of this study is controlling of spurious oscillations developing around discontinuous solutions of both linear and non linear wave equations or hyperbolic partial differential equations (PDEs). The equations include both first order and second order (wave) hyperbolic systems. In these systems even smooth initial conditions or smoothly varying source (load) terms could lead to discontinuous propagating solutions (fronts). For the first order hyperbolic PDEs the concept of central high resolution schemes is integrated with the multiresolution based adaptation to capture properly both discontinuous propagating fronts and effects of fine scale responses on those of larger scales in the multiscale manner. This integration leads to using central high resolution schemes on non uniform grids; however such simulation is unstable as the central schemes are originally developed to work properly on uniform cells grids. Hence the main concern is stable collaboration of central schemes and multiresoltion based cell adapters. Regarding central schemes the considered approaches are: 1) Second order central and central upwind schemes; 2) Third order central schemes; 3) Third and fourth order central weighted non oscillatory schemes (central WENO or CWENO); 4) Piece wise parabolic methods (PPMs) obtained with two different local stencils. For these methods corresponding (nonlinear) stability conditions are studied and modified as well. Based on these stability conditions several limiters are modified developed as follows: 1) Several second order limiters with total variation diminishing (TVD) feature 2) Second order uniformly high order accurate non oscillatory (UNO) limiters 3) Two third order nonlinear scaling limiters 4) Two new limiters for PPMs. Numerical results show that adaptive solvers lead to cost effective computations (e.g. in some 1 D problems number of adapted grid points are less than 200 points during simulations while in the uniform grid case to have the same accuracy using of 2049 points is essential). Also in some cases it is confirmed that fine scale responses have considerable effects on higher scales. In numerical simulation of nonlinear first order hyperbolic systems the two main concerns are: convergence and uniqueness. The former is important due to developing of the spurious oscillations the numerical dispersion and the numerical dissipation. Convergence in a numerical solution does not guarantee that it is the physical real one (the uniqueness feature). Indeed a nonlinear systems can converge to several numerical results (which mathematically all of them are true). In this work the convergence and uniqueness are directly studied on non uniform grids cells by the concepts of local numerical truncation error and numerical entropy production respectively. Also both of these concepts have been used for cell grid adaptations. So the performance of these concepts is also compared by the multiresolution based method. Several 1 D and 2 D numerical examples are examined to confirm the efficiency of the adaptive solver. Examples involve problems with convex and non convex fluxes. In the latter case due to developing of complex waves proper capturing of real answers needs more attention. For this purpose using of method adaptation seems to be essential (in parallel to the cell grid adaptation). This new type of adaptation is also performed in the framework of the multiresolution analysis. ...