학술논문

Divide-and-conquer solver in tensor-train format for $d$-dimensional time-space fractional diffusion equations.
Document Type
Journal
Author
Huang, Yun-Chi (PRC-SHTO) AMS Author Profile; Chou, Lot-Kei (PRC-UMAC-M) AMS Author Profile; Lei, Siu-Long (PRC-UMAC-M) AMS Author Profile
Source
Journal of Scientific Computing (J. Sci. Comput.) (20230101), 96, no.~1, Paper No 29, 35~pp. ISSN: 0885-7474 (print).eISSN: 1573-7691.
Subject
26 Real functions -- 26A Functions of one variable
  26A33 Fractional derivatives and integrals

41 Approximations and expansions
  41-04 Explicit machine computation and programs

41 Approximations and expansions -- 41A Approximations and expansions
  41A63 Multidimensional problems

65 Numerical analysis -- 65G Error analysis and interval analysis
  65G50 Roundoff error

65 Numerical analysis -- 65M Partial differential equations, initial value and time-dependent initial-boundary value problems
  65M12 Stability and convergence of numerical methods
Language
English
Abstract
In this work, the authors consider the $d$-dimensional time-space fractional partial differential equations, with $d\ge 3$, where the Caputo fractional derivative of order $\gamma$ is considered in a time variable, with $0 < \gamma < 1$, and the $\alpha_k$-order left and right Riemann-Liouville fractional derivatives are considered in the spatial variables, with ${1 < \alpha_k < 2}$. \par To solve the continuous problem, an alternating direction implicit method (ADI) is used to discretize in time on a uniform mesh. The Caputo derivative is approximated by a $L_1$ approximation, and the linear combination of Riemann-Liouville fractional derivatives are approximated by shifted Grünwald approximations. These types of approximations, for the fractional derivatives, in practice require a huge computational cost and storage. \par To reduce those difficulties, a tensor-train format is used, which is a good choice to obtain a faster solver in the context of the numerical resolution of fractional partial differential equations. In order to apply this technique, the ADI scheme is reformulated in an appropriate way, and therefore the divide-and-conquer solver is better defined. \par In my opinion, the main interest of the article is that an adequate method is used to considerably reduce the computational cost and the storage need of the numerical method used to solve an important type of partial differential equation. These problems are more difficult than classical problems, because fractional derivatives appear in both the time and the spatial variables. The new algorithm permits one to obtain good and efficient solutions for the type of problems considered, which are good models for many problems that appear in different areas related to applied mathematics.