<p>In this study, a two-dimensional axisymmetric model is established using OpenFOAM. The bubble interface is captured using the volume of fluid method, where the liquid’s viscoelastic behavior is characterized by the Giesekus model. The authors systematically investigated bubble formation in viscoelastic liquids, focusing on the effects of rheological parameters (such as zero-shear viscosity, relaxation time, and mobility factor) and operating conditions (gas injection velocity and orifice radius) on bubble evolution (such as aspect ratio, contact line diameter, and duration of each stage) and departure behaviors (including departure time and volume). By analyzing the evolution of bubble contact line diameter and bubble height, bubble formation process can be divided into five stages: (1) nucleation, (2) expansion, (3) elongation, (4) retraction, and (5) departure. The present results indicate that mobility factor exerts a comparatively minor impact on bubble evolution and departure. In contrast, with increasing relaxation time, the apparent viscosity and the elastic stress of liquid decrease significantly, leading to a marked change in bubble morphology. Specifically, the bubble tail gradually changes from a pointed shape to a concave shape, and the gas line also presents at the bubble tail. Both the contact line diameter and the maximum pre-departure bubble height decrease with increasing relaxation time, but they increase with increasing zero-shear viscosity, injection velocity, and orifice radius. Departure time increases only with zero-shear viscosity and decreases with all other parameters studied. Bubble departure volume decreases solely with increasing relaxation time but it grows with an increase in the other parameters.</p>