The Flux Reconstruction Method with Lax–Wendroff Type Temporal Discretization for Hyperbolic Conservation Laws

In this paper, we develop a Lax–Wendroff type time discretization method for high order Flux Reconstruction scheme to solve hyperbolic conservation laws. Through Cauchy–Kowalewski procedure, the resulting Lax–Wendroff Flux Reconstruction (LWFR) scheme is an alternative spatial–temporal coupling method to the popular Runge–Kutta Flux Reconstruction (RKFR) scheme. LWFR is a one-step explicit high order discontinuous finite element method and its discretization procedure is more compact and effective for certain problems than that of RKFR. Furthermore, aiming at accurate simulation of discontinuity, we propose a robust local artificial viscosity formulation of LWFR for the first time. A collection of successful numerical experiments show that LWFR can give essentially non-oscillatory and sharp solutions for discontinuity, and maintain designed order of accuracy for smooth regions, both in one-dimensional and two-dimensional Euler equations. In conclusion, LWFR scheme is cost effective and accuracy-preserving for certain problems and can be an attractive candidate to solve hyperbolic conservation laws.