Modeling discrete gas bubble formation and mobilization during subsurface heating of contaminated zones