An integral equation based numerical method for the forced heat equation on complex domains

by   Fredrik Fryklund, et al.

Integral equation based numerical methods are directly applicable to homogeneous elliptic PDEs, and offer the ability to solve these with high accuracy and speed on complex domains. In this paper, extensions to problems with inhomogeneous source terms and time dependent PDEs, such as the heat equation, have been introduced. One such approach for the heat equation is to first discretise in time, and in each time-step solve a so-called modified Helmholtz equation with a parameter depending on the time step size. The modified Helmholtz equation is then split into two parts: a homogeneous part solved with a boundary integral method and a particular part, where the solution is obtained by evaluating a volume potential over the inhomogeneous source term over a simple domain. In this work, we introduce two components which are critical for the success of this approach: a method to efficiently compute a high-regularity extension of a function outside the domain where it is defined, and a special quadrature method to accurately evaluate singular and nearly singular integrals in the integral formulation of the modified Helmholtz equation for all time step sizes.


page 17

page 21

page 22

page 29

page 31


An integral equation method for the advection-diffusion equation on time-dependent domains in the plane

Boundary integral methods are attractive for solving homogeneous linear ...

An Integral Equation Method for the Cahn-Hilliard Equation in the Wetting Problem

We present an integral equation approach to solve the Cahn-Hilliard equa...

An integral-like numerical approach for solving Burgers' equation

An integral-like approach established on spline polynomial interpolation...

Parallel algorithm for numerical solution of heat equation in complex cylindrical domain

In this article we present a parallel algorithm for simulation of the he...

Immersed Boundary Double Layer Method

The Immersed Boundary (IB) method of Peskin (J. Comput. Phys., 1977) is ...

An extended Krylov subspace method for decoding edge-based compressed images by homogeneous diffusion

The heat equation is often used in order to inpaint dropped data in inpa...

Efficient extraction of resonant states in systems with defects

We introduce a new numerical method to compute resonances induced by loc...

Please sign up or login with your details

Forgot password? Click here to reset