Abstract
In this paper a random walk on arbitrary rectangles (2D) and parallelepipeds (3D) algorithm is developed for solving transient anisotropic drift-diffusion-reaction equations. The method is meshless, both in space and time. The approach is based on a rigorous representation of the first passage time and exit point distributions for arbitrary rectangles and parallelepipeds. The probabilistic representation is then transformed to a form convenient for stochastic simulation. The method can be used to calculate fluxes to any desired part of the boundary, from arbitrary sources. A global version of the method we call here as a stochastic expansion from cell to cell (SECC) algorithm for calculating the whole solution field is suggested. Application of this method to solve a system of transport equations for electrons and holes in a semicoductor is discussed. This system consists of the continuity equations for particle densities and a Poisson equation for electrostatic potential. To validate the method we have derived a series of exact solutions of the drift-diffusion-reaction problem in a three-dimensional layer presented in the last section in details.
Funding source: Russian Science Foundation
Award Identifier / Grant number: 19-11-00019
Funding statement: Support of the Russian Science Foundation under Grant 19-11-00019 is gratefully acknowledged.
References
[1] M. Deaconu and A. Lejay, A random walk on rectangles algorithm, Methodol. Comput. Appl. Probab. 8 (2006), no. 1, 135–151. 10.1007/s11009-006-7292-3Search in Google Scholar
[2] L. Devroye, The series method for random variate generation and its application to the Kolmogorov–Smirnov distribution, Amer. J. Math. Management Sci. 1 (1981), no. 4, 375–379. 10.1080/01966324.1981.10737080Search in Google Scholar
[3] S. M. Ermakov and A. S. Sipin, The random walk on semi-spheres and applications for solving boundary value problems (in Russian), Vestnik St. Petersburg Univ. Math. 1 (2009), no. 3, 9–18. 10.3103/S1063454109030029Search in Google Scholar
[4] D. Grebenkov, Efficient Monte Carlo methods for simulating diffusion-reaction processes in complex systems, First-Passage Phenomena and Their Applications, World Scientific, Hackensack (2014), 571–595. 10.1142/9789814590297_0023Search in Google Scholar
[5] K. Itô and H. P. McKean, Jr., Diffusion Processes and Their Sample Paths, Grundlehren Math. Wiss. 125, Springer, Berlin, 1965. Search in Google Scholar
[6] A. Kireeva, K. K. Sabelfeld and S. Kireev, Synchronous Multi-particle cellular automaton model of diffusion with self-annihilation, Proceedings of the PACT-2019, to appear. 10.1007/978-3-030-25636-4_27Search in Google Scholar
[7] P. A. Markowich, C. A. Ringhofer and C. Schmeiser, Semiconductor Equations, Springer, Vienna, 1990. 10.1007/978-3-7091-6961-2Search in Google Scholar
[8] Y. Medvedev, Multi-particle cellular-automata models for diffusion simulation, Methods and Tools of Parallel Programming Multicomputers—MTPP 2010, Lecture Notes in Comput. Sci. 6083, Springer, Berlin (2010), 204–211. 10.1007/978-3-642-14822-4_23Search in Google Scholar
[9] S. Milewski, Combination of the meshless finite difference approach with the Monte Carlo random walk technique for solution of elliptic problems, Comput. Math. Appl. 76 (2018), no. 4, 854–876. 10.1016/j.camwa.2018.05.025Search in Google Scholar
[10] G. N. Milstein and M. V. Tretyakov, Simulation of a space-time bounded diffusion, Ann. Appl. Probab. 9 (1999), no. 3, 732–779. 10.1214/aoap/1029962812Search in Google Scholar
[11] A. Onischuk, A. Levykin, V. Strunun, M. Ushakova, R. Samoilova, K. K. Sabelfeld and V. Panfilov, Aggregate under homogeneous Silane thermal decomposition, J. Aerosol Sci. 31 (2000), no. 11, 1263–1281. 10.1016/S0021-8502(00)00031-8Search in Google Scholar
[12] A. D. Polyanin, Handbook of Linear Partial Differential Equations for Engineers and Scientists, Chapman & Hall/CRC, Boca Raton, 2002. 10.1201/9781420035322Search in Google Scholar
[13] A. P. Prudnikov, Y. A. Brychkov and O. I. Marichev, Integrals and Series. Vol. 1: Elementary Functions, “Nauka”, Moscow, 1981. Search in Google Scholar
[14] K. K. Sabelfeld, Random walk on spheres method for solving drift-diffusion problems, Monte Carlo Methods Appl. 22 (2016), no. 4, 265–275. 10.1515/mcma-2016-0118Search in Google Scholar
[15] K. K. Sabelfeld, A mesh free floating random walk method for solving diffusion imaging problems, Statist. Probab. Lett. 121 (2017), 6–11. 10.1016/j.spl.2016.10.006Search in Google Scholar
[16] K. K. Sabelfeld, Random walk on spheres algorithm for solving transient drift-diffusion-reaction problems, Monte Carlo Methods Appl. 23 (2017), no. 3, 189–212. 10.1515/mcma-2017-0113Search in Google Scholar
[17] K. K. Sabelfeld, Meshfree stochastic algorithms for systems of diffusion-advection-reaction equations and anisotropic diffusion flux calculations, Prob. Eng. Mech., to appear. Search in Google Scholar
[18] K. K. Sabelfeld and A. E. Kireeva, A meshless random walk on parallelepipeds algorithm for solving transient anisotropic diffusion-recombination equations and applications to cathodoluminescence imaging, Numer. Math., to appear. Search in Google Scholar
[19] C. Vamoş, N. Suciu and H. Vereecken, Generalized random walk algorithm for the numerical modeling of complex diffusion processes, J. Comput. Phys. 186 (2003), no. 2, 527–544. 10.1016/S0021-9991(03)00073-1Search in Google Scholar
[20] T. Vesala, U. Rannik, M. Leclerc, T. Foken and K. K. Sabelfeld, Flux and concentration footprints, Agricultural Forest Meteorol. 127 (2004), no. 3–4, 111–116. 10.1016/j.agrformet.2004.07.007Search in Google Scholar
© 2019 Walter de Gruyter GmbH, Berlin/Boston