uu.seUppsala University Publications
Change search
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • ieee
  • modern-language-association
  • 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
Well-posed and discretely stable perfectly matched layers for elastic wave equations in second order formulation
Uppsala University, Disciplinary Domain of Science and Technology, Mathematics and Computer Science, Department of Information Technology, Division of Scientific Computing. Uppsala University, Disciplinary Domain of Science and Technology, Mathematics and Computer Science, Department of Information Technology, Numerical Analysis.
Uppsala University, Disciplinary Domain of Science and Technology, Mathematics and Computer Science, Department of Information Technology, Division of Scientific Computing. Uppsala University, Disciplinary Domain of Science and Technology, Mathematics and Computer Science, Department of Information Technology, Numerical Analysis.
2010 (English)Report (Other academic)
Place, publisher, year, edition, pages
2010.
Series
Technical report / Department of Information Technology, Uppsala University, ISSN 1404-3203 ; 2010-004
National Category
Computational Mathematics
Identifiers
URN: urn:nbn:se:uu:diva-117577OAI: oai:DiVA.org:uu-117577DiVA: diva2:298003
Available from: 2010-02-19 Created: 2010-02-19 Last updated: 2014-06-16Bibliographically approved
In thesis
1. Perfectly matched layers for second order wave equations
Open this publication in new window or tab >>Perfectly matched layers for second order wave equations
2010 (English)Licentiate thesis, comprehensive summary (Other academic)
Abstract [en]

Numerical simulation of propagating waves in unbounded spatial domains is a challenge common to many branches of engineering and applied mathematics. Perfectly matched layers (PML) are a novel technique for simulating the absorption of waves in open domains. The equations modeling the dynamics of phenomena of interest are usually posed as differential equations (or integral equations) which must be solved at every time instant. In many application areas like general relativity, seismology and acoustics, the underlying equations are systems of second order hyperbolic partial differential equations. In numerical treatment of such problems, the equations are often rewritten as first order systems and are solved in this form. For this reason, many existing PML models have been developed for first order systems. In several studies, it has been reported that there are drawbacks with rewriting second order systems into first order systems before numerical solutions are obtained. While the theory and numerical methods for first order systems are well developed, numerical techniques to solve second order hyperbolic systems is an on-going research.

In the first part of this thesis, we construct PML equations for systems of second order hyperbolic partial differential equations in two space dimensions, focusing on the equations of linear elasto-dynamics. One advantage of this approach is that we can choose auxiliary variables such that the PML is strongly hyperbolic, thus strongly well-posed. The second is that it requires less auxiliary variables as compared to existing first order formulations. However, in continuum the stability of both first order and second order formulations are linearly equivalent. A turning point is in numerical approximations. We have found that if the so-called geometric stability condition is violated, approximating the first order PML with standard central differences leads to a high frequency instability for any given resolution. The second order discretization behaves much more stably. In the second order setting instability occurs only if unstable modes are well resolved.

The second part of this thesis discusses the construction of PML equations for the time-dependent Schrödinger equation. From mathematical perspective, the Schrödinger equation is unique, in the sense that it is only first order in time but second order in space. However, with slight modifications, we carry over our ideas from the hyperbolic systems to the Schrödinger equations and derive a set of asymptotically stable PML equations. The new model can be viewed as a modified complex absorbing potential (CAP). The PML model can easily be adapted to existing codes developed for CAP by accurately discretizing the auxiliary variables and appending them accordingly. Numerical experiments are presented illustrating the accuracy and absorption properties of the new PML model.

We are hopeful that the results obtained in this thesis will find useful applications in time-dependent wave scattering calculations.

Place, publisher, year, edition, pages
Uppsala University, 2010
Series
Information technology licentiate theses: Licentiate theses from the Department of Information Technology, ISSN 1404-5117 ; 2010-004
National Category
Computational Mathematics
Research subject
Scientific Computing with specialization in Numerical Analysis
Identifiers
urn:nbn:se:uu:diva-124538 (URN)
Supervisors
Available from: 2010-05-07 Created: 2010-05-04 Last updated: 2017-08-31Bibliographically approved

Open Access in DiVA

No full text

Other links

http://www.it.uu.se/research/publications/reports/2010-004/

Authority records BETA

Duru, KennethKreiss, Gunilla

Search in DiVA

By author/editor
Duru, KennethKreiss, Gunilla
By organisation
Division of Scientific ComputingNumerical Analysis
Computational Mathematics

Search outside of DiVA

GoogleGoogle Scholar

urn-nbn

Altmetric score

urn-nbn
Total: 525 hits
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • ieee
  • modern-language-association
  • 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