uu.seUppsala University Publications

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

Perfectly matched layers for second order 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();
}
}
2010 (English)Licentiate thesis, comprehensive summary (Other academic)
##### Abstract [en]

##### 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: urn:nbn:se:uu:diva-124538OAI: oai:DiVA.org:uu-124538DiVA, id: diva2:317660
#####

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

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

PrimeFaces.cw("AccordionPanel","widget_formSmash_j_idt1130",{id:"formSmash:j_idt1130",widgetVar:"widget_formSmash_j_idt1130",multiple:true}); Available from: 2010-05-07 Created: 2010-05-04 Last updated: 2017-08-31Bibliographically approved
##### List of papers

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.

1. Well-posed and discretely stable perfectly matched layers for elastic wave equations in second order formulation$(function(){PrimeFaces.cw("OverlayPanel","overlay298003",{id:"formSmash:j_idt1179:0:j_idt1183",widgetVar:"overlay298003",target:"formSmash:j_idt1179: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_idt1179:1:j_idt1183",widgetVar:"overlay516176",target:"formSmash:j_idt1179:1:partsLink",showEvent:"mousedown",hideEvent:"mousedown",showEffect:"blind",hideEffect:"fade",appendToBody:true});});

3. Stable perfectly matched layers for the Schrödinger equations$(function(){PrimeFaces.cw("OverlayPanel","overlay359730",{id:"formSmash:j_idt1179:2:j_idt1183",widgetVar:"overlay359730",target:"formSmash:j_idt1179:2:partsLink",showEvent:"mousedown",hideEvent:"mousedown",showEffect:"blind",hideEffect:"fade",appendToBody:true});});

urn-nbn$(function(){PrimeFaces.cw("Tooltip","widget_formSmash_j_idt1840",{id:"formSmash:j_idt1840",widgetVar:"widget_formSmash_j_idt1840",showEffect:"fade",hideEffect:"fade",showDelay:500,hideDelay:300,target:"formSmash:altmetricDiv"});});

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