In present study a mathematical model of bubble interaction with fluid flow based on the physical observations of the zero-shear-stress and no separation on the surface of gaseous clean bubble is introduced. The potential theory is employed to derive the model. This theory can be used easily in two-dimensional and also three dimensional axial symmetry analyses. Also the rotation of bubble in non-uniform flow is modeled by generating a vortex blob. This model explains that the bubbles near the solid boundaries generate more vorticity than when they are far from the wall. Bubble dynamics comprises the forces of drag, lift, added mass, pressure gradient and buoyancy. The basset force is neglected in simulation of bubble because this force is a historical force and acts just in a small period of time. A model is also used to consider the effect of wall in two directions of normal and parallel to the wall.