Rupture risk in abdominal aortic aneurysms: A realistic assessment of the explicit GPU approach.

Accurate estimation of peak wall stress (PWS) is the crux of biomechanically motivated rupture risk assessment for abdominal aortic aneurysms aimed to improve clinical outcomes. Such assessments often use the finite element (FE) method to obtain PWS, albeit at a high computational cost, motivating simplifications in material or element formulations. These simplifications, while useful, come at a cost of reliability and accuracy. We achieve research-standard accuracy and maintain clinically applicable speeds by using novel computational technologies. We present a solution using our custom finite element code based on graphics processing unit (GPU) technology that is able to account for added complexities involved with more physiologically relevant solutions, e.g. strong anisotropy and heterogeneity. We present solutions up to 17× faster relative to an established finite element code using state-of-the-art nonlinear, anisotropic and nearly-incompressible material descriptions. We show a realistic assessment of the explicit GPU FE approach by using complex problem geometry, biofidelic material law, double-precision floating point computation and full element integration. Due to the increased solution speed without loss of accuracy, shown on five clinical cases of abdominal aortic aneurysms, the method shows promise for clinical use in determining rupture risk of abdominal aortic aneurysms.

[1]  Eric Darve,et al.  Assembly of finite element methods on graphics processors , 2011 .

[2]  J Swedenborg,et al.  Survival in patients with abdominal aortic aneurysms. Comparison between operative and nonoperative management. , 1990, European journal of vascular surgery.

[3]  Jos Vander Sloten,et al.  GPU-Based Fast Finite Element Solution for Nonlinear Anisotropic Material Behavior and Comparison of Integration Strategies , 2016 .

[4]  Kathrin Abendroth,et al.  Nonlinear Finite Elements For Continua And Structures , 2016 .

[5]  Manuel Doblaré,et al.  Assessing the Use of the “Opening Angle Method” to Enforce Residual Stresses in Patient-Specific Arteries , 2007, Annals of Biomedical Engineering.

[6]  Karol Miller,et al.  An efficient hourglass control implementation for the uniform strain hexahedron using the Total Lagrangian formulation , 2007 .

[7]  Karol Miller,et al.  A Total Lagrangian based method for recovering the un-deformed configuration in finite elasticity , 2015 .

[8]  R. Courant,et al.  On the Partial Difference Equations, of Mathematical Physics , 2015 .

[9]  Enrique J. Gómez,et al.  High-quality conforming hexahedral meshes of patient-specific abdominal aortic aneurysms including their intraluminal thrombi , 2014, Medical & Biological Engineering & Computing.

[10]  Michael A. Puso,et al.  A stabilized nodally integrated tetrahedral , 2006 .

[11]  Philippe Pierre Pebay,et al.  The verdict geometric quality library. , 2006 .

[12]  Stephane Cotin,et al.  Efficient Nonlinear FEM for Soft Tissue Modelling and Its GPU Implementation within the Open Source Framework SOFA , 2008, ISBMS.

[13]  Ted Belytschko,et al.  A survey of numerical methods and computer programs for dynamic structural analysis , 1976 .

[14]  T. Prevost,et al.  Growth rates and risk of rupture of abdominal aortic aneurysms , 1998, The British journal of surgery.

[15]  Karol Miller,et al.  An adaptive dynamic relaxation method for solving nonlinear finite element problems. Application to brain shift estimation , 2011, International journal for numerical methods in biomedical engineering.

[16]  Eric F Darve,et al.  A new sparse matrix vector multiplication graphics processing unit algorithm designed for finite element problems , 2015 .

[17]  Manuel Doblaré,et al.  Numerical framework for patient‐specific computational modelling of vascular tissue , 2010 .

[18]  J Vander Sloten,et al.  Non-invasive, energy-based assessment of patient-specific material properties of arterial tissue , 2015, Biomechanics and modeling in mechanobiology.

[19]  K. Miller,et al.  Total Lagrangian explicit dynamics finite element algorithm for computing soft tissue deformation , 2006 .

[20]  Peter Huthwaite,et al.  Accelerated finite element elastodynamic simulations using the GPU , 2014, J. Comput. Phys..

[21]  Nicola Parolini,et al.  An explicit dynamics GPU structural solver for thin shell finite elements , 2015 .

[22]  Lianghao Han,et al.  Development of patient-specific biomechanical models for predicting large breast deformation , 2012, Physics in medicine and biology.

[23]  Elena S. Di Martino,et al.  Is There a Role for Biomechanical Engineering in Helping to Elucidate the Risk Profile of the Thoracic Aorta? , 2016, The Annals of thoracic surgery.

[24]  P Regitnig,et al.  Effects of age on the elastic properties of the intraluminal thrombus and the thrombus-covered wall in abdominal aortic aneurysms: biaxial extension behaviour and material modelling. , 2011, European journal of vascular and endovascular surgery : the official journal of the European Society for Vascular Surgery.

[25]  Karol Miller,et al.  Real-Time Nonlinear Finite Element Computations on GPU - Application to Neurosurgical Simulation. , 2010, Computer methods in applied mechanics and engineering.

[26]  R. Ogden,et al.  A New Constitutive Framework for Arterial Wall Mechanics and a Comparative Study of Material Models , 2000 .

[27]  Karol Miller,et al.  Computation of intra-operative brain shift using dynamic relaxation. , 2009, Computer methods in applied mechanics and engineering.

[28]  Elliot L Chaikof,et al.  The care of patients with an abdominal aortic aneurysm: the Society for Vascular Surgery practice guidelines. , 2009, Journal of vascular surgery.

[29]  Sébastien Ourselin,et al.  NiftySim: A GPU-based nonlinear finite element package for simulation of soft tissue biomechanics , 2014, International Journal of Computer Assisted Radiology and Surgery.

[30]  G. Moneta,et al.  Screening for Abdominal Aortic Aneurysm: A Best-Evidence Systematic Review for the U.S. Preventive Services Task Force , 2006 .

[31]  Sebastian Vogt,et al.  IN VIVO DETERMINATION OF ELASTIC PROPERTIES OF THE HUMAN AORTA BASED ON 4D ULTRASOUND DATA , 2012 .

[32]  Martin Auer,et al.  Reconstruction and Finite Element Mesh Generation of Abdominal Aortic Aneurysms From Computerized Tomography Angiography Data With Minimal User Interactions , 2010, IEEE Transactions on Medical Imaging.

[33]  Peter Regitnig,et al.  Layer-specific damage experiments and modeling of human thoracic and abdominal aortas with non-atherosclerotic intimal thickening. , 2012, Journal of the mechanical behavior of biomedical materials.

[34]  Thomas Sangild Sørensen,et al.  Exploring Parallel Algorithms for Volumetric Mass-Spring-Damper Models in CUDA , 2008, ISBMS.

[35]  J. D. Humphrey,et al.  A Microstructurally Motivated Model of Arterial Wall Mechanics with Mechanobiological Implications , 2013, Annals of Biomedical Engineering.

[36]  J. C. Schulz Finite element hourglassing control , 1985 .

[37]  Gerhard A Holzapfel,et al.  A method for incorporating three-dimensional residual stretches/stresses into patient-specific finite element simulations of arteries. , 2015, Journal of the mechanical behavior of biomedical materials.

[38]  Pascal Verdonck,et al.  Full-hexahedral structured meshing for image-based computational vascular modeling. , 2011, Medical engineering & physics.

[39]  Frank Lindseth,et al.  Medical image segmentation on GPUs - A comprehensive review , 2015, Medical Image Anal..

[40]  R. Ogden,et al.  Hyperelastic modelling of arterial layers with distributed collagen fibre orientations , 2006, Journal of The Royal Society Interface.

[41]  David Atkinson,et al.  On modelling of anisotropic viscoelasticity for soft tissue simulation: Numerical solution and GPU execution , 2009, Medical Image Anal..

[42]  Jos Vander Sloten,et al.  Analyzing the potential of GPGPUs for real-time explicit finite element analysis of soft tissue deformation using CUDA , 2015 .

[43]  Franz Maier,et al.  Human thoracic and abdominal aortic aneurysmal tissues: Damage experiments, statistical analysis and constitutive modeling. , 2015, Journal of the mechanical behavior of biomedical materials.

[44]  F. N. van de Vosse,et al.  Patient-specific initial wall stress in abdominal aortic aneurysms with a backward incremental method. , 2007, Journal of biomechanics.

[45]  Sébastien Ourselin,et al.  Fast Deformation Simulation of Breasts Using GPU-Based Dynamic Explicit Finite Element Method , 2010, Digital Mammography / IWDM.

[46]  Mark F Fillinger,et al.  In vivo analysis of mechanical wall stress and abdominal aortic aneurysm rupture risk. , 2002, Journal of vascular surgery.

[47]  Joris Degroote,et al.  A computational method to assess the in vivo stresses and unloaded configuration of patient-specific blood vessels , 2013, J. Comput. Appl. Math..

[48]  Evelyn P Whitlock,et al.  Screening for Abdominal Aortic Aneurysm: A Best-Evidence Systematic Review for the U.S. Preventive Services Task Force , 2005, Annals of Internal Medicine.

[49]  Gerhard Sommer,et al.  Determination of layer-specific mechanical properties of human coronary arteries with nonatherosclerotic intimal thickening and related constitutive modeling. , 2005, American journal of physiology. Heart and circulatory physiology.

[50]  Gerhard A. Holzapfel,et al.  A generalized prestressing algorithm for finite element simulations of preloaded geometries with application to the aorta , 2014, International journal for numerical methods in biomedical engineering.

[51]  J. Grandy Efficient computation of volume of hexahedral cells , 1997 .

[52]  T. Belytschko,et al.  A uniform strain hexahedron and quadrilateral with orthogonal hourglass control , 1981 .

[53]  J Swedenborg,et al.  Biomechanical rupture risk assessment of abdominal aortic aneurysms: model complexity versus predictability of finite element simulations. , 2010, European journal of vascular and endovascular surgery : the official journal of the European Society for Vascular Surgery.

[54]  Edoardo Mazza,et al.  On non-physical response in models for fiber-reinforced hyperelastic materials , 2010 .