Nyström method for BEM of the heat equation with moving boundaries

Abstract

A direct boundary integral equation method for the heat equation based on Nystrom discretization is proposed and analyzed. For problems with moving geometries, a weakly and strongly singular Green’s integral equation is formulated. Here the hypersingular integral operator, i.e., the normal trace of the double-layer potential, must be understood as a Hadamard finite part integral. The thermal layer potentials are regarded as generalized Abel integral operators in time and discretized with a singularity-corrected trapezoidal rule. The spatial discretization is a standard quadrature rule for smooth surface integrals. The discretized systems lead to an explicit time stepping scheme and is effective for solving the Dirichlet and Neumann boundary value problems based on both the weakly and/or strongly singular integral equations.