Related papers: A fast algorithm for the wave equation using time-…
Standard fractional models on manifolds often conflate geometric anisotropy with medium heterogeneity. In this Letter, we overcome this rigidity by deriving the fundamental solution for a weighted space-time fractional ultrahyperbolic…
In this article we present a new family of high order accurate Arbitrary Lagrangian-Eulerian one-step WENO finite volume schemes for the solution of stiff hyperbolic balance laws. High order accuracy in space is obtained with a standard…
Two fundamental difficulties are encountered in the numerical evaluation of time-dependent layer potentials. One is the quadratic cost of history dependence, which has been successfully addressed by splitting the potentials into two parts -…
We propose two easy-to-implement fast algorithms based on moment-matching to compute the nonlocal potential $\varphi(\textbf{x})=(U\ast \rho)(\textbf{x})$ on bounded domain, where the kernel $U$ is singular at the origin and the density…
This paper presents a solution to an initial value problem for the 1-dimensional wave equation on time scales through the application of a Fourier transform and its inverse via contour integrals. The time scale of the spatial dimension is…
This paper proposes a frequency-time hybrid solver for the time-dependent wave equation in two-dimensional interior spatial domains. The approach relies on four main elements, namely, 1) A multiple scattering strategy that decomposes a…
This paper, Part I in a two-part series, presents (i) A simple and highly efficient algorithm for evaluation of quasi-periodic Green functions, as well as (ii) An associated boundary-integral equation method for the numerical solution of…
The pointwise space-time behavior of the Green's function of the three-dimensional modified Vlasov-Poisson-Boltzmann system is studied in this paper. It is shown that the Green's function has a decomposition of the macroscopic diffusive…
Over the last two decades, several fast, robust, and high-order accurate methods have been developed for solving the Poisson equation in complicated geometry using potential theory. In this approach, rather than discretizing the partial…
In this paper, we present a fast boundary integral method accelerated by the fast multipole method (FMM) for acoustic wave scattering governed by the scalar Helmholtz equation in multi-layered two-dimensional media. Multiple scatterers are…
A high-performance shooting algorithm is developed to compute time-periodic solutions of the free-surface Euler equations with spectral accuracy in double and quadruple precision. The method is used to study resonance and its effect on…
The nonlinear Fourier transform, which is also known as the forward scattering transform, decomposes a periodic signal into nonlinearly interacting waves. In contrast to the common Fourier transform, these waves no longer have to be…
We develop a new algorithm for the Brownian dynamics of soft matter systems that evolves time by spatially correlated Monte Carlo moves. The algorithm uses vector wavelets as its basic moves and produces hydrodynamics in the low Reynolds…
Frequency domain Mie solutions to scattering from spheres have been used for a long time. However, deriving their transient analogue is a challenge as it involves an inverse Fourier transform of the spherical Hankel functions (and their…
Hawking radiation has become experimentally testable thanks to the many analogue systems which mimic the effects of the event horizon on wave propagation. These systems are typically dominated by dispersion, and give rise to a numerically…
This paper provides a new analytical method to obtain Green's functions of linear dispersive partial differential equations. The Euler-Bernoulli beam equation and the one-dimensional heat conduction equation (dissipation equation) under…
In calculating integral or discrete transforms, use has been made of fast algorithms for multiplying vectors by matrices whose elements are specified as values of special (Chebyshev, Legendre, Laguerre, etc.) functions. The currently…
We present a fast algorithm for the subset convolution problem: given functions f and g defined on the lattice of subsets of an n-element set N, compute their subset convolution f*g, defined for all S\subseteq N by (f * g)(S) = \sum_{T…
During the past three years, Wapenaar, Snieder, Broggini and others have developed an algorithm to compute the Green's function for any point inside a medium to points on the surface from measurements on that surface only. Their algorithm…
We introduce an efficient boundary-adapted spectral method for peridynamic diffusion problems with arbitrary boundary conditions. The spectral approach transforms the convolution integral in the peridynamic formulation into a multiplication…