Generating Initial Data in General Relativity using Adaptive Finite Element Methods
arXiv:0801.3142
Abstract
The conformal formulation of the Einstein constraint equations is first reviewed, and we then consider the design, analysis, and implementation of adaptive multilevel finite element-type numerical methods for the resulting coupled nonlinear elliptic system. We derive weak formulations of the coupled constraints, and review some new developments in the solution theory for the constraints in the cases of constant mean extrinsic curvature (CMC) data, near-CMC data, and arbitrarily prescribed mean extrinsic curvature data. We then outline some recent results on a priori and a posteriori error estimates for a broad class of Galerkin-type approximation methods for this system which includes techniques such as finite element, wavelet, and spectral methods. We then use these estimates to construct an adaptive finite element method (AFEM) for solving this system numerically, and outline some new convergence and optimality results. We then describe in some detail an implementation of the methods using the FETK software package, which is an adaptive multilevel finite element code designed to solve nonlinear elliptic and parabolic systems on Riemannian manifolds. We finish by describing a simplex mesh generation algorithm for compact binary objects, and then look at a detailed example showing the use of FETK for numerical solution of the constraints.
54 pages, 15 figures. Published as LSU CCT Technical Report 2008-09
References in corpus (12)
- The Dynamics of General Relativity
- Numerical Relativity: A review
- Excision boundary conditions for black hole initial data
- Einstein Constraints on Asymptotically Euclidean Manifolds
- Rough Solutions of the Einstein Constraint Equations on Compact Manifolds
- Strongly hyperbolic second order Einstein's evolution equations
- Gravitational Collapse of Gravitational Waves in 3D Numerical Relativity
- Far-from-constant mean curvature solutions of Einstein's constraint equations with positive Yamabe metrics
- Near-Constant Mean Curvature Solutions of the Einstein Constraint Equations with Non-Negative Yamabe Metrics
- Dynamical Gauge Conditions for the Einstein Evolution Equations
- Solving the Einstein constraint equations on multi-block triangulations using finite element methods
- Rough solutions of the Einstein constraint equations with nonconstant mean curvature