Solving the Einstein constraint equations on multi-block triangulations using finite element methods

DOI10.1088/0264-9381/26/14/145007MaRDI QIDQ3393265zbMATH OpenOpenAlexWikidataFDO

Authors Oleg Korobkin, Burak Aksoylu, Enrique Pazos, Manuel Tiglio, Michael Holst

Publication date 19 August 2009

Published in Classical and Quantum Gravity (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/0801.1823

Mathematics Subject Classification ID

Software, source code, etc. for problems pertaining to relativity and gravitational theory (83-04) Einstein's equations (general structure, canonical formalism, Cauchy problems) (83C05) Exact solutions to problems in general relativity and gravitational theory (83C15) Gravitational waves (83C35)

Abstract: In order to generate initial data for nonlinear relativistic simulations, one needs to solve the Einstein constraints, which can be cast into a coupled set of nonlinear elliptic equations. Here we present an approach for solving these equations on three-dimensional multi-block domains using finite element methods. We illustrate our approach on a simple example of Brill wave initial data, with the constraints reducing to a single linear elliptic equation for the conformal factor $psi$. We use quadratic Lagrange elements on semi-structured simplicial meshes, obtained by triangulation of multi-block grids. In the case of uniform refinement the scheme is superconvergent at most mesh vertices, due to local symmetry of the finite element basis with respect to local spatial inversions. We show that in the superconvergent case subsequent unstructured mesh refinements do not improve the quality of our initial data. As proof of concept that this approach is feasible for generating multi-block initial data in three dimensions, after constructing the initial data we evolve them in time using a high order finite-differencing multi-block approach and extract the gravitational waves from the numerical solution.

Recommendations

• An application of the multi-grid method to the construction of initial data for Brill waves
• THE INITIAL VALUE PROBLEM IN NUMERICAL RELATIVITY
• Solving the Einstein constraints in periodic spaces with a multigrid approach
• Toward computing gravitational initial data without elliptic solvers
• scientific article; zbMATH DE number 1183037

Cited in

(4)

• An operator-based local discontinuous Galerkin method compatible with the BSSN formulation of the Einstein equations
• Solving the Einstein constraints in periodic spaces with a multigrid approach
• Characteristic evolution and matching
• An application of the multi-grid method to the construction of initial data for Brill waves

Describes a project that uses

Uses Software

• Cactus
• MPI/MPICH
• LAM-MPI

This page was built for publication: Solving the Einstein constraint equations on multi-block triangulations using finite element methods

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q3393265)