Source author record

Thomas Hagstrom

Thomas Hagstrom appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA Numerical Analysis math-ph math.AP math.MP Mathematical Software nlin.PS

Catalog footprint

What is connected

8works

7topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

An Energy-Based Discontinuous Galerkin Method with Tame CFL Numbers for the Wave Equation

We extend and analyze the energy-based discontinuous Galerkin method for second order wave equations on staggered and structured meshes. By combining spatial staggering with local time-stepping near boundaries, the method overcomes the typical numerical stiffness associated with high order piecewise polynomial approximations. In one space dimension with periodic boundary conditions and suitably chosen numerical fluxes, we prove bounds on the spatial operators that establish stability for CFL numbers $c \frac {Δt}{h} < C$ independent of order when stability-enhanced explicit time-stepping schemes of matching order are used. For problems on bounded domains and in higher dimensions we demonstrate numerically that one can march explicitly with large time steps at high order temporal and spatial accuracy.

preprint2020arXiv

An energy-based discontinuous Galerkin method for semilinear wave equations

We generalize the energy-based discontinuous Galerkin method proposed in [SIAM J. Num. Anal., 53(6):2705-2726, 2015.] to second-order semilinear wave equations. A stability and convergence analysis is presented along with numerical experiments demonstrating optimal convergence for certain choices of the interelement fluxes. Applications to the sine-Gordon equation include simulations of breathers, kink, and anti-kink solitons.

preprint2016arXiv

GPU Acceleration of Hermite Methods for the Simulation of Wave Propagation

The Hermite methods of Goodrich, Hagstrom, and Lorenz (2006) use Hermite interpolation to construct high order numerical methods for hyperbolic initial value problems. The structure of the method has several favorable features for parallel computing. In this work, we propose algorithms that take advantage of the many-core architecture of Graphics Processing Units. The algorithm exploits the compact stencil of Hermite methods and uses data structures that allow for efficient data load and stores. Additionally the highly localized evolution operator of Hermite methods allows us to combine multi-stage time-stepping methods within the new algorithms incurring minimal accesses of global memory. Using a scalar linear wave equation, we study the algorithm by considering Hermite interpolation and evolution as individual kernels and alternatively combined them into a monolithic kernel. For both approaches we demonstrate strategies to increase performance. Our numerical experiments show that although a two kernel approach allows for better performance on the hardware, a monolithic kernel can offer a comparable time to solution with less global memory usage.

preprint2015arXiv

On the stability of time-domain integral equations for acoustic wave propagation

We give a principled approach for the selection of a boundary integral, retarded potential representation for the solution of scattering problems for the wave equation in an exterior domain.

preprint2015arXiv

Variations on Hermite methods for wave propagation

Hermite methods, as introduced by Goodrich et al., combine Hermite interpolation and staggered (dual) grids to produce stable high order accurate schemes for the solution of hyperbolic PDEs. We introduce three variations of this Hermite method which do not involve time evolution on dual grids. Computational evidence is presented regarding stability, high order convergence, and dispersion/dissipation properties for each new method. Hermite methods may also be coupled to discontinuous Galerkin (DG) methods for additional geometric flexibility. An example illustrates the simplification of this coupling of this coupling for the Hermite methods.

preprint2013arXiv

The solution of the scalar wave equation in the exterior of a sphere

We derive new, explicit representations for the solution to the scalar wave equation in the exterior of a sphere, subject to either Dirichlet or Robin boundary conditions. Our formula leads to a stable and high-order numerical scheme that permits the evaluation of the solution at an arbitrary target, without the use of a spatial grid and without numerical dispersion error. In the process, we correct some errors in the analytic literature concerning the asymptotic behavior of the logarithmic derivative of the spherical modified Hankel function. We illustrate the performance of the method with several numerical examples.

preprint2012arXiv

On generalized discrete PML optimized for propagative and evanescent waves

We suggest a unified spectrally matched optimal grid approach for finite-difference and finite-element approximation of the PML. The new approach allows to combine optimal discrete absorption for both evanescent and propagative waves.

preprint2010arXiv

Analysis of supercontinuum generation under general dispersion characteristics and beyond the slowly varying envelope approximation

The generation of broadband supercontinua (SC) in air-silica microstructured fibers results from a delicate balance of dispersion and nonlinearity. We analyze two models aimed at better understanding SC. In the first one, we characterize linear dispersion in the Fourier domain from the calculated group velocity dispersion (GVD) without using a Taylor approximation for the propagation constant. Results of our numerical simulations are in a good agreement with experiments. A novel relevant length scale, namely the length for shock formation is introduced and its role is discussed. The second part shows similar dynamics for a model that goes beyond the slowly varying approximation for optical pulse propagation.