A Feynman-Kac based numerical method for the exit time probability of a class of transport problems

JOURNAL OF COMPUTATIONAL PHYSICS

卷:444

Article

Yang, Minglei¹  Zhang, Guannan²  del-Castillo-Negrete, Diego¹  Stoyanov, Miroslav² 

[1] Oak Ridge Natl Lab, Fus Energy Div, Oak Ridge, TN USA

[2] Oak Ridge Natl Lab, Comp Sci & Math Div, Oak Ridge, TN 37830 USA

关键词: Feynman-Kac formula;    Fokker-Planck equation;    Stochastic differential equations;    First exit time;    Adjoint equations;    Transport;   

DOI  :  10.1016/j.jcp.2021.110564

来源: Elsevier

【 摘 要 】

The exit time probability, which gives the likelihood that an initial condition leaves a prescribed region of the phase space of a dynamical system at, or before, a given time, is arguably one of the most natural and important transport problems. Here we present an accurate and efficient numerical method for computing this probability for systems described by non-autonomous (time-dependent) stochastic differential equations (SDEs) or their equivalent Fokker-Planck partial differential equations. The method is based on the direct approximation of the Feynman-Kac formula that establishes a link between the adjoint Fokker-Planck equation and the forward SDE. The Feynman-Kac formula is approximated using the Gauss-Hermite quadrature rules and piecewise cubic Hermite interpolating polynomials, and a GPU accelerated matrix representation is used to compute the entire time evolution of the exit time probability using a single pass of the algorithm. The method is unconditionally stable, exhibits second order convergence in space, first order convergence in time, and it is straightforward to parallelize. Applications are presented to the advection diffusion of a passive tracer in a fluid flow exhibiting chaotic advection, and to the runaway acceleration of electrons in a plasma in the presence of an electric field, collisions, and radiation damping. Benchmarks against analytical solutions as well as comparisons with explicit and implicit finite difference standard methods for the adjoint Fokker-Planck equation are presented. (C) 2021 Elsevier Inc. All rights reserved.

【 授权许可】

Free