This study is directed at a rigorous characterization of the consistency and convergence of discontinuous finite element schemes formulated using Flux Reconstruction (FR). We show that the FR formulation is consistent for linear advection and converges to the exact solution for any scheme that is stable in the von Neumann sense. The numerical solution for a scheme of polynomial order P is composed of P + 1 eignemodes, of which, one and exactly one is 'physical' such that it exhibits the analytical dispersion behavior in the limit of asymptotic grid resolution. The remaining P modes are 'spurious' such that the fraction of energy received by them from the initial condition vanishes in the asymptotic limit. On grid refinement, the rate of convergence of the numerical solution is a function of time, starting from a short-time rate at t = 0 + , associated with interpolation, and asymptotically approaching a long-time rate as t → ∞, associated with numerical differentiation. Both these rates can be inferred directly from the eigensystem of the numerical derivative operator. We verify these analytical expectations using simple experiments in 1-D and 2-D.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.