Nonlinear conjugate gradient method for identifying Young’s modulus of the elasticity imaging inverse problem. (English) Zbl 07480128

Summary: Application of elasticity imaging inverse problem to identify Young’s modulus in the elasticity problems in human’s life is an interesting research area. In this study, we identify the modulus of elasticity for solving elasticity imaging inverse problem using a modified output least-squares method. Numerical convergence in the displacements of the direct problem for elasticity is investigated. To study the elasticity imaging inverse problem in an optimization framework, we utilize the sensitivity and adjoint problems to conceptualize a new model for computing the gradient of the minimizer. Discrete formulae in the model are then used to devise a scheme for an efficient computation gradient of the modified output least-squares objective function using the nonlinear conjugate gradient method. Numerical experiments demonstrate the effectiveness of the proposed technique.


65N20 Numerical methods for ill-posed problems for boundary value problems involving PDEs
65N30 Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs
65N21 Numerical methods for inverse problems for boundary value problems involving PDEs


OpenQSEI; Matlab
Full Text: DOI


[1] Smyl, D.; Antin, K-N; Liu, D., Coupled digital image correlation and quasi-static elasticity imaging of inhomogeneous orthotropic composite structures, Inverse Probl, 34, 12 (2018)
[2] Smyl, D.; Bossuyt, S.; Liu, D., OpenQSEI: a MATLAB package for quasi static elasticity imaging, SoftwareX, 9, 73-76 (2019)
[3] Raghavan, K.; Yagle, AE., Forward and inverse problems in elasticity imaging of soft tissues, IEEE Trans Nucl Sci, 41, 4, 1639-1648 (1994)
[4] Doyley, MM; Meaney, PM; Bamber, JC., Evaluation of an iterative reconstruction method for quantitative elastography, Phys Med Biol, 45, 6, 1521-1540 (2000)
[5] Constantinescu, A., On the identification of elastic moduli from displacement-force boundary measurements, Inverse Probl Eng, 1, 293-313 (1995)
[6] Dominguez, N.; Gibiat, V.; Esquerre, Y., Time domain topological gradient and time reversal analogy: an inverse method for ultrasonic target detection, Wave Motion, 42, 1, 31-52 (2005)
[7] Knopoff, DA; Fernández, D.; Torres, GA, Adjoint method for a tumor growth PDE-constrained optimization problem, Comput Math Appl, 66, 6, 1104-1119 (2013)
[8] Oberai, AA; Gokhale, NH; Feijóo, GR., Solution of inverse problems in elasticity imaging using the adjoint method, Inverse Probl, 19, 2, 297-313 (2003)
[9] Lozano, C., Discrete surprises in the computation of sensitivities from boundary integrals in the continuous adjoint approach to inviscid aerodynamic shape optimization, Comput Fluids, 56, 118-127 (2012)
[10] Jadamba, B.; Khan, A.; Rus, G., A new convex inversion framework for parameter identification in saddle point problems with an application to the elasticity imaging inverse problem of predicting tumor location, SIAM J Appl Math, 74, 5, 1486-1510 (2014)
[11] Parker, KJ; Huang, SR; Musulin, RA, Tissue response to mechanical vibrations for ‘sonoelasticity imaging’, Ultrasound Med Biol, 16, 3, 241-246 (1990)
[12] Fehrenbach, J.; Masmoudi, M.; Souchon, R., Detection of small inclusions by elastography, Inverse Probl, 22, 3, 1055-1069 (2006)
[13] Arnold, A.; Reichling, S.; Bruhns, OT, Efficient computation of the elastography inverse problem by combining variational mesh adaption and a clustering technique, Phys Med Biol, 55, 7, 2035-2056 (2010)
[14] Ji, L.; McLaughlin, J., Recovery of the Lamé parameter μ in biological tissues, Inverse Probl, 20, 1, 1-24 (2003)
[15] McLaughlin, JR; Yoon, J-R., Unique identifiability of elastic parameters from time-dependent interior displacement measurement, Inverse Probl, 20, 1, 25-45 (2003)
[16] Peltz, JJ; Cacuci, DG., Inverse predictive modeling of a spent fuel dissolver model, Nucl Sci Eng, 184, 1-15 (2016)
[17] Cacuci, DG., Inverse predictive modeling of radiation transport through optically thick media in the presence of counting uncertainties, Nucl Sci Eng, 186, 199-223 (2017)
[18] Aguilo, MA; Aquino, W.; Brigham, JC, An inverse problem approach for elasticity imaging through vibroacoustics, IEEE Trans Med Imaging, 29, 1012-1021 (2010)
[19] Jadamba, B.; Khan, AA; Oberai, AA, First-order and second-order adjoint methods for parameter identification problems with an application to the elasticity imaging inverse problem, Inverse Probl Sci Eng, 25, 1768-1787 (2017)
[20] Babaniyi, OA; Oberai, AA; Barbone, PE., Recovering vector displacement estimates in quasistatic elastography using sparse relaxation of the momentum equation, Inverse Probl Sci Eng, 25, 3, 326-362 (2017)
[21] Abdelhamid, T., Simultaneous identification of the spatio-temporal dependent heat transfer coefficient and spatially dependent heat flux using an MCGM in a parabolic system, J Comput Appl Math, 328, 9, 164-176 (2018)
[22] Abdelhamid, T.; Deng, X.; Chen, R., A new method for simultaneously reconstructing the space-time dependent Robin coefficient and heat flux in a parabolic system, Int J Numer Anal Model, 14, 6, 893-915 (2017)
[23] Abdelhamid, T.; Elsheikh, AH; Elazab, A., Simultaneous reconstruction of the time-dependent Robin coefficient and heat flux in heat conduction problems, Inverse Probl Sci Eng, 26, 9, 1231-1248 (2018)
[24] Ammari, H.; Garapon, P.; Jouve, F., Separation of scales in elasticity imaging: a numerical study, J Comput Math, 28, 354-370 (2010)
[25] Harrigana, TP; Konofagou, EE., Estimation of material elastic moduli in elastography: a local method, and an investigation of Poissons ratio sensitivity, J Biomech, 37, 1215-1221 (2004)
[26] Babuška, I.; Suri, M., Locking effects in the finite element approximation of elasticity problems, Numer Math, 62, 1, 439-463 (1992)
[27] Kwon, YW; Bang, H., The finite element method using MATLAB (2000), London: CRC Press
[28] Reddy, JN., An introduction to the finite element method, 2 (1993), New York: McGraw-Hill, New York
[29] Gockenbach, MS; Khan, AA., Identification of Lamé parameters in linear elasticity: a fixed point approach, J Indust Manag Optim, 1, 4, 487-497 (2005)
[30] Jadamba, B.; Khan, AA; Raciti, F., On the inverse problem of identifying Lamé coefficients in linear elasticity, Comput Math Appl, 56, 2, 431-443 (2008)
[31] Cai, Z.; Starke, G., Least-squares methods for linear elasticity, SIAM J Numer Anal, 42, 826-842 (2004)
[32] Houston, DSP; Wihler, TP., An hp-adaptive mixed discontinuous Galerkin FEM for nearly incompressible linear elasticity, Comput Methods Appl Mech Eng, 195, 3224-3246 (2006)
[33] Brezzi, F.; Fortin, M., Mixed and hybrid finite element methods (2012), New York: Springer-Verlag, New York
[34] Alifanov, OM., Inverse heat transfer problems (2012), Berlin: Springer Science & Business Media
[35] Jarny, Y.; Ozisik, M.; Bardon, J., A general optimization method using adjoint equation for solving multidimensional inverse heat conduction, Int J Heat Mass Transfer, 34, 11, 2911-2919 (1991)
[36] Alifanov, OM., Solution of an inverse problem of heat conduction by iteration methods, J Eng Phys Thermophys, 26, 4, 471-476 (1974)
[37] Hansen, P.; O’Leary, D., The use of the L-curve in the regularization of discrete ill-posed problems, SIAM J Sci Comput, 14, 6, 1487-1503 (1993)
[38] Abdelhamid, T.; Elsheikh, AH; Omisore, OM, Reconstruction of the heat transfer coefficients and heat fluxes in heat conduction problems, Math Comput Simul, 187, 134-154 (2021)
[39] Abdelhamid, T.; Chen, R.; Alam, MM., Prediction of multi-parameters in the inverse heat conduction problems, J Phys: Conf Ser, 1707 (2020)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.