The main aim of this research is focused on determining the velocity and particle density profiles across the flame propagation of microlycopodium dust particles. In this model, it is tried to incorporate the forces acting on the particles such as thermophoretic, gravitational, and buoyancy in the Lagrangian equation of motion. For this purpose, it is considered that the flame structure has four zones (i.e., preheat, vaporization, reaction, and postflame zones) and the temperature profile, as the unknown parameter in the thermophoretic force, is extracted from this model. Consequently, employing the Lagrangian equation with the known elements results in the velocity distribution versus the forefront of the combustion region. Satisfactory agreement is achieved between the present model and previously published experiments. It is concluded that the maximum particle concentration and velocity are gained on the flame front with the gradual decrease in the distance away from this location.