Two numerical models for studying the dynamics of formation and rise of single bubbles in high viscosity ionic liquids were implemented using the Level-Set method. The models describe two stages of bubble dynamics separately: bubble formation at the inlet nozzle and bubble displacement across the column. The models were experimentally validated through a laboratory scale bubble column using water-glycerol mixtures and two Imidazolium-type ionic liquids. The results indicated that the models are consistent with the experimental tests, particularly for Reynolds below 5; outside this range, the models tend to underestimate the bubble terminal velocity, which can be explained by the effect of the high velocity and pressure gradients close to the gas-liquid interface. Furthermore, the models predicted the velocity and pressure fields near the bubble surface before and after detachment.