Abstract
We propose a new nonconforming finite element algorithm to approximate the solution to the elliptic problem involving the fractional Laplacian. We first derive an integral representation of the bilinear form corresponding to the variational problem. The numerical approximation of the action of the corresponding stiffness matrix consists of three steps: (1) apply a sinc quadrature scheme to approximate the integral representation by a finite sum where each term involves the solution of an elliptic partial differential equation defined on the entire space, (2) truncate each elliptic problem to a bounded domain, (3) use the finite element method for the space approximation on each truncated domain. The consistency error analysis for the three steps is discussed together with the numerical implementation of the entire algorithm. The results of computations are given illustrating the error behavior in terms of the mesh size of the physical domain, the domain truncation parameter and the quadrature spacing parameter.
Similar content being viewed by others
References
Acosta, G., Borthagaray, J.P.: A fractional Laplace equation: regularity of solutions and finite element approximations. SIAM J. Numer. Anal. 55(2), 472–495 (2017). https://doi.org/10.1137/15M1033952
Auscher, P., Hofmann, S., Lewis, J.L., Tchamitchian, P.: Extrapolation of Carleson measures and the analyticity of Kato’s square-root operators. Acta Math. 187(2), 161–190 (2001). https://doi.org/10.1007/BF02392615
Bacuta, C.: Interpolation Between Subspaces of Hilbert Spaces and Applications to Shift Theorems for Elliptic Boundary Value Problems and Finite Element Methods. ProQuest LLC, Ann Arbor, MI (2000). http://gateway.proquest.com/openurl?url_ver=Z39.88-2004&rft_val_fmt=info:ofi/fmt:kev:mtx:dissertation&res_dat=xri:pqdiss&rft_dat=xri:pqdiss:9994203. Thesis (Ph.D.)—Texas A&M University
Bangerth, W., Hartmann, R., Kanschat, G.: deal.II—a general-purpose object-oriented finite element library. ACM Trans. Math. Softw. 33(4), 24 (2007). https://doi.org/10.1145/1268776.1268779
Biler, P., Karch, G., Woyczyński, W.A.: Critical nonlinearity exponent and self-similar asymptotics for Lévy conservation laws. Ann. Inst. H. Poincaré Anal. Non Linéaire 18(5), 613–637 (2001). https://doi.org/10.1016/S0294-1449(01)00080-4
Bonito, A., Lei, W., Pasciak, J.E.: On sinc quadrature approximations of fractional powers of regularly accretive operators. J. Numer. Math. (2018). https://doi.org/10.1515/jnma-2017-0116
Bonito, A., Lei, W., Pasciak, J.E.: The approximation of parabolic equations involving fractional powers of elliptic operators. J. Comput. Appl. Math. 315, 32–48 (2017). https://doi.org/10.1016/j.cam.2016.10.016
Bonito, A., Pasciak, J.E.: Numerical approximation of fractional powers of elliptic operators. Math. Comput. 84(295), 2083–2110 (2015). https://doi.org/10.1090/S0025-5718-2015-02937-8
Bonito, A., Pasciak, J.E.: Numerical approximation of fractional powers of regularly accretive operators. IMA J. Numer. Anal. 37(3), 1245–1273 (2017). https://doi.org/10.1093/imanum/drw042
Borthagaray, J.P., Del Pezzo, L.M., Martínez, S.: Finite element approximation for the fractional eigenvalue problem. J. Sci. Comput. 77(1), 308–329 (2018). https://doi.org/10.1007/s10915-018-0710-1
Bramble, J.H., Pasciak, J.E., Vassilevski, P.S.: Computational scales of Sobolev norms with application to preconditioning. Math. Comput. 69(230), 463–480 (2000). https://doi.org/10.1090/S0025-5718-99-01106-0
Bramble, J.H., Xu, J.: Some estimates for a weighted \(L^2\) projection. Math. Comput. 56(194), 463–476 (1991). https://doi.org/10.2307/2008391
Bramble, J.H., Zhang, X.: The analysis of multigrid methods. In: Ciarlet, P.G., Lions, J.L. (eds.) Handbook of Numerical Analysis, vol. VII, pp. 173–415. North-Holland, Amsterdam (2000)
Chandler-Wilde, S.N., Hewett, D.P., Moiola, A.: Interpolation of Hilbert and Sobolev spaces: quantitative estimates and counterexamples. Mathematika 61(2), 414–443 (2015). https://doi.org/10.1112/S0025579314000278
Ciarlet, P.G.: The finite element method for elliptic problems. In: Classics in Applied Mathematics, vol. 40. Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA (2002). https://doi.org/10.1137/1.9780898719208. Reprint of the 1978 original [North-Holland, Amsterdam; MR0520174 (58 #25001)]
Constantin, P.: Energy spectrum of quasigeostrophic turbulence. Phys. Rev. Lett. 89(18), 184,501 (2002)
Constantin, P., Majda, A.J., Tabak, E.: Formation of strong fronts in the \(2\)-D quasigeostrophic thermal active scalar. Nonlinearity 7(6), 1495–1533 (1994)
Cont, R., Tankov, P.: Financial Modelling with Jump Processes. Chapman & Hall/CRC Financial Mathematics Series. Chapman & Hall/CRC, Boca Raton (2004)
Córdoba, A., Córdoba, D.: A maximum principle applied to quasi-geostrophic equations. Commun. Math. Phys. 249(3), 511–528 (2004). https://doi.org/10.1007/s00220-004-1055-1
Davis, T.A.: Umfpack Version 5.2.0 User Guide. University of Florida, Gainesville (2007)
D’Elia, M., Gunzburger, M.: The fractional Laplacian operator on bounded domains as a special case of the nonlocal diffusion operator. Comput. Math. Appl. 66(7), 1245–1260 (2013). https://doi.org/10.1016/j.camwa.2013.07.022
Demengel, F., Demengel, G.: Functional spaces for the theory of elliptic partial differential equations. Universitext. Springer, London; EDP Sciences, Les Ulis (2012). https://doi.org/10.1007/978-1-4471-2807-6. Translated from the 2007 French original by Reinie Erné
Droniou, J.: A numerical method for fractal conservation laws. Math. Comput. 79(269), 95–124 (2010). https://doi.org/10.1090/S0025-5718-09-02293-5
Du, Q., Gunzburger, M., Lehoucq, R.B., Zhou, K.: Analysis and approximation of nonlocal diffusion problems with volume constraints. SIAM Rev. 54(4), 667–696 (2012). https://doi.org/10.1137/110833294
Dupont, T., Scott, R.: Polynomial approximation of functions in Sobolev spaces. Math. Comput. 34(150), 441–463 (1980). https://doi.org/10.2307/2006095
Faermann, B.: Localization of the Aronszajn–Slobodeckij norm and application to adaptive boundary element methods. II. The three-dimensional case. Numer. Math. 92(3), 467–499 (2002). https://doi.org/10.1007/s002110100319
Gatto, P., Hesthaven, J.S.: Numerical approximation of the fractional Laplacian via \(hp\)-finite elements, with an application to image denoising. J. Sci. Comput. 65(1), 249–270 (2015). https://doi.org/10.1007/s10915-014-9959-1
Grisvard, P.: Elliptic problems in nonsmooth domains. In: Classics in Applied Mathematics, vol. 69. Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA (2011). Reprint of the 1985 original [MR0775683], With a foreword by Susanne C. Brenner
Grubb, G.: Fractional Laplacians on domains, a development of Hörmander’s theory of \(\mu \)-transmission pseudodifferential operators. Adv. Math. 268, 478–528 (2015). https://doi.org/10.1016/j.aim.2014.09.018
Held, I.M., Pierrehumbert, R.T., Garner, S.T., Swanson, K.L.: Surface quasi-geostrophic dynamics. J. Fluid Mech. 282, 1–20 (1995). https://doi.org/10.1017/S0022112095000012
Jerison, D., Kenig, C.E.: The inhomogeneous Dirichlet problem in Lipschitz domains. J. Funct. Anal. 130(1), 161–219 (1995). https://doi.org/10.1006/jfan.1995.1067
Kato, T.: Fractional powers of dissipative operators. J. Math. Soc. Jpn. 13, 246–274 (1961). https://doi.org/10.2969/jmsj/01330246
Kilbas, A.A., Srivastava, H.M., Trujillo, J.J.: Theory and Applications of Fractional Differential Equations, North-Holland Mathematics Studies, vol. 204. Elsevier Science B.V, Amsterdam (2006)
Levendorskiĭ, S.Z.: Pricing of the American put under Lévy processes. Int. J. Theor. Appl. Finance 7(3), 303–335 (2004). https://doi.org/10.1142/S0219024904002463
Lions, J.L., Magenes, E.: Non-Homogeneous Boundary Value Problems and Applications, vol. II. Springer, New York-Heidelberg (1972). Translated from the French by P. Kenneth, Die Grundlehren der mathematischen Wissenschaften, Band 182
Lund, J., Bowers, K.L.: Sinc Methods for Quadrature and Differential Equations. Society for Industrial and Applied Mathematics (SIAM), Philadelphia (1992). https://doi.org/10.1137/1.9781611971637
McLean, W.: Strongly Elliptic Systems and Boundary Integral Equations. Cambridge University Press, Cambridge (2000)
Mueller, C.: The heat equation with Lévy noise. Stoch. Process. Appl. 74(1), 67–82 (1998). https://doi.org/10.1016/S0304-4149(97)00120-8
Olver, F., Lozier, D., Boisvert, R., Clark, C.: NIST digital library of mathematical functions. Online companion to [65]: http://dlmf.nist.gov (2010)
Pham, H.: Optimal stopping, free boundary, and American option in a jump-diffusion model. Appl. Math. Optim. 35(2), 145–164 (1997). https://doi.org/10.1007/s002459900042
Ros-Oton, X., Serra, J.: The Dirichlet problem for the fractional Laplacian: regularity up to the boundary. J. Math. Pures Appl. (9) 101(3), 275–302 (2014). https://doi.org/10.1016/j.matpur.2013.06.003
Sauter, S.A., Schwab, C.: Boundary Element Methods. Springer Series in Computational Mathematics, vol. 39. Springer, Berlin (2011). https://doi.org/10.1007/978-3-540-68093-2. Translated and expanded from the 2004 German original
Scott, L.R., Zhang, S.: Finite element interpolation of nonsmooth functions satisfying boundary conditions. Math. Comput. 54(190), 483–493 (1990). https://doi.org/10.2307/2008497
Stinga, P.R., Torrea, J.L.: Extension problem and Harnack’s inequality for some fractional operators. Commun. Partial Differ. Equ. 35(11), 2092–2122 (2010). https://doi.org/10.1080/03605301003735680
Zaslavsky, G.M.: Chaos, fractional kinetics, and anomalous transport. Phys. Rep. 371(6), 461–580 (2002). https://doi.org/10.1016/S0370-1573(02)00331-9
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Andrea Bonito and Wenyu Lei were supported in part by the National Science Foundation through Grant DMS-1254618 while the Wenyu Lei and Joseph E. Pasciak were supported in part by the National Science Foundation through Grant DMS-1216551.
A. Proof of Lemma 7.6
A. Proof of Lemma 7.6
The proof of Lemma 7.6 requires the following auxiliary localization result. We refer to [26] for a similar result in two dimensional space.
Lemma A.1
For \(r\in (0,1/2)\), let v be in \(H^r(D)\) and \({\widetilde{v}}\) denote the extension by zero of u to \({{{\mathbb {R}}}^d}\). There exists a constant C independent of h such that
with a constant C independent of h.
Proof
Let \(\widetilde{\mathcal {T}}_h(D)\) be any quasi-uniform mesh [(satisfying (63) and (64)] which extends \({\mathcal {T}_h(D)}\) beyond a unit size neighborhood of D. Fix \(\delta >0\) and for \(\tau \in {\mathcal {T}_h(D)}\) set
Let
and let \(\widetilde{\mathcal {T}}_h(D^\delta _h)\) denote the set of \(\tau \in \widetilde{\mathcal {T}}_h(D)\) contained in \(D_h^\delta \). Finally, for \(\tau \in \widetilde{\mathcal {T}}_h(D^\delta _h){\setminus } {\mathcal {T}_h(D)}\), set
Fix \(v\in H^r(D)\). Since \({\widetilde{v}}\) vanishes outside of \(D_h^\delta \),
The second integral above is bounded by
Expanding the first integral gives
Applying the arithmetic-geometric mean inequality gives
As in (97),
Now,
Using this and Fubini’s Theorem gives
Thus \(J_4\le 4 J_5\) and is bounded by the right hand side of (100).
For \(J_3\), we clearly have
For any element \(\tau ' \in \widetilde{\mathcal {T}}_h(D^\delta _h)\), let \(v_{\tau '}\) denote \({\widetilde{v}}\) restricted to \(\tau '\) and extended by zero outside. As \(r\in (0,1/2)\), \(v_{\tau '} \in H^r({{{\mathbb {R}}}^d})\) and satisfies
The constant C above only depends on Lipschitz constants associated with \(\tau '\) (see [22, 28]), which in turn only depend on the constants appearing in (63). We use the triangle inequality to get
and hence a Cauchy–Schwarz inequality implies that
with \(N_\tau \) denoting the number of elements in \({\widetilde{\tau }}\). As the mesh is quasi-uniform, \(N_\tau \) can be bounded independently of h. In addition, the mesh quasi-uniformity condition also implies that each \(\tau ' \in {\mathcal {T}_h(D)}\) is contained in a most a fixed number (independent of h) of \({\widetilde{\tau }}\) (with \(\tau \in \widetilde{\mathcal {T}}_h(D^\delta _h)\)). Thus,
Combining the estimates for \(J_2,J_3\) and \(J_4\) completes the proof of the lemma. \(\square \)
Proof of Lemma 7.6
In this proof, C denotes a generic constant independent of h and j defined later. The inequality (4.1) of [43] guarantees that for \(\tau \in {{\mathcal {T}}}_h\), we have
for any linear polynomial p and \(v\in H^{1}(S_\tau )\). Here \(S_\tau \) denotes the union of \(\tau ^\prime \in {{\mathcal {T}}}_h\) with \(\tau \cap \tau ^\prime \ne \emptyset \).
Now, we map \(\tau \) to the reference element using an affine transformation. The mapping takes \(S_\tau \) to \({\widehat{S}}_\tau \). Our aim is to take advantage of the averaged Taylor polynomial constructed in [25], which requires the domain to be star-shaped with respect to a ball (of uniform diameter). The patch \({\widehat{S}}_\tau \) may not satisfy this property. However, it can be written as the (overlapping) union of domains \({\widehat{D}}_j\) with each \( {\widehat{D}}_j\) consisting of the union of pairs of elements of \({\widehat{S}}_\tau \) sharing a common face. These \({\widehat{D}}_j\) are star-shaped with respect to balls of diameter depending on the shape regularity constant of the subdivision, which is uniform thanks to (63). Hence, the averaged Taylor polynomial \({{\mathcal {Q}}}_j\) constructed in [25] satisfies (see Theorem 6.1 of [25]), for all \(v\in H^{\beta }({\widehat{S}}_j)\),
Taking \(\Vert \cdot \Vert _{{\widehat{D}}_j}\) to be \(\Vert \cdot \Vert _{L^2({\widehat{D}}_j)} \) or \(\Vert \cdot \Vert _{H^1({\widehat{D}}_j)}\) and \(|\cdot |_{{\widehat{D}}_j} = |\cdot |_{H^\beta ({\widehat{D}}_j)}\) in Theorem 7.1 of [25] implies that (102) holds with \({\widehat{D}}_j\) replaced by \({\widehat{S}}_j\). This, (101) and a Bramble-Hilbert argument implies that for \(v\in H^\beta (D)\cap H^1_0(D)\),
Inequality (88) follows from (103) and interpolation.
We cannot use Theorem 7.1 of [25] to derive (87) because of the non-locality of the norm \(|\cdot |_{H^\beta (D)}\). Instead, we apply Lemma A.1, (103), and the fact that \(|\pi _h^{sz}v|_{H^\beta (\tau )}=0\) to obtain, for \(v\in H^\beta (D)\cap H^1_0(D)\),
The norms in (87) can be replaced by \(\Vert \cdot \Vert _{H^\beta (D)}\) and hence (87) follows from (103) and (104). \(\square \)
Rights and permissions
About this article
Cite this article
Bonito, A., Lei, W. & Pasciak, J.E. Numerical approximation of the integral fractional Laplacian. Numer. Math. 142, 235–278 (2019). https://doi.org/10.1007/s00211-019-01025-x
Received:
Revised:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00211-019-01025-x