A stochastic method for solving time-fractional differential equations

DOI10.1016/J.CAMWA.2024.02.020arXiv2303.15458MaRDI QIDQ6202618FDOQ6202618

Authors:

Publication date: 26 March 2024

Published in: Computers & Mathematics with Applications (Search for Journal in Brave)

Abstract: We present a stochastic method for efficiently computing the solution of time-fractional partial differential equations (fPDEs) that model anomalous diffusion problems of the subdiffusive type. After discretizing the fPDE in space, the ensuing system of fractional linear equations is solved resorting to a Monte Carlo evaluation of the corresponding Mittag-Leffler matrix function. This is accomplished through the approximation of the expected value of a suitable multiplicative functional of a stochastic process, which consists of a Markov chain whose sojourn times in every state are Mittag-Leffler distributed. The resulting algorithm is able to calculate the solution at conveniently chosen points in the domain with high efficiency. In addition, we present how to generalize this algorithm in order to compute the complete solution. For several large-scale numerical problems, our method showed remarkable performance in both shared-memory and distributed-memory systems, achieving nearly perfect scalability up to 16,384 CPU cores.

Full work available at URL: https://arxiv.org/abs/2303.15458

zbMATH Keywords

Mittag-Leffler functionMonte Carlo methodparallel algorithmsanomalous diffusionhigh performance computingmatrix functionstime-fractional differential equations

Mathematics Subject Classification ID

Probability theory and stochastic processes (60-XX) Numerical analysis (65-XX)

Cites Work

• Numerical evaluation of two and three parameter Mittag-Leffler functions
• Functions of Matrices
• Title not available (Why is that?)
• Title not available (Why is that?)
• Fractional Poisson process
• Applications of fractional calculus in physics
• Title not available (Why is that?)
• Title not available (Why is that?)
• Advances in Fractional Calculus
• The random walk's guide to anomalous diffusion: A fractional dynamics approach
• On Mittag-Leffler functions and related distributions
• Title not available (Why is that?)
• Accuracy and Stability of Numerical Algorithms
• Title not available (Why is that?)
• Random numbers from the tails of probability distributions using the transformation method
• Fractional Partial Differential Equations and Their Numerical Solutions
• A Schur-Parlett Algorithm for Computing Matrix Functions
• Fractional moment estimation of Linnik and Mittag-Leffler parameters
• Title not available (Why is that?)
• On the convergence of Krylov subspace methods for matrix Mittag-Leffler functions
• Exponential mixture representation of geometric stable distributions
• Computational science and engineering.
• Numerical algorithm for calculating the generalized Mittag-Leffler function
• Monte Carlo Methods for Applied Scientists
• Computation of the generalized Mittag-Leffler function and its inverse in the complex plane
• Subdiffusive discrete time random walks via Monte Carlo and subordination
• Title not available (Why is that?)
• A new \textit{walk on equations} Monte Carlo method for solving systems of linear algebraic equations
• Convergence analysis of Markov chain Monte Carlo linear solvers using Ulam-von Neumann algorithm
• Mittag-Leffler functions, related topics and applications
• Computing the matrix Mittag-Leffler function with applications to fractional calculus
• Computer simulation of geometric stable distributions
• Parallel resolvent Monte Carlo algorithms for linear algebra problems
• Analysis of Monte Carlo accelerated iterative methods for sparse linear systems.
• Elements of random walk and diffusion processes
• Fractional calculus view of complexity. Tomorrow's science
• What is the fractional Laplacian? A comparative review with new results
• Fractional diffusion on the human proteome as an alternative to the multi-organ damage of SARS-CoV-2
• Mittag-Leffler functions and their applications in network science
• A highly parallel algorithm for computing the action of a matrix exponential on a vector based on a multilevel Monte Carlo method
• A Monte Carlo method for computing the action of a matrix exponential on a vector
• Efficient generation of exponential and normal deviates

Cited In (1)

• A Fast Monte Carlo Algorithm for Evaluating Matrix Functions with Application in Complex Networks

This page was built for publication: A stochastic method for solving time-fractional differential equations

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6202618)