The coating of surfaces by plasma spraying is an important manufacturing process with many industrial applications. In the last several decades, numerical modeling of plasma spraying processes has met with considerable attention [1,2,3]. That is in order to well understand the complex phenomena the plasma spray involves, for economic constraints and to well predict the plasma-inflight-particles exchanges since this affects directly the coating formability and microstructure. This study deals with the investigation of plasma jets using an axisymmetric LB thermal model. Plasma jets have been very successful in many applications (such as spraying, cutting, welding,…). The excellent choice of high performance plasma gases and spraying materials has been the subject of several experimental and numerical efforts. An excellent choice will be the response of efficient numerical studies and the results of experimental tests. Plasma jets are high temperature flows (>8000K). Therefore, all diffusion parameters involved in conservation equations are temperature dependent. In the following, we present a plasma jet investigation in an axisymmetric LBM (Jian’s model ). In the context of our knowledge, it is the first attempt to tackle this field by using the LBM. Further reading on solution procedure, the model implementation and assumptions may be found in [5,6].