scholarly journals Greedy maximin distance sampling based model order reduction of prestressed and parametrized abdominal aortic aneurysms

Author(s):  
Alexander Schein ◽  
Michael W. Gee

AbstractThis work proposes a framework for projection-based model order reduction (MOR) of computational models aiming at a mechanical analysis of abdominal aortic aneurysms (AAAs). The underlying full-order model (FOM) is patient-specific, stationary and nonlinear. The quantities of interest are the von Mises stress and the von Mises strain field in the AAA wall, which result from loading the structure to the level of diastolic blood pressure at a fixed, imaged geometry (prestressing stage) and subsequent loading to the level of systolic blood pressure with associated deformation of the structure (deformation stage). Prestressing is performed with the modified updated Lagrangian formulation (MULF) approach. The proposed framework aims at a reduction of the computational cost in a many-query context resulting from model uncertainties in two material and one geometric parameter. We apply projection-based MOR to the MULF prestressing stage, which has not been presented to date. Additionally, we propose a reduced-order basis construction technique combining the concept of subspace angles and greedy maximin distance sampling. To further achieve computational speedup, the reduced-order model (ROM) is equipped with the energy-conserving mesh sampling and weighting hyper reduction method. Accuracy of the ROM is numerically tested in terms of the quantities of interest within given bounds of the parameter domain and performance of the proposed ROM in the many-query context is demonstrated by comparing ROM and FOM statistics built from Monte Carlo sampling for three different patient-specific AAAs.

2010 ◽  
Vol 132 (7) ◽  
Author(s):  
Clark A. Meyer ◽  
Carine Guivier-Curien ◽  
James E. Moore

How much and how the thrombus supports the wall of an abdominal aortic aneurysm (AAA) is unclear. While some previous studies have indicated that thrombus lacks the mechanical integrity to support much load compared with the aneurysm wall, others have shown that removing thrombus in computational AAA models drastically changes aneurysm wall stress. Histopathological studies have shown that thrombus properties vary through the thickness and it can be porous. The goal of this study is to explore the variations in thrombus properties, including the ability to isolate pressure from the aneurysm wall, incomplete attachment, and their effects on aneurysm wall stress, an important parameter in determining risk for rupture. An analytical model comprised of cylinders and two patient specific models were constructed with pressurization boundary conditions applied at the lumen or the thrombus/aneurysm wall interface (to simulate complete transmission of pressure through porous thrombus). Aneurysm wall stress was also calculated in the absence of thrombus. The potential importance of partial thrombus attachment was also analyzed. Pressurizing at either surface (lumen versus interface) made little difference to mean von Mises aneurysm wall stress values with thrombus completely attached (3.1% analytic, 1.2% patient specific) while thrombus presence reduced mean von Mises stress considerably (79% analytic, 40–46% patient specific) in comparison to models without it. Peak von Mises stresses were similarly influenced with pressurization surface differing slightly (3.1% analytic, 1.4% patient specific) and reductions in stress by thrombus presence (80% analytic, 28–37% patient specific). The case of partial thrombus attachment was investigated using a cylindrical model in which there was no attachment between the thrombus and aneurysm wall in a small area (10 deg). Applying pressure at the lumen resulted in a similar stress field to fully attached thrombus, whereas applying pressure at the interface resulted in a 42% increase in peak aneurysm wall stress. Taken together, these results show that the thrombus can have a wall stress reducing role even if it does not shield the aneurysm wall from direct pressurization—as long as the thrombus is fully attached to the aneurysm wall. Furthermore, the potential for porous thrombus to transmit pressure to the interface can result in a considerable increase in aneurysm wall stress in cases of partial attachment. In the search for models capable of accurately assessing the risk for rupture, the nature of the thrombus and its attachment to the aneurysm wall must be carefully assessed.


2019 ◽  
Vol 170 ◽  
pp. 95-106 ◽  
Author(s):  
Nathan Lauzeral ◽  
Domenico Borzacchiello ◽  
Michael Kugler ◽  
Daniel George ◽  
Yves Rémond ◽  
...  

2013 ◽  
Vol 136 (1) ◽  
Author(s):  
Sergio Ruiz de Galarreta ◽  
Aitor Cazón ◽  
Raúl Antón ◽  
Ender A. Finol

The goal of this work is to develop a framework for manufacturing nonuniform wall thickness replicas of abdominal aortic aneurysms (AAAs). The methodology was based on the use of computed tomography (CT) images for virtual modeling, additive manufacturing for the initial physical replica, and a vacuum casting process and range of polyurethane resins for the final rubberlike phantom. The average wall thickness of the resulting AAA phantom was compared with the average thickness of the corresponding patient-specific virtual model, obtaining an average dimensional mismatch of 180 μm (11.14%). The material characterization of the artery was determined from uniaxial tensile tests as various combinations of polyurethane resins were chosen due to their similarity with ex vivo AAA mechanical behavior in the physiological stress configuration. The proposed methodology yields AAA phantoms with nonuniform wall thickness using a fast and low-cost process. These replicas may be used in benchtop experiments to validate deformations obtained with numerical simulations using finite element analysis, or to validate optical methods developed to image ex vivo arterial deformations during pressure-inflation testing.


Author(s):  
Amirhossein Arzani ◽  
Shawn C. Shadden

Abdominal aortic aneurysms (AAA) are characterized by disturbed flow patterns, low and oscillatory wall shear stress with high gradients, increased particle residence time, and mild turbulence. Diameter is the most common metric for rupture prediction, although this metric can be unreliable. We hypothesize that understanding the flow topology and mixing inside AAA could provide useful insight into mechanisms of aneurysm growth. AAA morphology has high variability, as with AAA hemodynamics, and therefore we consider patient-specific analyses over several small to medium sized AAAs. Vortical patterns dominate AAA hemodynamics and traditional analyses based on the Eulerian fields (e.g. velocity) fail to convey the complex flow structures. The computation of finite-time Lyapunov exponent (FTLE) fields and underlying Lagrangian coherent structures (LCS) help reveal a Lagrangian template for quantifying the flow [1].


Author(s):  
David M. Pierce ◽  
Thomas E. Fastl ◽  
Hannah Weisbecker ◽  
Gerhard A. Holzapfel ◽  
Borja Rodriguez-Vila ◽  
...  

Through progress in medical imaging, image analysis and finite element (FE) meshing tools it is now possible to extract patient-specific geometries from medical images of, e.g., abdominal aortic aneurysms (AAAs), and thus to study clinically relevant problems via FE simulations. Medical imaging is most often performed in vivo, and hence the reconstructed model geometry in the problem of interest will represent the in vivo state, e.g., the AAA at physiological blood pressure. However, classical continuum mechanics and FE methods assume that constitutive models and the corresponding simulations start from an unloaded, stress-free reference condition.


2010 ◽  
Vol 43 (7) ◽  
pp. 1408-1416 ◽  
Author(s):  
Barry J. Doyle ◽  
Aidan J. Cloonan ◽  
Michael T. Walsh ◽  
David A. Vorp ◽  
Timothy M. McGloughlin

2021 ◽  
Vol 7 ◽  
Author(s):  
Jasamine Coles-Black ◽  
Damien Bolton ◽  
Jason Chuen

Introduction: 3D printed patient-specific vascular phantoms provide superior anatomical insights for simulating complex endovascular procedures. Currently, lack of exposure to the technology poses a barrier for adoption. We offer an accessible, low-cost guide to producing vascular anatomical models using routine CT angiography, open source software packages and a variety of 3D printing technologies.Methods: Although applicable to all vascular territories, we illustrate our methodology using Abdominal Aortic Aneurysms (AAAs) due to the strong interest in this area. CT aortograms acquired as part of routine care were converted to representative patient-specific 3D models, and then printed using a variety of 3D printing technologies to assess their material suitability as aortic phantoms. Depending on the technology, phantoms cost $20–$1,000 and were produced in 12–48 h. This technique was used to generate hollow 3D printed thoracoabdominal aortas visible under fluoroscopy.Results: 3D printed AAA phantoms were a valuable addition to standard CT angiogram reconstructions in the simulation of complex cases, such as short or very angulated necks, or for positioning fenestrations in juxtarenal aneurysms. Hollow flexible models were particularly useful for device selection and in planning of fenestrated EVAR. In addition, these models have demonstrated utility other settings, such as patient education and engagement, and trainee and anatomical education. Further study is required to establish a material with optimal cost, haptic and fluoroscopic fidelity.Conclusion: We share our experiences and methodology for developing inexpensive 3D printed vascular phantoms which despite material limitations, successfully mimic the procedural challenges encountered during live endovascular surgery. As the technology continues to improve, 3D printed vascular phantoms have the potential to disrupt how endovascular procedures are planned and taught.


2019 ◽  
Vol 24 (1) ◽  
pp. 30 ◽  
Author(s):  
Shadi Alameddin ◽  
Amélie Fau ◽  
David Néron ◽  
Pierre Ladevèze ◽  
Udo Nackenhorst

The solution of structural problems with nonlinear material behaviour in a model order reduction framework is investigated in this paper. In such a framework, greedy algorithms or adaptive strategies are interesting as they adjust the reduced order basis (ROB) to the problem of interest. However, these greedy strategies may lead to an excessive increase in the size of the ROB, i.e., the solution is no more represented in its optimal low-dimensional expansion. Here, an optimised strategy is proposed to maintain, at each step of the greedy algorithm, the lowest dimension of a Proper Generalized Decomposition (PGD) basis using a randomised Singular Value Decomposition (SVD) algorithm. Comparing to conventional approaches such as Gram–Schmidt orthonormalisation or deterministic SVD, it is shown to be very efficient both in terms of numerical cost and optimality of the ROB. Examples with different mesh densities are investigated to demonstrate the numerical efficiency of the presented method.


2018 ◽  
Vol 226 ◽  
pp. 04036
Author(s):  
Yuriy M. Manatskov ◽  
Torsten Bertram ◽  
Danil V. Shaykhutdinov ◽  
Nikolay I. Gorbatenko

Complex dynamic linear systems of equations are solved by numerical iterative methods, which need much computation and are timeconsuming ones, and the optimization stage requires repeated solution of these equation systems that increases the time on development. To shorten the computation time, various methods can be applied, among them preliminary (estimated) calculation or oversimple models calculation, however, while testing and optimizing the full model is used. Reduced order models are very popular in solving this problem. The main idea of a reduced order model is to find a simplified model that may reflect the required properties of the original model as accurately as possible. There are many methods for the model order reduction, which have their advantages and disadvantages. In this article, a method based on Krylov subspaces and SVD methods is considered. A numerical experiments is given.


Sign in / Sign up

Export Citation Format

Share Document