This paper is devoted to the presentation of a modeling approach of DC arc plasma flows by implementation of the CFD PHOENICS code. The equations of mass, momentum and energy, of elliptic type, are discretized using a finite/control volumes method. The closure of the equations set is obtained with the standard k-ε turbulence model. Finally, the algebraic equations set is solved by means of the SIMPLE algorithm. 3D calculations are performed from the point of the cold gas injection into the torch to the plasma ejection into the surrounding atmosphere, taking into account some currently used spray parameters (i.e., plasma gas flow rate, power, etc.). Results are discussed and compared with experimental values taken from the literature.