Mass Lumping for Constant-Density Acoustics

dc.contributor.authorSymes, William W.
dc.contributor.authorTerentyev, Igor
dc.date.accessioned2018-06-19T17:45:05Z
dc.date.available2018-06-19T17:45:05Z
dc.date.issued2009-04
dc.date.noteApril 2009
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.
dc.format.extent18 pp
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>.
dc.identifier.digitalTR09-06
dc.identifier.urihttps://hdl.handle.net/1911/102111
dc.language.isoeng
dc.titleMass Lumping for Constant-Density Acoustics
dc.typeTechnical report
dc.type.dcmiText
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