INITIAL-BOUNDARY VALUE PROBLEM FOR THE HEAT EQUATION-A STOCHASTIC ALGORITHM

• 成果类型: Article
• 署名作者: Deaconu, Madalina; Herrmann, Samuel
• 署名单位: Inria; Universite de Lorraine; Centre National de la Recherche Scientifique (CNRS); CNRS - National Institute for Mathematical Sciences (INSMI); Universite Bourgogne Europe; Centre National de la Recherche Scientifique (CNRS); CNRS - National Institute for Mathematical Sciences (INSMI)
• 刊物名称: ANNALS OF APPLIED PROBABILITY
• ISSN/ISSBN: 1050-5164
• DOI: 10.1214/17-AAP1348
• 发表日期: 2018
• 页码: 1943-1976
• 关键词: bessel processes; moving spheres; walk; PROPERTY; time

摘要

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.