When spraying ceramic particles with a low thermal conductivity such as zirconia using Ar-H 2 direct-current (d.c.) plasma jets where the heat transfer is important, heat propagation phenomena take place with the propagation of melting, evaporation or even solidification fronts. Most models neglect these heat propagation phenomena assuming the particle as a lumped media. This work is aimed at developing a model coupling the effect of heat propagation with the particle dynamic within plasma jets. It uses an adaptative grid in which the coordinates of the phase change fronts are fixed. It allows minimizing the calculation costs (approximately 10 seconds on PC under windows XP against 1hour with an enthalpy model). Such calculations are illustrated for dense and porous agglomerated zirconia as well as iron particles which evaporation in an Ar- H 2 (25 vol %) plasma is important.