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

by   Fredrik Fryklund, et al.

Boundary integral methods are attractive for solving homogeneous linear constant coefficient elliptic partial differential equations on complex geometries, since they can offer accurate solutions with a computational cost that is linear or close to linear in the number of discretization points on the boundary of the domain. However, these numerical methods are not straightforward to apply to time-dependent equations, which often arise in science and engineering. We address this problem with an integral equation-based solver for the advection-diffusion equation on moving and deforming geometries in two space dimensions. In this method, an adaptive high-order accurate time-stepping scheme based on semi-implicit spectral deferred correction is applied. One time-step then involves solving a sequence of non-homogeneous modified Helmholtz equations, a method known as elliptic marching. Our solution methodology utilizes several recently developed methods, including special purpose quadrature, a function extension technique and a spectral Ewald method for the modified Helmholtz kernel. Special care is also taken to handle the time-dependent geometries. The numerical method is tested through several numerical examples to demonstrate robustness, flexibility and accuracy


page 26

page 27

page 33


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

Integral equation based numerical methods are directly applicable to hom...

Matrix-oriented FEM formulation for stationary and time-dependent PDEs on x-normal domains

When numerical solution of elliptic and parabolic partial differential e...

ITVOLT: An Iterative Solver for Volterra Integral Equations with Application to the Time-Dependent Schrödinger Equation

We present a novel iterative method for solving Volterra integral equati...

High-order gas-kinetic scheme for radiation hydrodynamics in equilibrium-diffusion limit

In this paper, a high-order gas-kinetic scheme is developed for the equa...

Fast and High-order Accuracy Numerical Methods for Time-Dependent Nonlocal Problems in R^2

In this paper, we study the Crank-Nicolson method for temporal dimension...

The Smooth Forcing Extension Method: A High-Order Technique for Solving Elliptic Equations on Complex Domains

High-order numerical methods for solving elliptic equations over arbitra...

Spectral solver for Cauchy problems in polar coordinates using discrete Hankel transforms

We introduce a Fourier-Bessel-based spectral solver for Cauchy problems ...

Please sign up or login with your details

Forgot password? Click here to reset