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: 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 models were consistent with the experimental tests for Reynolds numbers < 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 gasliquid interface. The models also predicted the velocity and pressure fields near the bubble surface before and after detachment.