The understanding of flow problems, and finding their solution, has been important for most of human history, from the design of aqueducts to boats and airplanes. The use of physical miniature models and wind tunnels were, and still are, useful tools for design, but with the development of computers, an increasingly large part of the design process is assisted by computational fluid dynamics (CFD).

Many industrial CFD codes have their origins in the 1980s and 1990s, when the low order finite volume method (FVM) was prevalent. Discontinuous Galerkin methods (DGM) have, since the turn of the century, been seen as the successor of these methods, since it is potentially of arbitrarily high order. In its lowest order form DGM is equivalent to FVM. However, many existing codes are not compatible with standard DGM and would need a complete rewrite to obtain the advantages of the higher order.

This thesis shows how to extend existing vertex-centered and edge-based FVM codes to higher order, using a special kind of DGM discretization, which is different from the standard cell-centered type. Two model problems are examined to show the necessary data structures that need to be constructed, the order of accuracy for the method, and the use of an hp-adaptation scheme to resolve a developing shock. Then the method is further developed to solve the steady Euler equations, within the existing industrial Edge code, using acceleration techniques such as local time stepping and multigrid.

With the ever increasing need for more efficient and accurate solvers and algorithms in CFD, the modified DGM presented in this thesis could be used to help and accelerate the adoption of high order methods in industry.