uu.seUppsala University Publications

CiteExport$(function(){PrimeFaces.cw("TieredMenu","widget_formSmash_upper_j_idt158",{id:"formSmash:upper:j_idt158",widgetVar:"widget_formSmash_upper_j_idt158",autoDisplay:true,overlay:true,my:"left top",at:"left bottom",trigger:"formSmash:upper:exportLink",triggerEvent:"click"});}); $(function(){PrimeFaces.cw("OverlayPanel","widget_formSmash_upper_j_idt159_j_idt161",{id:"formSmash:upper:j_idt159:j_idt161",widgetVar:"widget_formSmash_upper_j_idt159_j_idt161",target:"formSmash:upper:j_idt159:permLink",showEffect:"blind",hideEffect:"fade",my:"right top",at:"right bottom",showCloseIcon:true});});

Perfectly Matched Layers and High Order Difference Methods for Wave EquationsPrimeFaces.cw("AccordionPanel","widget_formSmash_some",{id:"formSmash:some",widgetVar:"widget_formSmash_some",multiple:true}); PrimeFaces.cw("AccordionPanel","widget_formSmash_all",{id:"formSmash:all",widgetVar:"widget_formSmash_all",multiple:true});
function selectAll()
{
var panelSome = $(PrimeFaces.escapeClientId("formSmash:some"));
var panelAll = $(PrimeFaces.escapeClientId("formSmash:all"));
panelAll.toggle();
toggleList(panelSome.get(0).childNodes, panelAll);
toggleList(panelAll.get(0).childNodes, panelAll);
}
/*Toggling the list of authorPanel nodes according to the toggling of the closeable second panel */
function toggleList(childList, panel)
{
var panelWasOpen = (panel.get(0).style.display == 'none');
// console.log('panel was open ' + panelWasOpen);
for (var c = 0; c < childList.length; c++) {
if (childList[c].classList.contains('authorPanel')) {
clickNode(panelWasOpen, childList[c]);
}
}
}
/*nodes have styleClass ui-corner-top if they are expanded and ui-corner-all if they are collapsed */
function clickNode(collapse, child)
{
if (collapse && child.classList.contains('ui-corner-top')) {
// console.log('collapse');
child.click();
}
if (!collapse && child.classList.contains('ui-corner-all')) {
// console.log('expand');
child.click();
}
}
PrimeFaces.cw("AccordionPanel","widget_formSmash_responsibleOrgs",{id:"formSmash:responsibleOrgs",widgetVar:"widget_formSmash_responsibleOrgs",multiple:true}); 2012 (English)Doctoral thesis, comprehensive summary (Other academic)
##### Abstract [en]

##### Place, publisher, year, edition, pages

Uppsala: Acta Universitatis Upsaliensis, 2012. , 47 p.
##### Series

Digital Comprehensive Summaries of Uppsala Dissertations from the Faculty of Science and Technology, ISSN 1651-6214 ; 931
##### Keyword [en]

Elastic waves, Surface waves, Perfectly matched layers, High order difference methods, Stability, Summation-by-parts operators, Boundary treatments
##### National Category

Computational Mathematics
##### Research subject

Scientific Computing with specialization in Numerical Analysis
##### Identifiers

URN: urn:nbn:se:uu:diva-173009ISBN: 978-91-554-8365-4 OAI: oai:DiVA.org:uu-173009DiVA: diva2:516613
##### Public defence

2012-06-08, Room 2446, Polacksbacken, Lägerhyddsvägen 2D, Uppsala, 10:00 (English)
##### Opponent

PrimeFaces.cw("AccordionPanel","widget_formSmash_j_idt467",{id:"formSmash:j_idt467",widgetVar:"widget_formSmash_j_idt467",multiple:true});
##### Supervisors

PrimeFaces.cw("AccordionPanel","widget_formSmash_j_idt473",{id:"formSmash:j_idt473",widgetVar:"widget_formSmash_j_idt473",multiple:true});
#####

PrimeFaces.cw("AccordionPanel","widget_formSmash_j_idt481",{id:"formSmash:j_idt481",widgetVar:"widget_formSmash_j_idt481",multiple:true});
##### Funder

Swedish Research Council, VR 2009-5852
Available from: 2012-05-14 Created: 2012-04-17 Last updated: 2012-10-05Bibliographically approved
##### List of papers

The perfectly matched layer (PML) is a novel technique to simulate the absorption of waves in unbounded domains. The underlying equations are often a system of second order hyperbolic partial differential equations. In the numerical treatment, second order systems are often rewritten and solved as first order systems. There are several benefits with solving the equations in second order formulation, though. However, while the theory and numerical methods for first order hyperbolic systems are well developed, numerical techniques to solve second order hyperbolic systems are less complete.

We construct a strongly well-posed PML for second order systems in two space dimensions, focusing on the equations of linear elasto-dynamics. In the continuous setting, the stability of both first order and second order formulations are linearly equivalent. 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 at most resolutions. In the second order setting growth occurs only if growing modes are well resolved. We determine the number of grid points that can be used in the PML to ensure a discretely stable PML, for several anisotropic elastic materials.

We study the stability of the PML for problems where physical boundaries are important. First, we consider the PML in a waveguide governed by the scalar wave equation. To ensure the accuracy and the stability of the discrete PML, we derived a set of equivalent boundary conditions. Second, we consider the PML for second order symmetric hyperbolic systems on a half-plane. For a class of stable boundary conditions, we derive transformed boundary conditions and prove the stability of the corresponding half-plane problem. Third, we extend the stability analysis to rectangular elastic waveguides, and demonstrate the stability of the discrete PML.

Building on high order summation-by-parts operators, we derive high order accurate and strictly stable finite difference approximations for second order time-dependent hyperbolic systems on bounded domains. Natural and mixed boundary conditions are imposed weakly using the simultaneous approximation term method. Dirichlet boundary conditions are imposed strongly by injection. By constructing continuous strict energy estimates and analogous discrete strict energy estimates, we prove strict stability.

1. A well-posed and discretely stable perfectly matched layer for elastic wave equations in second order formulation$(function(){PrimeFaces.cw("OverlayPanel","overlay480335",{id:"formSmash:j_idt518:0:j_idt522",widgetVar:"overlay480335",target:"formSmash:j_idt518:0:partsLink",showEvent:"mousedown",hideEvent:"mousedown",showEffect:"blind",hideEffect:"fade",appendToBody:true});});

2. Discrete stability of perfectly matched layers for anisotropic wave equations in first and second order formulation$(function(){PrimeFaces.cw("OverlayPanel","overlay516176",{id:"formSmash:j_idt518:1:j_idt522",widgetVar:"overlay516176",target:"formSmash:j_idt518:1:partsLink",showEvent:"mousedown",hideEvent:"mousedown",showEffect:"blind",hideEffect:"fade",appendToBody:true});});

3. On the accuracy and stability of the perfectly matched layer in transient waveguides$(function(){PrimeFaces.cw("OverlayPanel","overlay516160",{id:"formSmash:j_idt518:2:j_idt522",widgetVar:"overlay516160",target:"formSmash:j_idt518:2:partsLink",showEvent:"mousedown",hideEvent:"mousedown",showEffect:"blind",hideEffect:"fade",appendToBody:true});});

4. Boundary waves and stability of the perfectly matched layer$(function(){PrimeFaces.cw("OverlayPanel","overlay515967",{id:"formSmash:j_idt518:3:j_idt522",widgetVar:"overlay515967",target:"formSmash:j_idt518:3:partsLink",showEvent:"mousedown",hideEvent:"mousedown",showEffect:"blind",hideEffect:"fade",appendToBody:true});});

5. Numerical interaction of boundary waves with perfectly matched layers in elastic waveguides$(function(){PrimeFaces.cw("OverlayPanel","overlay515971",{id:"formSmash:j_idt518:4:j_idt522",widgetVar:"overlay515971",target:"formSmash:j_idt518:4:partsLink",showEvent:"mousedown",hideEvent:"mousedown",showEffect:"blind",hideEffect:"fade",appendToBody:true});});

6. Stable and high-order accurate boundary treatments for the elastic wave equation on second-order form$(function(){PrimeFaces.cw("OverlayPanel","overlay516183",{id:"formSmash:j_idt518:5:j_idt522",widgetVar:"overlay516183",target:"formSmash:j_idt518:5:partsLink",showEvent:"mousedown",hideEvent:"mousedown",showEffect:"blind",hideEffect:"fade",appendToBody:true});});

CiteExport$(function(){PrimeFaces.cw("TieredMenu","widget_formSmash_lower_j_idt1290",{id:"formSmash:lower:j_idt1290",widgetVar:"widget_formSmash_lower_j_idt1290",autoDisplay:true,overlay:true,my:"left top",at:"left bottom",trigger:"formSmash:lower:exportLink",triggerEvent:"click"});}); $(function(){PrimeFaces.cw("OverlayPanel","widget_formSmash_lower_j_idt1293_j_idt1296",{id:"formSmash:lower:j_idt1293:j_idt1296",widgetVar:"widget_formSmash_lower_j_idt1293_j_idt1296",target:"formSmash:lower:j_idt1293:permLink",showEffect:"blind",hideEffect:"fade",my:"right top",at:"right bottom",showCloseIcon:true});});