Startseite Laplace inversion for the solution of an abstract heat equation without the forward transform of the source term
Artikel
Lizenziert
Nicht lizenziert Erfordert eine Authentifizierung

Laplace inversion for the solution of an abstract heat equation without the forward transform of the source term

  • Shu-Lin Wu EMAIL logo
Veröffentlicht/Copyright: 29. Juni 2016

Abstract

We consider the discretization in time of inhomogeneous parabolic equations, using the technique of Laplace inversion along a contour located in the complex left half-plane which, after transformation to a finite interval, is then evaluated to high accuracy by a contour quadrature rule. This reduces the problem to a finite set of elliptic equations with complex coefficients, which may be solved in parallel. A serious problem is how to treat the source term f(t), because at each quadrature node along the contour we need its Laplace forward transform, which unfortunately is often unavailable. In this paper, we propose a new contour quadrature which does not require direct use of the Laplace forward transform of f(t). Compared to the existing contour quadratures, error analysis shows that the new quadrature possesses competitive asymptotic order of accuracy and numerical results show that when regularity of the initial term and/or differentiability of f(t) is not satisfied, the new quadrature is more accurate.

JEL Classification: 65R20; 45L05; 65L20

Funding statement: This work was supported by the NSF of China (11301362, 11371157, 61573010), the project of Sichuan University of Science and Engineering (2015LX01), the China Postdoctoral Foundation(2015M580777, 2016T- 90841), the project of Technology and Education of Sichuan Province (2014 JQ0035, 15A0220).

Acknowledgement

The author is very grateful to the anonymous referees for the careful reading of a preliminary version of the manuscript and their valuable suggestions and comments, which greatly improved the quality of this paper.

References

[1] R. Bellman, R. E. Kalaba, and J. A. Lockett, Numerical Inversion of the Laplace Transform: Applications to Biology, Economics, Engineering, and Physics, Elsevier, New York, 1966.Suche in Google Scholar

[2] C. A. Brebbia, J. C. F. Telles, and L. Wrobel, Boundary Element Techniques: Theory and Applications in Engineering, Springer, Berlin–New York, 1984.10.1007/978-3-642-48860-3Suche in Google Scholar

[3] A. M. Cohen, Numerical Methods for Laplace Transform Inversion, Springer, New York, 2007.Suche in Google Scholar

[4] D. G. Duffy, Transform Methods for Solving Partial Differential Equations 2nd edition, CRC, Boca Raton–London–New York–Washington, D.C., 2004, pp. 191–220.10.1201/9781420035148Suche in Google Scholar

[5] B. Dingfelder and J. A. C. Weideman, An improved Talbot method for numerical Laplace transform inversion, Numer. Algor., 68 (2015), 167–183.10.1007/s11075-014-9895-zSuche in Google Scholar

[6] J. Douglas Jr., J. E. Santos, D. Sheen, and L. S. Bennethum, Frequency domain treatment of one-dimensional scalar waves, Math. Models Methods Appl. Sci., 3 (1993), 171–194.10.1142/S0218202593000102Suche in Google Scholar

[7] J. Douglas Jr., J. E. Santos, and D. Sheen, Approximation of scalar waves in the space frequency domain, Math. Models Methods Appl. Sci. 4 (1994), 509–531.10.1142/S0218202594000297Suche in Google Scholar

[8] I. P. Gavrilyuk and V. L. Makarov, Exponentially convergent algorithms for the operator exponential with applications to inhomogeneous problems in Banach spaces, SIAM J. Numer. Anal., 43 (2005), 2144–2171.10.1137/040611045Suche in Google Scholar

[9] N. Hale and J. A. C. Weideman, Contour integral solution of elliptic PDEs in cylindrical domains, SIAM J. Sci. Comput., 37 (2015), A2630–A2655.10.1137/15M1032764Suche in Google Scholar

[10] K. J. In’t Hout and J. A. C. Weideman, A contour integral method for the Black–Scholes and Heston equations, SIAM J. Sci. Comput., 33 (2011), 763–785.10.1137/090776081Suche in Google Scholar

[11] K. L. Kuhlman, Review of inverse Laplace transform algorithms for Laplace-space numerical approaches, Numer. Algorithm., 63 (2013), 339-355.10.1007/s11075-012-9625-3Suche in Google Scholar

[12] K. Kwon and D. Sheen, A parallel method for the numerical solution of integro–differential equation with positive memory, Comput. Meth. Appl. Mech. Engrg., 192 (2003), 4641–4658.10.1016/S0045-7825(03)00442-0Suche in Google Scholar

[13] P. O. Kano, M. Brio, and J. V. Moloney, Application of Weeks method for the numerical inversion of the Laplace transform to the matrix exponential, Comm. Math. Sci., 3 (2005), 335–372.10.4310/CMS.2005.v3.n3.a4Suche in Google Scholar

[14] C. Lanczos, Applied Analysis, Prentice-Hall, Englewood Cliffs, N.J., 1956.Suche in Google Scholar

[15] M. López-Fernández, C. Palencia, and A. Schädle, A spectral order method for inverting sectorial Laplace transforms, SIAM J. Numer. Anal., 44 (2006), 1332–1350.10.1137/050629653Suche in Google Scholar

[16] W. Liao and Y. Yan, Singly diagonally implicit Runge–Kutta method for time-dependent reaction diffusion equation, Numer. Meth. Part. Diff. Equ., 27 (2011), 1423–1441.10.1002/num.20589Suche in Google Scholar

[17] W. Mclean, I. H. Sloan, and V. Thomée, Time discretization via Laplace transformation of an integro-differential equation of parabolic type, Numer. Math., 102 (2006), 497–522.10.1007/s00211-005-0657-7Suche in Google Scholar

[18] W. Mclean and V. Thomée, Time discretization of an evolution equation via Laplace transformation, IMA J. Numer. Anal., 24 (2004), 439–463.10.1093/imanum/24.3.439Suche in Google Scholar

[19] W. Mclean and V. Thomée, Numerical solution via Laplace transformation of a fractional order evolution equation, J. Integr. Equ. Appl., 22 (2010), 57–94.10.1216/JIE-2010-22-1-57Suche in Google Scholar

[20] W. Mclean and V. Thomée, Maximum-norm error analysis of a numerical solution via Laplace transformation and quadrature of a fractional-order evolution equation, IMA J. Numer. Anal., 30 (2010), 208–230.10.1093/imanum/drp004Suche in Google Scholar

[21] E. Morales-Casique and S. P. Neuman, Laplace-transform finite element solution of nonlocal and localized stochastic moment equations of transport, Commun. comput. phys., 6 (2009), 131–16110.4208/cicp.2009.v6.p131Suche in Google Scholar

[22] R. A. Schapery, Two simple approximate methods of Laplace transform inversion for viscoelastic stress analysis, Technical Report GALCIT SM 61-23, California Institute of Technology, 1961.Suche in Google Scholar

[23] D. Sheen, I. H. Sloan, and V. Thomée, A parallel method for time discretization of parabolic problems based on contour integral representation and quadrature, Math. Comp., 69 (1999), 177–195.10.1090/S0025-5718-99-01098-4Suche in Google Scholar

[24] D. Sheen, I. H. Sloan, and V. Thomée, A parallel method for time discretization of parabolic equations based on Laplace transformation and quadrature, IMA J. Numer. Anal., 23 (2003), 269–299.10.1093/imanum/23.2.269Suche in Google Scholar

[25] E. A. Sudicky, The Laplace transform Galerkin technique: A time-continuous finite element theory and application to mass transport in groundwater, Water Resour. Res., 25 (1989), 1833–1846.10.1029/WR025i008p01833Suche in Google Scholar

[26] A. Talbot, The accurate numerical inversion of Laplace transforms, J. Inst. Math. Appl., 23 (1979), 97–120.10.1093/imamat/23.1.97Suche in Google Scholar

[27] L. N. Trefethen and J. A. C. Weideman, The exponentially convergent trapezoidal rule, SIAM Rev., 56 (2014), 385–458.10.1137/130932132Suche in Google Scholar

[28] J. A. C. Weideman, Improved contour integral methods for parabolic PDEs, IMA J. Numer. Anal., 30 (2010), 334–350.10.1093/imanum/drn074Suche in Google Scholar

[29] J. A. C. Weideman and L. N. Trefethen, Parabolic and hyperbolic contours for computing the Bromwich integral, Math. Comp., 76 (2007), 1341–1356.10.1090/S0025-5718-07-01945-XSuche in Google Scholar

Received: 2016-2-22
Accepted: 2016-6-24
Published Online: 2016-6-29
Published in Print: 2017-9-26

© 2016 by Walter de Gruyter Berlin/Boston

Heruntergeladen am 14.10.2025 von https://www.degruyterbrill.com/document/doi/10.1515/jnma-2016-1014/html?lang=de
Button zum nach oben scrollen