Communications on Applied Mathematics and Computation ›› 2021, Vol. 3 ›› Issue (1): 61-90.doi: 10.1007/s42967-020-00070-w

• ORIGINAL PAPER • 上一篇    

A Unified Petrov–Galerkin Spectral Method and Fast Solver for Distributed-Order Partial Differential Equations

Mehdi Samiee1,2, Ehsan Kharazmi1,2, Mark M. Meerschaert3, Mohsen Zayernouri1,3   

  1. 1 Department of Mechanical Engineering, Michigan State University, East Lansing, MI 48824, USA;
    2 Department of Computational Mathematics, Science and Engineering, Michigan State University, East Lansing, MI 48824, USA;
    3 Department of Probability and Statistics, Michigan State University, East Lansing, MI 48824, USA
  • 收稿日期:2019-07-25 修回日期:2020-03-26 发布日期:2021-03-15
  • 通讯作者: Mohsen Zayernouri, zayern@msu.edu E-mail:zayern@msu.edu
  • 基金资助:
    This work was supported by the AFOSR Young Investigator Program (YIP) award (FA9550-17-1-0150), the MURI/ARO (W911NF-15-1-0562), the National Science Foundation Award (DMS-1923201), and the ARO Young Investigator Program Award (W911NF-19-1-0444)

A Unified Petrov–Galerkin Spectral Method and Fast Solver for Distributed-Order Partial Differential Equations

Mehdi Samiee1,2, Ehsan Kharazmi1,2, Mark M. Meerschaert3, Mohsen Zayernouri1,3   

  1. 1 Department of Mechanical Engineering, Michigan State University, East Lansing, MI 48824, USA;
    2 Department of Computational Mathematics, Science and Engineering, Michigan State University, East Lansing, MI 48824, USA;
    3 Department of Probability and Statistics, Michigan State University, East Lansing, MI 48824, USA
  • Received:2019-07-25 Revised:2020-03-26 Published:2021-03-15
  • Contact: Mohsen Zayernouri, zayern@msu.edu E-mail:zayern@msu.edu
  • Supported by:
    This work was supported by the AFOSR Young Investigator Program (YIP) award (FA9550-17-1-0150), the MURI/ARO (W911NF-15-1-0562), the National Science Foundation Award (DMS-1923201), and the ARO Young Investigator Program Award (W911NF-19-1-0444)

摘要: Fractional calculus and fractional-order modeling provide effective tools for modeling and simulation of anomalous diffusion with power-law scalings. In complex multi-fractal anomalous transport phenomena, distributed-order partial differential equations appear as tractable mathematical models, where the underlying derivative orders are distributed over a range of values, hence taking into account a wide range of multi-physics from ultraslow-to-standard-to-superdiffusion/wave dynamics. We develop a unified, fast, and stable Petrov–Galerkin spectral method for such models by employing Jacobi poly-fractonomials and Legendre polynomials as temporal and spatial basis/test functions, respectively. By defining the proper underlying distributed Sobolev spaces and their equivalent norms, we rigorously prove the well-posedness of the weak formulation, and thereby, we carry out the corresponding stability and error analysis. We finally provide several numerical simulations to study the performance and convergence of proposed scheme.

关键词: Distributed Sobolev space, Well-posedness analysis, Discrete inf-sup condition, Spectral convergence, Jacobi poly-fractonomials, Legendre polynomials

Abstract: Fractional calculus and fractional-order modeling provide effective tools for modeling and simulation of anomalous diffusion with power-law scalings. In complex multi-fractal anomalous transport phenomena, distributed-order partial differential equations appear as tractable mathematical models, where the underlying derivative orders are distributed over a range of values, hence taking into account a wide range of multi-physics from ultraslow-to-standard-to-superdiffusion/wave dynamics. We develop a unified, fast, and stable Petrov–Galerkin spectral method for such models by employing Jacobi poly-fractonomials and Legendre polynomials as temporal and spatial basis/test functions, respectively. By defining the proper underlying distributed Sobolev spaces and their equivalent norms, we rigorously prove the well-posedness of the weak formulation, and thereby, we carry out the corresponding stability and error analysis. We finally provide several numerical simulations to study the performance and convergence of proposed scheme.

Key words: Distributed Sobolev space, Well-posedness analysis, Discrete inf-sup condition, Spectral convergence, Jacobi poly-fractonomials, Legendre polynomials

中图分类号: