In this study, higher-order flux reconstruction methods are used to simulate the flow around a medium and large-radius VFE-2 delta wing using PyFR. The objective of the study was to validate the ability of higher-order methods to accurately predict blunt leading edge separation and adequately resolve vortex dynamics and to compare this approach to standard Reynolds-Averaged Navier-Stokes (RANS) approaches. Wall-resolved implicit Large Eddy Simulations with fourth-order spatial and temporal accuracy were validated against experimental data for a medium-radius wing at a Reynolds number of 1 · 106 with very good agreement. Higher-order simulations of a large-radius wing at a Reynolds number of 40,000 were then compared to RANS simulations with two turbulence models. Noticeable differences in the computed flow structures and turbulence quantities were observed between the three methods along with minor differences in the predicted surface pressure distribution.