JOREK [1] simulations of mitigated axisymmetric (2D) ITER disruptions with vertical plasma motion (VDE) and runaway electrons (RE) will be presented. An RE fluid model [2] including the effects of partially-ionized impurities [3,4] is used, and the surrounding conducting structures are included via STARWALL coupling. Both hot and cold VDE scenarios are examined with focus on the effects of different possible current profile flattening, efficacy of Neon second injection, and impurity flushout via massive deuterium second-injection onto the RE beam. Further plans for full 3D simulations of RE beam termination will also be briefly discussed. [1] Hoelzl et. al. NF 2021 [2] Bandaru et al. Phys. Rev. E. 2019 [3] Hesslow et al. PPCF 2018 [4] Hesslow et al. NF 2019