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
A New Absorbing Layer for Elastic Waves
KTH, School of Computer Science and Communication (CSC), Numerical Analysis and Computer Science, NADA.
KTH, School of Computer Science and Communication (CSC), Numerical Analysis and Computer Science, NADA.
2006 (English)In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 215, no 2, 642-660 p.Article in journal (Refereed) Published
Abstract [en]

A new perfectly matched layer (PML) for the simulation of elastic waves in anisotropic media on an unbounded domain is constructed. Theoretical and numerical results, showing that the stability properties of the present layer are better than previously suggested layers, are presented. In addition, the layer can be formulated with fewer auxiliary variables than the split-field PML.

Place, publisher, year, edition, pages
2006. Vol. 215, no 2, 642-660 p.
Keyword [en]
Absorbing layers, Elastic waves, Perfectly matched layers, PML
National Category
Computational Mathematics
Identifiers
URN: urn:nbn:se:kth:diva-7435DOI: 10.1016/j.jcp.2005.11.006ISI: 000237458200016Scopus ID: 2-s2.0-33645909260OAI: oai:DiVA.org:kth-7435DiVA: diva2:12461
Note
QC 20100830. Uppdaterad från Submitted till Published 20100830.Available from: 2005-10-14 Created: 2005-10-14 Last updated: 2017-12-14Bibliographically approved
In thesis
1. Absorbing Layers and Non-Reflecting Boundary Conditions for Wave Propagation Problems
Open this publication in new window or tab >>Absorbing Layers and Non-Reflecting Boundary Conditions for Wave Propagation Problems
2005 (English)Doctoral thesis, comprehensive summary (Other scientific)
Abstract [en]

The presence of wave motion is the defining feature in many fields of application,such as electro-magnetics, seismics, acoustics, aerodynamics,oceanography and optics. In these fields, accurate numerical simulation of wave phenomena is important for the enhanced understanding of basic phenomenon, but also in design and development of various engineering applications.

In general, numerical simulations must be confined to truncated domains, much smaller than the physical space were the wave phenomena takes place. To truncate the physical space, artificial boundaries, and corresponding boundary conditions, are introduced. There are four main classes of methods that can be used to truncate problems on unbounded or large domains: boundary integral methods, infinite element methods, non-reflecting boundary condition methods and absorbing layer methods.

In this thesis, we consider different aspects of non-reflecting boundary conditions and absorbing layers. In paper I, we construct discretely non-reflecting boundary conditions for a high order centered finite difference scheme. This is done by separating the numerical solution into spurious and physical waves, using the discrete dispersion relation.

In paper II-IV, we focus on the perfectly matched layer method, which is a particular absorbing layer method. An open issue is whether stable perfectly matched layers can be constructed for a general hyperbolic system.

In paper II, we present a stable perfectly matched layer formulation for 2 x 2 symmetric hyperbolic systems in (2 + 1) dimensions. We also show how to choose the layer parameters as functions of the coefficient matrices to guarantee stability.

In paper III, we construct a new perfectly matched layer for the simulation of elastic waves in an anisotropic media. We present theoretical and numerical results, showing that the stability properties of the present layer are better than previously suggested layers.

In paper IV, we develop general tools for constructing PMLs for first order hyperbolic systems. We present a model with many parameters which is applicable to all hyperbolic systems, and which we prove is well-posed and perfectly matched. We also use an automatic method, derived in paper V, for analyzing the stability of the model and establishing energy inequalities. We illustrate our techniques with applications to Maxwell s equations, the linearized Euler equations, as well as arbitrary 2 x 2 systems in (2 + 1) dimensions.

In paper V, we use the method of Sturm sequences for bounding the real parts of roots of polynomials, to construct an automatic method for checking Petrowsky well-posedness of a general Cauchy problem. We prove that this method can be adapted to automatically symmetrize any well-posed problem, producing an energy estimate involving only local quantities.

Place, publisher, year, edition, pages
Stockholm: KTH, 2005. iv, 19 p.
Series
Trita-NA, ISSN 0348-2952 ; 2005:34
Keyword
PML, perfectly matched layer
National Category
Computational Mathematics
Identifiers
urn:nbn:se:kth:diva-448 (URN)91-7178-174-9 (ISBN)
Public defence
2005-10-28, Salongen, KTHB, Osquars backe 31, KTH, Stockholm, 10:15
Opponent
Supervisors
Note
QC 20100830Available from: 2005-10-14 Created: 2005-10-14 Last updated: 2010-08-30Bibliographically approved

Open Access in DiVA

No full text

Other links

Publisher's full textScopusScienceDirect

Search in DiVA

By author/editor
Appelö, DanielKreiss, Gunilla
By organisation
Numerical Analysis and Computer Science, NADA
In the same journal
Journal of Computational Physics
Computational Mathematics

Search outside of DiVA

GoogleGoogle Scholar

doi
urn-nbn

Altmetric score

doi
urn-nbn
Total: 101 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