Abstract
In simulation of cardiovascular processes and diseases patient-specific model parameters, such as constitutive properties, are usually not easy to obtain. Instead of using population mean values to perform “patient-specific” simulations, thereby neglecting the inter- and intra-patient variations present in these parameters, these uncertainties have to be considered in the computational assessment. However, due to limited computational resources and several shortcomings of traditional uncertainty quantification approaches, parametric uncertainties, modeled as random fields, have not yet been considered in patient-specific, nonlinear, large-scale, and complex biomechanical applications. Hence, the purpose of this study is twofold. First, we present an uncertainty quantification framework based on multi-fidelity sampling and Bayesian formulations. The key feature of the presented method is the ability to rigorously exploit and incorporate information from an approximate, low fidelity model. Most importantly, response statistics of the corresponding high fidelity model can be computed accurately even if the low fidelity model provides only a very poor approximation. The approach merely requires that the low fidelity model and the corresponding high fidelity model share a similar stochastic structure, i.e., dependence on the random input. This results in a tremendous flexibility in choice of the approximate model. The flexibility and capabilities of the framework are demonstrated by performing uncertainty quantification using two patient-specific, large-scale, nonlinear finite element models of abdominal aortic aneurysms. One constitutive parameter of the aneurysmatic arterial wall is modeled as a univariate three-dimensional, non-Gaussian random field, thereby taking into account inter-patient as well as intra-patient variations of this parameter. We use direct Monte Carlo to evaluate the proposed method and found excellent agreement with this reference solution. Additionally, the employed approach results in a tremendous reduction of computational costs, rendering uncertainty quantification with complex patient-specific nonlinear biomechanical models practical for the first time. Second, we also analyze the impact of the uncertainty in the input parameter on mechanical quantities typically related to abdominal aortic aneurysm rupture potential such as von Mises stress, von Mises strain and strain energy. Thus, providing first estimates on the variability of these mechanical quantities due to an uncertain constitutive parameter, and revealing the potential error made by assuming population averaged mean values in patient-specific simulations of abdominal aortic aneurysms. Moreover, the influence of correlation length of the random field is investigated in a parameter study using MC.
Similar content being viewed by others
References
Astrand H, Stalhand J, Karlsson J, Karlsson M, Sonesson B, Lanne T (2011) In vivo estimation of the contribution of elastin and collagen to the mechanical properties in the human abdominal aorta: effect of age and sex. J Appl Physiol 110(1):176–187
Au S, Beck J (2001) Estimation of small failure probabilities in high dimensions by subset simulation. Probabilist Eng Mech 16(4):263–277
Bocchini P, Deodatis G (2008) Critical review and latest developments of a class of simulation algorithms for strongly non-Gaussian random fields. Probabilist Eng Mech 23(4):393–407
Bourgund U, Bucher CG (1986) Importance sampling procedure using design points (ISPUD)—a user’s manual. Institute of Engineering Mechanics, University of Innsbruck
Charmpis DC, Schuëller GI, Pellissetti MF (2007) The need for linking micromechanics of materials with stochastic finite elements: a challenge for materials science. Comput Mater Sci 41(1):27–37
Chen P, Quarteroni A, Rozza G (2013) Simulation based uncertainty quantification of human arterial network hemodynamics. Int J Numer Method Biomed Eng 29(6):698–721
Cho H, Venturi D, Karniadakis GE (2013) Karhunen-Loève expansion for multi-correlated stochastic processes. Probabilist Eng Mech 34:157–167
Cliffe KA, Giles MB, Scheichl R, Teckentrup AL (2011) Multilevel Monte Carlo methods and applications to elliptic PDEs with random coefficients. Comput Vis Sci 14(1):3–15
Del Moral P, Doucet A, Jasra A (2006) Sequential Monte Carlo samplers. J R Stat Soc Ser B (Statistical Methodology) 68(3):411–436
Doucet A, de Freitas N (2001) Sequential Monte Carlo methods in practice. Springer, Berlin
Eldred MS (2009) Recent advances in non-intrusive polynomial chaos and stochastic collocation methods for uncertainty analysis and design. In: 50th AIAA/ASME/ASCE/AHS/ASC structures, structural dynamics, and materials conference, pp 2009–2274
Field RV Jr, Grigoriu M (2012) A method for the efficient construction and sampling of vector-valued translation random fields. Probabilist Eng Mech 29:79–91
Gasser TC, Holzapfel GA (2002) A rate-independent elastoplastic constitutive model for biological fiber-reinforced composites at finite strains: continuum basis, algorithmic formulation and finite element implementation. Comput Mech 29(4–5):340–360
Gasser TC, Görgülü G, Folkesson M, Swedenborg J (2008) Failure properties of intraluminal thrombus in abdominal aortic aneurysm under static and pulsating mechanical loads. J Vasc Surg 48(1):179–188
Gasser TC, Auer M, Labruto F, Swedenborg J, Roy J (2010) Biomechanical rupture risk assessment of abdominal aortic aneurysms: model complexity versus predictability of finite element simulations. Eur J Vasc Endovasc Surg 40(2):176–185
Gee M, Förster C, Wall W (2010) A computational strategy for prestressing patient-specific biomechanical problems under finite deformation. Int J Numer Method Biomed Eng 26(1):52–72
Gee MW, Reeps C, Eckstein HH, Wall WA (2009) Prestressing in finite deformation abdominal aortic aneurysm simulation. J Biomech 42(11):1732–1739
Ghanem R, Spanos P (2003) Stochastic finite elements: a spectral approach. Dover, New York
Grigoriu M (1995) Applied non-Gaussian processes. Prentice Hall, Englewood Cliffs
Hurtado J, Barbat A (1998) Monte Carlo techniques in computational stochastic mechanics. Arch Comput Method E 5(1):3–29
Kemp F (2003) An introduction to sequential Monte Carlo methods. J R Stat Soc Ser D (The Statistician) 52(4):694–695
Kennedy MC, O’Hagan A (2000) Predicting the output from a complex computer code when fast approximations are available. Biometrika 87(1):1–13
Koutsourelakis P (2009) Accurate uncertainty quantification using inaccurate models. SIAM J Sci Comput 31(5):3274–3300
Ma B, Lu J, Harbaugh RE, Raghavan ML (2007) Nonlinear anisotropic stress analysis of anatomically realistic cerebral aneurysms. J Biomech Eng 129(1):88
Maier A, Gee MW, Reeps C, Pongratz J, Eckstein HH, Wall WA (2010) A comparison of diameter, wall stress, and rupture potential index for abdominal aortic aneurysm rupture risk prediction. Ann Biomed Eng 38(10):3124–3134
Marini G, Maier A, Reeps C, Eckstein HH, Wall WA, Gee MW (2011) A continuum description of the damage process in the arterial wall of abdominal aortic aneurysms. Int J Numer Method Biomed Eng 28(1):87–99
McKay MD, Beckman RJ, Conover WJ (1979) Comparison of three methods for selecting values of input variables in the analysis of output from a computer code. Technometrics 21(2):239–245
Metropolis N, Ulam S (1949) The Monte Carlo method. J Am Stat Assoc 44(247):335–341
Miller K, Lu J (2013) On the prospect of patient-specific biomechanics without patient-specific properties of tissues. J Mech Behav Biomed Mater 27:154–166
Moore JA, Steinman DA (1997) Computational blood flow modelling: errors associated with reconstructing finite element models from magnetic resonance images. J Biomech 31(2):179–184
Neal RM (2001) Annealed importance sampling. Stat Comput 11(2):125–139
Olsson A, Sandberg G (2002) Latin hypercube sampling for stochastic finite element analysis. J Eng Mech 128:121
Panayirci HM, Schuëller GI (2011) On the capabilities of the polynomial chaos expansion method within SFE analysis—an overview. Arch Comput Method E 18(1):43–55
Pellissetti M, Schuëller GI (2009) Scalable uncertainty and reliability analysis by integration of advanced Monte Carlo simulation and generic finite element solvers. Comput Struct 87(13–14):930–947
Peña E, Calvo B, Martínez MA, Doblaré M (2008) On finite-strain damage of viscoelastic-fibred materials. Application to soft biological tissues. Int J Numer Meth Eng 74(7):1198–1218
Popescu R, Deodatis G, Prevost J (1998) Simulation of homogeneous nonGaussian stochastic vector fields. Probabilist Eng Mech 13(1):1–13
Pradlwarter H, Schuëller GI, Koutsourelakis P, Charmpis D (2007) Application of line sampling simulation method to reliability benchmark problems. Struct Saf 29(3):208–221
Raghavan M, Vorp D (2000) Toward a biomechanical tool to evaluate rupture potential of abdominal aortic aneurysm: identification of a finite strain constitutive model and evaluation of its applicability. J Biomech 33(4):475–482
Raghavan M, Kratzberg J, Castro de Tolosa E, Hanaoka M, Walker P, da Silva E (2006) Regional distribution of wall thickness and failure properties of human abdominal aortic aneurysm. J Biomech 39(16):3010–3016
Raghavan ML, Webster MW, Vorp DA (1996) Ex vivo biomechanical behavior of abdominal aortic aneurysm: assessment using a new mathematical model. Ann Biomed Eng 24(5):573–582
Raghavan ML, Hanaoka MM, Kratzberg JA, Higuchi MdL, Da Silva ES (2011) Biomechanical failure properties and microstructural content of ruptured and unruptured abdominal aortic aneurysms. J Biomech 44(13):2501–2507
Reeps C, Gee M, Maier A, Gurdan M, Eckstein HH, Wall WA (2010) The impact of model assumptions on results of computational mechanics in abdominal aortic aneurysm. J Vasc Surg 51(3):679–688
Reeps C, Maier A, Pelisek J, Härtl F, Grabher-Meier V, Wall WA, Essler M, Eckstein HH, Gee MW (2012) Measuring and modeling patient-specific distributions of material properties in abdominal aortic aneurysm wall. Biomech Model Mechanobiol 12(4):1–17
Roccabianca S, Figueroa CA, Tellides G, Humphrey JD (2014) Quantification of regional differences in aortic stiffness in the aging human. J Mech Behav Biomed Mater 29:618–634
Rodriguez JF, Alastrue V, Doblare M (2008) Finite element implementation of a stochastic three dimensional finite-strain damage model for fibrous soft tissue. Comput Method Appl M 197(9):946–958
Sankaran S, Marsden AL (2011) A stochastic collocation method for uncertainty quantification and propagation in cardiovascular simulations. J Biomech Eng 133(3):031,001
Schulze-Bauer CA, Holzapfel GA (2003) Determination of constitutive equations for human arteries from clinical data. J Biomech 36(2):165–169
Shields MD, Deodatis G, Bocchini P (2011) A simple and efficient methodology to approximate a general non-Gaussian stationary stochastic process by a translation process. Probabilist Eng Mech 26(4):511–519
Shinozuka M, Deodatis G (1996) Simulation of multi-dimensional Gaussian stochastic fields by spectral representation. Appl Mech Rev 49(1):29–53
Shinozuka M, Jan C (1972) Digital simulation of random processes and its applications. J Sound Vib 25(1):111–128
Simo J (1987) On a fully three-dimensional finite-strain viscoelastic damage model: formulation and computational aspects. Comput Methods Appl Mech Eng 60(2):153–173
Spanos PD, Zeldin BA (1998) Monte Carlo treatment of random fields: a broad perspective. Appl Mech Rev 51(2):219–237
Speelman L, Bosboom EMH, Schurink GWH, Buth J, Breeuwer M, Jacobs MJ, van de Vosse FN (2009) Initial stress and nonlinear material behavior in patient-specific AAA wall stress analysis. J Biomech 42(11):1713–1719
Stefanou G (2009) The stochastic finite element method: past, present and future. Comput Method Appl M 198(9):1031–1051
Stefanou G, Papadrakakis M (2007) Assessment of spectral representation and Karhunen-Loève expansion methods for the simulation of Gaussian stochastic fields. Comput Method Appl M 196(21):2465–2477
Thubrikar J, Labrosse M, Robicsek F, Al-Soudi J, Fowler B (2001) Mechanical properties of abdominal aortic aneurysm wall. J Med Eng Technol 25(4):133–142
Vallabhaneni SR, Gilling-Smith GL, How TV, Carter SD, Brennan JA, Harris PL (2004) Heterogeneity of tensile strength and matrix metalloproteinase activity in the wall of abdominal aortic aneurysms. J Endovasc Ther 11(4):494–502
Vande Geest J, Di Martino E, Bohra A, Makaroun M, Vorp DA (2006a) A biomechanics-based rupture potential index for abdominal aortic aneurysm risk assessment: demonstrative application. Ann NY Acad Sci 1085(1):11–21
Vande Geest J, Sacks M, Vorp D (2006b) The effects of aneurysm on the biaxial mechanical behavior of human abdominal aorta. J Biomech 39(7):1324–1334
Vande Geest J, Wang D, Wisniewski S, Makaroun M, Vorp D (2006c) Towards a noninvasive method for determination of patient-specific wall strength distribution in abdominal aortic aneurysms. Ann Biomed Eng 34(7):1098–1106
Vanmarcke E (2010) Random fields: analysis and synthesis. World Scientific, Singapore
Vappou J, Luo J, Konofagou E (2010) Pulse wave imaging for noninvasive and quantitative measurement of arterial stiffness in vivo. Am J Hypertens 23(4):393–398
Volokh K, Vorp D (2008) A model of growth and rupture of abdominal aortic aneurysm. J Biomech 41(5):1015–1021
Volokh KY (2010) Comparison of biomechanical failure criteria for abdominal aortic aneurysm. J Biomech 43(10):2032–2034
Vořechovský M (2008) Simulation of simply cross correlated random fields by series expansion methods. Struct Saf 30(4):337–363
Wall WA, Gee MW (2014) BACI: a parallel multiphysics simulation environment. Tech. rep., Institute for Computational Mechanics, Technische Universität München
Xiu D (2009) Fast numerical methods for stochastic computations: a review. Commun Comput Phys 5(2–4):242–272
Xiu D (2010) Numerical methods for stochastic computations: a spectral method approach. Princeton University Press, Princeton
Xiu D, Sherwin SJ (2007) Parametric uncertainty analysis of pulse wave propagation in a model of a human arterial network. J Comput Phys 226(2):1385–1407
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Biehler, J., Gee, M.W. & Wall, W.A. Towards efficient uncertainty quantification in complex and large-scale biomechanical problems based on a Bayesian multi-fidelity scheme. Biomech Model Mechanobiol 14, 489–513 (2015). https://doi.org/10.1007/s10237-014-0618-0
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10237-014-0618-0