Unpublished conference/Abstract (Scientific congresses, symposiums and conference proceedings)
Hyperelastic Elastography in a Large-Scale Bayesian Inversion Setting
Hale, Jack; Farrel, Patrick E.; Bordas, Stéphane
201513th US National Congress on Computational Mechanics


Full Text
Publisher postprint (2.61 kB)

All documents in ORBilu are protected by a user license.

Send to


Keywords :
Bayesian; adjoints; uncertainty; hyperelastic; noise; elastography
Abstract :
[en] We consider the problem of recovering the material parameters of a hyperelastic material [1] in the Bayesian inversion setting. In the Bayesian setting we can extract the statistics associated with various sources of uncertainty, including noise, experimental deficiencies and incomplete observations of the domain. This will allow medical practitioners to make superior diagnosis decisions when presented with a quantitative measure of uncertainty in the recovered parameters. On the assumption of a Gaussian additive noise model, a Gaussian prior and a linear forward model, the posterior distribution of the material parameters given the observations will also be Gaussian. To ensure that the assumption of a linear forward model is valid, and that the posterior is approximated sufficiently well by a Gaussian distribution, we place a limit on the strain regime in which our current methodology applies. We are developing MCMC methods for exploring the non-Gaussian statistics of the posterior distribution. In the linear case, the covariance matrix of the posterior distribution is then characterised by the inverse of the Hessian of the objective functional evaluated at its minimiser. To extract statistical information from the large and dense Hessian we perform a low-rank approximation of the Hessian [2]. The eigenvectors associated with the lowest eigenvalues are the directions in parameter space that are least constrained by the observations. We implement this work within the dolfin-adjoint [3] software package. We derive the MPI-parallel finite element discretisation of the forward, adjoint (1st and 2nd order), and tangent linear models using the high-level differentiation tools available within the FEniCS project. We show results demonstrating the effects of partial observations and poor experimental design on the reliability of the recovered parameters. [1] N. H. Gokhale, P. E. Barbone, and A. A. Oberai, “Solution of the nonlinear elasticity imaging inverse problem: the compressible case,” Inverse Problems, 10.1088/0266-5611/24/4/045010 [2] H. P. Flath, L. C. Wilcox, V. Akçelik, J. Hill, B. van Bloemen Waanders, and O. Ghattas, “Fast Algorithms for Bayesian Uncertainty Quantification in Large-Scale Linear Inverse Problems Based on Low-Rank Partial Hessian Approximations,” SIAM J. Sci. Comput., 10.1137/090780717 [3] P. Farrell, D. Ham, S. Funke, and M. Rognes, “Automated Derivation of the Adjoint of High-Level Transient Finite Element Programs,” SIAM J. Sci. Comput., 10.1137/120873558
Disciplines :
Materials science & engineering
Author, co-author :
Hale, Jack  ;  University of Luxembourg > Faculty of Science, Technology and Communication (FSTC) > Engineering Research Unit
Farrel, Patrick E.;  Oxford University > Mathematical Institute
Bordas, Stéphane ;  University of Luxembourg > Faculty of Science, Technology and Communication (FSTC) > Engineering Research Unit
External co-authors :
Language :
Title :
Hyperelastic Elastography in a Large-Scale Bayesian Inversion Setting
Publication date :
July 2015
Number of pages :
Event name :
13th US National Congress on Computational Mechanics
Event organizer :
United States Association for Computational Mechanics
Event place :
San Diego, United States
Event date :
from 26-07-2015 to 30-07-2015
Audience :
Focus Area :
Computational Sciences
European Projects :
FP7 - 279578 - REALTCUT - Towards real time multiscale simulation of cutting in non-linear materials with applications to surgical simulation and computer guided surgery
Funders :
FNR - Fonds National de la Recherche [LU]
EPSRC - Engineering and Physical Sciences Research Council [GB]
CE - Commission Européenne [BE]
Available on ORBilu :
since 30 March 2015


Number of views
396 (18 by Unilu)
Number of downloads
32 (2 by Unilu)


Similar publications

Contact ORBilu