Discontinuous Galerkin (DG) methods combine features in finite element methods (weak formulation, finite dimensional solution and test function spaces) and in finite volume methods (numerical fluxes, nonlinear limiters) and are particularly suitable for simulating convection dominated problems, such as linear and nonlinear waves including shock waves. In this article we will give a brief survey of DG methods, emphasizing their applications in computational fluid dynamics (CFD). We will discuss essential ingredients and properties of DG methods, and will also give a few examples of recent developments of DG methods which have facilitated their applications in CFD.

An important property for finite difference schemes designed on curvilinear meshes is the exact preservation of free-stream solutions. This property is difficult to fulfill for high order conservative essentially non-oscillatory (WENO) finite difference schemes. In this paper we explore an alternative flux formulation for such finite difference schemes which can preserve free-stream solutions, based on the numerical technique for the metric terms, which can be applied to this alternative flux formulation but is difficult to be applied to the standard finite difference formulation. Free-stream and vortex preservation properties are investigated, and comparison with standard finite difference WENO schemes is made. Theoretical derivation and numerical results show that the finite difference WENO schemes based on the alternative flux formulation can preserve free-stream and vortex solutions on both stationary and dynamically generalized coordinate systems, hence giving much better performance than the standard finite difference WENO schemes for such problems.

We are interested in front propagation problems in the presence of obstacles. We extend a previous work (Bokanowski, Cheng and Shu 09}), to propose a simple and direct discontinuous Galerkin (DG) method adapted to such front propagation problems. We follow the formulation of Bokanowski et al. (Bokanowski_Forcadel_Zidani_2010), leading to a level set formulation driven by $\min(u_t + H(x,\nabla u), u-g(x))=0$, where $g(x)$ is an obstacle function. The DG scheme is motivated by the variational formulation when the Hamiltonian $H$ is a linear function of $\nabla u$, corresponding to linear convection problems in presence of obstacles. The scheme is then generalized to nonlinear equations, written in an explicit form. Stability analysis are performed for the linear case with Euler forward, a Heun scheme and a Runge-Kutta third order time discretization using the technique proposed in Zhang and Shu 2010. Several numerical examples are provided to demonstrate the robustness of the method. Finally, a narrow band approach is considered in order to reduce the computational cost.