Modeling the merger of binary black holes requires numerical simulations. Solving the equations of general relativity on the computer, however, presents unique challenges. I will describe some of these challenges and will discuss how these problems were solved. I will then review some recent simulations at the coalescence and merger of binary black holes, and will discuss some surprising and unexpected astrophysical phenomena that these simulations revealed.