Change search
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • harvard1
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf
Finite Difference Schemes for Second Order Systems Describing Black Holes
KTH, School of Computer Science and Communication (CSC), Numerical Analysis and Computer Science, NADA.
Albert Einstein Institute, Max Planck Gesellschaft.
Albert Einstein Institute, Max Planck Gesellschaft.
KTH, School of Computer Science and Communication (CSC), Numerical Analysis and Computer Science, NADA.
Show others and affiliations
2006 (English)In: Physical Review D. Particles and fields, ISSN 0556-2821, E-ISSN 1089-4918, Vol. 73, no 12, 124008-1-124008-14 p.Article in journal (Refereed) Published
Abstract [en]

In the harmonic description of general relativity, the principal part of Einstein's equations reduces to 10 curved space wave equations for the components of the space-time metric. We present theorems regarding the stability of several evolution-boundary algorithms for such equations when treated in second order differential form. The theorems apply to a model black hole space-time consisting of a spacelike inner boundary excising the singularity, a timelike outer boundary and a horizon in between. These algorithms are implemented as stable, convergent numerical codes and their performance is compared in a 2-dimensional excision problem.

Place, publisher, year, edition, pages
2006. Vol. 73, no 12, 124008-1-124008-14 p.
Keyword [en]
wave-equation, numerical relativity, general-relativity, approximations, stability
National Category
Astronomy, Astrophysics and Cosmology Computational Mathematics
Identifiers
URN: urn:nbn:se:kth:diva-8300DOI: 10.1103/PhysRevD.73.124008ISI: 000238698900048Scopus ID: 2-s2.0-33744907537OAI: oai:DiVA.org:kth-8300DiVA: diva2:13584
Note
QC 20100830Available from: 2008-04-30 Created: 2008-04-30 Last updated: 2010-08-30Bibliographically approved
In thesis
1. Topics in Analysis and Computation of Linear Wave Propagation
Open this publication in new window or tab >>Topics in Analysis and Computation of Linear Wave Propagation
2008 (English)Doctoral thesis, comprehensive summary (Other scientific)
Abstract [en]

This thesis concerns the analysis and numerical simulation of wave propagation problems described by systems of linear hyperbolic partial differential equations.

A major challenge in wave propagation problems is numerical simulation of high frequency waves. When the wavelength is very small compared to the overall size of the computational domain, we encounter a multiscale problem. Examples include the forward and the inverse seismic wave propagation, radiation and scattering problems in computational electromagnetics and underwater acoustics. In direct numerical simulations, the accuracy of the approximate solution is determined by the number of grid points or elements per wavelength. The computational cost to maintain constant accuracy grows algebraically with the frequency, and for sufficiently high frequency, direct numerical simulations are no longer feasible. Other numerical methods are therefore needed. Asymptotic methods, for instance, are good approximations for very high frequency waves. They are based on constructing asymptotic expansions of the solution. The accuracy increases with increasing frequency for a fixed computational cost. Most asymptotic techniques rely on geometrical optics equations with frequency independent unknowns. There are however two deficiencies in the geometrical optics solution. First, it does not include diffraction effects. Secondly, it breaks down at caustics. Geometrical theory of diffraction provides a technique for adding diffraction effects to the geometrical optics approximation by introducing diffracted rays. In papers 1 and 2 we present a numerical algorithm for computing an important type of diffracted rays known as creeping rays. Another asymptotic model which is valid also at caustics is based on Gaussian beams. In papers 3 and 4, we present an error analysis of Gaussian beams approximation and develop a new numerical algorithm for computing Gaussian beams, respectively.

Another challenge in computation of wave propagation problems arises when the system of equations consists of second order hyperbolic equations involving mixed space-time derivatives. Examples include the harmonic formulation of Einstein’s equations and wave equations governing elasticity and acoustics. The classic computational treatment of such second order hyperbolic systems has been based on reducing the systems to first order differential forms. This treatment has however the disadvantage of introducing auxiliary variables with their associated constraints and boundary conditions. In paper 5, we treat the problem in the second order differential form, which has advantages for both computational efficiency and accuracy over the first order formulation.

Finally, paper 6 concerns the concept of well-posedness for a class of linear hyperbolic initial boundary value problems which are not boundary stable. The well-posedness is well established for boundary stable hyperbolic systems for which we can obtain sharp estimates of the solution including estimates at boundaries. There are, however, problems which are not boundary stable but are well-posed in a weaker sense, i.e., the problems for which an energy estimate can be obtained in the interior of the domain but not on the boundaries. We analyze a model problem of this type. Possible applications arise in elastic wave equations and Maxwell’s equations describing glancing and surface waves.

Place, publisher, year, edition, pages
Stockholm: KTH, 2008. xiii, 36 p.
Series
Trita-CSC-A, ISSN 1653-5723 ; 2008:07
National Category
Computational Mathematics
Identifiers
urn:nbn:se:kth:diva-4715 (URN)978-91-7178-961-7 (ISBN)
Public defence
2008-05-20, Sal D2, KTH, Lindstedtsvägen 5, Stockholm, 10:15
Opponent
Supervisors
Note
QC 20100830Available from: 2008-04-30 Created: 2008-04-30 Last updated: 2010-08-30Bibliographically approved

Open Access in DiVA

No full text

Other links

Publisher's full textScopus

Search in DiVA

By author/editor
Motamed, MohammadKreiss, Heinz-Otto
By organisation
Numerical Analysis and Computer Science, NADA
In the same journal
Physical Review D. Particles and fields
Astronomy, Astrophysics and CosmologyComputational Mathematics

Search outside of DiVA

GoogleGoogle Scholar

doi
urn-nbn

Altmetric score

doi
urn-nbn
Total: 112 hits
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • harvard1
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf