Related papers: A fourth-order compact time-splitting Fourier pseu…
Spectral deferred correction (SDC) methods are an attractive approach to iteratively computing collocation solutions to an ODE by performing so-called sweeps with a low-order time stepping method. SDC allows to easily construct high order…
We present new splitting methods designed for the numerical integration of near-integrable Hamiltonian systems, and in particular for planetary N-body problems, when one is interested in very accurate results over a large time span. We…
In this paper, we present a class of high-order and efficient compact difference schemes for nonlinear convection diffusion equations, which can preserve both bounds and mass. For the one-dimensional problem, we first introduce a high-order…
This paper presents a novel, efficient, high-order accurate, and stable spectral element-based model for computing the complete three-dimensional linear radiation and diffraction problem for floating offshore structures. We present a…
In this paper, we discuss the nonlinear stability and convergence of a fully discrete Fourier pseudospectral method coupled with a specially designed second order time-stepping for the numerical solution of the "good" Boussinesq equation.…
In this paper, we consider a second-order scalar auxiliary variable (SAV) Fourier spectral method to solve the nonlinear fractional generalized wave equation. Unconditional energy conservation or dissipation properties of the fully discrete…
Four-component Dirac Hartree--Fock is an accurate mean-field method for treating molecular systems where relativistic effects are important. However, the computational cost and complexity of the two-electron interaction makes this method…
We present a practical algorithm based on symplectic splitting methods to integrate numerically in time the Schr\"odinger equation. When discretized in space, the Schr\"odinger equation can be recast as a classical Hamiltonian system…
The present work proposes a second-order time splitting scheme for a linear dispersive equation with a variable advection coefficient subject to transparent boundary conditions. For its spatial discretization, a dual Petrov--Galerkin method…
Dispersive time-decay estimates are proved for a one-parameter family of one-dimensional Dirac Hamiltonians with dislocations; these are operators which interpolate between two phase-shifted massive Dirac Hamiltonians at $x=+\infty$ and…
We compare the long-time error bounds and spatial resolution of finite difference methods with different spatial discretizations for the Dirac equation with small electromagnetic potentials characterized by $\varepsilon \in (0, 1]$ a…
The Allen-Cahn equation is solved numerically by operator splitting Fourier spectral methods. The basic idea of the operator splitting method is to decompose the original problem into sub-equations and compose the approximate solution of…
In this paper, we present a new SDC scheme for solving semi-explicit DAEs with the ability to be parallelized in which only the differential equations are numerically integrated is presented. In Shu et al. (2007) it was shown that SDC for…
The Dirac equation may be thought as originating from a theory of five-dimensional (5D) space-time. We define a special 5D Clifford algebra and introduce a spin-1/2 constraint equation to describe null propagation in a 5D space-time…
The problem of scattering of neutral fermions in two-dimensional space-time is approached with a pseudoscalar potential step in the Dirac equation. Some unexpected aspects of the solutions beyond the absence of Klein\'{}s paradox are…
We formulate the Dirac equation for a massive neutral spin-half particle on a rotating black hole spacetime, and we consider its (quasi)bound states: gravitationally-trapped modes which are regular across the future event horizon. These…
Numerical relativity has traditionally been pursued via finite differencing. Here we explore pseudospectral collocation (PSC) as an alternative to finite differencing, focusing particularly on the solution of the Hamiltonian constraint (an…
A Fourier transform method is introduced for a class of hybrid time-frequency methods that solve the acoustic scattering problem in regimes where the solution exhibits both highly oscillatory behavior and slow decay in time. This extends…
We approximate the solution for the time dependent Schr\"odinger equation (TDSE) in two steps. We first use a pseudo-spectral collocation method that uses samples of functions on rank-1 or rank-r lattice points with unitary Fourier…
We present a novel numerical method and algorithm for the solution of the 3D axially symmetric time-dependent Schr\"odinger equation in cylindrical coordinates, involving singular Coulomb potential terms besides a smooth time-dependent…