Spectral, Tensor and Domain Decomposition Methods for Fractional PDEs

Tianyi Shi; Harbir Antil; Drew P. Kouri

doi:10.1515/cmam-2021-0118

Artikel

Spectral, Tensor and Domain Decomposition Methods for Fractional PDEs

Tianyi Shi , Harbir Antil und Drew P. Kouri

Veröffentlicht/Copyright: 12. Februar 2022

Veröffentlicht von

Veröffentlichen auch Sie bei De Gruyter Brill

Manuskript einreichen Informationen für Autor*innen Erkunden Sie dieses Fachgebiet

Aus der Zeitschrift Computational Methods in Applied Mathematics Band 22 Heft 4

Abstract

Fractional PDEs have recently found several geophysics and imaging science applications due to their nonlocal nature and their flexibility in capturing sharp transitions across interfaces. However, this nonlocality makes it challenging to design efficient solvers for such problems. In this paper, we introduce a spectral method based on an ultraspherical polynomial discretization of the Caffarelli–Silvestre extension to solve such PDEs on rectangular and disk domains. We solve the discretized problem using tensor equation solvers and thus can solve higher-dimensional PDEs. In addition, we introduce both serial and parallel domain decomposition solvers. We demonstrate the numerical performance of our methods on a 3D fractional elliptic PDE on a cube as well as an application to optimization problems with fractional PDE constraints.

Keywords: Ultraspherical polynomial; Fractional PDEs; Tensor Methods; Domain Decomposition; 3D Problems; Optimal Control

MSC 2010: 49J20; 49K20; 35S15; 65R20; 65N30

Dedicated to our Friend Francisco-Javier Sayas

Funding source: Air Force Office of Scientific Research

Award Identifier / Grant number: FA9550-19-1-0036

Award Identifier / Grant number: F4FGA09135G001

Funding source: National Science Foundation

Award Identifier / Grant number: DMS-2110263

Award Identifier / Grant number: DMS-1913004

Funding source: National Nuclear Security Administration

Award Identifier / Grant number: DE-NA-0003525

Funding statement: This work is partially supported by the Air Force Office of Scientific Research under award numbers FA9550-19-1-0036 and F4FGA09135G001 and NSF grants DMS-2110263 and DMS-1913004. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology and Engineering Solutions of Sandia, LLC., a wholly owned subsidiary of Honeywell International, Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under grant DE-NA-0003525. This paper describes objective technical results and analysis. Any subjective views or opinions that might be expressed in the paper do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

A Visualization of Matrices in Ultraspherical Discretization

In this appendix, we provide the readers with some details of the matrices we use during discretization in Sections 3 and 4.

𝐷 is a diagonal matrix representing second derivative of ( 1 - ρ 2 ) ⁢ C ~ ( 3 / 2 ) ⁢ ( ρ ) , with diagonal elements D j , j = - ( j ⁢ ( j + 3 ) + 2 ) .
𝑀 is a symmetric penta-diagonal matrix with 0 super- and subdiagonals representing multiplication of 1 - ρ 2 in C ~ ( 3 / 2 ) basis, with elements
M j , j = 2 ⁢ ( j + 1 ) ⁢ ( j + 2 ) ( 2 ⁢ j + 1 ) ⁢ ( 2 ⁢ j + 5 ) , M j , j + 1 = 0 , M j , j + 2 = - 1 ( 2 ⁢ j + 3 ) ⁢ ( 2 ⁢ j + 5 ) ⁢ ( j + 4 ) ! ⁢ ( 2 ⁢ j + 3 ) j ! ⁢ ( 2 ⁢ j + 7 ) .
B 1 and B 2 are two banded matrices representing multiplications in different ultraspherical polynomial basis, where the bandwidth is determined by the degree of the polynomial approximation in the respective basis. In addition, B 1 can be shown to be a Toeplitz-plus-Hankel-plus-rank-1 operator [26], and both B 1 and B 2 satisfy three-term recurrence [30, Chapter 6].
D 2 represents second derivative of Chebyshev basis, with the form
D 2 = [ 0 0 2 3 4 ⋱ ] .
Here, elements of D 2 on the diagonal and first upper top diagonal are all 0.
S 2 , 0 converts Chebyshev basis to C ( 2 ) basis, and can be calculated by S 2 , 0 = S 1 ⁢ S 0 , where
S 0 = [ 1 0 - 1 / 2 1 / 2 0 - 1 / 2 1 / 2 0 ⋱ 1 / 2 ⋱ ⋱ ] , S 1 = [ 1 0 - 1 / 3 1 / 2 0 - 1 / 4 1 / 3 0 ⋱ 1 / 4 ⋱ ⋱ ] .

References

[1] M. Ainsworth and C. Glusa, Hybrid finite element–spectral method for the fractional Laplacian: Approximation theory and efficient solver, SIAM J. Sci. Comput. 40 (2018), no. 4, A2383–A2405. 10.1137/17M1144696Suche in Google Scholar

[2] M. Ainsworth and Z. Mao, Analysis and approximation of a fractional Cahn–Hilliard equation, SIAM J. Numer. Anal. 55 (2017), no. 4, 1689–1718. 10.1137/16M1075302Suche in Google Scholar

[3] H. Antil and S. Bartels, Spectral approximation of fractional PDEs in image processing and phase field modeling, Comput. Methods Appl. Math. 17 (2017), no. 4, 661–678. 10.1515/cmam-2017-0039Suche in Google Scholar

[4] H. Antil, Z. W. Di and R. Khatri, Bilevel optimization, deep learning and fractional Laplacian regularization with applications in tomography, Inverse Problems 36 (2020), no. 6, Article ID 064001. 10.1088/1361-6420/ab80d7Suche in Google Scholar

[5] H. Antil, J. Pfefferer and S. Rogovs, Fractional operators with inhomogeneous boundary conditions: Analysis, control, and discretization, Commun. Math. Sci. 16 (2018), no. 5, 1395–1426. 10.4310/CMS.2018.v16.n5.a11Suche in Google Scholar

[6] I. Babuška and B. Guo, The h, p and hp version of the finite element method; basis theory and applications, Adv. Eng. Softw. 15 (1992), no. 3–4, 159–174. 10.1016/0965-9978(92)90097-YSuche in Google Scholar

[7] L. Banjai, J. M. Melenk, R. H. Nochetto, E. Otárola, A. J. Salgado and C. Schwab, Tensor FEM for spectral fractional diffusion, Found. Comput. Math. 19 (2019), no. 4, 901–962. 10.1007/s10208-018-9402-3Suche in Google Scholar

[8] R. H. Bartels and G. W. Stewart, Solution of the matrix equation AX+ XB= C, Commun. ACM 15 (1972), no. 9, 820–826. 10.1145/361573.361582Suche in Google Scholar

[9] P. Benner, R.-C. Li and N. Truhar, On the ADI method for Sylvester equations, J. Comput. Appl. Math. 233 (2009), no. 4, 1035–1045. 10.1016/j.cam.2009.08.108Suche in Google Scholar

[10] A. Bonito and J. E. Pasciak, Numerical approximation of fractional powers of elliptic operators, Math. Comp. 84 (2015), no. 295, 2083–2110. 10.1090/S0025-5718-2015-02937-8Suche in Google Scholar

[11] L. Caffarelli and L. Silvestre, An extension problem related to the fractional Laplacian, Comm. Partial Differential Equations 32 (2007), no. 7–9, 1245–1260. 10.1080/03605300600987306Suche in Google Scholar

[12] S. Chen and J. Shen, An efficient and accurate numerical method for the spectral fractional Laplacian equation, J. Sci. Comput. 82 (2020), no. 1, Paper No. 17. 10.1007/s10915-019-01122-xSuche in Google Scholar

[13] P. Constantin and J. Wu, Behavior of solutions of 2D quasi-geostrophic equations, SIAM J. Math. Anal. 30 (1999), no. 5, 937–948. 10.1137/S0036141098337333Suche in Google Scholar

[14] A. de Pablo, F. Quirós, A. Rodríguez and J. L. Vázquez, A fractional porous medium equation, Adv. Math. 226 (2011), no. 2, 1378–1409. 10.1016/j.aim.2010.07.017Suche in Google Scholar

[15] E. Di Nezza, G. Palatucci and E. Valdinoci, Hitchhiker’s guide to the fractional Sobolev spaces, Bull. Sci. Math. 136 (2012), no. 5, 521–573. 10.1016/j.bulsci.2011.12.004Suche in Google Scholar

[16] D. Fortunato, N. Hale and A. Townsend, The ultraspherical spectral element method, J. Comput. Phys. 436 (2021), Paper No. 110087. 10.1016/j.jcp.2020.110087Suche in Google Scholar

[17] D. Fortunato and A. Townsend, Fast Poisson solvers for spectral methods, IMA J. Numer. Anal. 40 (2020), no. 3, 1994–2018. 10.1093/imanum/drz034Suche in Google Scholar

[18] N. Hale and L. N. Trefethen, Chebfun and numerical quadrature, Sci. China Math. 55 (2012), no. 9, 1749–1760. 10.1007/s11425-012-4474-zSuche in Google Scholar

[19] T. G. Kolda and B. W. Bader, Tensor decompositions and applications, SIAM Rev. 51 (2009), no. 3, 455–500. 10.1137/07070111XSuche in Google Scholar

[20] N. Laskin, Fractional quantum mechanics and Lévy path integrals, Phys. Lett. A 268 (2000), no. 4–6, 298–305. 10.1016/S0375-9601(00)00201-2Suche in Google Scholar

[21] A. Lischke, G. Pang, M. Gulian and et al., What is the fractional Laplacian? A comparative review with new results, J. Comput. Phys. 404 (2020), Article ID 109009. 10.1016/j.jcp.2019.109009Suche in Google Scholar

[22] P. Martinsson, The hierarchical Poincaré–Steklov (HPS) solver for elliptic PDEs: A tutorial, preprint (2015), https://arxiv.org/abs/1506.01308. Suche in Google Scholar

[23] D. Meidner, J. Pfefferer, K. Schürholz and B. Vexler, h ⁢ p -finite elements for fractional diffusion, SIAM J. Numer. Anal. 56 (2018), no. 4, 2345–2374. 10.1137/17M1135517Suche in Google Scholar

[24] R. H. Nochetto, E. Otárola and A. J. Salgado, A PDE approach to fractional diffusion in general domains: A priori error analysis, Found. Comput. Math. 15 (2015), no. 3, 733–791. 10.1007/s10208-014-9208-xSuche in Google Scholar

[25] F. W. J. Olver, D. W. Lozier, R. F. Boisvert and C. W. Clark, NIST Handbook of Mathematical Functions, Cambridge University, Cambridge, 2010. Suche in Google Scholar

[26] S. Olver and A. Townsend, A fast and well-conditioned spectral method, SIAM Rev. 55 (2013), no. 3, 462–489. 10.1137/120865458Suche in Google Scholar

[27] B. Schmitt, B. N. Khoromskij, V. Khoromskaia and V. Schulz, Tensor method for optimal control problems constrained by fractional three-dimensional elliptic operator with variable coefficients, Numer. Linear Algebra Appl. 2021 (2021), Article ID e2404. 10.1002/nla.2404Suche in Google Scholar

[28] F. Song, C. Xu and G. E. Karniadakis, Computing fractional Laplacians on complex-geometry domains: Algorithms and simulations, SIAM J. Sci. Comput. 39 (2017), no. 4, A1320–A1344. 10.1137/16M1078197Suche in Google Scholar

[29] P. R. Stinga and J. L. Torrea, Extension problem and Harnack’s inequality for some fractional operators, Comm. Partial Differential Equations 35 (2010), no. 11, 2092–2122. 10.1080/03605301003735680Suche in Google Scholar

[30] A. Townsend, Computing with functions in two dimensions, PhD thesis, University of Oxford, 2014. Suche in Google Scholar

[31] A. Townsend and S. Olver, The automatic solution of partial differential equations using a global spectral method, J. Comput. Phys. 299 (2015), 106–123. 10.1016/j.jcp.2015.06.031Suche in Google Scholar

[32] C. J. Weiss, B. G. van Bloemen Waanders and H. Antil, Fractional operators applied to geophysical electromagnetics, Geophys. J. Int. 220 (2020), no. 2, 1242–1259. 10.1093/gji/ggz516Suche in Google Scholar

[33] H. Wilber, A. Townsend and G. B. Wright, Computing with functions in spherical and polar geometries II. The disk, SIAM J. Sci. Comput. 39 (2017), no. 3, C238–C262. 10.1137/16M1070207Suche in Google Scholar

[34] M. Zayernouri, M. Ainsworth and G. E. Karniadakis, A unified Petrov–Galerkin spectral method for fractional PDEs, Comput. Methods Appl. Mech. Engrg. 283 (2015), 1545–1569. 10.1016/j.cma.2014.10.051Suche in Google Scholar

Received: 2021-06-24

Revised: 2021-12-02

Accepted: 2022-01-04

Published Online: 2022-02-12

Published in Print: 2022-10-01

Sie haben derzeit keinen Zugang zu diesem Inhalt.

Artikel in diesem Heft

https://doi.org/10.1515/cmam-2021-0118

Schlagwörter für diesen Artikel

Ultraspherical polynomial; Fractional PDEs; Tensor Methods; Domain Decomposition; 3D Problems; Optimal Control