A mathematical model that combines turbulent phenomena of particle-laden gas jets and chemical reactions has been developed to predict gas mixing and solid reaction in a two-dimensional cylindrical system. A simplified gas-solid reaction model and isothermality were introduced to simplify model predictions. The gas field was viewed from an Eulerian frame work. The k-ε model was used for the prediction of turbulent behavior. The SIMPLER algorithm was also used for the correction of pressure gradient in the gas-phase momentum equation with the continuity equation. Lagrangian treatment of the particles was performed, pre-senting the particle field as a series of the trajectories. Particle dispersion in the turbulent gas field is modeled through the particle number density. The application of the developed mathematical model can be found in the atomization process, flash-smelting process, carbothermal reduction process, and the oxidation process of inorganic materials for the production of pigments.