This paper deals mainly with the numerical simulation on edge localized mode (ELM) pacing by pellet injection that is useful for fuelling and control of plasma profiles to achieve enhanced tokamak operations. The fuelling and pellet-induced ELMs are simulated with a 1.5-dimensional core transport code, which includes a neutral gas shielding model and a grad-B drift model for pellet deposition in H-mode tokamak plasmas. Fuelling and ELM pacing experiments by pellet injections at JET are introduced as a current experimental approach. For the description of ELM triggering by pellet injection based on ideal ballooning mode criteria, three possible models are suggested and discussed on their ELM characteristics, respectively: (i) the density enhanced ELMs in the post-pellet phase, (ii) the modification of the surface averaged pressure profiles in a transport time scale and (iii) the local increase in the pressure (density and/or temperature) gradients perturbed by pellets. Among them, the pellet-induced density perturbation model is adopted, in practice, to carry out an ELM pacing simulation in preparation for future experiments in KSTAR. The numerical simulation shows that the artificially induced ELM by pellets releases the reduced energy bursts, compared with spontaneous ELMs. The energy loss per burst by the pellet-induced ELM turns out to be much smaller than that by the spontaneous ELM as the pellet injection frequency becomes higher in ELM pacing. Based on the simulation results showing good agreement with the general ELM characteristics observed in pellet pacing experiments, the ELM pacing by pellet injection is very promising for mitigating the ELM energy bursts to the divertor by controlling the injection frequency.