Fast High-Order Integral Equation Methods for Solving Boundary Value Problems of Two Dimensional Heat Equation in Complex Geometry
Document Type
Article
Publication Date
5-15-2019
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.
Identifier
85056744518 (Scopus)
Publication Title
Journal of Scientific Computing
External Full Text Location
https://doi.org/10.1007/s10915-018-0872-x
ISSN
08857474
First Page
787
Last Page
808
Issue
2
Volume
79
Grant
1720405
Fund Ref
National Science Foundation
Recommended Citation
Wang, Shaobo; Jiang, Shidong; and Wang, Jing, "Fast High-Order Integral Equation Methods for Solving Boundary Value Problems of Two Dimensional Heat Equation in Complex Geometry" (2019). Faculty Publications. 7591.
https://digitalcommons.njit.edu/fac_pubs/7591
