  Lax-Wendroff Flux Reconstruction (LWFR) is a single-stage, high order,
quadrature free method for solving hyperbolic conservation laws. We perform a
cell average decomposition of the LWFR scheme that is similar to the one used
in the admissibility preserving framework of Zhang and Shu (2010). By
performing a flux limiting of the time averaged numerical flux, the
decomposition is used to obtain an admissibility preserving LWFR scheme. The
admissibility preservation framework is further extended to a newly proposed
extension of LWFR scheme for conservation laws with source terms. This is the
first extension of the high order LW scheme that can handle source terms. The
admissibility and accuracy are verified by numerical experiments on the Ten
Moment equations of Livermore et al.
