Regimes of operation in tokamaks that are devoid of large ELMs have to be better understood to extrapolate their applicability to reactor-relevant devices. This paper describes non-linear extended MHD simulations that use an experimental equilibrium from an EDA H-mode in ASDEX Upgrade. Linear ideal MHD analysis indicates that the operational point lies slightly inside of the stable region. The non-linear simulations with the visco-resistive extended MHD code, JOREK, sustain non-axisymmetric perturbations that are linearly most unstable with toroidal mode numbers of n = \{6 \dots 9\}, but non-linearly higher and lower n become driven and the low-n become dominant. The poloidal mode velocity during the linear phase is found to correspond to the expected velocity for resistive ballooning modes. The perturbations that exist in the simulations have somewhat smaller poloidal wavenumbers (k_{\theta} \sim 0.1 to 0.5 cm^{-1} ) than the experimental expectations for the quasi-coherent mode in EDA, and cause non-negligible transport in both the heat and particle channels. In the transition from linear to non-linear phase, the mode frequency chirps down from approximately 35 kHz to 13 kHz, which corresponds approximately to the lower end of frequencies that are typically observed in EDA H-modes in ASDEX Upgrade.