The paper reports numerical simulation results of a direct current (DC) suspension plasma spray with an axial injection system. In the numerical modelling, a two-way coupled Eulerian-Lagrangian approach was employed to simulate plasma flow and suspension behavior. As effects of two-way interaction, momentum transfer and energy transfer were chosen. The plasma spray was assumed to have a rotationally symmetric, two-dimensional shape around the injection axis of suspensions (the central axis), and therefore the axisymmetric two-dimensional approximation was adopted in the numerical modelling. Working gas was pure argon and was supplied from both the anode-side and the cathode-side ports. A total volume flow rate of the working gas from anode and cathode was set to 46.5 slm. A feed rate of suspension was parametrically set to 0, 15, 25, or 35 g/min. Numerical results indicate that the temperature of a plasma hot-core region and the velocity of a plasma jet around the central axis drop more with increasing feed rate of suspension mainly because of a decrease in Joule heating with increasing it. The numerical results also suggest that the increase in feed rate of suspension leads to practically lengthening flight distance of suspensions required for completing evaporation process of disperse medium.