5 Numerical Hydrodynamics on Null 4 Cauchy-Characteristic Matching4.7 3D Cauchy-Characteristic Matching

4.8 The Binary Black Hole Inner Boundary 

It is clear that the three-dimensional inspiral and coalescence of black holes challenges the limits of present computational know-how. CCM offers a new approach for excising an interior trapped region which might provide the enhanced flexibility to solve this problem. In a binary system, there are major computational advantages in posing the Cauchy evolution in a frame which is co-rotating with the orbiting black holes. Such a description seems necessary in order to keep the numerical grid from being intrinsically twisted. In this co-orbiting description, the Cauchy evolution requires an inner boundary condition inside the black holes and also an outer boundary condition on a worldtube outside of which the grid rotation is likely to be superluminal. An outgoing characteristic code can routinely handle such superluminal gauge flows in the exterior [30]. Thus, successful implementation of CCM would solve the exterior boundary problem for this co-orbiting description.

CCM also has the potential to handle the two black holes inside the Cauchy region. As described earlier with respect to Fig.  5, an ingoing characteristic code can evolve a moving black hole with long term stability [73, 70]. This means that CCM might also be able to provide the inner boundary condition for Cauchy evolution once stable matching has been accomplished. In this approach, the interior boundary of the Cauchy evolution is located outside the apparent horizon and matched to a characteristic evolution based upon ingoing null cones. The inner boundary for the characteristic evolution is a trapped or marginally trapped surface, whose interior is excised from the evolution.

In addition to restricting the Cauchy evolution to the region outside the black holes, this strategy offers several other advantages. Although finding a marginally trapped surface on the ingoing null hypersurfaces remains an elliptic problem, there is a natural radial coordinate system tex2html_wrap_inline2365 to facilitate its solution. Motion of the black hole through the grid reduces to a one-dimensional radial problem, leaving the angular grid intact and thus reducing the computational complexity of excising the inner singular region. (The angular coordinates can even rotate relative to the Cauchy coordinates in order to accommodate spinning black holes.) The chief danger in this approach is that a caustic might be encountered on the ingoing null hypersurface before entering the trapped region. This is a gauge problem whose solution lies in choosing the right location and geometry of the surface across which the Cauchy and characteristic evolutions are matched. There is a great deal of flexibility here because the characteristic initial data can be posed without constraints.

This global strategy is tailor-made to treat two black holes in the co-orbiting gauge, as illustrated in Fig.  6 . Two disjoint characteristic evolutions based upon ingoing null cones are matched across worldtubes to a central Cauchy region. The interior boundary of each of these interior characteristic regions border a trapped surface. At the outer boundary of the Cauchy region, a matched characteristic evolution based upon outgoing null hypersurfaces propagates the radiation to infinity.


Click on thumbnail to view image

Figure 6: CCM for binary black holes, portrayed in a co-rotating frame. The Cauchy evolution is matched across two inner worldtubes tex2html_wrap_inline1679 and tex2html_wrap_inline1681 to two ingoing null evolutions whose inner boundaries excise the individual black holes. The outer Cauchy boundary is matched across the worldtube tex2html_wrap_inline1683 to an outgoing null evolution extending to tex2html_wrap_inline1655 .

Present characteristic and Cauchy codes can handle the individual pieces of this problem. Their unification appears to offer the best chance for simulating the inspiral and merger of two black holes. The CCM module is in place and calibrated for accuracy. The one missing ingredient is its long term stability, which would make future reviews of this subject very exciting.

5 Numerical Hydrodynamics on Null 4 Cauchy-Characteristic Matching4.7 3D Cauchy-Characteristic Matching

image Characteristic Evolution and Matching
Jeffrey Winicour
© Max-Planck-Gesellschaft. ISSN 1433-8351
Problems/Comments to livrev@aei-potsdam.mpg.de