Eulerian–Lagrangian method for three-dimensional thermal reacting flow with application to coal gasifiers