Predictive Models with Patient Specific Material Properties for the Biomechanical Behavior of Ascending Thoracic Aneurysms

Significance of Dynamic Axial Stretching on Estimating Biomechanical Behavior and Properties of the Human Ascending Aorta

Contrary to most vessels, the ascending thoracic aorta (ATA) not only distends but also elongates in the axial direction. The purpose of this study is to investigate the biomechanical behavior of the ascending thoracic aorta (ATA) in response to dynamic axial stretching during the cardiac cycle. In addition, the implications of neglecting this dynamic axial stretching when estimating the constitutive model parameters of the ATA are investigated. The investigations were performed through in silico simulations by assuming a Gasser-Ogden-Holzapfel (GOH) constitutive model representative of ATA tissue material. The GOH model parameters were obtained from biaxial tests performed on four human ATA tissues in a previous study. Pressure-diameter curves were simulated as synthetic data to assess the effect of neglecting dynamic axial stretching on estimating constitutive model parameters. Our findings reveal a significant increase in axial stress (~ 16%) and stored strain energy (~ 18%) in the vessel when dynamic axial stretching is considered, as opposed to assuming a fixed axial stretch. All but one artery showed increased volume compliance while considering a dynamic axial stretching condition. Furthermore, we observe a notable difference in the estimated constitutive model parameters when dynamic axial stretching of the ATA is neglected, compared to the ground truth model parameters. These results underscore the critical importance of accounting for axial deformations when conducting in vivo biomechanical characterization of the ascending thoracic aorta.

Towards the mechanical characterisation of unruptured intracranial aneurysms: Numerical modelling of interactions between a deformation device and the aneurysm wall

Intracranial aneurysm is a critical pathology related to the arterial wall deterioration. This work is an essential aspect of a large scale project aimed at providing clinicians with a non-invasive patient-specific decision support tool regarding the rupture risk assessment. A machine learning algorithm links the aneurysm shape observed and a database of UIA clinical images associated with in vivo wall mechanical properties and rupture characterisation. The database constitution is derived from a device prototype coupled with medical imaging. It provides the mechanical characterisation of the aneurysm from the wall deformation obtained by inverse analysis based on the variation of luminal volume. Before performing in vivo tests of the device on small animals, a numerical model was built to quantify the device's impact on the aneurysm wall under natural blood flow conditions. As the clinician will never be able to precisely situate the device, several locations were considered. In preparation for the inverse analysis procedure, artery material laws of increasing complexity were studied (linear elastic, hyper elastic Fung-like). Considering all the device locations and material laws, the device induced relative displacements to the Systole peak (worst case scenario with the highest mechanical stimulus linked to the blood flow) ranging from 375 μm to 1.28 mm. The variation of luminal volume associated with the displacements was between 0.95 % and 4.3 % compared to the initial Systole volume of the aneurysm. Significant increase of the relative displacements and volume variations were found with the study of different cardiac cycle moments between the blood flow alone and the device application. For forthcoming animal model studies, Spectral Photon CT Counting, with a minimum spatial resolution of 250 μm, was selected as the clinical imaging technique. Based on this preliminary study, the displacements and associated volume variations (baseline for inverse analyse), should be observable and exploitable.

In vivo mechanical characterization of arterial wall using an inverse analysis procedure: application on an animal model of intracranial aneurysm

Intracranial aneurysm is a pathology related to the deterioration of the arterial wall. This work is an essential part of a large-scale project aimed at providing clinicians with a non-invasive patient-specific decision support tool to facilitate the rupture risk assessment. It will lean on the link between the aneurysm shape clinically observed and a database derived from the in vivo mechanical characterization of aneurysms. To supply this database, a deformation device prototype of the arterial wall was developed. Its use coupled with medical imaging (spectral photon-counting computed tomography providing a spatial resolution down to 250 μm) is used to determine the in vivo mechanical properties of the wall based on the inverse analysis of the quantification of the wall deformation observed experimentally. This study presents the in vivo application of this original procedure to an animal model of aneurysm. The mechanical properties of the aneurysm wall identified were consistent with the literature, and the errors between the numerical and experimental results were less than 10%. Based on these parameters, this study allows the assessment of the aneurysm stress state for a known solicitation and points towards the definition of a rupture criterion.

On the Material Constitutive Behavior of the Aortic Root in Patients with Transcatheter Aortic Valve Implantation

BACKGROUND

Transcatheter aortic valve implantation (TAVI) is a minimally invasive procedure used to treat patients with severe aortic valve stenosis. However, there is limited knowledge on the material properties of the aortic root in TAVI patients, and this can impact the credibility of computer simulations. This study aimed to develop a non-invasive inverse approach for estimating reliable material constituents for the aortic root and calcified valve leaflets in patients undergoing TAVI.

METHODS

The identification of material parameters is based on the simultaneous minimization of two cost functions, which define the difference between model predictions and cardiac-gated CT measurements of the aortic wall and valve orifice area. Validation of the inverse analysis output was performed comparing the numerical predictions with actual CT shapes and post-TAVI measures of implanted device diameter.

RESULTS

A good agreement of the peak systolic shape of the aortic wall was found between simulations and imaging, with similarity index in the range in the range of 83.7% to 91.5% for n.20 patients. Not any statistical difference was observed between predictions and CT measures of orifice area for the stenotic aortic valve. After TAVI simulations, the measurements of SAPIEN 3 Ultra (S3) device diameter were in agreement with those from post-TAVI angio-CT imaging. A sensitivity analysis demonstrated a modest impact on the S3 diameters when altering the elastic material property of the aortic wall in the range of inverse analysis solution.

CONCLUSIONS

Overall, this study demonstrates the feasibility and potential benefits of using non-invasive imaging techniques and computational modeling to estimate material properties in patients undergoing TAVI.

Image-Based Finite Element Modeling Approach for Characterizing In Vivo Mechanical Properties of Human Arteries

Mechanical properties of the arterial walls could provide meaningful information for the diagnosis, management and treatment of cardiovascular diseases. Classically, various experimental approaches were conducted on dissected arterial tissues to obtain their stress-stretch relationship, which has limited value clinically. Therefore, there is a pressing need to obtain biomechanical behaviors of these vascular tissues in vivo for personalized treatment. This paper reviews the methods to quantify arterial mechanical properties in vivo. Among these methods, we emphasize a novel approach using image-based finite element models to iteratively determine the material properties of the arterial tissues. This approach has been successfully applied to arterial walls in various vascular beds. The mechanical properties obtained from the in vivo approach were compared to those from ex vivo experimental studies to investigate whether any discrepancy in material properties exists for both approaches. Arterial tissue stiffness values from in vivo studies generally were in the same magnitude as those from ex vivo studies, but with lower average values. Some methodological issues, including solution uniqueness and robustness; method validation; and model assumptions and limitations were discussed. Clinical applications of this approach were also addressed to highlight their potential in translation from research tools to cardiovascular disease management.

Ascending aorta mechanics and dimensions in aortopathy – from science to application

The ascending aorta has a unique microstructure and biomechanical properties that allow it to absorb energy during systole and return energy during diastole (Windkessel effect). Derangements in aortic architecture can result in changes to biomechanics and inefficiencies in function. Ultimately biomechanical failure may occur resulting in aortic dissection or rupture. By measuring aortic biomechanics with either in vivo or ex vivo methods, one may be able to predict tissue failure in patients with aortic disease such as aneurysms. An understanding of the biomechanical changes that lead to these tissue-level failures may help guide therapy, disease surveillance, surgical intervention, and aid in the development of new treatments for this deadly condition.

A population-specific material model for sagittal craniosynostosis to predict surgical shape outcomes

Sagittal craniosynostosis consists of premature fusion (ossification) of the sagittal suture during infancy, resulting in head deformity and brain growth restriction. Spring-assisted cranioplasty (SAC) entails skull incisions to free the fused suture and insertion of two springs (metallic distractors) to promote cranial reshaping. Although safe and effective, SAC outcomes remain uncertain. We aimed hereby to obtain and validate a skull material model for SAC outcome prediction. Computed tomography data relative to 18 patients were processed to simulate surgical cuts and spring location. A rescaling model for age matching was created using retrospective data and validated. Design of experiments was used to assess the effect of different material property parameters on the model output. Subsequent material optimization-using retrospective clinical spring measurements-was performed for nine patients. A population-derived material model was obtained and applied to the whole population. Results showed that bone Young's modulus and relaxation modulus had the largest effect on the model predictions: the use of the population-derived material model had a negligible effect on improving the prediction of on-table opening while significantly improved the prediction of spring kinematics at follow-up. The model was validated using on-table 3D scans for nine patients: the predicted head shape approximated within 2 mm the 3D scan model in 80% of the surface points, in 8 out of 9 patients. The accuracy and reliability of the developed computational model of SAC were increased using population data: this tool is now ready for prospective clinical application.

Identification of in vivo nonlinear anisotropic mechanical properties of ascending thoracic aortic aneurysm from patient-specific CT scans

Accurate identification of in vivo nonlinear, anisotropic mechanical properties of the aortic wall of individual patients remains to be one of the critical challenges in the field of cardiovascular biomechanics. Since only the physiologically loaded states of the aorta are given from in vivo clinical images, inverse approaches, which take into account of the unloaded configuration, are needed for in vivo material parameter identification. Existing inverse methods are computationally expensive, which take days to weeks to complete for a single patient, inhibiting fast feedback for clinicians. Moreover, the current inverse methods have only been evaluated using synthetic data. In this study, we improved our recently developed multi-resolution direct search (MRDS) approach and the computation time cost was reduced to 1~2 hours. Using the improved MRDS approach, we estimated in vivo aortic tissue elastic properties of two ascending thoracic aortic aneurysm (ATAA) patients from pre-operative gated CT scans. For comparison, corresponding surgically-resected aortic wall tissue samples were obtained and subjected to planar biaxial tests. Relatively close matches were achieved for the in vivo-identified and ex vivo-fitted stress-stretch responses. It is hoped that further development of this inverse approach can enable an accurate identification of the in vivo material parameters from in vivo image data.

Computational modeling of bicuspid aortopathy: Towards personalized risk strategies

This paper describes current advances on the application of in-silico for the understanding of bicuspid aortopathy and future perspectives of this technology on routine clinical care. This includes the impact that artificial intelligence can provide to develop computer-based clinical decision support system and that wearable sensors can offer to remotely monitor high-risk bicuspid aortic valve (BAV) patients. First, we discussed the benefit of computational modeling by providing tangible examples of in-silico software products based on computational fluid-dynamic (CFD) and finite-element method (FEM) that are currently transforming the way we diagnose and treat cardiovascular diseases. Then, we presented recent findings on computational hemodynamic and structural mechanics of BAV to highlight the potentiality of patient-specific metrics (not-based on aortic size) to support the clinical-decision making process of BAV-associated aneurysms. Examples of BAV-related personalized healthcare solutions are illustrated.

Computational predictions of damage propagation preceding dissection of ascending thoracic aortic aneurysms

Dissections of ascending thoracic aortic aneurysms (ATAAs) cause significant morbidity and mortality worldwide. They occur when a tear in the intima-media of the aorta permits the penetration of the blood and the subsequent delamination and separation of the wall in 2 layers, forming a false channel. To predict computationally the risk of tear formation, stress analyses should be performed layer-specifically and they should consider internal or residual stresses that exist in the tissue. In the present paper, we propose a novel layer-specific damage model based on the constrained mixture theory, which intrinsically takes into account these internal stresses and can predict appropriately the tear formation. The model is implemented in finite-element commercial software Abaqus coupled with user material subroutine. Its capability is tested by applying it to the simulation of different exemplary situations, going from in vitro bulge inflation experiments on aortic samples to in vivo overpressurizing of patient-specific ATAAs. The simulations reveal that damage correctly starts from the intimal layer (luminal side) and propagates across the media as a tear but never hits the adventitia. This scenario is typically the first stage of development of an acute dissection, which is predicted for pressures of about 2.5 times the diastolic pressure by the model after calibrating the parameters against experimental data performed on collected ATAA samples. Further validations on a larger cohort of patients should hopefully confirm the potential of the model in predicting patient-specific damage evolution and possible risk of dissection during aneurysm growth for clinical applications.

Estimation of in vivo mechanical properties of the aortic wall: A multi-resolution direct search approach

The patient-specific biomechanical analysis of the aorta requires in vivo mechanical properties of individual patients. Existing approaches for estimating in vivo material properties often demand high computational cost and mesh correspondence of the aortic wall between different cardiac phases. In this paper, we propose a novel multi-resolution direct search (MRDS) approach for estimation of the nonlinear, anisotropic constitutive parameters of the aortic wall. Based on the finite element (FE) updating scheme, the MRDS approach consists of the following three steps: (1) representing constitutive parameters with multiple resolutions using principal component analysis (PCA), (2) building links between the discretized PCA spaces at different resolutions, and (3) searching the PCA spaces in a 'coarse to fine' fashion following the links. The estimation of material parameters is achieved by minimizing a node-to-surface error function, which does not need mesh correspondence. The method was validated through a numerical experiment by using the in vivo data from a patient with ascending thoracic aortic aneurysm (ATAA), the results show that the number of FE iterations was significantly reduced compared to previous methods. The approach was also applied to the in vivo CT data from an aged healthy human patient, and using the estimated material parameters, the FE-computed geometry was well matched with the image-derived geometry. This novel MRDS approach may facilitate the personalized biomechanical analysis of aortic tissues, such as the rupture risk analysis of ATAA, which requires fast feedback to clinicians.

A new inverse method for estimation of in vivo mechanical properties of the aortic wall

The aortic wall is always loaded in vivo, which makes it challenging to estimate the material parameters of its nonlinear, anisotropic constitutive equation from in vivo image data. Previous approaches largely relied on either computationally expensive finite element models or simplifications of the geometry or material models. In this study, we investigated a new inverse method based on aortic wall stress computation. This approach consists of the following two steps: (1) computing an "almost true" stress field from the in vivo geometries and loading conditions, (2) building an objective function based on the "almost true" stress fields, constitutive equations and deformation relations, and estimating the material parameters by minimizing the objective function. The method was validated through numerical experiments by using the in vivo data from four ascending aortic aneurysm (AsAA) patients. The results demonstrated that the method is computationally efficient. This novel approach may facilitate the personalized biomechanical analysis of aortic tissues in clinical applications, such as in the rupture risk analysis of ascending aortic aneurysms.

Biaxial rupture properties of ascending thoracic aortic aneurysms

Although hundreds of samples obtained from ascending thoracic aortic aneurysms (ATAA) of patients undergoing elective surgical repair have already been characterized biomechanically, their rupture properties were always derived from uniaxial tensile tests. Due to their bulge shape, ATAAs are stretched biaxially in vivo. In order to understand the biaxial rupture of ATAAs, our group developed a novel methodology based on bulge inflation and full-field optical measurements. The objective of the current paper is threefold. Firstly, we will review the failure properties (maximum stress, maximum stretch) obtained by bulge inflation testing on a cohort of 31 patients and compare them with failure properties obtained by uniaxial tension in a previously published study. Secondly, we will investigate the relationship between the failure properties and the age of patients, showing that patients below 55years of age display significantly higher strength. Thirdly, we will define a rupture risk based on the extensibility of the tissue and we will show that this rupture risk is strongly correlated with the physiological elastic modulus of the tissue independently of the age, ATAA diameter or the aortic valve phenotype of the patient.

STATEMENT OF SIGNIFICANCE

Despite their medical importance, rupture properties of ascending thoracic aortic aneurysms (ATAA) subjected to biaxial tension were inexistent in the literature. In order to address this lack, our group developed a novel methodology based on bulge inflation and full-field optical measurements. Here we report rupture properties obtained with this methodology on 31 patients. It is shown for the first time that rupture occurs when the stretch applied to ATAAs reaches the maximum extensibility of the tissue and that this maximum extensibility correlates strongly with the elastic properties. The outcome is a better detection of at-risk individuals for elective surgical repair.

From Finite Element Meshes to Clouds of Points: A Review of Methods for Generation of Computational Biomechanics Models for Patient-Specific Applications

It has been envisaged that advances in computing and engineering technologies could extend surgeons' ability to plan and carry out surgical interventions more accurately and with less trauma. The progress in this area depends crucially on the ability to create robustly and rapidly patient-specific biomechanical models. We focus on methods for generation of patient-specific computational grids used for solving partial differential equations governing the mechanics of the body organs. We review state-of-the-art in this area and provide suggestions for future research. To provide a complete picture of the field of patient-specific model generation, we also discuss methods for identifying and assigning patient-specific material properties of tissues and boundary conditions.