Related papers: Spectral methods for the wave equation in second-o…
Second-order gravitational self-force theory has recently led to the breakthrough calculation of ``first post-adiabatic'' (1PA) compact-binary waveforms [Phys. Rev. Lett. 130, 241402 (2023)]. The computations underlying those waveforms…
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…
New implicit and implicit-explicit time-stepping methods for the wave equation in second-order form are described with application to two and three-dimensional problems discretized on overset grids. The implicit schemes are single step,…
We present a parametric family of semi-implicit second order accurate numerical methods for non-conservative and conservative advection equation for which the numerical solutions can be obtained in a fixed number of forward and backward…
This study employs spectral methods to capture the behaviour of wave equation with dispersive-nonlinearity. We describe the evolution of hump initial data and track the conservation of the mass and energy functionals. The…
Spectral methods are now common in the solution of ordinary differential eigenvalue problems in a wide variety of fields, such as in the computation of black hole quasinormal modes. Most of these spectral codes are based on standard…
A new development of the ``monodromy transform'' method for analysis of hyperbolic as well as elliptic integrable reductions of Einstein equations is presented. Compatibility conditions for some alternative representations of the…
While not generally a conservation law, any symmetry of the equations of motion implies a useful reduction of any second-order equationto a first-order equation between invariants, whose solutions (first integrals) can then be integrated by…
Here we present a numerical method for finding non-hydrostatic coastal-trapped wave and instability solutions to the non-hydrostatic Boussinesq equations in the presence of a background flow and complicated coastal topography. We use…
We propose some finite element schemes to solve a class of fourth-order nonlinear PDEs, which include the vector-valued Landau--Lifshitz--Baryakhtar equation, the Swift--Hohenberg equation, and various Cahn--Hilliard-type equations with…
Spectral methods are well suited for solving hydrodynamic problems in which the self-gravity of the flow needs to be considered. Because Poisson's equation is linear, the numerical solution for the gravitational potential for each…
In this work numerical methods for solving Einstein's equations are developed and applied to the study of inhomogeneous cosmological models. A two-dimensional computer code is described which implements two advanced numerical methods:…
We present a numerical scheme for determining hyperboloidal initial data sets for the conformal field equations by using pseudo-spectral methods. This problem is split into two parts. The first step is the determination of a suitable…
This paper generalizes the earlier work on the energy-based discontinuous Galerkin method for second-order wave equations to fourth-order semilinear wave equations. We first rewrite the problem into a system with a second-order spatial…
We modify the well-known interior penalty finite element discretization method so that it allows for element-by-element assembly. This is possible due to the introduction of additional unknowns associated with the interfaces between…
An efficient numerical method to compute solitary wave solutions to the free surface Euler equations is reported. It is based on the conformal mapping technique combined with an efficient Fourier pseudo-spectral method. The resulting…
The ground state energy of a many-electron system can be approximated by an variational approach in which the total energy of the system is minimized with respect to one and two-body reduced density matrices (RDM) instead of many-electron…
A detailed account is given on approximation schemes to the Einstein theory of general relativity where the iteration starts from the Newton theory of gravity. Two different coordinate conditions are used to represent the Einstein field…
We present a high order, Fourier penalty method for the Maxwell's equations in the vicinity of perfect electric conductor boundary conditions. The approach relies on extending the smooth non-periodic domain of the equations to a periodic…
This paper develops an explicit spectral representation for solutions of a one-dimensional linear wave equation with a constant time delay. The model is considered on a bounded interval with non-homogeneous Dirichlet boundary data and a…