Numerical simulations of gas puff imaging using a multi-component model of the neutral–plasma interaction in the tokamak boundary