This work presents a methodology to address the problems of bubbles and drops evolving in incompressible, viscous flows due to the effect of surface tension. It is based on the combination of a recently developed immersed interface method to resolve discontinuities, with the level set (LS) method to reproduce the evolving interfaces. The paramount feature of this immersed interface method is the use of Lagrange interpolation enclosing grid points positioned in the vicinity of the interface and few exceptional grid points positioned on the interface. Different problems are considered to assert the accurateness of the proposed methodology, involving both simple and complex interface geometries. Precisely, the following problems are addressed: circular flow with a fixed interface, the dispersion of capillary waves, initially circular, ellipse, star shaped bubbles oscillating to an equilibrium state and circular drops deforming in shear flows. The transient evolution of bubbles/drops in terms of their shapes, pressure profiles, velocity vectors, deformation ratios of major and minor axis is analyzed to observe the effect of surface tension. The proposed methodology is seen to recover the exact numerical equilibrium between the surface tension and pressure gradient in the vicinity of complex interface geometries as well, while recreating the flow physics with an adequate level of accuracy with well representation of overall trends. Moreover, the numerical results yield a good level of agreement with the reference data. Keywords immersed interface level set (LS) surface tension multiphase Navier-Stokes equations shear flows HOC