Vapour shielding at solid walls is expected for a future fusion reactor in case of large heat load due to transient events. In order to estimate the dissipation of the incoming heat flux by the vapour shielding, a simulation code called PIXY, using a weighted particle model, is developed. First, a validation study is undertaken. Vapour shielding of an aluminium‐coated tungsten specimen was experimentally observed in a plasma gun device. Results are reproduced by the PIXY simulation. Then, the code is applied to the transient edge‐localized mode (ELM) heat flux for the reactor condition. It was found that the ion–neutral collision, as well as radiative cooling, is a dominant phenomenon for the heat flux dissipation in both cases. As the maximum heat flux increases, oscillating phenomena appear remarkably, which have been observed in other experimental studies on vapour shielding.