Using the nonlinear 3D MHD code JOREK with reduced MHD equations (visco-resistive MHD), we have successfully simulated a recent n??????=??????2 resonant magnetic perturbation (RMP)-driven edge localized mode (ELM) suppression in KSTAR. We have found that such ELM suppression has been attributable not only to the degraded pedestal but also to the direct coupling between the peeling-ballooning mode (PBM) and RMP-driven plasma response. Notably, the pedestal pressure gradient is reduced as the radial transport is enhanced because of the formation of the stochastic layer and increased convection fluxes due to tearing and the kink-peeling mode driven by RMPs. The increased transport in the stochastic layer is due to the parallel transport across...