Abstract:We show how the use of the normal projection of the Einstein tensor as a set of boundary conditions relates to the propagation of the constraints, for two representations of the Einstein equations with vanishing shift vector: the ADM formulation, which is ill posed, and the Einstein-Christoffel formulation, which is symmetric hyperbolic. Essentially, the components of the normal projection of the Einstein tensor that act as non-trivial boundary conditions are linear combinations of the evolution equations with the constraints that are not preserved at the boundary, in both cases. In the process, the relationship of the normal projection of the Einstein tensor to the recently introduced ``constraint-preserving'' boundary conditions becomes apparent.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to correctly handle the problem of constraint propagation in the initial - boundary value problem of Einstein's equations in numerical simulations. Specifically, the paper explores how to use the projection of the Einstein tensor on the boundary as a boundary condition to ensure that the constraint conditions are satisfied at the boundary, thereby reducing the phenomenon of constraint violation in numerical simulations.
### Background and Motivation
In numerical relativity, solving the initial - boundary value problem of the vacuum Einstein equations is a common task. The characteristic of this problem is the existence of four additional differential equations, which usually appear in the form of constraints, that is, they do not contain the time derivatives of the basic variables. These four constraint equations represent that the projection of the Einstein tensor along the normal of the time slice is zero, that is, \( G_{ab}n^a = 0 \).
However, in actual numerical simulations, even if the initial data satisfy these constraint conditions, as time evolves, the constraint conditions often gradually become invalid. Early numerical simulations usually solved this problem by forcing the constraint conditions at each time step, but this method is inefficient and difficult to implement. In recent years, researchers are more inclined to perform unconstrained simulations, that is, ignoring the constraint equations and only stopping the simulation when the growth of constraint violation is too large.
### The Core Problem of the Paper
The core problem of the paper is to explore how to use the projection of the Einstein tensor on the boundary \( G_{ab}e^a \) as a boundary condition to ensure that the constraint conditions are satisfied at the boundary. Specifically, the author proves that under two different representations of Einstein's equations, namely the ADM representation and the Einstein - Christoffel representation, the projection of the Einstein tensor on the boundary can be regarded as non - trivial boundary conditions, and these boundary conditions are closely related to the propagation of the constraint conditions.
### Main Findings
1. **ADM Representation**:
- In the ADM representation, although the evolution equations are ill - conditioned, the propagation of the constraint conditions is well - conditioned. The author proves that a part of the projection of the Einstein tensor on the boundary \( G_{ab}e^a \) can be regarded as non - trivial boundary conditions, and these conditions are closely related to the propagation of the constraint conditions.
- Specifically, \( G_{x}^t \) and \( G_{x}^x \) are related to the Hamiltonian constraint \( C \) and the momentum constraint \( C_x \), while \( G_{x}^y \) and \( G_{x}^z \) are related to the evolution equations.
2. **Einstein - Christoffel Representation**:
- In the Einstein - Christoffel representation, the evolution equations are strongly hyperbolic, so the entire system is well - conditioned. The author also proves that the projection of the Einstein tensor on the boundary \( G_{ab}e^a \) can be regarded as non - trivial boundary conditions.
- Different from the ADM representation, in the Einstein - Christoffel representation, there are three non - trivial boundary conditions, corresponding to three non - preserved constraint conditions respectively.
### Conclusion
Through the above analysis, the author proves that the projection of the Einstein tensor on the boundary can be used as an effective boundary condition, which helps to reduce the phenomenon of constraint violation in numerical simulations. This finding is of great significance for improving the stability and accuracy of numerical relativity simulations.