Quick search Find article
Quick search
Find article

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

Oleg Korobkin1,2, Burak Aksoylu2,3, Michael Holst4, Enrique Pazos1,2 and Manuel Tiglio5,6

Show affiliations


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 ψ. 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.


PACS

04.20.Ex Initial value problem, existence and uniqueness of solutions

04.30.-w Gravitational waves

MSC

83C05 Einstein's equations (general structure, canonical formalism, Cauchy problems)

83C35 Gravitational waves

Subjects

Gravitation and cosmology

Dates

Issue 14 (21 July 2009)

Received 5 November 2008, in final form 8 April 2009

Published 25 June 2009



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

    Oleg Korobkin et al 2009 Class. Quantum Grav. 26 145007

  2. Moment inversion problem for piecewise D-finite functions

    Dmitry Batenkov 2009 Inverse Problems 25 105001

  3. Connection between structural characteristics of glasses and their crystalline counterparts

    S N Taraskin 2007 J. Phys.: Condens. Matter 19 455215

  4. Variational calculus with constraints on general algebroids

    Katarzyna Grabowska and Janusz Grabowski 2008 J. Phys. A: Math. Theor. 41 175204

  5. Interaction and coalescence of nanovoids and dynamic fracture in silica glass: multimillion-to-billion atom molecular dynamics simulations

    Ken-ichi Nomura et al 2009 J. Phys. D: Appl. Phys. 42 214011

  6. Operational status of TAMA300 with the seismic attenuation system (SAS)

    R Takahashi et al 2008 Class. Quantum Grav. 25 114036

  7. Baxterization of GLq(2) and its application to the Liouville model and some other models on a lattice

    Andrei G Bytsko 2008 J. Phys. A: Math. Theor. 41 194003

  8. Non-equilibrium relaxation and near-arrest dynamics in colloidal suspensions

    M Medina-Noyola and Pedro Ramírez-González 2009 J. Phys.: Condens. Matter 21 504103

  9. Fully connected network of superconducting qubits in a cavity

    Dimitris I Tsomokos et al 2008 New J. Phys. 10 113020

  10. Low-temperature growth of single-wall carbon nanotubes

    S-M Kim et al 2007 Nanotechnology 18 495606

View by subject




Export








Please login to access our web services, or create an account if you don't yet have one.

You must have cookies enabled in your web browser to be able to login.

Username
Password

Forgotten your password? Get a new one here.