In this thesis, through numerical simulations, we investigate the transport of heavy impurities in plasmas where instabilities are suppressed by applied external 3D magnetic fields. As the basis of the model, we build upon the nonlinear extended magnetohydrodynamics (MHD) code JOREK as it can self-consistently simulate the plasma response to magnetic perturbations and simulations of the plasma up to the first wall and divertor. To correctly assess the production of impurities, we need to predict the power fluxes and plasma-wall interactions impacted by MHD processes during operation, including an accurate description of the MHD processes themselves, as well as the edge plasma and plasma-wall interaction processes. Therefore, we have developed a full-f, Monte-Carlo, kinetic neutral model two-way coupled to the plasma model, including recycling of plasma and neutrals at the plasma-facing components. To assess the transport of impurities, we utilize the same particle framework in JOREK to describe the full-orbit transport and the full distribution (full-f) of the impurity species. We have introduced a collision operator that reproduces collisional neoclassical transport in the case of well-defined flux surfaces.
