Adjoint-based sensitivity methods are powerful design tools for engineers who use computational fluid dynamics. In recent years, these engineers have started to use scale-resolving simulations like large-eddy simulations (LES) and direct numerical simulations (DNS), which resolve more scales in complex flows with unsteady separation and jets than Reynolds-averaged Navier-Stokes (RANS) methods. However, the conventional adjoint method computes large, unusable sensitivities for scale-resolving simulations, which unlike RANS simulations exhibit the chaotic dynamics inherent in turbulent flows. Sensitivity analysis based on least-squares shadowing (LSS) avoids the issues encountered by conventional adjoint methods, but has a high computational cost even for relatively small simulations. The following talk discusses a new, more computationally efficient formulation of LSS and its application to turbulent flows simulated with Eddy, a discontinuous-Galkerin spectral-element-method LES/DNS solver. First, the new LSS formulation, called “non-intrusive” LSS, is outlined, followed by a cost analysis of the method. Results are presented for the minimal flow unit, a turbulent channel flow with a limited streamwise and spanwise domain.