Abstract:We present a fast and accurate algorithm to solve Poisson problems in complex geometries, using regular Cartesian grids. We consider a variety of configurations, including Poisson problems with interfaces across which the solution is discontinuous (of the type arising in multi-fluid flows). The algorithm is based on a combination of the Correction Function Method (CFM) and Boundary Integral Methods (BIM). Interface and boundary conditions can be treated in a fast and accurate manner using boundary integral equations, and the associated BIM. Unfortunately, BIM can be costly when the solution is needed everywhere in a grid, e.g. fluid flow problems. We use the CFM to circumvent this issue. The solution from the BIM is used to rewrite the problem as a series of Poisson problems in rectangular domains—which requires the BIM solution at interfaces/boundaries only. These Poisson problems involve discontinuities at interfaces, of the type that the CFM can handle. Hence we use the CFM to solve them (to high order of accuracy) with finite differences and a Fast Fourier Transform based fast Poisson solver. We present 2-D examples of the algorithm applied to Poisson problems involving complex geometries, including cases in which the solution is discontinuous. We show that the algorithm produces solutions that converge with either 3rd or 4th order of accuracy, depending on the type of boundary condition and solution discontinuity.

Solution of 1D Poisson Equation with Neumann-Dirichlet and Dirichlet-Neumann Boundary Conditions, Using the Finite Difference Method

High Accuracy Non-Equidistant Method for Singular Perturbation Reaction-Diffusion Problem

An Approach to Numerically Solving the Poisson Equation

A Compact Exponential Scheme for Solving 1D Unsteady Convection-Diffusion Equation with Neumann Boundary Conditions

A direct method for the solution of poisson's equation with neumann boundary conditions on a staggered grid of arbitrary size

Fast Algorithm for the Three-Dimensional Poisson Equation in Infinite Domains

Neural-Network-Assisted Finite Difference Discretization for Numerical Solution of Partial Differential Equations

On direct methods for solving Poisson's equation

Compact Difference Schemes for Heat Equation with Neumann Boundary Conditions (II)

A finite difference scheme for two-dimensional singularly perturbed convection-diffusion problem with discontinuous source term

High order solution of Poisson problems with piecewise constant coefficients and interface jumps

A numerical study on solving a fractional time-space diffusion equation via the finite difference method

Solution of the mixed boundary problem for the Poisson equation on two-dimensional irregular domains

A locally second order symmetric method for discontinuous solution of Poisson's equation on uniform cartesian grids

Efficient finite element method for 2D singularly perturbed parabolic convection diffusion problems with discontinuous source term

An efficient neural-network and finite-difference hybrid method for elliptic interface problems with applications

Solution of the Dirichlet problem by a finite difference analog of the boundary integral equation

Stability analysis and numerical results for some schemes discretising 2D nonconstant coefficient advection–diffusion equations

Numerical method for a boundary value problem for a linear singularly perturbed parabolic delay differential equation of reaction-diffusion type with discontinuous source term

A fast numerical method for the conductivity of heterogeneous media with Dirichlet boundary conditions based on discrete sine–cosine transforms

A fast modified finite difference method for time fractional diffusion equations with weakly singular solution