Endovascular device design in the future: transformation from trial and error to computational design

Computational evaluation of inferior vena cava filters through computational fluid dynamics methods

Numerical simulation is growing in its importance toward the design, testing and evaluation of medical devices. Computational fluid dynamics and finite element analysis allow improved calculation of stress, heat transfer, and flow to better understand the medical device environment. Current research focuses not only on improving medical devices, but also on improving the computational tools themselves. As methods and computer technology allow for faster simulation times, iterations and trials can be performed faster to collect more data. Given the adverse events associated with long-term inferior vena cava (IVC) filter placement, IVC filter design and device evaluation are of paramount importance. This work reviews computational methods used to develop, test, and improve IVC filters to ultimately serve the needs of the patient.

Endograft Specific Haemodynamics After Endovascular Aneurysm Repair: Flow Characteristics of Four Stent Graft Systems

OBJECTIVES

The implication of haemodynamics in the occurrence of complications after endovascular aneurysm repair (EVAR) has been raised in the literature. Different aortic stent graft configurations may lead to different haemodynamic properties. The current study deals with the post-operative haemodynamic variability between four stent graft systems with different structure, material, and type of fixation.

METHODS

Computed tomography data of 32 patients were used, equally distributed among the four endograft groups, namely the AFX, Endurant, Excluder, and Nellix. Velocity, wall shear stress (WSS), and helicity statistics were calculated, in regions around the flow division where disturbances are expected. The haemodynamic data were compared between and within the groups.

RESULTS

The morphology of AAAs pre-operatively did not vary significantly among the four groups. Before the flow division, lowest velocity was observed in Endurant cases and highest in Nellix cases. Endurant induced the lowest peak WSS and Nellix the highest (p = .03). The helicity levels were low in AFX and Nellix cases and high in Endurant and Excluder cases. After the flow division, the trend in the results was preserved. Nellix induced the highest velocity and WSS, followed closely by Excluder and AFX. There was a significant increase of helicity before and after flow division in AFX (p <0.001, R² = 0.09) and Nellix (p <0.001) cases.

CONCLUSIONS

It has been shown that different types of endografts induce variable haemodynamic conditions around the flow division. The parallel limb structure, featured by Nellix, seems to induce favourable flow conditions in terms of velocity and WSS, while helical flow before the flow division is suppressed. High WSS is generally considered to be a desirable flow characteristic in endovascular devices, whereas helicity extremes (very low or high) are potentially a negative sign. Endurant, with the stiffer material and the short neck structure, was associated with the lowest blood velocity and WSS values but preserved high helicity levels. The AFX and Excluder, which include the same material, induced similar haemodynamic conditions.

Haemodynamic performance of AFX and Nellix endografts: a computational fluid dynamics study

OBJECTIVES

The objective of this study is to analyse the flow conditions in the AFX and Nellix endografts (EGs) accounting for their postimplantation configuration in patients with an endovascular aneurysm repair-treated abdominal aortic aneurysm.

METHODS

We reconstructed post-endovascular aneurysm repair computed tomography scans of patients treated with an AFX or Nellix EG creating post-implantation EG models. We examined 16 patients, 8 in each group. The blood flow properties were obtained by computational fluid dynamics simulations and were subsequently compared with physiological infrarenal blood flow properties measured in 5 healthy subjects. Specifically, pressure drop, maximum velocity and wall shear stress were measured at peak systole and mean helicity at mid-diastole.

RESULTS

Our statistical analyses showed that the haemodynamic properties in both control regions did not vary statistically after the implantation of either the AFX or the Nellix EG, except for helicity that was significantly lower in the abdominal part of the Nellix EG compared with the expected physiological measurement. Regardless of the overall blood flow restoration, it is important to note that low pressure drop was detected along the limbs of the AFX and suppressed blood helical motion was detected at the entrance of the Nellix device.

CONCLUSIONS

It is observed from the results that the AFX EG has achieved absolute restoration of blood flow after endovascular aneurysm repair, although the development of secondary flow in the upper part of the EG and the low pressure drop in its limbs should be acknowledged. The Nellix EG also seems to be haemodynamically efficient. However, the suppression of helical flow before blood enters the device might raise concerns about its clinical application.

Abdominal aortic aneurysm endovascular repair: profiling post-implantation morphometry and hemodynamics with image-based computational fluid dynamics

Endovascular aneurysm repair (EVAR) has disseminated rapidly as an alternative to open surgical repair for the treatment of abdominal aortic aneurysms (AAAs), because of its reduced invasiveness, low mortality and morbidity rate. The effectiveness of the endovascular devices used in EVAR is always at question as postoperative adverse events can lead to re-intervention or to a possible fatal scenario for the circulatory system. Motivated by the assessment of the risks related to thrombus formation, here the impact of two different commercial endovascular grafts on local hemodynamics is explored through 20 image-based computational hemodynamic models of EVAR-treated patients (N=10 per each endograft model). Hemodynamic features, susceptible to promote thrombus formation, such as flow separation and recirculation, are quantitatively assessed and compared with the local hemodynamics established in image-based infrarenal abdominal aortic models of healthy subjects (N=10). The hemodynamic analysis is complemented by a geometrical characterization of the EVAR-induced reshaping of the infrarenal abdominal aortic vascular region. The findings of this study indicate that: (1) the clinically observed propensity to thrombus formation in devices used in EVAR strategies can be explained in terms of local hemodynamics by means of image-based computational hemodynamics approach; (2) reportedly pro-thrombotic hemodynamic structures are strongly correlated with the geometry of the aortoiliac tract postoperatively. In perspective, our study suggests that future clinical follow up studies could include a geometric analysis of the region of the implant, monitoring shape variations that can lead to hemodynamic disturbances of clinical significance.

A framework for designing patient-specific bioprosthetic heart valves using immersogeometric fluid-structure interaction analysis

Numerous studies have suggested that medical image derived computational mechanics models could be developed to reduce mortality and morbidity due to cardiovascular diseases by allowing for patient-specific surgical planning and customized medical device design. In this work, we present a novel framework for designing prosthetic heart valves using a parametric design platform and immersogeometric fluid-structure interaction (FSI) analysis. We parameterize the leaflet geometry using several key design parameters. This allows for generating various perturbations of the leaflet design for the patient-specific aortic root reconstructed from the medical image data. Each design is analyzed using our hybrid arbitrary Lagrangian-Eulerian/immersogeometric FSI methodology, which allows us to efficiently simulate the coupling of the deforming aortic root, the parametrically designed prosthetic valves, and the surrounding blood flow under physiological conditions. A parametric study is performed to investigate the influence of the geometry on heart valve performance, indicated by the effective orifice area and the coaptation area. Finally, the FSI simulation result of a design that balances effective orifice area and coaptation area reasonably well is compared with patient-specific phase contrast magnetic resonance imaging data to demonstrate the qualitative similarity of the flow patterns in the ascending aorta.

Hemodynamic and Anatomic Predictors of Renovisceral Stent-Graft Occlusion Following Chimney Endovascular Repair of Juxtarenal Aortic Aneurysms

PURPOSE

To identify anatomic and hemodynamic changes associated with impending visceral chimney stent-graft occlusion after endovascular aneurysm repair (EVAR) with the chimney technique (chEVAR).

METHODS

A retrospective evaluation was performed of computed tomography scans from 41 patients who underwent juxtarenal chEVAR from 2008 to 2012 to identify stent-grafts demonstrating conformational changes following initial placement. Six subjects (mean age 74 years; 3 men) were selected for detailed reconstruction and computational hemodynamic analysis; 4 had at least 1 occluded chimney stent-graft. This subset of repairs was systematically analyzed to define the anatomic and hemodynamic impact of these changes and identify signature patterns associated with impending renovisceral stent-graft occlusion. Spatial and temporal analyses of cross-sectional area, centerline angle, intraluminal pressure, and wall shear stress (WSS) were performed within the superior mesenteric and renal artery chimney grafts used for repair.

RESULTS

Conformational changes in the chimney stent-grafts and associated perturbations, in both local WSS and pressure, were responsible for the 5 occlusions in the 13 stented branches. Anatomic and hemodynamic signatures leading to occlusion were identified within 1 month postoperatively, with a lumen area <14 mm² (p=0.04), systolic pressure gradient >25 Pa/mm (p=0.03), and systolic WSS >45 Pa (p=0.03) associated with future chimney stent-graft occlusion.

CONCLUSION

Chimney stent-grafts at increased risk for occlusion demonstrated anatomic and hemodynamic signatures within 1 month of juxtarenal chEVAR. Analysis of these parameters in the early postoperative period may be useful for identifying and remediating these high-risk stent-grafts.

Hemodynamic Profile of Two Aortic Endografts Accounting for Their Postimplantation Position

Endovascular aneurysm repair (EVAR) is a clinically effective technique for treating anatomically eligible abdominal aortic aneurysms (AAAs), involving the deployment of an endograft (EG) that is designed to prevent blood leakage in the aneurysmal sac. While most EGs have equivalent operating principles, the hemodynamic environment established by different EGs is not necessarily the same. So, to unveil the post-EVAR hemodynamic properties, we need an EG-specific computational approach that currently lacks from the literature. Endurant and Excluder are two EGs with similar pre-installation designs. We assumed that the flow conditions in the particular EGs do not vary significantly. The hypothesis was tested combining image reconstructions, computational fluid dynamics (CFD), and statistics, taking into account the postimplantation position of the EGs. Ten patients with Endurant EGs and ten patients with Excluder EGs were included in this study. The two groups were matched with respect to the preoperative morphological characteristics of the AAAs. The EG models are derived from image reconstructions of postoperative computed tomography scans. Wall shear stress (WSS), displacement force, velocity, and helicity were calculated in regions of interest within the EG structures, i.e., the main body, the upper and lower part of the limbs. Excluder generated higher WSS compared to Endurant, especially on the lower part of the limbs (p = 0.001). Spatial fluctuations of WSS were observed on the upper part of the Excluder limbs. Higher blood velocity was induced by Excluder in all the regions of interest (p = 0.04, p = 0.01, and p = 0.004). Focal points of secondary flow were detected in the main body of Endurant and the limbs of Excluder. The displacement force acting on the lower part of the Excluder limbs was stronger compared to the Endurant one (p = 0.03). The results showed that two similar EGs implanted in similar AAAs can induce significantly different flow properties. The delineation of the hemodynamic features associated with the various commercially available EGs could further promote the personalization of treatment offered to aneurysmal patients and inspire ideas for the improvement of EG designs in the future.

An in Vitro Twist Fatigue Test of Fabric Stent-Grafts Supported by Z-Stents vs. Ringed Stents

Whereas buckling can cause type III endoleaks, long-term twisting of a stent-graft was investigated here as a mechanism leading to type V endoleak or endotension. Two experimental device designs supported with Z-stents having strut angles of 35° or 45° were compared to a ringed control under accelerated twisting. Damage to each device was assessed and compared after different durations of twisting, with focus on damage that may allow leakage. Stent-grafts with 35° Z-stents had the most severe distortion and damage to the graft fabric. The 45° Z-stents caused less fabric damage. However, consistent stretching was still seen around the holes for sutures, which attach the stents to the graft fabric. Larger holes may become channels for fluid percolation through the wall. The ringed stent-graft had the least damage observed. Stent apexes with sharp angles appear to be responsible for major damage to the fabrics. Device manufacturers should consider stent apex angle when designing stent-grafts, and ensure their devices are resistant to twisting.

Contemporary Role of Computational Analysis in Endovascular Treatment for Thoracic Aortic Disease

In the past decade, thoracic endovascular aortic repair (TEVAR) has become the primary treatment option in descending aneurysm and dissection. The clinical outcome of this minimally invasive technique is strictly related to an appropriate patient/stent graft selection, hemodynamic interactions, and operator skills. In this context, a quantitative assessment of the biomechanical stress induced in the aortic wall due to the stent graft may support the planning of the procedure. Different techniques of medical imaging, like computed tomography or magnetic resonance imaging, can be used to evaluate dynamics in the thoracic aorta. Such information can also be combined with dedicated patient-specific computer-based simulations, to provide a further insight into the biomechanical aspects. In clinical practice, computational analysis might show the development of aortic disease, such as the aortic wall segments which experience higher stress in places where rupture and dissection may occur. In aortic dissections, the intimal tear is usually located at the level of the sino-tubular junction and/or at the origin of the left subclavian artery. Besides, computational models may potentially be used preoperatively to predict stent graft behavior, virtually testing the optimal stent graft sizing, deployment, and conformability, in order to provide the best endovascular treatment. The present study reviews the current literature regarding the use of computational tools for TEVAR biomechanics, highlighting their potential clinical applications.

Evaluating Design of Abdominal Aortic Aneurysm Endografts in a Patient-Specific Model Using Computational Fluid Dynamics

Computer modeling of blood flow in patient-specific anatomies can be a powerful tool for evaluating the design of implantable medical devices. We assessed three different endograft designs, which are implantable devices commonly used to treat patients with abdominal aortic aneurysms (AAAs). Once implanted, the endograft may shift within the patient’s aorta allowing blood to flow into the aneurismal sac. One potential cause for this movement is the pulsatile force experienced by the endograft over the cardiac cycle. We used contrast-enhanced computed tomography angiography (CTA) data from four patients with diagnosed AAAs to build patient-specific models using 3D segmentation. For each of the four patients, we constructed a baseline model from the patient’s preoperative CTA data. In addition, geometries characterizing three distinct endograft designs were created, differing by where each device bifurcated into two limbs (proximal bifurcation, mid bifurcation, and distal bifurcation). Computational fluid dynamics (CFD) was used to simulate blood flow, utilizing patient-specific boundary conditions. Pressures, flows, and displacement forces on the endograft surface were calculated. The curvature and surface area of each device was quantified for all patients. The magnitude of the total displacement force on each device ranged from 2.43 N to 8.68 N for the four patients examined. Within each of the four patient anatomies, the total displacement force was similar (varying at least by 0.12 N and at most by 1.43 N), although there were some differences in the direction of component forces. Proximal bifurcation and distal bifurcation geometries consistently generated the smallest and largest displacement forces, respectively, with forces observed in the mid bifurcation design falling in between the two devices. The smallest curvature corresponded to the smallest total displacement force, and higher curvature values generally corresponded to higher magnitudes of displacement force. The same trend was seen for the surface area of each device, with lower surface areas resulting in lower displacement forces and vise versa. The patient with the highest blood pressure displayed the highest magnitudes of displacement force. The data indicate that curvature, device surface area, and patient blood pressure impact the magnitude of displacement force acting on the device. Endograft design may influence the displacement force experienced by an implanted endograft, with the proximal bifurcation design showing a small advantage for minimizing the displacement force on endografts.

Increased artery wall stress post-stenting leads to greater intimal thickening

Since the first human procedure in the late 1980s, vascular stent implantation has been accepted as a standard form of treatment of atherosclerosis. Despite their tremendous success, these medical devices are not without their problems, as excessive neointimal hyperplasia can result in the formation of a new blockage (restenosis). Clinical data suggest that stent design is a key factor in the development of restenosis. Additionally, computational studies indicate that the biomechanical environment is strongly dependent on the geometrical configuration of the stent, and, therefore, possibly involved in the development of restenosis. We hypothesize that stents that induce higher stresses on the artery wall lead to a more aggressive pathobiologic response, as determined by the amount of neointimal hyperplasia. The aim of this investigation was to examine the role of solid biomechanics in the development of restenosis. A combination of computational modeling techniques and in vivo analysis were employed to investigate the pathobiologic response to two stent designs that impose greater or lesser levels of stress on the artery wall. Stent designs were implanted in a porcine model (pigs) for approximately 28 days and novel integrative pathology techniques (quantitative micro-computed tomography, histomorphometry) were utilized to quantify the pathobiologic response. Concomitantly, computational methods were used to quantify the mechanical loads that the two stents place on the artery. Results reveal a strong correlation between the computed stress values induced on the artery wall and the pathobiologic response; the stent that subjected the artery to the higher stresses had significantly more neointimal thickening at stent struts (high-stress stent: 0.197±0.020 mm vs low-stress stent: 0.071±0.016 mm). Therefore, we conclude that the pathobiologic differences are a direct result of the solid biomechanical environment, confirming the hypothesis that stents that impose higher wall stresses will provoke a more aggressive pathobiological response.

Computer-Aided Patient-Specific Coronary Artery Graft Design Improvements Using CFD Coupled Shape Optimizer

This study aims to (i) demonstrate the efficacy of a new surgical planning framework for complex cardiovascular reconstructions, (ii) develop a computational fluid dynamics (CFD) coupled multi-dimensional shape optimization method to aid patient-specific coronary artery by-pass graft (CABG) design and, (iii) compare the hemodynamic efficiency of the sequential CABG, i.e., raising a daughter parallel branch from the parent CABG in patient-specific 3D settings. Hemodynamic efficiency of patient-specific complete revascularization scenarios for right coronary artery (RCA), left anterior descending artery (LAD), and left circumflex artery (LCX) bypasses were investigated in comparison to the stenosis condition. Multivariate 2D constraint optimization was applied on the left internal mammary artery (LIMA) graft, which was parameterized based on actual surgical settings extracted from 2D CT slices. The objective function was set to minimize the local variation of wall shear stress (WSS) and other hemodynamic indices (energy dissipation, flow deviation angle, average WSS, and vorticity) that correlate with performance of the graft and risk of re-stenosis at the anastomosis zone. Once the optimized 2D graft shape was obtained, it was translated to 3D using an in-house "sketch-based" interactive anatomical editing tool. The final graft design was evaluated using an experimentally validated second-order non-Newtonian CFD solver incorporating resistance based outlet boundary conditions. 3D patient-specific simulations for the healthy coronary anatomy produced realistic coronary flows. All revascularization techniques restored coronary perfusions to the healthy baseline. Multi-scale evaluation of the optimized LIMA graft enabled significant wall shear stress gradient (WSSG) relief (~34%). In comparison to original LIMA graft, sequential graft also lowered the WSSG by 15% proximal to LAD and diagonal bifurcation. The proposed sketch-based surgical planning paradigm evaluated the selected coronary bypass surgery procedures based on acute hemodynamic readjustments of aorta-CA flow. This methodology may provide a rational to aid surgical decision making in time-critical, patient-specific CA bypass operations before in vivo execution.