Show simple item record

dc.contributor.authorMarkakis, Cen_US
dc.contributor.authorBarack, Len_US
dc.description9 pages, 5 figures
dc.description9 pages, 5 figuresen_US
dc.description9 pages, 5 figuresen_US
dc.description.abstractHigh order finite-difference or spectral methods are typically problematic in approximating a function with a jump discontinuity. Some common remedies come with a cost in accuracy near discontinuities, or in computational cost, or in complexity of implementation. However, for certain classes of problems involving piecewise analytic functions, the jump in the function and its derivatives are known or easy to compute. We show that high-order or spectral accuracy can then be recovered by simply adding to the Lagrange interpolation formula a linear combination of the jumps. Discretizations developed for smooth problems are thus easily extended to nonsmooth problems. Furthermore, in the context of one-dimensional finite-difference or pseudospectral discretizations, numerical integration and differentiation amount to matrix multiplication. We construct the matrices for such operations, in the presence of known discontinuities, by operating on the corrected Lagrange formula. In a method-of-lines framework, this provides a simple and efficient way to obtain solutions with moving discontinuities to evolution partial differential equations.en_US
dc.titleHigh-order difference and pseudospectral methods for discontinuous problemsen_US
dc.rights.holder© The Author(s) 2014
pubs.notesNo embargoen_US

Files in this item


There are no files associated with this item.

This item appears in the following Collection(s)

Show simple item record