The flux reconstruction (FR) method for modeling high Reynolds’ (Re) number flows is studied. First, we show that stability is not trivial for under-resolved cases and split forms are a promising technique to address this issue. We use kinetic energy preserving split forms and show their stability, as well as demonstrate the utility of FR, in doing implicit large eddy simulations (ILES) using such split forms. We then introduce wall modeling using an algebraic wall function approach and show that with this we can capture flow characteristics well. Finally, we simulate a high Re case of flow over a multi-element airfoil.