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

From MaRDI portal
Publication:2000499



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.


The authors develop a new method for the evaluation of layer heat potentials in two dimensions. Using an exponential change of variables, the phenomenon of ``geometrically induced stiffness, which prevents the robust application of high-order Gauss-Jacobi-type quadrature rules is overcome. In the hybrid scheme presented by the authors, local asymptotic approximation is combined with Gauss-Legendre quadrature in the transformed time variable. The corresponding spatial boundary integral operators involve only Gaussian kernels, permitting the application 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 straightforward and should allow layer heat potentials to become flexible and powerful tools for modeling diffusion processes.





Describes a project that uses

Uses Software






This page was built for publication: Hybrid asymptotic/numerical methods for the evaluation of layer heat potentials in two dimensions

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2000499)