Drop Coalescence Simulations Using Level Sets Coupled with Boundary Integral Methods