Particle-laden electrostatic discharges are widely used in research as well as in industrial applications. The mechanism of particle lifting in these highly unsteady flows are not well understood. Modifications are made to the FLASH research code to simulate this multiphase system using the Lagrangian-Eulerian (LE) approach. Various Lagrangian submodels are implemented in the code to study the particle phase evolution in the system. Results generated from the newly developed code are compared to the experiments of the dispersion of glass beads in an electrostatic discharge of energy 2 J performed at the New Jersey Institute of Technology (NJIT). The simulation results are qualitatively in good agreement with the experiments. Furthermore, the prominent mechanism of particle lifting was identified to be the Saffman lifting force generated on the particles due to velocity gradients in the gas flow field. Results suggest that the newly developed code is flexible in modeling particle transportation on larger length scales such as flow generated by “fireball” plasmas.