Mass Lumping for Constant-Density Acoustics

dc.contributor.authorSymes, William W.en_US
dc.contributor.authorTerentyev, Igoren_US
dc.date.accessioned2018-06-19T17:45:05Zen_US
dc.date.available2018-06-19T17:45:05Zen_US
dc.date.issued2009-04en_US
dc.date.noteApril 2009en_US
dc.description.abstractConforming Galerkin discretization of the constant-density acoustic wave equation provides optimal order convergence even in the presence of very rough coefficients, provided that the time dependence of the data (right-hand side) is minimally smooth. Such discretizations avoid the well-known first-order error ("stairstep diffraction") phenomenon produced by standard finite difference methods. On the other hand, Galerkin methods in themselves are inefficient for high frequency wave simulation, due to the implicit nature of the time step system. Mass lumping renders the time step explicit, and provides an avenue for efficient time-stepping of time-dependent problems with conforming finite element spatial discretization. Typical justifications for mass lumping use quadrature error estimates which do not hold for nonsmooth coefficients. In this paper, we show that the mass-lumped semidiscrete system for the constant-density acoustic wave equation with rectangular multilinear elements exhibits optimal order convergence even when the coefficient (bulk modulus) is merely bounded and measurable, provided that the right-hand side possesses some smoothness in time. We illustrate the theory with numerical examples involving discontinuous, non-grid-aligned bulk moduli, in which the coefficient averaging implicit in mass lumping eliminates the stairstep diffraction effect.en_US
dc.format.extent18 ppen_US
dc.identifier.citationSymes, William W. and Terentyev, Igor. "Mass Lumping for Constant-Density Acoustics." (2009) <a href="https://hdl.handle.net/1911/102111">https://hdl.handle.net/1911/102111</a>.en_US
dc.identifier.digitalTR09-06en_US
dc.identifier.urihttps://hdl.handle.net/1911/102111en_US
dc.language.isoengen_US
dc.titleMass Lumping for Constant-Density Acousticsen_US
dc.typeTechnical reporten_US
dc.type.dcmiTexten_US
Files
Original bundle
Now showing 1 - 1 of 1
Loading...
Thumbnail Image
Name:
TR09-06.pdf
Size:
279.37 KB
Format:
Adobe Portable Document Format