In this work we use the IP-Glasma+MUSIC+UrQMD framework to systematically study a wide range of hadronic flow observables at 2.76 TeV. In addition to the single particle spectra and anisotropic flow coefficients $v_n$ previously studied in cite{1609.02958}, we consider event-plane correlations, non-linear response coefficients $chi_{npq}$, and event shape engineering. Taken together, these observables provide a wealth of insight into the collective behavior of the QGP and initial state fluctuations. They shed light on flow fluctuations, flow at fixed system size but different initial geometries, as well as the non-linear hydrodynamic response to the initial state spatial eccentricities. By synthesizing this information we can gain further insight into the transport properties of the QGP as well as the fluctuation spectrum of the initial state.