TY - JOUR

T1 - Stochastic simulation algorithms for solving narrow escape diffusion problems by introducing a drift to the target

AU - Sabelfeld, Karl

PY - 2020/6/1

Y1 - 2020/6/1

N2 - We suggest in this paper a new stochastic simulation algorithm for solving narrow escape problems governed by drift-diffusion-reaction equations of high dimension. The developed method drastically improves the efficiency of the diffusion trajectory tracking algorithm by introducing an artificial drift directed to the target position. The method is especially appropriate to solve narrow escape problems for domains of very long extension in one direction which is the case in many practical problems. We present simulation results for a diffusion transport problem of calculation of cathodoluminescence intensity, a diffusion flux of excitons to a threading dislocation, and the electron beam induced current in a semiconductor. The diffusion tracking algorithm is based on the random walk on spheres process. The method is meshless both in space and time, and is well applied to solve high-dimensional problems in complicated domains. The algorithms are based on tracking the trajectories of the diffusing particles exactly in accordance with the probabilistic distributions derived from the explicit representation of the relevant Green functions. They can be conveniently used not only for the solutions, but also for a direct calculation of fluxes to any part of the boundary without calculating the whole solution in the domain.

AB - We suggest in this paper a new stochastic simulation algorithm for solving narrow escape problems governed by drift-diffusion-reaction equations of high dimension. The developed method drastically improves the efficiency of the diffusion trajectory tracking algorithm by introducing an artificial drift directed to the target position. The method is especially appropriate to solve narrow escape problems for domains of very long extension in one direction which is the case in many practical problems. We present simulation results for a diffusion transport problem of calculation of cathodoluminescence intensity, a diffusion flux of excitons to a threading dislocation, and the electron beam induced current in a semiconductor. The diffusion tracking algorithm is based on the random walk on spheres process. The method is meshless both in space and time, and is well applied to solve high-dimensional problems in complicated domains. The algorithms are based on tracking the trajectories of the diffusing particles exactly in accordance with the probabilistic distributions derived from the explicit representation of the relevant Green functions. They can be conveniently used not only for the solutions, but also for a direct calculation of fluxes to any part of the boundary without calculating the whole solution in the domain.

KW - Cathodoluminescence imaging

KW - Drift-diffusion-reaction equations

KW - First passage time

KW - Narrow escape problems

KW - Random walk on spheres

KW - PROBABILITY

KW - RANDOM-WALK

KW - WINDOWS

KW - FLUX

KW - EQUATION

UR - http://www.scopus.com/inward/record.url?scp=85081978698&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2020.109406

DO - 10.1016/j.jcp.2020.109406

M3 - Article

AN - SCOPUS:85081978698

VL - 410

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

M1 - 109406

ER -