Numerical modeling of reactive bubbly flows based on Euler-Lagrange approach