Numerical simulations have been carried out to investigate the formation and motion of single bubble in liquids using volume-of-fluid (VOF) method using the software platform of FLUENT 6.3. Transient conservation mass and momentum equations with considering the effects of surface tension and gravitational force were solved by the pressure implicit splitting operator (PISO) algorithm to simulate the behavior of gas-liquid interface movements in the VOF method. The simulation results of bubble formation and characteristics were in reasonable agreement with experimental observations and available literature results. Effects of fluid physical properties, operation conditions such as orifice diameter on bubble behavior, detachment time, bubble formation frequency and bubble diameter were numerically studied. The simulations showed that bubble size and bubble detachment times are linear functions of surface tension and decrease exponentially with the increase in liquid density. In contrast, only a small influence of the fluid viscosity on bubble size and detachment time was observed. Bubble collapse at a free surface simulation with VOF method was also investigated.