kth.sePublikationer
Ändra sökning
RefereraExporteraLänk till posten
Permanent länk

Direktlänk
Referera
Referensformat
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Annat format
Fler format
Språk
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Annat språk
Fler språk
Utmatningsformat
  • html
  • text
  • asciidoc
  • rtf
Absorbing Layers and Non-Reflecting Boundary Conditions for Wave Propagation Problems
KTH, Skolan för datavetenskap och kommunikation (CSC), Numerisk Analys och Datalogi, NADA.
2005 (Engelska)Doktorsavhandling, sammanläggning (Övrigt vetenskapligt)
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.

Ort, förlag, år, upplaga, sidor
Stockholm: KTH , 2005. , s. iv, 19
Serie
Trita-NA, ISSN 0348-2952 ; 2005:34
Nyckelord [en]
PML, perfectly matched layer
Nationell ämneskategori
Beräkningsmatematik
Identifikatorer
URN: urn:nbn:se:kth:diva-448ISBN: 91-7178-174-9 (tryckt)OAI: oai:DiVA.org:kth-448DiVA, id: diva2:12464
Disputation
2005-10-28, Salongen, KTHB, Osquars backe 31, KTH, Stockholm, 10:15
Opponent
Handledare
Anmärkning
QC 20100830Tillgänglig från: 2005-10-14 Skapad: 2005-10-14 Senast uppdaterad: 2022-06-26Bibliografiskt granskad
Delarbeten
1. Discretely nonreflecting boundary conditions for higher order centered schemes for wave equations
Öppna denna publikation i ny flik eller fönster >>Discretely nonreflecting boundary conditions for higher order centered schemes for wave equations
2003 (Engelska)Ingår i: Proceedings of the WAVES-2003 conference, Berlin: Springer Verlag , 2003, s. 130-135Kapitel i bok, del av antologi (Övrigt vetenskapligt)
Abstract [en]

Using the framework introduced by Rawley and Colonius [2] we construct a nonreflecting boundary condition for the one-way wave equation spatially discretized with a fourth order centered difference scheme. The boundary condition, which can be extended to arbitrary order accuracy, is shown to be well posed. Numerical simulations have been performed showing promising results.

Ort, förlag, år, upplaga, sidor
Berlin: Springer Verlag, 2003
Nationell ämneskategori
Beräkningsmatematik
Identifikatorer
urn:nbn:se:kth:diva-7433 (URN)3-540-40127-X (ISBN)
Anmärkning
QC 20100830. Konferens: : 6th International Conference on Mathematical and Numerical Aspects of Wave Propagation (WAVES 2003) JYVASKYLA, FINLAND, JUN 30-JUL 04, 2003.Tillgänglig från: 2005-10-14 Skapad: 2005-10-14 Senast uppdaterad: 2022-06-26Bibliografiskt granskad
2. Construction of stable PMLs for general 2 x 2 symmetric hyperbolic systems
Öppna denna publikation i ny flik eller fönster >>Construction of stable PMLs for general 2 x 2 symmetric hyperbolic systems
2004 (Engelska)Ingår i: Proceedings of the HYP2004 conference, 2004, s. 1-8Konferensbidrag, Publicerat paper (Refereegranskat)
Abstract [en]

The perfectly matched layer (PML) has emerged as animportant tool for accurately solving certain hyperbolic systems onunbounded domains. An open issue is whether stable PMLs can beconstructed in general. In this work we consider the specializationof our general PML formulation to 2 × 2 symmetric hyperbolicsystems in 2 + 1 dimensions. We show how to choose the layerparameters as functions of the coefficient matrices to guaranteestability.

Nationell ämneskategori
Beräkningsmatematik
Identifikatorer
urn:nbn:se:kth:diva-7434 (URN)
Konferens
The HYP2004 conference, September 13-17 2004, Osaka, Japan
Anmärkning
QC 20100830Tillgänglig från: 2005-10-14 Skapad: 2005-10-14 Senast uppdaterad: 2022-06-26Bibliografiskt granskad
3. A New Absorbing Layer for Elastic Waves
Öppna denna publikation i ny flik eller fönster >>A New Absorbing Layer for Elastic Waves
2006 (Engelska)Ingår i: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716, Vol. 215, nr 2, s. 642-660Artikel i tidskrift (Refereegranskat) 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.

Nyckelord
Absorbing layers, Elastic waves, Perfectly matched layers, PML
Nationell ämneskategori
Beräkningsmatematik
Identifikatorer
urn:nbn:se:kth:diva-7435 (URN)10.1016/j.jcp.2005.11.006 (DOI)000237458200016 ()2-s2.0-33645909260 (Scopus ID)
Anmärkning
QC 20100830. Uppdaterad från Submitted till Published 20100830.Tillgänglig från: 2005-10-14 Skapad: 2005-10-14 Senast uppdaterad: 2022-06-26Bibliografiskt granskad
4. Perfectly matched layers for hyperbolic systems: General formulation, well-posedness and stability
Öppna denna publikation i ny flik eller fönster >>Perfectly matched layers for hyperbolic systems: General formulation, well-posedness and stability
2006 (Engelska)Ingår i: SIAM Journal on Applied Mathematics, ISSN 0036-1399, E-ISSN 1095-712X, Vol. 67, nr 1, s. 1-23Artikel i tidskrift (Refereegranskat) Published
Abstract [en]

Since its introduction the perfectly matched layer (PML) has proven to be an accurate and robust method for domain truncation in computational electromagnetics. However, the mathematical analysis of PMLs has been limited to special cases. In particular, the basic question of whether or not a stable PML exists for arbitrary wave propagation problems remains unanswered. In this work 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 introduce an automatic method 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, and arbitrary 2 x 2 systems in (2 + 1) dimensions.

Nyckelord
Perfectly matched layers, Stability, Computational methods, Euler equations, Mathematical models, Maxwell equations, Parameter estimation, Wave propagation, Computational electromagnetics, Hyperbolic systems, Mathematical analysis, Perfectly matched layers, Magnetoelectric effects
Nationell ämneskategori
Matematik
Identifikatorer
urn:nbn:se:kth:diva-7436 (URN)10.1137/050639107 (DOI)000243279300001 ()2-s2.0-33847744194 (Scopus ID)
Anmärkning
QC 20100830. Uppdaterad från Submitted till Published 20100830.Tillgänglig från: 2005-10-14 Skapad: 2005-10-14 Senast uppdaterad: 2022-06-26Bibliografiskt granskad
5. Automatic symmetrization and energy estimates using local operators for partial differential equations
Öppna denna publikation i ny flik eller fönster >>Automatic symmetrization and energy estimates using local operators for partial differential equations
2007 (Engelska)Ingår i: Communications in Partial Differential Equations, ISSN 0360-5302, E-ISSN 1532-4133, Vol. 32, nr 7, s. 1129-1145Artikel i tidskrift (Refereegranskat) Published
Abstract [en]

We develop a method for automatically symmetrizing Petrowsky well-posed Cauchy problems for constant coefficient linear partial differential equations. The method is rooted in the Sturm sequence technique for establishing the location of the roots of a complex polynomial and can be automated using standard symbolic computation tools. In the special case of homogeneous strictly hyperbolic scalar equations, we show that the resulting estimates are strong enough to control all principal order derivatives and thus can be used in place of the Leray energies. We also illustrate the method by applying it to various problems of mixed type.

Nyckelord
Cauchy problem, Energy estimates, Sturm sequences, Well-posedness
Nationell ämneskategori
Matematik
Identifikatorer
urn:nbn:se:kth:diva-7437 (URN)10.1080/03605300600854258 (DOI)000250012900005 ()2-s2.0-34547762609 (Scopus ID)
Anmärkning
QC 20100830. Tidigare titel: On symmetrization and energy estimates using local operators for partial differential equations. Titel ändrad samt uppdaterad från Manuskript till Artikel i tidskrift 20100830Tillgänglig från: 2005-10-14 Skapad: 2005-10-14 Senast uppdaterad: 2022-06-26Bibliografiskt granskad

Open Access i DiVA

fulltext(1453 kB)3167 nedladdningar
Filinformation
Filnamn FULLTEXT01.pdfFilstorlek 1453 kBChecksumma SHA-1
c09fb15b160fb8e437d315e235adeb7b7829ed79aff7e4150c9cd24c8134419640a4f7cb
Typ fulltextMimetyp application/pdf

Sök vidare i DiVA

Av författaren/redaktören
Appelö, Daniel
Av organisationen
Numerisk Analys och Datalogi, NADA
Beräkningsmatematik

Sök vidare utanför DiVA

GoogleGoogle Scholar
Totalt: 3167 nedladdningar
Antalet nedladdningar är summan av nedladdningar för alla fulltexter. Det kan inkludera t.ex tidigare versioner som nu inte längre är tillgängliga.

isbn
urn-nbn

Altmetricpoäng

isbn
urn-nbn
Totalt: 1897 träffar
RefereraExporteraLänk till posten
Permanent länk

Direktlänk
Referera
Referensformat
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Annat format
Fler format
Språk
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Annat språk
Fler språk
Utmatningsformat
  • html
  • text
  • asciidoc
  • rtf