Related papers: A fast multi-resolution lattice Green's function m…
We develop efficient and high-order accurate finite difference methods for elliptic partial differential equations in complex geometry in the Difference Potentials framework. The main novelty of the developed schemes is the use of local…
Fast Fourier transform (FFT) based methods have turned out to be an effective computational approach for numerical homogenisation. In particular, Fourier-Galerkin methods are computational methods for partial differential equations that are…
High-order numerical methods for solving elliptic equations over arbitrary domains typically require specialized machinery, such as high-quality conforming grids for finite elements method, and quadrature rules for boundary integral…
We present a fast direct solver for boundary integral equations on complex surfaces in three dimensions using an extension of the recently introduced recursive strong skeletonization scheme. For problems that are not highly oscillatory, our…
This paper presents a windowed Green function (WGF) method for the numerical solution of problems of elastic scattering by "locally-rough surfaces" (i.e., local perturbations of a half space), under either Dirichlet or Neumann boundary…
We introduce a finite element method for numerical upscaling of second order elliptic equations with highly heterogeneous coefficients. The method is based on a mixed formulation of the problem and the concepts of the domain decomposition…
To numerically solve a generic elliptic equation on two-dimensional domains with rectangular Cartesian grids, we propose a cut-cell geometric multigrid method that features (1) general algorithmic steps that apply to two-dimensional…
We propose an unfitted finite element method for numerically solving the time-harmonic Maxwell equations on a smooth domain. The model problem involves a Lagrangian multiplier to relax the divergence constraint of the vector unknown. The…
In this study, we address the challenge of obtaining a Green's function operator for linear partial differential equations (PDEs). The Green's function is well-sought after due to its ability to directly map inputs to solutions, bypassing…
We present an Augmented Lagrangian formulation and its real-space implementation for non-periodic orbital-free Density Functional Theory (OF-DFT) calculations. In particular, we rewrite the constrained minimization problem of OF-DFT as a…
We introduce a generalized finite difference method for solving a large range of fully nonlinear elliptic partial differential equations in three dimensions. Methods are based on Cartesian grids, augmented by additional points carefully…
This paper presents a new finite difference method, called {\varphi}-FD, inspired by the {\phi}-FEM approach for solving elliptic partial differential equations (PDEs) on general geometries. The proposed method uses Cartesian grids,…
This article presents an $O(N\log N)$ method for numerical solution of Maxwell's equations for dielectric scatterers using a 3D boundary integral equation (BIE) method. The underlying BIE method used is based on a hybrid…
Many boundary element integral equation kernels are based on the Green's functions of the Laplace and Helmholtz equations in three dimensions. These include, for example, the Laplace, Helmholtz, elasticity, Stokes, and Maxwell's equations.…
In this paper, we develop two fast implicit difference schemes for solving a class of variable-coefficient time-space fractional diffusion equations with integral fractional Laplacian (IFL). The proposed schemes utilize the graded $L1$…
The Laplace transform method for solving of a wide class of initial value problems for fractional differential equations is introduced. The method is based on the Laplace transform of the Mittag-Leffler function in two parameters. To extend…
A Lagrangian numerical scheme for solving nonlinear degenerate Fokker-Planck equations in space dimensions $d\ge2$ is presented. It applies to a large class of nonlinear diffusion equations, whose dynamics are driven by internal energies…
Localized collocation methods based on radial basis functions (RBFs) for elliptic problems appear to be non-robust in the presence of Neumann boundary conditions. In this paper we overcome this issue by formulating the RBF-generated finite…
Efficient computation of lattice defect geometries such as point defects, dislocations, disconnections, grain boundaries, interfaces and free surfaces requires accurate coupling of displacements near the defect to the long-range elastic…
We study the lattice Green's function (LGF) of the screened Poisson equation on a two-dimensional rectangular lattice. This LGF arises in numerical analysis, random walks, solid-state physics, and other fields. Its defining characteristic…