The major issue addressed in this thesis is the development of a Mach-uniform method that uses a staggered scheme on planar unstructured grids with triangular cells, that has a superlinear rate of convergence. In our staggered grid, the normal velocity components are associated with cell faces, whereas scalar variables are assigned to cell centers. In order to obtain a superlinearly convergent scheme, it is necessary to interpolate vector fields from staggered components with sufficient accuracy in a robust way. Several methods to achieve this objective are investigated. Like all higher order methods, the superlinearly convergent scheme obtained suffers from unphysical oscillations near discontinuities such as shocks. This difficulty is ove...