Front tracking provides sharp resolution of wave fronts through the active
tracking of interfaces between distinct materials. A major challenge to thi
s method is to handle changes in the interface topology. We describe two al
gorithms, implemented in the front tracking code FronTier, to model dynamic
changes in three-dimensional interfaces. The two methods can be combined t
o give a hybrid method that is superior to each individual method. The succ
ess of these algorithms is shown by simulations of Rayleigh Taylor instabil
ity, which is an interfacial instability driven by an acceleration directed
across a material interface. Our numerical results are validated by compar
ing the numerical computation of the velocity of a single rising bubble wit
h an analytic model for the bubble velocity.