Hemodynamic analysis of a novel stent graft design with slit perforations in thoracic aortic aneurysm

Numerical analysis of blood flow in a branched modular stent-graft for aneurysms covering all zones of the aortic arch

Due to the anatomical complexity of the aortic arch for the development of stent-grafts for total repair, this region remains without a validated and routinely used endovascular option. In this work, a modular stent-graft for aneurysms that covers all aortic arch zones, proposed by us and previously structurally evaluated, was evaluated from the point of view of haemodynamics using fluid-structural numerical simulations. Blood was assumed to be non-Newtonian shear-thinning using the Carreau model, and the arterial wall was assumed to be anisotropic hyperelastic using the Holzapfel model. Nitinol and expanded polytetrafluoroethylene (PTFE-e) were used as materials for the stents and the graft, respectively. Nitinol was modelled as a superelastic material with shape memory by the Auricchio model, and PTFE-e was modelled as an isotropic linear elastic material. To validate the numerical model, a silicone model representative of the aneurysmal aorta was subjected to tests on an experimental bench representative of the circulatory system. The numerical results showed that the stent-graft restored flow behaviour, making it less oscillatory, but increasing the strain rate, turbulence kinetic energy, and viscosity compared to the pathological case. Taking the mean of the entire cycle, the increase in turbulence kinetic energy was 198.82% in the brachiocephalic trunk, 144.63% in the left common carotid artery and 284.03% in the left subclavian artery after stent-graft implantation. Based on wall shear stress parameters, it was possible to identify that the internal branches of the stent-graft and the stent-graft fixation sites in the artery were the most favourable regions for the deposition and accumulation of thrombus. In these regions, the oscillating shear index reached the maximum value of 0.5 and the time-averaged wall shear stress was close to zero, which led the relative residence time to reach values above 15 Pa-1. The stent-graft was able to preserve flow in the supra-aortic branches.

Intracranial bypass for giant aneurysms treatment assessed by computational fluid dynamics (CFD) analysis

Unruptured giant intracranial aneurysms (GIA) are those with diameters of 25 mm or greater. As aneurysm size is correlated with rupture risk, GIA natural history is poor. Parent artery occlusion or trapping plus bypass revascularization should be considered to encourage intra-aneurysmal thrombosis when other treatment options are contraindicated. The mechanistic background of these methods is poorly studied. Thus, we assessed the potential of computational fluid dynamics (CFD) and fluid-structure interaction (FSI) analyses for clinical use in the preoperative stage. A CFD investigation in three patient-specific flexible models of whole arterial brain circulation was performed. A C6 ICA segment GIA model was created based on CT angiography. Two models were then constructed that simulated a virtual bypass in combination with proximal GIA occlusion, but with differing middle cerebral artery (MCA) recipient vessels for the anastomosis. FSI and CFD investigations were performed in three models to assess changes in flow pattern and haemodynamic parameters alternations (wall shear stress (WSS), oscillatory shear index (OSI), maximal time averaged WSS (TAWSS), and pressure). General flow splitting across the entire domain was affected by virtual bypass procedures, and any deficiency was partially compensated by a specific configuration of the circle of Willis. Following the implementation of bypass procedures, a reduction in haemodynamic parameters was observed within the aneurysm in both cases under analysis. In the case of the temporal MCA branch bypass, the decreases in the studied parameters were slightly greater than in the frontal MCA branch bypass. The reduction in the magnitude of the chosen area-averaged parameters (averaged over the aneurysm wall surface) was as follows: WSS 35.7%, OSI 19.0%, TAWSS 94.7%, and pressure 24.2%. FSI CFD investigation based on patient-specific anatomy models with subsequent stimulation of virtual proximal aneurysm occlusion in conjunction with bypass showed that this method creates a pro-thrombotic favourable environment whilst reducing intra-aneurysmal pressure leading to shrinking. MCA branch recipient selection for optimum haemodynamic conditions should be evaluated individually in the preoperative stage.

Study on the hemodynamic effects of different pulsatile working modes of a rotary blood pump using a microfluidic platform that realizes in vitro cell culture effectively

Rotary blood pumps (RBPs) operating at a constant speed generate non-physiologic blood pressure and flow rate, which can cause endothelial dysfunction, leading to adverse clinical events in peripheral blood vessels and other organs. Notably, pulsatile working modes of the RBP can increase vascular pulsatility to improve arterial endothelial function. However, the laws and related mechanisms of differentially regulating arterial endothelial function under different pulsatile working modes are still unclear. This knowledge gap hinders the optimal selection of the RBP working modes. To address these issues, this study developed a multi-element in vitro endothelial cell culture system (ECCS), which could realize in vitro cell culture effectively and accurately reproduce blood pressure, shear stress, and circumferential strain in the arterial endothelial microenvironment. Performance of this proposed ECCS was validated with numerical simulation and flow experiments. Subsequently, this study investigated the effects of four different pulsation frequency modes that change once every 1-4-fold cardiac cycles (80, 40, 80/3, and 20 cycles per min, respectively) of the RBP on the expression of nitric oxide (NO) and reactive oxygen species (ROS) in endothelial cells. Results indicated that the 2-fold and 3-fold cardiac cycles significantly increased the production of NO and prevented the excessive generation of ROS, potentially minimizing the occurrence of endothelial dysfunction and related adverse events during the RBP support, and were consistent with animal study findings. In general, this study may provide a scientific basis for the optimal selection of the RBP working modes and potential treatment options for heart failure.

Correct Closure of the Left Atrial Appendage Reduces Stagnant Blood Flow and the Risk of Thrombus Formation: A Proof-of-Concept Experimental Study Using 4D Flow Magnetic Resonance Imaging

OBJECTIVE

The study was conducted to investigate the effect of correct occlusion of the left atrial appendage (LAA) on intracardiac blood flow and thrombus formation in patients with atrial fibrillation (AF) using four-dimensional (4D) flow magnetic resonance imaging (MRI) and three-dimensional (3D)-printed phantoms.

MATERIALS AND METHODS

Three life-sized 3D-printed left atrium (LA) phantoms, including a pre-occlusion (i.e., before the occlusion procedure) model and correctly and incorrectly occluded post-procedural models, were constructed based on cardiac computed tomography images from an 86-year-old male with long-standing persistent AF. A custom-made closed-loop flow circuit was set up, and pulsatile simulated pulmonary venous flow was delivered by a pump. 4D flow MRI was performed using a 3T scanner, and the images were analyzed using MATLAB-based software (R2020b; Mathworks). Flow metrics associated with blood stasis and thrombogenicity, such as the volume of stasis defined by the velocity threshold (|V̅| < 3 cm/s), surface-and-time-averaged wall shear stress (WSS), and endothelial cell activation potential (ECAP), were analyzed and compared among the three LA phantom models.

RESULTS

Different spatial distributions, orientations, and magnitudes of LA flow were directly visualized within the three LA phantoms using 4D flow MRI. The time-averaged volume and its ratio to the corresponding entire volume of LA flow stasis were consistently reduced in the correctly occluded model (70.82 mL and 39.0%, respectively), followed by the incorrectly occluded (73.17 mL and 39.0%, respectively) and pre-occlusion (79.11 mL and 39.7%, respectively) models. The surface-and-time-averaged WSS and ECAP were also lowest in the correctly occluded model (0.048 Pa and 4.004 Pa^-1 , respectively), followed by the incorrectly occluded (0.059 Pa and 4.792 Pa^-1 , respectively) and pre-occlusion (0.072 Pa and 5.861 Pa^-1 , respectively) models.

CONCLUSION

These findings suggest that a correctly occluded LAA leads to the greatest reduction in LA flow stasis and thrombogenicity, presenting a tentative procedural goal to maximize clinical benefits in patients with AF.

Abnormal flow pattern of low wall shear stress and high oscillatory shear index in spontaneous vertebral artery dissection with vertebral artery hypoplasia

Introduction

Spontaneous vertebral artery dissection (sVAD) might tend to develop in vertebral artery hypoplasia (VAH) with hemodynamic dysfunction and it is crucial to assess hemodynamics in sVAD with VAH to investigate this hypothesis. This retrospective study aimed to quantify hemodynamic parameters in patients with sVAD with VAH.

Methods

Patients who had suffered ischemic stroke due to an sVAD of VAH were enrolled in this retrospective study. The geometries of 14 patients (28 vessels) were reconstructed using Mimics and Geomagic Studio software from CT angiography (CTA). ANSYS ICEM and ANSYS FLUENT were utilized for mesh generation, set boundary conditions, solve governing equations, and perform numerical simulations. Slices were obtained at the upstream area, dissection or midstream area and downstream area of each VA. The blood flow patterns were visualized through instantaneous streamline and pressure at peak systole and late diastole. The hemodynamic parameters included pressure, velocity, time-averaged blood flow, time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), endothelial cell action potential (ECAP), relative residence time (RRT) and time-averaged nitric oxide production rate (TAR_NO).

Results

Significant focal increased velocity was present in the dissection area of steno-occlusive sVAD with VAH compared to other nondissected areas (0.910 m/s vs. 0.449 vs. 0.566, p < 0.001), while focal slow flow velocity was observed in the dissection area of aneurysmal dilatative sVAD with VAH according to velocity streamlines. Steno-occlusive sVAD with VAH arteries had a lower time-averaged blood flow (0.499 cm³/s vs. 2.268, p < 0.001), lower TAWSS (1.115 Pa vs. 2.437, p = 0.001), higher OSI (0.248 vs. 0.173, p = 0.006), higher ECAP (0.328 Pa^-1 vs. 0.094, p = 0.002), higher RRT (3.519 Pa^-1 vs. 1.044, p = 0.001) and deceased TAR_NO (104.014 nM/s vs. 158.195, p < 0.001) than the contralateral VAs.

Conclusion

Steno-occlusive sVAD with VAH patients had abnormal blood flow patterns of focal increased velocity, low time-averaged blood flow, low TAWSS, high OSI, high ECAP, high RRT and decreased TAR_NO. These results provide a good basis for further investigation of sVAD hemodynamics and support the applicability of the CFD method in testing the hemodynamic hypothesis of sVAD. More detailed hemodynamic conditions with different stages of sVAD are warranted in the future.

Hemodynamic numerical simulation of aortic arch modular inner branched stent-graft in eight early patients from the first-in-human case series

Background

The modular inner branched stent-graft (MIBSG) (WeFlow-Arch™) is an emerging device for challenging aortic arch pathologies. Hemodynamic numerical simulation is conducive to predicting long-term outcomes as well as optimizing the stent-graft design.

Objective

This study aims to analyze the hemodynamic characteristics of the MIBSG devices based on numerical simulation analyses.

Methods

From June 2019 to June 2021, MIBSGs were utilized in eight cases. Numerical simulation analyses of branch perfusion and indicators including the time-averaged wall shear stress, oscillatory shear index, and relative residence time were performed.

Results

Lesions involved Zone 1 (n = 2), Zone 2 (n = 4), and Zone 3 (n = 2). Branched stent-grafts were deployed in the innominate artery and left common carotid artery (n = 5) or in the innominate artery and left subclavian artery (n = 3). The hemodynamic change in common was increased perfusion in the descending aorta and left common carotid artery. Half of the patients had increased cerebral perfusion of 8.7% at most, and the other half of the patients showed a reduction of 5.3% or less. Case 3 was considered to have acquired the greatest improvement in hemodynamic features.

Conclusion

The MIBSG showed improved hemodynamic features in most cases. The design of the MIBSG could be partly modified to acquire better hemodynamic performance.

Review of the Development of Hemodynamic Modeling Techniques to Capture Flow Behavior in Arteries Affected by Aneurysm, Atherosclerosis, and Stenting

Cardiovascular diseases (CVDs) are the leading cause of death in the developed world. CVD can include atherosclerosis, aneurysm, dissection, or occlusion of the main arteries. Many CVDs are caused by unhealthy hemodynamics. Some CVDs can be treated with the implantation of stents and stent grafts. Investigations have been carried out to understand the effects of stents and stent grafts have on arteries and the hemodynamic changes post-treatment. Numerous studies on stent hemodynamics have been carried out using computational fluid dynamics (CFD) which has yielded significant insight into the effect of stent mesh design on near-wall blood flow and improving hemodynamics. Particle image velocimetry (PIV) has also been used to capture behavior of fluids that mimic physiological hemodynamics. However, PIV studies have largely been restricted to unstented models or intra-aneurysmal flow rather than peri or distal stent flow behaviors. PIV has been used both as a standalone measurement method and as a comparison to validate the CFD studies. This article reviews the successes and limitations of CFD and PIV-based modeling methods used to investigate the hemodynamic effects of stents. The review includes an overview of physiology and relevant mechanics of arteries as well as consideration of boundary conditions and the working fluids used to simulate blood for each modeling method along with the benefits and limitations introduced.

A Finite Element Based Analysis of a Hemodynamics Efficient Flow Stent Suitable for Different Abdominal Aneurysm Shapes

This research aimed to examine the impact of a proposed flow stent (PFS) on different abdominal artery shapes. For that purpose, a finite element-based model using the computational fluid dynamics (CFD) method is developed. The effect of PFS intervention on the hemodynamic efficiency is estimated by all of the significant criteria used for the evaluation of aneurysm occlusion and possible rupture; the flow velocity, pressure, wall shear stress (WSS), and WSS-related indices. Results showed that PFS intervention preserves the effects of high flow rate and decreases irregular flow recirculation in the sac of the aneurysm. The flow velocity decreases inside the aneurysm sac in the range of 55% to 80%. The time-averaged wall shear stress (TAWSS) was reduced from 42% to 53% by FPS deployment. The simulation results implies that PFS could heal an aneurysm efficiently with a mechanism that causes the development of thrombus and ultimately leads to aneurysm resorption.

Haemodynamic Analysis of Branched Endografts for Complex Aortic Arch Repair

This study aims to investigate the haemodynamic response induced by implantation of a double-branched endograft used in thoracic endovascular aortic repair (TEVAR) of the aortic arch. Anatomically realistic models were reconstructed from CT images obtained from patients who underwent TEVAR using the RelayPlus double-branched endograft implanted in the aortic arch. Two cases (Patient 1, Patient 2) were included here, both patients presented with type A aortic dissection before TEVAR. To examine the influence of inner tunnel branch diameters on localised flow patterns, three tunnel branch diameters were tested using the geometric model reconstructed for Patient 1. Pulsatile blood flow through the models was simulated by numerically solving the Navier–Stokes equations along with a transitional flow model. The physiological boundary conditions were imposed at the model inlet and outlets, while the wall was assumed to be rigid. Our simulation results showed that the double-branched endograft allowed for the sufficient perfusion of blood to the supra-aortic branches and restored flow patterns expected in normal aortas. The diameter of tunnel branches in the device plays a crucial role in the development of flow downstream of the branches and thus must be selected carefully based on the overall geometry of the vessel. Given the importance of wall shear stress in vascular remodelling and thrombus formation, longitudinal studies should be performed in the future in order to elucidate the role of tunnel branch diameters in long-term patency of the supra-aortic branches following TEVAR with the double-branched endograft.

On the Potential Self-Amplification of Aneurysms Due to Tissue Degradation and Blood Flow Revealed From FSI Simulations

Tissue degradation plays a crucial role in the formation and rupture of aneurysms. Using numerical computer simulations, we study the combined effects of blood flow and tissue degradation on intra-aneurysm hemodynamics. Our computational analysis reveals that the degradation-induced changes of the time-averaged wall shear stress (TAWSS) and oscillatory shear index (OSI) within the aneurysm dome are inversely correlated. Importantly, their correlation is enhanced in the process of tissue degradation. Regions with a low TAWSS and a high OSI experience still lower TAWSS and higher OSI during degradation. Furthermore, we observed that degradation leads to an increase of the endothelial cell activation potential index, in particular, at places experiencing low wall shear stress. These findings are robust and occur for different geometries, degradation intensities, heart rates and pressures. We interpret these findings in the context of recent literature and argue that the degradation-induced hemodynamic changes may lead to a self-amplification of the flow-induced progressive damage of the aneurysmal wall.

Aortic arch aneurysm repair - Unsteady hemodynamics and perfusion at different heart rates

The aortic arch aneurysm is a complex disease that requires branching of one or more aortic arch vessels and can be fatal if left untreated. In this in vitro study, we examine the effect of the treatment approach on the unsteady hemodynamics and blood perfusion to the upper vessel's in models of an aortic arch aneurysm, and of the three common repair approaches: open-chest surgical repair, chimney, and hybrid approach. A particle image velocimetry method was used to quantify the unsteady hemodynamics in the four models simulated in a mock circulatory loop, to evaluate unsteady hemodynamic parameters and measure perfusion to the brain and the upper body. According to the findings, in terms of perfusion to the brain and upper body, the surgery model has the highest flow rate comparing to the other models in most heart-rate conditions. It also shows oscillatory parameters in the upper vessels which in normal arteries are correlated with a better arterial function. Between the two endovascular procedures, the hybrid model exhibits slightly better hemodynamic characteristics than the chimney model, with lower shear stresses and more oscillatory flow and WSS in the upper vessels. The hybrid model had lower perfusion flow rates to upper vessels during rest conditions (90BPM). However, unlike the other models, perfusion in the hybrid model increased with heart rate, thus at 135 BPM, it results in flow rate to upper vessels similar to that of the chimney model. The results of this study may shed light on future endograft' design and placement techniques.

Effect of Geometric Accuracy at the Proximal Landing Zone on Simulation Results for Thoracic Endovascular Repair Patients

PURPOSE

Existing hemodynamic studies on aortic dissection after thoracic endovascular aortic repair (TEVAR) apply geometric simplifications. This study aims to evaluate the necessity of more accurate geometries at the proximal landing zone in computational fluid dynamic (CFD) studies.

METHODS

Three patient-specific 3D aortic dissection models with different geometric accuracies at the proximal landing zone were manually fabricated for CFD simulations: (i) model 1 without the stent graft (SG), (ii) model 2 with the metal stent, and (iii) model 3 with the SG. The flow distribution, flow pattern, and wall shear stress (WSS)-related indicators in these three models were compared.

RESULTS

The flow distributions were quite similar for the three models, with a maximum absolute difference of 0.27% at the left suclavian artery (LSA) between models 1 and 3 because of partial coverage. A more chaotic flow pattern was observed at the proximal landing zone in model 3, with significant regional differences in the WSS-related indicator distributions. The upstream and downstream WSS-related indicator distributions were quite similar for the three models.

CONCLUSIONS

The flow pattern and hemodynamic parameter distributions were affected by the geometric accuracy only in a small region near the proximal landing zone. The flow split was hardly affected by the LSA partial coverage, indicating that the coverage may have slight effects on short-term blood perfusion. However, this conclusion needs to be verified in future studies with larger sample sizes.

The effect of the entry and re-entry size in the aortic dissection: a two-way fluid-structure interaction simulation

Aortic dissection (AD) is one of the most catastrophic cardiovascular diseases. AD occurs when a layer inside the aorta is disrupted and gives rise to the formation of a true lumen and a false lumen. These lumens can be connected through tears in the intimal flap which are known as entries. Despite being known for about two centuries, the effects of many factors on the morbidity and mortality of this disease are still unknown. As the blood interaction with the aorta is crucial in the severity and the progression of the aortic dissection, a biomechanical approach is chosen to investigate the influence of different morphologies on the severity of this disease. Using the finite element method (FEM) and the fluid-structure interaction (FSI) approach, we have evaluated the blood flow characteristics along the diseased aorta, in conjunction with the deformation of the aortic wall. In this study, an idealized geometry of a dissected descending aorta (type B) with two entries has been studied. The values for the diameter of the entry tear were chosen to be 5 mm and 10 mm. Therefore, a total of four conditions were investigated. According to our results, the retrograde flow through the proximal tear is dependent on the size of the distal re-entry and vice versa. Our results revealed that when both entry and re-entry tears are 10 mm in diameter, the flow passes through the true and false lumens with smaller resistance, resulting in a smaller flutter of the intimal flap, and therefore more stable intimal flap. Major oscillation frequencies of 2.5 Hz and 7.4 Hz were observed for the oscillation of the intimal flap, and amplitudes of the waves with higher frequencies were negligible.