{"title":"Random walk on spheres method for solving anisotropic transient diffusion problems and flux calculations","authors":"Irina Shalimova, Karl Sabelfeld","doi":"10.1515/mcma-2023-2022","DOIUrl":null,"url":null,"abstract":"Abstract The Random Walk on Spheres (RWS) algorithm for solving anisotropic transient diffusion problems based on a stochastic learning procedure for calculation of the exit position of the anisotropic diffusion process on a sphere is developed. Direct generalization of the Random Walk on Spheres method to anisotropic diffusion equations is not possible, therefore, we have numerically calculated the probability density for the exit position on a sphere. The first passage time is then represented explicitly. The method can easily be implemented to solve diffusion problems with spatially varying diffusion coefficients for complicated three-dimensional domains. Particle tracking algorithm is highly efficient for calculation of fluxes to boundaries. We apply the developed algorithm for solving an exciton transport in a semiconductor material with a threading dislocation where the measured functions are the exciton fluxes to the semiconductor’s substrate and on the dislocation surface.","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":"7 3","pages":"0"},"PeriodicalIF":0.8000,"publicationDate":"2023-11-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Monte Carlo Methods and Applications","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/10.1515/mcma-2023-2022","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q3","JCRName":"STATISTICS & PROBABILITY","Score":null,"Total":0}
引用次数: 0
Abstract
Abstract The Random Walk on Spheres (RWS) algorithm for solving anisotropic transient diffusion problems based on a stochastic learning procedure for calculation of the exit position of the anisotropic diffusion process on a sphere is developed. Direct generalization of the Random Walk on Spheres method to anisotropic diffusion equations is not possible, therefore, we have numerically calculated the probability density for the exit position on a sphere. The first passage time is then represented explicitly. The method can easily be implemented to solve diffusion problems with spatially varying diffusion coefficients for complicated three-dimensional domains. Particle tracking algorithm is highly efficient for calculation of fluxes to boundaries. We apply the developed algorithm for solving an exciton transport in a semiconductor material with a threading dislocation where the measured functions are the exciton fluxes to the semiconductor’s substrate and on the dislocation surface.