Related papers: Method of Lines Transpose: A Fast Implicit Wave Pr…
We describe a fourth-order accurate finite-difference time-domain scheme for solving dispersive Maxwell's equations with nonlinear multi-level carrier kinetics models. The scheme is based on an efficient single-step three time-level…
This paper is devoted to the efficient numerical solution of the Helmholtz equation in a two- or three-dimensional rectangular domain with an absorbing boundary condition (ABC). The Helmholtz problem is discretized by standard bilinear and…
In this paper we study the stability of explicit finite difference discretizations of linear advection-diffusion equations (ADE) with arbitrary order of accuracy in the context of method of lines. The analysis first focuses on the stability…
We deal with an initial-boundary value problem for the multidimensional acoustic wave equation, with the variable speed of sound. For a three-level semi-explicit in time higher-order vector compact scheme, we prove stability and derive 4th…
A spectral method is developed for the direct solution of linear ordinary differential equations with variable coefficients. The method leads to matrices which are almost banded, and a numerical solver is presented that takes O(m^2n)…
We present a novel computational methodology for solving the scalar nonlinear Helmholtz equation (NLH) that governs the propagation of laser light in Kerr dielectrics. The methodology addresses two well-known challenges in nonlinear optics:…
In this paper, a high-order exponential scheme is developed to solve the 1D unsteady convection-diffusion equation with Neumann boundary conditions. The present method applies fourth-order compact exponential difference scheme in spatial…
Given a Hilbert space, we investigate the well-posedness of the Cauchy problem for the wave equation for operators with discrete non-negative spectrum acting on it. We consider the cases when the time-dependent propagation speed is regular,…
We present an implicit-explicit (IMEX) scheme for semilinear wave equations with strong damping. By treating the nonlinear, nonstiff term explicitly and the linear, stiff part implicitly, we obtain a method which is not only unconditionally…
Recently, a nonlinear Poisson equation has been introduced to model nonlinear and nonlocal hyperpolarization effects in electrostatic solute-solvent interaction for biomolecular solvation analysis. Due to a strong nonlinearity associated…
In numerical linear algebra, considerable effort has been devoted to obtaining faster algorithms for linear systems whose underlying matrices exhibit structural properties. A prominent success story is the method of generalized nested…
This work presents the discontinuous Galerkin discretization of the consistent splitting scheme proposed by Liu [J. Liu, J. Comp. Phys., 228(19), 2009]. The method enforces the divergence-free constraint implicitly, removing…
The dispersive character of the Hall-MHD solutions, in particular the whistler waves, is a strong restriction to numerical treatments of this system. Numerical stability demands a time step dependence of the form $\Delta t\propto (\Delta…
This paper deals with the construction and analysis of two integrators for (semi-linear) second-order partial differential-algebraic equations of semi-explicit type. More precisely, we consider an implicit-explicit Crank-Nicolson scheme as…
We propose a boundary integral formulation for the dynamic problem of electromagnetic scattering and transmission by homogeneous dielectric obstacles. In the spirit of Costabel and Stephan, we use the transmission conditions to reduce the…
We formulate a new family of high order on-surface radiation conditions to approximate the outgoing solution to the Helmholtz equation in exterior domains. Motivated by the pseudo-differential expansion of the Dirichlet-to-Neumann operator…
Numerical discretization of the large-scale Maxwell's equations leads to an ill-conditioned linear system that is challenging to solve. The key requirement for successive solutions of this linear system is to choose an efficient solver. In…
In this paper, we propose high order numerical methods to solve a 2D advection diffusion equation, in the highly oscillatory regime. We use an integrator strategy that allows the construction of arbitrary high-order schemes {leading} to an…
We introduce a fast direct solver for variable-coefficient elliptic partial differential equations on surfaces based on the hierarchical Poincar\'e-Steklov method. The method takes as input an unstructured, high-order quadrilateral mesh of…
The famous Fisher-KPP reaction diffusion model combines linear diffusion with the typical Fisher-KPP reaction term, and appears in a number of relevant applications. It is remarkable as a mathematical model since, in the case of linear…