Correlation between micro and macrostructural biaxial behavior of ascending thoracic aneurysm: a novel experimental technique

Numerical investigation on circular and elliptical bulge tests for inverse soft tissue characterization

The acquisition of insights concerning the mechanobiology of aneurysmatic aortic tissues is an important field of investigation. The complete characterization of aneurysm mechanical behaviour can be carried out by biaxial experimental tests on ex vivo specimens. In literature, several works proposed bulge inflation tests as a valid method to analyse aneurysmatic tissue. Bulge test data processing requires the adoption of digital image correlation and inverse analysis approaches to estimate strain and stress distributions, respectively. In this context, however, the accuracy of inverse analysis method has not been evaluated yet. This aspect appears particularly interesting given the anisotropic behaviour of the soft tissue and the possibility to adopt different die geometries. The goal of this study is to provide an accuracy characterization of the inverse analysis applied to the bulge test technique using a numerical approach. In particular, different cases of bulge inflation were simulated in a finite element environment as a reference. To investigate the effect of tissue anisotropic degree and bulge die geometries (circular and elliptical), different input parameters were considered to obtain multiple test cases. The specimen deformed shapes, resulting from the reference finite element simulations, were then analysed through an inverse analysis approach to produce an estimation of stress distributions. The estimated stresses were, at last, compared with the values from the reference finite element simulations. The results demonstrated that the circular die geometry produces a satisfactory estimation accuracy only under certain conditions of material quasi-isotropy. On the other hand, the choice of an elliptical bulge die was proven to be more suitable for the analysis of anisotropic tissues.

A review on the biomechanical behaviour of the aorta

Large aortic aneurysm and acute and chronic aortic dissection are pathologies of the aorta requiring surgery. Recent advances in medical intervention have improved patient outcomes; however, a clear understanding of the mechanisms leading to aortic failure and, hence, a better understanding of failure risk, is still missing. Biomechanical analysis of the aorta could provide insights into the development and progression of aortic abnormalities, giving clinicians a powerful tool in risk stratification. The complexity of the aortic system presents significant challenges for a biomechanical study and requires various approaches to analyse the aorta. To address this, here we present a holistic review of the biomechanical studies of the aorta by categorising articles into four broad approaches, namely theoretical, in vivo, experimental and combined investigations. Experimental studies that focus on identifying mechanical properties of the aortic tissue are also included. By reviewing the literature and discussing drawbacks, limitations and future challenges in each area, we hope to present a more complete picture of the state-of-the-art of aortic biomechanics to stimulate research on critical topics. Combining experimental modalities and computational approaches could lead to more comprehensive results in risk prediction for the aortic system.

Uncertainty Quantification in the In Vivo Image-Based Estimation of Local Elastic Properties of Vascular Walls

Introduction: Patient-specific computational models are a powerful tool for planning cardiovascular interventions. However, the in vivo patient-specific mechanical properties of vessels represent a major source of uncertainty. In this study, we investigated the effect of uncertainty in the elastic module (E) on a Fluid–Structure Interaction (FSI) model of a patient-specific aorta. Methods: The image-based χ-method was used to compute the initial E value of the vascular wall. The uncertainty quantification was carried out using the generalized Polynomial Chaos (gPC) expansion technique. The stochastic analysis was based on four deterministic simulations considering four quadrature points. A deviation of about ±20% on the estimation of the E value was assumed. Results: The influence of the uncertain E parameter was evaluated along the cardiac cycle on area and flow variations extracted from five cross-sections of the aortic FSI model. Results of stochastic analysis showed the impact of E in the ascending aorta while an insignificant effect was observed in the descending tract. Conclusions: This study demonstrated the importance of the image-based methodology for inferring E, highlighting the feasibility of retrieving useful additional data and enhancing the reliability of in silico models in clinical practice.

An image-based approach for the estimation of arterial local stiffness in vivo

The analysis of mechanobiology of arterial tissues remains an important topic of research for cardiovascular pathologies evaluation. In the current state of the art, the gold standard to characterize the tissue mechanical behavior is represented by experimental tests, requiring the harvesting of ex-vivo specimens. In recent years though, image-based techniques for the in vivo estimation of arterial tissue stiffness were presented. The aim of this study is to define a new approach to provide local distribution of arterial stiffness, estimated as the linearized Young's Modulus, based on the knowledge of in vivo patient-specific imaging data. In particular, the strain and stress are estimated with sectional contour length ratios and a Laplace hypothesis/inverse engineering approach, respectively, and then used to calculate the Young's Modulus. After describing the method, this was validated by using a set of Finite Element simulations as input. In particular, idealized cylinder and elbow shapes plus a single patient-specific geometry were simulated. Different stiffness distributions were tested for the simulated patient-specific case. After the validation from Finite Element data, the method was then applied to patient-specific ECG-gated Computed Tomography data by also introducing a mesh morphing approach to map the aortic surface along the cardiac phases. The validation process revealed satisfactory results. In the simulated patient-specific case, root mean square percentage errors below 10% for the homogeneous distribution and below 20% for proximal/distal distribution of stiffness. The method was then successfully used on the three ECG-gated patient-specific cases. The resulting distributions of stiffness exhibited significant heterogeneity, nevertheless the resulting Young's moduli were always contained within the 1-3 MPa range, which is in line with literature.

Introduction of a Novel Image-Based and Non-Invasive Method for the Estimation of Local Elastic Properties of Great Vessels

Background: In the context of a growing demand for the use of in silico models to meet clinical requests, image-based methods play a crucial role. In this study, we present a parametric equation able to estimate the elasticity of vessel walls, non-invasively and indirectly, from information uniquely retrievable from imaging. Methods: A custom equation was iteratively refined and tuned from the simulations of a wide range of different vessel models, leading to the definition of an indirect method able to estimate the elastic modulus E of a vessel wall. To test the effectiveness of the predictive capability to infer the E value, two models with increasing complexity were used: a U-shaped vessel and a patient-specific aorta. Results: The original formulation was demonstrated to deviate from the ground truth, with a difference of 89.6%. However, the adoption of our proposed equation was found to significantly increase the reliability of the estimated E value for a vessel wall, with a mean percentage error of 9.3% with respect to the reference values. Conclusion: This study provides a strong basis for the definition of a method able to estimate local mechanical information of vessels from data easily retrievable from imaging, thus potentially increasing the reliability of in silico cardiovascular models.

On the Role and Effects of Uncertainties in Cardiovascular in silico Analyses

The assessment of cardiovascular hemodynamics with computational techniques is establishing its fundamental contribution within the world of modern clinics. Great research interest was focused on the aortic vessel. The study of aortic flow, pressure, and stresses is at the basis of the understanding of complex pathologies such as aneurysms. Nevertheless, the computational approaches are still affected by sources of errors and uncertainties. These phenomena occur at different levels of the computational analysis, and they also strongly depend on the type of approach adopted. With the current study, the effect of error sources was characterized for an aortic case. In particular, the geometry of a patient-specific aorta structure was segmented at different phases of a cardiac cycle to be adopted in a computational analysis. Different levels of surface smoothing were imposed to define their influence on the numerical results. After this, three different simulation methods were imposed on the same geometry: a rigid wall computational fluid dynamics (CFD), a moving-wall CFD based on radial basis functions (RBF) CFD, and a fluid-structure interaction (FSI) simulation. The differences of the implemented methods were defined in terms of wall shear stress (WSS) analysis. In particular, for all the cases reported, the systolic WSS and the time-averaged WSS (TAWSS) were defined.

Fully-Coupled FSI Computational Analyses in the Ascending Thoracic Aorta Using Patient-Specific Conditions and Anisotropic Material Properties

Computational hemodynamics has become increasingly important within the context of precision medicine, providing major insight in cardiovascular pathologies. However, finding appropriate compromise between speed and accuracy remains challenging in computational hemodynamics for an extensive use in decision making. For example, in the ascending thoracic aorta, interactions between the blood and the aortic wall must be taken into account for the sake of accuracy, but these fluid structure interactions (FSI) induce significant computational costs, especially when the tissue exhibits a hyperelastic and anisotropic response. The objective of the current study is to use the Small On Large (SOL) theory to linearize the anisotropic hyperelastic behavior in order to propose a reduced-order model for FSI simulations of the aorta. The SOL method is tested for fully-coupled FSI simulations in a patient-specific aortic geometry presenting an Ascending Thoracic Aortic Aneurysm (aTAA). The same model is also simulated with a fully-coupled FSI with non-linear material behavior, without SOL linearization. Eventually, the results and computational times with and without the SOL are compared. The SOL approach is demonstrated to provide a significant reduction of computational costs for FSI analysis in the aTAA, and the results in terms of stress state distribution are comparable. The method is implemented in ANSYS and will be further evaluated for clinical applications.