Related papers: A generalized optimal fourth-order finite differen…
We develop a stable finite difference method for the elastic wave equation in bounded media, where the material properties can be discontinuous at curved interfaces. The governing equation is discretized in second order form by a fourth or…
This paper introduces a numerical approach to solve singularly perturbed convection diffusion boundary value problems for second-order ordinary differential equations that feature a small positive parameter {\epsilon} multiplying the…
We present a novel approach for high-order accurate numerical differentiation on unstructured meshes of quadrilateral elements. To differentiate a given function, an auxiliary function with greater smoothness properties is defined which…
For scattering problems of time-harmonic waves, the boundary integral equation (BIE) methods are highly competitive, since they are formulated on lower-dimension boundaries or interfaces, and can automatically satisfy outgoing radiation…
We investigate error bounds for numerical solutions of divergence structure linear elliptic PDEs on compact manifolds without boundary. Our focus is on a class of monotone finite difference approximations, which provide a strong form of…
Two-fluid relativistic plasma flow equations combine the equations of relativistic hydrodynamics with Maxwell's equations for electromagnetic fields, which involve divergence constraints for the magnetic and electric fields. When developing…
Computing the dispersion relation for two-dimensional photonic crystals is a notoriously challenging task: It involves solving parameterized Helmholtz eigenvalue problems with high-contrast coefficients. To resolve the challenge, we propose…
Efficient frequency-domain Full Waveform Inversion (FWI) of long-offset/wide-azimuth node data can be designed with a few discrete frequencies. However, 3D frequency-domain seismic modeling remains challenging since it requires solving a…
The nonlinear Helmholtz equation (NLH) models the propagation of electromagnetic waves in Kerr media, and describes a range of important phenomena in nonlinear optics and in other areas. In our previous work, we developed a fourth order…
A novel fourth-order finite difference formula coupling the Crank-Nicolson explicit linearized method is proposed to solve Riesz space fractional nonlinear reaction-diffusion equations in two dimensions. Theoretically, under the Lipschitz…
The Helmholtz equation is related to seismic exploration, sonar, antennas, and medical imaging applications. It is one of the most challenging problems to solve in terms of accuracy and convergence due to the scalability issues of the…
We consider fourth order singularly perturbed boundary value problems with two small parameters, and the approximation of their solution by the $hp$ version of the Finite Element Method on the {\emph{Spectral Boundary Layer}} mesh from…
In this paper we consider the numerical solution of the Hamiltonian wave equation in two spatial dimension. We use the Mimetic Finite Difference (MFD) method to approximate the continuous problem combined with a symplectic integration in…
We present strongly stable semi-discrete finite difference approximations to the quarter space problem (x>0, t>0) for the first order in time, second order in space wave equation with a shift term. We consider space-like (pure outflow) and…
We consider the numerical solution of partial differential equations with coefficients that are strongly heterogeneous in space. We provide an overview of higher-order localized orthogonal decomposition (LOD) methods for the elliptic…
Due to the lack of corresponding analysis on appropriate mapping operator between two grids, high-order two-grid difference algorithms are rarely studied. In this paper, we firstly discuss the boundedness of a local bi-cubic Lagrange…
Algorithmic and architecture-oriented optimizations are essential for achieving performance worthy of anticipated energy-austere exascale systems. In this paper, we present an extreme scale FMM-accelerated boundary integral equation solver…
In recent work, Li et al.\ (Comm.\ Math.\ Sci., 7:81-107, 2009) developed a diffuse-domain method (DDM) for solving partial differential equations in complex, dynamic geometries with Dirichlet, Neumann, and Robin boundary conditions. The…
A fast method is proposed for solving the high frequency Helmholtz equation. The building block of the new fast method is an overlapping source transfer domain decomposition method for layered medium, which is an extension of the source…
The main task in this paper is to prove that the perfectly matched layers (PML) method converges exponentially with respect to the PML parameter, for scattering problems with periodic surfaces. In [5], a linear convergence is proved for the…