The rising behavior of single bubble is investigated in 0.134Eo 8.588; 0.82We 4.83 and 0.5Ta 3.6 by means of CLSVOF method and CSF model. The new empirical relations are obtained by comparing experimental data, simulated data and the empirical relations of scholars. The three formulas can describe the bubble deformation more accurate. The new empirical relations provide a lower boundary of scope of application and provided reference to bubble research. Introduction Aerated flow[1] is a main flow phenomenon in high speed flow, the key and difficult problem is how to simulate and measure its complex bubble flow. Due to different bubble deformation and existence periodic length, Bubble formation, growth and movement process[2][3]show a different geometric structures and complicated hydrodynamic characteristics[4][5] Therefore, a thorough study of bubble movement is very important to extend bubble application field. The experiments and theoretical research have been carried out by scholars in the middle of last century, some have tried to correlate the aspect ratio E[6][7]as a function of a dimensionless parameter which could group the experimental data using Eo[8][9], We and Ta[10][11]. Ruzicka[12] studied bubble behavior by asymptotic expansion method, and analyzed the small deformation during the rising process. Raymond[13] investigated the rising behavior of bubble in liquid with different viscosity by model tests, and described bubble shape and aspect ratio using dimensionless parameters. Tomiyama[14] found that terminal velocity and the ratio of length to width is inversely proportional to the relationship. Churchill[15] proposed that a parameter which can describe the bubble deformation played an important academic role. Taylor[16] described bubble deformation using the relation of We, Eo and E. In order to make the research has universal significance, we select a group of characteristic quantity and handle variables by dimensionless method. The empirical expressions of dimensionless parameters are as follows: 2 ( ) f b gd Eo      ; 2 fU d We    ; a E b  ; 0.23 Ta=Re Mo  ; 4