Linearized solutions of the Einstein equations within a Bondi-Sachs framework, and implications for boundary conditions in numerical simulations

We linearize the Einstein equations when the metric is Bondi-Sachs, when the background is Schwarzschild or Minkowski, and when there is a matter source in the form of a thin shell whose density varies with time and angular position. By performing an eigenfunction decomposition, we reduce the problem to a system of linear ordinary differential equations which we are able to solve. The solutions are relevant to the characteristic formulation of numerical relativity: (a) as exact solutions against which computations of gravitational radiation can be compared; and (b) in formulating boundary conditions on the $r=2M$ Schwarzschild horizon.