Abstract: We present a novel framework based on a continuous fluid simulator for general simulation of realistic bubbles, with which we can handle as many significant dynamic bubble effects as possible. To capture nature of the very thin liquid film of bubbles, we have developed a regional level set method allowing multi-manifold interface tracking. The regional level set method is based on the definitions of regional distance and its five operators, which makes it very easy to implement. We can reconstruct an implicit surface of liquid film with arbitrary thickness from the representation of regional level set functions. To overcome the numerical instability caused by surface tension, we exploit a new semi-implicit surface tension model which is unconditionally stable and makes the simulation of surface tension dominated phenomena much more efficient. An approximated film thickness evolution model is proposed to control the bubble's lifecycle. All these new techniques combine into a general framework that can produce various realistic dynamic effects of bubbles.