The Bondi problem revisited: A spectral domain decomposition code

Abstract

In this paper, we present a simple domain decomposition code based on the Galerkin-Collocation method to integrate the field equations of the Bondi problem. The algorithm is stable, exhibits exponential convergence when considering the Bondi formula as an error measure, and is computationally economical. We have incorporated features of both Galerkin and Collocation methods along with the establishment of two nonoverlapping subdomains. We have further applied the code to show the decay of the Bondi mass in the nonlinear regime and its power-law late-time decay. Another application is the determination of the waveforms at the future null infinity connected with distinct initial data.