Bubble formation from an orifice submerged in quiescent polyacrylamide aqueous solution was investigated numerically with a sharp-interface coupled level-set/volume-of-fluid method based on the rheological characteristics of the fluid. In both non-Newtonian fluids and Newtonian fluids, the numerical approach was able to capture accurately the deformation of the bubble surface, validated by comparison with experimental results. The effects of orifice diameter, solution mass concentration, and gas flow rate on bubble volume and aspect ratio were evaluated. Both the instantaneous and detached volume decrease with the orifice diameter but increase with mass concentration and gas flow rate. The aspect ratio at the departing point tends to rise with the orifice diameter and mass concentration and falls with the gas flow rate.