In this work, a numerical simulation model was developed to investigate the complex heat and mass transfer process inside the molten pool during plasma surfacing. The temperature distribution and fluid flow were employed to depict the evolution of the molten pool. The findings revealed that the maximum temperature is located on the surface of the cladding layer corresponding to the position of heat source, and the temperature distribution follows a Gaussian distribution along the scanning direction. Thermal accumulation is obvious due to the continuous energy input in the initial stage, as a result, the maximum temperature and fluid flow velocity gradually increase with the deposition process. The maximum temperature and fluid flow velocity at 4.0 s are 1893 K and 0.281 m/s in case 110 A, respectively. The predicted shapes and dimensions are consistent with the results of the deposition experiments, with a maximum error of no more than 15 %. In order to investigate the impact of WC particles on the solidification microstructure of nickel composite layers, corresponding plasma surfacing experiment was implemented. The solidification characteristics of the microstructure follow planar-cellular-columnar-equiaxed crystals in the bonding layer and WC particles-columnar-equiaxed dendritic crystals in hard layer, respectively. Furthermore, the Ni dendritic near the retained WC particles transformed into columnar crystals due to temperature gradient in the local area.