#### Vol. 1, No. 4, 2019

 Recent Issues Volume 4, Issue 1 Volume 3, Issue 4 Volume 3, Issue 3 Volume 3, Issue 2 Volume 3, Issue 1 Volume 2, Issue 4 Volume 2, Issue 4 Volume 2, Issue 3 Volume 2, Issue 2 Volume 2, Issue 1 Volume 1, Issue 4 Volume 1, Issue 3 Volume 1, Issue 2 Volume 1, Issue 1
 The Journal About the Journal Editorial Board Subscriptions Submission Guidelines Submission Form Policies for Authors Ethics Statement ISSN (electronic): 2578-5885 ISSN (print): 2578-5893 Author Index To Appear Other MSP Journals
Explicit unconditionally stable methods for the heat equation via potential theory

### Alex Barnett, Charles L. Epstein, Leslie Greengard, Shidong Jiang and Jun Wang

Vol. 1 (2019), No. 4, 709–742
DOI: 10.2140/paa.2019.1.709
##### Abstract

We study the stability properties of explicit marching schemes for second-kind Volterra integral equations that arise when solving boundary value problems for the heat equation by means of potential theory. It is well known that explicit finite-difference or finite-element schemes for the heat equation are stable only if the time step $\Delta t$ is of the order $\mathsc{O}\left(\Delta {x}^{2}\right)$, where $\Delta x$ is the finest spatial grid spacing. In contrast, for the Dirichlet and Neumann problems on the unit ball in all dimensions $d\ge 1$, we show that the simplest Volterra marching scheme, i.e., the forward Euler scheme, is unconditionally stable. Our proof is based on an explicit spectral radius bound of the marching matrix, leading to an estimate that an ${L}^{2}$-norm of the solution to the integral equation is bounded by ${c}_{d}{T}^{d∕2}$ times the norm of the right-hand side. For the Robin problem on the half-space in any dimension, with constant Robin (heat transfer) coefficient $\kappa$, we exhibit a constant $C$ such that the forward Euler scheme is stable if $\Delta t, independent of any spatial discretization. This relies on new lower bounds on the spectrum of real symmetric Toeplitz matrices defined by convex sequences. Finally, we show that the forward Euler scheme is unconditionally stable for the Dirichlet problem on any smooth convex domain in any dimension, in the ${L}^{\infty }$-norm.

##### Keywords
heat equation, Abel equation, forward Euler scheme, Volterra integral equation, stability analysis, Toeplitz matrix, convex sequence, modified Bessel function of the first kind
Primary: 00A05