Hybrid asymptotic/numerical methods for the evaluation of layer heat potentials in two dimensions

Research output: Contribution to journalArticle

Abstract

We present a hybrid asymptotic/numerical method for the accurate computation of single- and double-layer heat potentials in two dimensions. It has been shown in previous work that simple quadrature schemes suffer from a phenomenon called “geometrically induced stiffness,” meaning that formally high-order accurate methods require excessively small time steps before the rapid convergence rate is observed. This can be overcome by analytic integration in time, requiring the evaluation of a collection of spatial boundary integral operators with non-physical, weakly singular kernels. In our hybrid scheme, we combine a local asymptotic approximation with the evaluation of a few boundary integral operators involving only Gaussian kernels, which are easily accelerated by a new version of the fast Gauss transform. This new scheme is robust, avoids geometrically induced stiffness, and is easy to use in the presence of moving geometries. Its extension to three dimensions is natural and straightforward, and should permit layer heat potentials to become flexible and powerful tools for modeling diffusion processes.

Original languageEnglish (US)
JournalAdvances in Computational Mathematics
DOIs
StateAccepted/In press - Jan 1 2018

Fingerprint

Asymptotic Methods
Numerical methods
Two Dimensions
Heat
Boundary Integral
Numerical Methods
Stiffness
Integral Operator
Evaluation
Gauss Transform
Weakly Singular Kernel
Gaussian Kernel
Local Approximation
High-order Methods
Asymptotic Approximation
Quadrature
Diffusion Process
Geometry
Three-dimension
Convergence Rate

Keywords

  • Gauss transform
  • Geometrically induced stiffness
  • Hybrid asymptotic/numerical method

ASJC Scopus subject areas

  • Computational Mathematics
  • Applied Mathematics

Cite this

@article{74cdbacffa46468bacd5f731ac05e10c,
title = "Hybrid asymptotic/numerical methods for the evaluation of layer heat potentials in two dimensions",
abstract = "We present a hybrid asymptotic/numerical method for the accurate computation of single- and double-layer heat potentials in two dimensions. It has been shown in previous work that simple quadrature schemes suffer from a phenomenon called “geometrically induced stiffness,” meaning that formally high-order accurate methods require excessively small time steps before the rapid convergence rate is observed. This can be overcome by analytic integration in time, requiring the evaluation of a collection of spatial boundary integral operators with non-physical, weakly singular kernels. In our hybrid scheme, we combine a local asymptotic approximation with the evaluation of a few boundary integral operators involving only Gaussian kernels, which are easily accelerated by a new version of the fast Gauss transform. This new scheme is robust, avoids geometrically induced stiffness, and is easy to use in the presence of moving geometries. Its extension to three dimensions is natural and straightforward, and should permit layer heat potentials to become flexible and powerful tools for modeling diffusion processes.",
keywords = "Gauss transform, Geometrically induced stiffness, Hybrid asymptotic/numerical method",
author = "Jun Wang and Leslie Greengard",
year = "2018",
month = "1",
day = "1",
doi = "10.1007/s10444-018-9641-5",
language = "English (US)",
journal = "Advances in Computational Mathematics",
issn = "1019-7168",
publisher = "Springer Netherlands",

}

TY - JOUR

T1 - Hybrid asymptotic/numerical methods for the evaluation of layer heat potentials in two dimensions

AU - Wang, Jun

AU - Greengard, Leslie

PY - 2018/1/1

Y1 - 2018/1/1

N2 - We present a hybrid asymptotic/numerical method for the accurate computation of single- and double-layer heat potentials in two dimensions. It has been shown in previous work that simple quadrature schemes suffer from a phenomenon called “geometrically induced stiffness,” meaning that formally high-order accurate methods require excessively small time steps before the rapid convergence rate is observed. This can be overcome by analytic integration in time, requiring the evaluation of a collection of spatial boundary integral operators with non-physical, weakly singular kernels. In our hybrid scheme, we combine a local asymptotic approximation with the evaluation of a few boundary integral operators involving only Gaussian kernels, which are easily accelerated by a new version of the fast Gauss transform. This new scheme is robust, avoids geometrically induced stiffness, and is easy to use in the presence of moving geometries. Its extension to three dimensions is natural and straightforward, and should permit layer heat potentials to become flexible and powerful tools for modeling diffusion processes.

AB - We present a hybrid asymptotic/numerical method for the accurate computation of single- and double-layer heat potentials in two dimensions. It has been shown in previous work that simple quadrature schemes suffer from a phenomenon called “geometrically induced stiffness,” meaning that formally high-order accurate methods require excessively small time steps before the rapid convergence rate is observed. This can be overcome by analytic integration in time, requiring the evaluation of a collection of spatial boundary integral operators with non-physical, weakly singular kernels. In our hybrid scheme, we combine a local asymptotic approximation with the evaluation of a few boundary integral operators involving only Gaussian kernels, which are easily accelerated by a new version of the fast Gauss transform. This new scheme is robust, avoids geometrically induced stiffness, and is easy to use in the presence of moving geometries. Its extension to three dimensions is natural and straightforward, and should permit layer heat potentials to become flexible and powerful tools for modeling diffusion processes.

KW - Gauss transform

KW - Geometrically induced stiffness

KW - Hybrid asymptotic/numerical method

UR - http://www.scopus.com/inward/record.url?scp=85055934281&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85055934281&partnerID=8YFLogxK

U2 - 10.1007/s10444-018-9641-5

DO - 10.1007/s10444-018-9641-5

M3 - Article

AN - SCOPUS:85055934281

JO - Advances in Computational Mathematics

JF - Advances in Computational Mathematics

SN - 1019-7168

ER -