Initial-boundary value problem for the heat equation.A stochastic algorithm

The initial-boundary value problem for the heat equation is solved by using an algorithm based on a random walk on heat balls. Even if it represents a sophisticated generalization of the Walk on Spheres (WOS) algorithm introduced to solve the Dirichlet problem for Laplace.s equation, its implementation is rather easy. The construction of this algorithm can be considered as a natural consequence of previous works the authors completed on the hitting time approximation for Bessel processes and Brownian motion [Ann. Appl. Probab. 23 (2013) 2259.2289, Math. Comput. Simulation 135 (2017) 28.38, Bernoulli 23 (2017) 3744.3771]. A similar procedure was introduced previously in the paper [Random Processes for Classical Equations of Mathematical Physics (1989) Kluwer Academic]. The definition of the random walk is based on a particular mean value formula for the heat equation. We present here a probabilistic view of this formula. The aim of the paper is to prove convergence results for this algorithm and to illustrate them by numerical examples. These examples permit to emphasize the efficiency and accuracy of the algorithm.