Related papers: A fully discrete Calderon Calculus for two dimensi…
We present a fast direct solver for the volume scattering problem of the Helmholtz equation. The algorithm is faster than existing methods. Moreover, discretization for our method is much simpler and more accurate than that for finite…
We show that any second order linear ordinary diffrential equation with constant coefficients (including the damped and undumped harmonic oscillator equation) admits an exact discretization, i.e., there exists a difference equation whose…
We investigate thin-slit diffraction problems for two-dimensional lattice waves. The peculiar structure allows us to consider the problems on the semi-infinite triangular lattice, consequently, we study Dirichlet problems for the…
We present two semidiscretizations of the Camassa-Holm equation in periodic domains based on variational formulations and energy conservation. The first is a periodic version of an existing conservative multipeakon method on the real line,…
Time-harmonic solutions to the wave equation can be computed in the frequency or in the time domain. In the frequency domain, one solves a discretized Helmholtz equation, while in the time domain, the periodic solutions to a discretized…
We study a class of fully-discrete schemes for the numerical approximation of solutions of stochastic Cahn--Hilliard equations with cubic nonlinearity and driven by additive noise. The spatial (resp. temporal) discretization is performed…
A novel discontinuous Galerkin (DG) method is developed to solve time-dependent bi-harmonic type equations involving fourth derivatives in one and multiple space dimensions. We present the spatial DG discretization based on a mixed…
Efficient and energy stable high order time marching schemes are very important but not easy to construct for the study of nonlinear phase dynamics. In this paper, we propose and study two linearly stabilized second order semi-implicit…
We consider the application of the WaveHoltz iteration to time-harmonic elastic wave equations with energy conserving boundary conditions. The original WaveHoltz iteration for acoustic Helmholtz problems is a fixed-point iteration that…
In this paper, we consider a semiclassical version of the fractional Klein-Gordon equation on the lattice, $h{\mathbb{Z}}^n.$ Contrary to the Euclidean case that was considered in [2], the discrete fractional Klein-Gordon equation is…
In this paper we propose an explicit fully discrete scheme to numerically solve the stochastic Allen-Cahn equation. The spatial discretization is done by a spectral Galerkin method, followed by the temporal discretization by a tamed…
The paper describes a formulation of discrete scalar wave propagation in an inhomogeneous medium by the use of elementary processes obeying a discrete Huygens' principle and satisfying fundamental symmetries such as time-reversal,…
We consider the discretization of electromagnetic wave propagation problems by a discontinuous Galerkin Method based on Trefftz polynomials. This method fits into an abstract framework for space-time discontinuous Galerkin methods for which…
We develop a new dispersion minimizing compact finite difference scheme for the Helmholtz equation in 2 and 3 dimensions. The scheme is based on a newly developed ray theory for difference equations. A discrete Helmholtz operator and a…
The complex Helmholtz equation $(\Delta + k^2)u=f$ (where $k\in{\mathbb R},u(\cdot),f(\cdot)\in{\mathbb C}$) is a mainstay of computational wave simulation. Despite its apparent simplicity, efficient numerical methods are challenging to…
In this paper, we use an implicit two-derivative deferred correction time discretization approach and combine it with a spatial discretization of the discontinuous Galerkin spectral element method to solve (non-)linear PDEs. The resulting…
We revisit second-order-in-time space-time discretizations of the linear and semilinear wave equations by establishing precise equivalences with first-order-in-time formulations. Focusing on schemes using continuous piecewise-polynomial…
Gravitational perturbations in a Kerr black hole background can not be decomposed into simple tensor harmonics in the time domain. Here, we make the mode decomposition only in the azimuthal direction and discuss the resulting…
We develop a convergence theory of space-time discretizations for the linear, 2nd-order wave equation in polygonal domains $\Omega\subset\mathbb{R}^2$, possibly occupied by piecewise homogeneous media with different propagation speeds.…
The simulation of high-dimensional problems with manageable computational resource represents a long-standing challenge. In a series of our recent work [25, 17, 18, 24], a class of sparse grid DG methods has been formulated for solving…