Monte Carlo method for parabolic equations involving fractional Laplacian
From MaRDI portal
Publication:6415271
DOI10.1515/MCMA-2022-2129arXiv2210.15192MaRDI QIDQ6415271
Publication date: 27 October 2022
Abstract: We apply the Monte Carlo method to solving the Dirichlet problem of linear parabolic equations with fractional Laplacian. This method exploit- s the idea of weak approximation of related stochastic differential equations driven by the symmetric stable L'evy process with jumps. We utilize the jump- adapted scheme to approximate L'evy process which gives exact exit time to the boundary. When the solution has low regularity, we establish a numeri- cal scheme by removing the small jumps of the L'evy 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.
Monte Carlo methods (65C05) Computational methods for stochastic equations (aspects of stochastic analysis) (60H35) Numerical solutions to stochastic differential and integral equations (65C30) Fractional partial differential equations (35R11)
This page was built for publication: Monte Carlo method for parabolic equations involving fractional Laplacian
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6415271)