Fast High-Order Integral Equation Methods for Solving Boundary Value Problems of Two Dimensional Heat Equation in Complex Geometry

Shaobo Wang, Shidong Jiang, Jing Wang

Research output: Contribution to journalArticlepeer-review

2 Scopus citations

Abstract

Efficient high-order integral equation methods have been developed for solving boundary value problems of the heat equation in complex geometry in two dimensions. First, the classical heat potential theory is applied to convert such problems to Volterra integral equations of the second kind via the heat layer potentials, where the unknowns are only on the space–time boundary. However, the heat layer potentials contain convolution integrals in both space and time whose direct evaluation requires O(NS2NT2) work and O(N S N T ) storage, where N S is the total number of discretization points on the spatial boundary and N T is the total number of time steps. In order to evaluate the heat layer potentials accurately and efficiently, they are split into two parts—the local part containing the temporal integration from t- δ to t and the history part containing the temporal integration from 0 to t- δ. The local part can be dealt with efficiently using conventional fast multipole type algorithms. For problems with complex stationary geometry, efficient separated sum-of-exponentials approximations are constructed for the heat kernel and used for the evaluation of the history part. Here all local and history kernels are compressed only once. The resulting algorithm is very efficient with quasilinear complexity in both space and time for both interior and exterior problems. For problems with complex moving geometry, the spectral Fourier approximation is applied for the heat kernel and nonuniform FFT is used to speed up the evaluation of the history part of heat layer potentials. The performance of both algorithms is demonstrated with several numerical examples.

Original languageEnglish (US)
Pages (from-to)787-808
Number of pages22
JournalJournal of Scientific Computing
Volume79
Issue number2
DOIs
StatePublished - May 15 2019

All Science Journal Classification (ASJC) codes

  • Software
  • Theoretical Computer Science
  • Numerical Analysis
  • Engineering(all)
  • Computational Theory and Mathematics
  • Computational Mathematics
  • Applied Mathematics

Keywords

  • Heat equation
  • Heat kernels
  • High-order methods
  • Integral equation methods
  • Nonuniform FFT
  • Sum-of-exponentials approximation

Fingerprint

Dive into the research topics of 'Fast High-Order Integral Equation Methods for Solving Boundary Value Problems of Two Dimensional Heat Equation in Complex Geometry'. Together they form a unique fingerprint.

Cite this