I have implemented an ADER-Discontinuous Galerkin scheme for the resolution of linear systems of conservation laws of the type of $\partial_t U + A \partial_x U + B \partial_y U=0 $ and observed that the CFL condition is very restrictive. In the bibliography, an upper bound for the time step $\Delta t \leq \frac{h}{d(2N+1)\lambda_{max}}$ can be found, where $h$ is the cell size, $d$ is the number of dimensions and $N$ is the max degree of the polynomials.Is there any way to circumvent this issue? I had been working with WENO-ADER Finite volume schemes and the CFL restrictions were much more relaxed. For instance, for a 5-th order scheme, a CFL lower than 0.04 must be imposed when using DG while CFL=0.4 can still be used in a WENO-ADER FV scheme.Why using DG schemes rather than ADER-FV, for instance, in computational aeroacustics (linearized Euler equations) or similar applications (gas dynamics, shallow water, magnetohydrodynamics)? Is the overall computational cost of the scheme similar than that of the ADER-FV, in spite of the much lower time step?Thoughts and suggestions for this are welcome. 
CFL condition in Discontinuous Galerkin schemes
fluid dynamics;finite volume;discontinuous galerkin;numerical modelling
The restrictive CFL of DG schemes typically comes from the combination of high order accuracy and a compact stencil (see this reference for example).  The CFL depends on bounding the variational form in terms of the $L^2$ norm of the solution, which depends on derivative and traces of polynomials.  Bounds for each of these quantities (using Bernstein or Markov brothers inequalities and discrete trace inequalities) give constants which depend inversely on $h$ and quadratically on the order $N$, resulting in an overall CFL of $O(h/N^2)$.FYI - I've seen the CFL you mention referenced before, but I can't recall where it's proven. I'd like to know how they avoid quadratic dependence on $N$ in their bound. Finite difference and WENO schemes (as well as B-spline based finite element methods on periodic meshes) have looser CFL conditions because the constants in the analogous bounds grow more slowly in $N$. This is in turn because the stencil size tends to increases with order $N$, which reduces some of these issues. DG methods are more expensive, but they can deal easily with unstructured meshes and can be implemented efficiently. There are high order versions of WENO (or similar reconstructions) for unstructured grids, though these can introduce additional mathematical or implementational complications.