{"title":"包含分数阶拉普拉斯式的抛物方程的蒙特卡罗方法","authors":"Caiyu Jiao, Changpin Li","doi":"10.1515/mcma-2022-2129","DOIUrl":null,"url":null,"abstract":"Abstract We apply the Monte Carlo method to solving the Dirichlet problem of linear parabolic equations with fractional Laplacian. This method exploits the idea of weak approximation of related stochastic differential equations driven by the symmetric stable Lévy process with jumps. We utilize the jump-adapted scheme to approximate Lévy process which gives exact exit time to the boundary. When the solution has low regularity, we establish a numerical scheme by removing the small jumps of the Lévy process and then show the convergence order. When the solution has higher regularity, we build up a higher-order numerical scheme by replacing small jumps with a simple process and then display the higher convergence order. Finally, numerical experiments including ten- and one hundred-dimensional cases are presented, which confirm the theoretical estimates and show the numerical efficiency of the proposed schemes for high-dimensional parabolic equations.","PeriodicalId":46576,"journal":{"name":"Monte Carlo Methods and Applications","volume":"29 1","pages":"33 - 53"},"PeriodicalIF":0.8000,"publicationDate":"2022-10-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"Monte Carlo method for parabolic equations involving fractional Laplacian\",\"authors\":\"Caiyu Jiao, Changpin Li\",\"doi\":\"10.1515/mcma-2022-2129\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"Abstract We apply the Monte Carlo method to solving the Dirichlet problem of linear parabolic equations with fractional Laplacian. This method exploits the idea of weak approximation of related stochastic differential equations driven by the symmetric stable Lévy process with jumps. We utilize the jump-adapted scheme to approximate Lévy process which gives exact exit time to the boundary. When the solution has low regularity, we establish a numerical scheme by removing the small jumps of the Lévy process and then show the convergence order. When the solution has higher regularity, we build up a higher-order numerical scheme by replacing small jumps with a simple process and then display the higher convergence order. Finally, numerical experiments including ten- and one hundred-dimensional cases are presented, which confirm the theoretical estimates and show the numerical efficiency of the proposed schemes for high-dimensional parabolic equations.\",\"PeriodicalId\":46576,\"journal\":{\"name\":\"Monte Carlo Methods and Applications\",\"volume\":\"29 1\",\"pages\":\"33 - 53\"},\"PeriodicalIF\":0.8000,\"publicationDate\":\"2022-10-27\",\"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-2022-2129\",\"RegionNum\":0,\"RegionCategory\":null,\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"Q3\",\"JCRName\":\"STATISTICS & PROBABILITY\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"Monte Carlo Methods and Applications","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/10.1515/mcma-2022-2129","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q3","JCRName":"STATISTICS & PROBABILITY","Score":null,"Total":0}
Monte Carlo method for parabolic equations involving fractional Laplacian
Abstract We apply the Monte Carlo method to solving the Dirichlet problem of linear parabolic equations with fractional Laplacian. This method exploits the idea of weak approximation of related stochastic differential equations driven by the symmetric stable Lévy process with jumps. We utilize the jump-adapted scheme to approximate Lévy process which gives exact exit time to the boundary. When the solution has low regularity, we establish a numerical scheme by removing the small jumps of the Lévy process and then show the convergence order. When the solution has higher regularity, we build up a higher-order numerical scheme by replacing small jumps with a simple process and then display the higher convergence order. Finally, numerical experiments including ten- and one hundred-dimensional cases are presented, which confirm the theoretical estimates and show the numerical efficiency of the proposed schemes for high-dimensional parabolic equations.