Abstract
The inverse eigenvalue problem for a weighted Helmholtz equation is investigated. Based on the finite spectral data, the density function is estimated. The inverse problem is formulated as a least squared functional with respect to the density function, with a \(L^2\) regularity term. The continuity of the eigenpairs with respect to the density is proved. Mathematical properties of the continuous and the discrete optimization problems are established. A conjugate gradient algorithm is proposed. Numerical results for 1D and 2D inverse eigenvalue problem of the weighted Helmholtz equation are presented to illustrate the effectiveness and efficiency of the proposed algorithm.
Similar content being viewed by others
Data Availability
The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.
References
Chu, M.T.: Inverse eigenvalue problems. SIAM Rev. 40(1), 1–39 (1998)
McLaughlin, J.R.: Analytical methods for recovering coefficients in differential equations from spectral data. SIAM Rev. 28(1), 53–72 (1986)
Nachman, A., Sylvester, J., Uhlmann, G.: An n-dimensional Borg–Levinson theorem. Commun. Math. Phys. 115(4), 595–605 (1988)
Rundell, W., Sacks, P.E.: Reconstruction techniques for classical inverse Sturm–Liouville problems. Math. Comput. 58(197), 161–183 (1992)
Barcilon, V.: A two-dimensional inverse eigenvalue problem. Inverse Prob. 6(1), 11–20 (1990)
McCarthy, C.M.: The inverse eigenvalue problem for a weighted Helmholtz equation. Appl. Anal. 77(1–2), 77–96 (2001)
McCarthy, C.M.: Recovery of a density from the eigenvalues of a nonhomogeneous membrane. In: Proceedings of the Third International Conference on Inverse Problems in Engineering: Theory and Practice, Port Ludlow, Washington (1999)
Andrew, A.L.: Numerical solution of inverse Sturm–Liouville problems. ANZIAM J. 45, 326–337 (2004)
Drignei, M.C.: Constructibility of an \({\rm L}_{\mathbb{R}}^2(0, a)\) solution to an inverse Sturm–Liouville problem using three dirichlet spectra. Inverse Prob. 26(2), 025003 (2009)
Jiang, X., Li, X., Xu, X.: Numerical algorithms for inverse Sturm–Liouville problems. Numer. Algorithms 89, 1278–1309 (2022)
Gao, Q., Huang, Z., Cheng, X.: A finite difference method for an inverse Sturm–Liouville problem in impedance form. Numer. Algorithms 70(3), 669–690 (2015)
Zhang, W., Du, Z., Sun, G., Guo, X.: A level set approach for damage identification of continuum structures based on dynamic responses. J. Sound Vib. 386, 100–115 (2017)
Zhang, Z., Dai, X., Chen, W.: A piecewise constant level set method for damage identification of continuum structures based on natural frequencies. Struct. Multidiscip. Optim. 60(6), 2477–2491 (2019)
Lee, U., Shin, J.: A frequency response function-based structural damage identification method. Comput. Struct. 80(2), 117–132 (2002)
Babuška, I., Osborn, J.: Eigenvalue problems. In: Finite Element Methods (Part 1). Handbook of Numerical Analysis, vol. 2, pp. 641–787. Elsevier, New York (1991)
Sun, J., Zhou, A.: Finite Element Methods for Eigenvalue Problems. CRC Press, New York (2016)
Babuška, I., Osborn, J.E.: Estimates for the errors in eigenvalue and eigenvector approximation by Galerkin methods, with particular attention to the case of multiple eigenvalues. SIAM J. Numer. Anal. 24(6), 1249–1276 (1987)
Babuška, I., Osborn, J.E.: Finite element-Galerkin approximation of the eigenvalues and eigenvectors of selfadjoint problems. Math. Comput. 52(186), 275–297 (1989)
Gilbert, A.D., Graham, I.G., Kuo, F.Y., Scheichl, R., Sloan, I.H.: Analysis of quasi-monte Carlo methods for elliptic eigenvalue problems with stochastic coefficients. Numer. Math. 142(4), 863–915 (2019)
Ciarlet, P.G.: The Finite Element Method for Elliptic Problems. Society for Industrial and Applied Mathematics, Paris (2002)
Jin, B., Zou, J.: Numerical estimation of the robin coefficient in a stationary diffusion equation. IMA J. Numer. Anal. 30(3), 677–701 (2010)
Fletcher, R., Reeves, C.M.: Function minimization by conjugate gradients. Comput. J. 7(2), 149–154 (1964)
Funding
This research is supported by Natural Science Foundation of Zhejiang Province, China (No. LY21A010011) and the European Union’s Horizon 2020 Research and Innovation Programme under the Marie Sklodowska-Cutie Grand Agreement (No. 823731 CONMECH).
Author information
Authors and Affiliations
Contributions
All authors contributed to the study conception and algorithm design. The theoretical part is derived mainly by ZZ and XC. The numerical part is mainly performed by XG. The first draft of the manuscript was written by XG and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Conflict of interest
The authors have no relevant financial or non-financial interests to disclose.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Zhang, Z., Gao, X. & Cheng, X. Numerical Estimation of the Inverse Eigenvalue Problem for a Weighted Helmholtz Equation. J Sci Comput 96, 16 (2023). https://doi.org/10.1007/s10915-023-02242-1
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1007/s10915-023-02242-1