Related papers: Sparse spectral-tau method for the three-dimension…
Polynomial spectral methods produce fast, accurate, and flexible solvers for broad ranges of PDEs with one bounded dimension, where the incorporation of general boundary conditions is well understood. However, automating extensions to…
We present a multidomain spectral approach for Fuchsian ordinary differential equations in the particular case of the hypergeometric equation. Our hybrid approach uses Frobenius' method and Moebius transformations in the vicinity of each of…
We develop a sparse spectral method for a class of fractional differential equations, posed on $\mathbb{R}$, in one dimension. These equations can include sqrt-Laplacian, Hilbert, derivative and identity terms. The numerical method utilizes…
We consider the 2+1 and 3+1 scalar wave equations reduced via a helical Killing field, respectively referred to as the 2-dimensional and 3-dimensional helically reduced wave equation (HRWE). The HRWE serves as the fundamental model for the…
We show that a pseudospectral representation of the wavefunction using multiple spatial domains of variable size yields a highly accurate, yet efficient method to solve the time-dependent Schr\"odinger equation. The overall spatial domain…
We present a new solver for coupled nonlinear elliptic partial differential equations (PDEs). The solver is based on pseudo-spectral collocation with domain decomposition and can handle one- to three-dimensional problems. It has three…
This paper is concerned with the numerical solution of the third kind Volterra integral equations with non-smooth solutions based on the recursive approach of the spectral Tau method. To this end, a new set of the fractional version of…
A domain decomposition method for the solution of general variable-coefficient elliptic partial differential equations on regular domains is introduced. The method is based on tessellating the domain into overlapping thin slabs or shells,…
We develop a high-order, explicit method for acoustic scattering in three space dimensions based on a combined-field time-domain integral equation. The spatial discretization, of Nystr\"om type, uses Gaussian quadrature on panels combined…
There are a large number of methods for solving under-determined linear inverse problem. Many of them have very high time complexity for large datasets. We propose a new method called Two-Stage Sparse Representation (TSSR) to tackle this…
A 3-D inverse medium problem in the frequency domain is considered. Another name for this problem is Coefficient Inverse Problem. The goal is to reconstruct spatially distributed dielectric constants from scattering data. Potential…
This paper provides an efficient recursive approach of the spectral Tau method to approximate the solution of system of generalized Abel-Volterra integral equations. In this regards, we first investigate the existence, uniqueness as well as…
Two-dimensional terahertz spectroscopy (2DTS) is a low-frequency analogue of two-dimensional optical spectroscopy that is rapidly maturing as a probe of a wide variety of condensed matter systems. However, a persistent problem of 2DTS is…
We present a novel computational scheme to solve acoustic wave transmission problems over composite scatterers, i.e. penetrable obstacles possessing junctions or triple points. Our continuous problem is cast as a multiple traces time-domain…
We construct a pseudospectral method for the solution of time-dependent, non-linear partial differential equations on a three-dimensional spherical shell. The problem we address is the treatment of tensor fields on the sphere. As a test…
Second-order self-force calculations will be critical for modelling extreme-mass-ratio inspirals, and they are now known to have high accuracy even for binaries with mass ratios $\sim 1:10$. Many of the challenges facing these calculations…
We introduce and analyse a sparse spectral method for the solution of Volterra integral equations using bivariate orthogonal polynomials on a triangle domain. The sparsity of the Volterra operator on a weighted Jacobi basis is used to…
The wave equation for vectors and symmetric tensors in spherical coordinates is studied under the divergence-free constraint. We describe a numerical method, based on the spectral decomposition of vector/tensor components onto spherical…
We present a scalable 2D Galerkin spectral element method solution to the linearized potential flow radiation problem for wave induced forcing of a floating offshore structure. The pseudo-impulsive formulation of the problem is solved in…
Global spectral methods offer the potential to compute solutions of partial differential equations numerically to very high accuracy. In this work, we develop a novel global spectral method for linear partial differential equations on cubes…