uu.seUppsala University Publications

References$(function(){PrimeFaces.cw("TieredMenu","widget_formSmash_upper_j_idt145",{id:"formSmash:upper:j_idt145",widgetVar:"widget_formSmash_upper_j_idt145",autoDisplay:true,overlay:true,my:"left top",at:"left bottom",trigger:"formSmash:upper:referencesLink",triggerEvent:"click"});}); $(function(){PrimeFaces.cw("OverlayPanel","widget_formSmash_upper_j_idt146_j_idt148",{id:"formSmash:upper:j_idt146:j_idt148",widgetVar:"widget_formSmash_upper_j_idt146_j_idt148",target:"formSmash:upper:j_idt146:permLink",showEffect:"blind",hideEffect:"fade",my:"right top",at:"right bottom",showCloseIcon:true});});

High Order Symmetric Finite Difference Schemes for the Acoustic Wave EquationPrimeFaces.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}); 2003 (English)Doctoral thesis, monograph (Other academic)
##### Abstract [en]

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

Uppsala: Acta Universitatis Upsaliensis , 2003. , 130 p.
##### Series

Uppsala Dissertations from the Faculty of Science and Technology, ISSN 1104-2516 ; 46
##### Keyword [en]

finite element method, mass lumping, high order, finite difference method, wave equation
##### National Category

Computational Mathematics
##### Research subject

Numerical Analysis
##### Identifiers

URN: urn:nbn:se:uu:diva-3511ISBN: 91-554-5685-5OAI: oai:DiVA.org:uu-3511DiVA: diva2:163066
##### Public defence

2003-09-26, Room 1311, Polacksbacken, Uppsala University, Uppsala, 10:15 (English)
##### Opponent

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

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

PrimeFaces.cw("AccordionPanel","widget_formSmash_j_idt387",{id:"formSmash:j_idt387",widgetVar:"widget_formSmash_j_idt387",multiple:true});
Available from: 2003-09-05 Created: 2003-09-05 Last updated: 2011-10-27Bibliographically approved

When someone speaks, or when a piano is played, sound waves are spread through the room. Physically the propagation of sound is governed by the acoustic wave equation. This thesis deals with the problem of finding solutions to this equation by numerical methods. These methods utilize the tremendous power of the modern computer to perform arithmetic operations.

The two most common methods for the solution of a partial differential equation like the wave equation are the finite element (FEM) and the finite difference method (FDM). Traditionally the FEM is the method of choice for stationary problems in structural mechanics, while the FDM is often preferred for time dependent problems. In this study we combine the spatial discretization of a specific FEM on a uniform grid with a time discretization from a FDM. The new scheme is explicit, symmetric and has fourth order of accuracy in space and time for Neumann and/or Dirichlet boundary conditions on a rectangle.

The Helmholtz equation is employed for the analysis of the spatial discretization. Concerning the time dependence, we focus on problems for which an explicit time marching scheme is the most adequate choice from a physical point of view.

When the FEM is applied to a time-dependent problem, a mass matrix appears in the front of the time evolution operator. To avoid the costly solution of a sparse linear system for each time step, the mass matrix is approximated by a diagonal matrix (mass lumping). For a particular basis of piecewise cubic polynomials we show that mass lumping with the rowsum preserves the order of accuracy for problems with Neumann boundary conditions. For Dirichlet boundary conditions mass lumping must be done differently in the vicinity of the boundaries to attain full accuracy. The theory of FDM is used for analysis of accuracy and stability.

We study the performance of the new scheme in the presence of media discontinuities. It is found that the continuity of pressure and particle velocity is fulfilled to the first order, even though the discontinuities are not aligned with the cartesian grid.

Numerical examples are presented which illustrate the outstanding performance of the new scheme.

References$(function(){PrimeFaces.cw("TieredMenu","widget_formSmash_lower_j_idt1080",{id:"formSmash:lower:j_idt1080",widgetVar:"widget_formSmash_lower_j_idt1080",autoDisplay:true,overlay:true,my:"left top",at:"left bottom",trigger:"formSmash:lower:referencesLink",triggerEvent:"click"});}); $(function(){PrimeFaces.cw("OverlayPanel","widget_formSmash_lower_j_idt1081_j_idt1083",{id:"formSmash:lower:j_idt1081:j_idt1083",widgetVar:"widget_formSmash_lower_j_idt1081_j_idt1083",target:"formSmash:lower:j_idt1081:permLink",showEffect:"blind",hideEffect:"fade",my:"right top",at:"right bottom",showCloseIcon:true});});