Multi-messenger astronomy promises to revolutionize our understanding of the universe by providing dramatically contrasting views of the same objects. Black holes are the most fascinating and ubiquitous objects in the universe. Unlike stellar-mass black holes, it is generally expected that supermassive black holes lurking at the center of most galaxies will merge in a gas-rich environment, because galaxies merge from time to time. In this talk, I will review the history of simulation efforts to model binary black holes in their astrophysical environments and also present some exciting new results in the context of magneto-hydro- dynamics simulation.