Computational fluid dynamics: hemodynamic changes in abdominal aortic aneurysm after stent-graft implantation

Fabrication of Soft Transparent Patient-Specific Vascular Models with Stereolithographic 3D printing and Thiol-Based Photopolymerizable Coatings

An ideal vascular phantom should be anatomically accurate, have mechanical properties as close as possible to the tissue, and be sufficiently transparent for ease of visualization. However, materials that enable the convergence of these characteristics have remained elusive. The fabrication of patient-specific vascular phantoms with high anatomical fidelity, optical transparency, and mechanical properties close to those of vascular tissue is reported. These final properties are achieved by 3D printing patient-specific vascular models with commercial elastomeric acrylic-based resins before coating them with thiol-based photopolymerizable resins. Ternary thiol-ene-acrylate chemistry is found optimal. A PETMP/allyl glycerol ether (AGE)/polyethylene glycol diacrylate (PEGDA) coating with a 30/70% AGE/PEGDA ratio applied on a flexible resin yielded elastic modulus, UTS, and elongation of 3.41 MPa, 1.76 MPa, and 63.2%, respectively, in range with the human aortic wall. The PETMP/AGE/PEGDA coating doubled the optical transmission from 40% to 80%, approaching 88% of the benchmark silicone-based elastomer. Higher transparency correlates with a decrease in surface roughness from 2000 to 90 nm after coating. Coated 3D-printed anatomical replicas are showcased for pre-procedural planning and medical training with good radio-opacity and echogenicity. Thiol-click chemistry coatings, as a surface treatment for elastomeric stereolithographic 3D-printed objects, address inherent limitations of photopolymer-based additive manufacturing.

Changes in Internal Thoracic Artery Blood Flow According to the Degree of Stenosis of the Anterior Descending Branch of the Left Coronary Artery

PURPOSE

Computational fluid dynamics has enabled the evaluation of coronary flow reserve. The purpose of this study was to clarify the hemodynamic variation and reserve potential of the left internal thoracic artery (LITA).

METHODS

Four patients were selected on the basis of various native coronary stenosis patterns and graft design. The wall shear stress and oscillatory shear index were measured, and one patient was selected. Next, we created three hypothetical lesions with 75%, 90%, and 99% stenosis in front of the graft anastomosis, and compared the changes in LITA blood flow and coronary flow distribution.

RESULTS

In the 75% to 90% stenosis model, blood flow was significantly higher in the native coronary flow proximal to the coronary artery bypass anastomosis regardless of time phase. In the 99% stenosis model, blood flow from the LITA was significantly dominant compared to native coronary flow at the proximal site of anastomosis. The range of LITA flow variability was the largest at 99% stenosis, with a difference of 70 ml/min.

CONCLUSION

The 99% stenosis model showed the highest LITA flow. The range of LITA flow variability is large, suggesting that it may vary according to the rate of native coronary stenosis.

The pulsatile 3D-Hemodynamics in a doubly afflicted human descending abdominal artery with iliac branching

The study of patient-specific human arterial flow dynamics is well known to face challenges like a) apt geometric modelling, b) bifurcation zone meshing, and c) capturing the hemodynamic prone to variations with multiple disease complications. Due to aneurysms and stenosis in the same arterial network, the blood flow dynamics get affected, which needs to be explored. This study develops a new protocol for accurate geometric modelling, bifurcation zone meshing and numerically investigates the arterial network with abdominal aortic aneurysms (AAA) and right internal iliac stenosis (RIIAS). A realistic arterial model is reconstructed from the computed tomography (CT) data of a human subject. To understand the combined effect of the aneurysm and aortoiliac occlusive diseases in a patient, an arterial network with AAA, RIIAS, multiple branches tapering, and curvature has been considered. Clinically significant pulsatile blood flow simulations have been carried out to trace the alteration in the flow dynamics with multiple pathological complications under consideration. The transient blood flow dynamics are investigated via wall shear stress, wall pressure, velocity contour, streamlines, vorticity, and swirling strength. During the systolic deceleration phase, the rhythmic nested rapid secondary oscillatory WSS, adverse pressure gradients, high WSS, and high WP bands are noticed. Also, the above studies will help researchers, clinicians, and doctors understand the influence of morphological changes on hemodynamics in cardiovascular studies.

Nanoparticle-Assisted Diagnosis and Treatment for Abdominal Aortic Aneurysm

An abdominal aortic aneurysm (AAA) is a localized dilatation of the aorta related to the regional weakening of the wall structure, resulting in substantial morbidity and mortality with the aortic ruptures as complications. Ruptured AAA is a dramatic catastrophe, and aortic emergencies constitute one of the leading causes of acute death in older adults. AAA management has been centered on surgical repair of larger aneurysms to mitigate the risks of rupture, and curative early diagnosis and effective pharmacological treatments for this condition are still lacking. Nanoscience provided a possibility of more targeted imaging and drug delivery system. Multifunctional nanoparticles (NPs) may be modified with ligands or biomembranes to target agents' delivery to the lesion site, thus reducing systemic toxicity. Furthermore, NPs can improve drug solubility, circulation time, bioavailability, and efficacy after systemic administration. The varied judiciously engineered nano-biomaterials can exist stably in the blood vessels for a long time without being taken up by cells. Here, in this review, we focused on the NP application in the imaging and treatment of AAA. We hope to make an overview of NP-assisted diagnoses and therapy in AAA and discussed the potential of NP-assisted treatment.

Numerical Simulation of Blood Flows in Patient-specific Abdominal Aorta with Primary Organs

The abdominal aorta is the largest artery in the abdominal cavity that supplies blood flows to vital organs through the complex visceral arterial branches, including the celiac trunk (the liver, stomach, spleen, etc.), the renal arteries (the kidneys) and the superior and inferior mesenteric arteries (the small and large intestine, pancreas, etc.). An accurate simulation of blood flows in this network of arteries is important for the understanding of the hemodynamics in various organs of healthy and diseased patients, but the computational cost is very high. As a result, most researchers choose to focus on a portion of the artery or use a low-dimensional approximation of the artery. In the present work, we introduce a parallel algorithm for the modeling of pulsatile flows in the abdominal aorta with branches to the primary organs, and an organ-based two-level method for calculating the resistances for the outflow boundary conditions. With this highly parallel approach, the simulation of the blood flow for a cardiac cycle of the anatomically detailed aorta can be obtained within a few hours, and the blood distribution to organs including liver, spleen and kidneys are also computed with certain accuracy. Moreover, we discuss the significant hemodynamic differences resulted from the influence of the peripheral branches. In addition, we examine the accuracy of the results with respect to the mesh size and time-step size and show the high parallel scalability of the proposed algorithm with up to 3000 processor cores.

Impact of Thoracic Endografting on the Hemodynamics of the Native Aorta: Pre- and Postoperative Assessments of Wall Shear Stress and Vorticity Using Computational Fluid Dynamics

PURPOSE

To quantify the hemodynamic consequences of thoracic endovascular aortic repair (TEVAR) by comparing the preoperative and postoperative wall shear stress (WSS) and vorticity profiles on computational fluid dynamics (CFD) simulations.

MATERIALS AND METHODS

The pre- and postoperative computed tomography (CT) scans from 20 consecutive patients (median age 69 years, range 20-87) treated for different thoracic aortic pathologies (11 aneurysms, 5 false aneurysms, 3 penetrating ulcers, and 1 traumatic aortic rupture) were segmented to construct patient-specific CFD models using a meshless code. The simulations were run over the cardiac cycle, and the WSS and vorticity values measured at the proximal and distal landing zones were compared.

RESULTS

The CFD runs provided 4-dimensional simulations of blood flow in all patients. WSS and vorticity profiles at the proximal landing zone (located in zones 0-3 in 15 patients) varied in 18 and 20 of the cases, respectively; WSS was increased in 11 cases and the vorticity in 9. Pre- and postoperative WSS median values were 4.19 and 4.90 Pa, respectively. Vorticity median values were 40.38 and 39.17 Hz, respectively.

CONCLUSION

TEVAR induces functional alterations in the native thoracic aorta, though the prognostic significance of these changes is still unknown. CFD appears to be a valuable tool to explore aortic hemodynamics, and its application in a larger series would help define a predictive role for these hemodynamic assessments.

Image-Based Displacements Analysis and Computational Blood Dynamics after Endovascular Aneurysm Repair

BACKGROUND

To examine intraheartbeat displacements (IHD) and geometrical changes of endografts for abdominal aortic aneurysm repair over the course of years, defined as follow-up displacements (FUD), and to correlate them with computational fluid dynamics (CFD). Despite the widespread use of endovascular aneurysm repair (EVAR), we still know little about endograft behavior after deployment.

METHODS

Two cases, treated with either expanded polytetrafluoroethylene on a nitinol stent frame (PI) or with woven polyester fabric sutured to a stainless-steel Z stent skeleton (PII), were submitted to dynamic computed tomography angiography at 1, 12, and 60 months after implantation. After segmentation, IHD were computed as displacements of the reconstructed surface with respect to the diastolic instant. Similarly, FUD were studied using imaging techniques that align temporal successive segmentations. In addition, numerical simulations for blood dynamics were performed to compute viscous forces, specifically wall shear stress and time-averaged wall shear stress (TAWSS).

RESULTS

IHD analysis showed slight translations without deformation for the PI endograft with respect to the stiffer stainless-steel endograft behavior of PII. FUD showed in PI motion of the metallic struts mainly focused on the distal main body of the endograft and in the zone overlapping with iliac branches. In PII, we observed a huge FUD in the middle and inferior-anterior regions of the main body. CFD analysis revealed changes of velocity patterns associated with remodeling of the iliac zone for PI and of the main body region for PII, where flow impinges the lumen wall and progressively induces deformation of the endograft wires. Measurement of TAWSS demonstrated flow disturbances in the enlarged region correlated with displacement analysis.

CONCLUSIONS

Image-based displacement analysis associated with CFD enabled very subtle evaluations of endograft behavior on different temporal scales. This kind of study could be helpful both for physicians, forecasting evolution during the life span of the endograft, and manufacturers, giving them useful information about endograft implant performance and design.

The effects of clinically-derived parametric data uncertainty in patient-specific coronary simulations with deformable walls

Cardiovascular simulations are increasingly used for noninvasive diagnosis of cardiovascular disease, to guide treatment decisions, and in the design of medical devices. Quantitative assessment of the variability of simulation outputs due to input uncertainty is a key step toward further integration of cardiovascular simulations in the clinical workflow. In this study, we present uncertainty quantification in computational models of the coronary circulation to investigate the effect of uncertain parameters, including coronary pressure waveform, intramyocardial pressure, morphometry exponent, and the vascular wall Young's modulus. We employ a left coronary artery model with deformable vessel walls, simulated via an Arbitrary-Lagrangian-Eulerian framework for fluid-structure interaction, with a prescribed inlet pressure and open-loop lumped parameter network outlet boundary conditions. Stochastic modeling of the uncertain inputs is determined from intra-coronary catheterization data or gathered from the literature. Uncertainty propagation is performed using several approaches including Monte Carlo, Quasi Monte Carlo sampling, stochastic collocation, and multi-wavelet stochastic expansion. Variabilities in the quantities of interest, including branch pressure, flow, wall shear stress, and wall deformation are assessed. We find that uncertainty in inlet pressures and intramyocardial pressures significantly affect all resulting QoIs, while uncertainty in elastic modulus only affects the mechanical response of the vascular wall. Variability in the morphometry exponent used to distribute the total downstream vascular resistance to the single outlets, has little effect on coronary hemodynamics or wall mechanics. Finally, we compare convergence behaviors of statistics of QoIs using several uncertainty propagation methods on three model benchmark problems and the left coronary simulations. From the simulation results, we conclude that the multi-wavelet stochastic expansion shows superior accuracy and performance against Quasi Monte Carlo and stochastic collocation methods.

Treatment of High Surgical Risk Thoracoabdominal Aortic Aneurysms with Stent Graft and Multilayer Bare Stents Joint Technique: Mid-Long-Term Clinical Results

BACKGROUND

This study aims to present the performance data on stent-graft and multilayer bare stents (MBS) joint technique in the treatment of high-risk thoracoabdominal aortic aneurysm (TAAA).

METHODS

From May 2012 to December 2015, 8 selective TAAA cases (ages 46-75 years) ineligible for surgical repair underwent the stent-graft and MBS joint procedure, and were closely followed up for a median of 32 months (range 14-58). Using computed tomography images, the aneurysm size, luminal blood flow diameter, and the covered visceral branches were analyzed.

RESULTS

Technical success was achieved in all patients (100%, 8/8). Twenty-four visceral branches were covered by MBS in total. There was no complication or death during hospital stay. During follow-up period, no death or complication occurred. Aneurysm shrinkage (maximum diameter decrease ≥5 mm) was observed in 7 patients. No aneurysm expansion was observed. Total aneurysm sac thrombosis was observed in all patients. The majority of covered side branches (23/24) were successfully preserved. No visceral ischemia or bleeding complications was observed during follow-up.

CONCLUSIONS

Total endovascular repair of TAAA using stent-graft and MBS joint technique may be a safe and effective alternative in high surgical risk patients. More approving clinical evidences about the safety and efficacy of this procedure are anticipated.

Computational study on hemodynamic changes in patient-specific proximal neck angulation of abdominal aortic aneurysm with time-varying velocity

Aneurysms are considered as a critical cardiovascular disease worldwide when they rupture. The clinical understanding of geometrical impact on the flow behaviour and biomechanics of abdominal aortic aneurysm (AAA) is progressively developing. Proximal neck angulations of AAAs are believed to influence the hemodynamic changes and wall shear stress (WSS) within AAAs. Our aim was to perform pulsatile simulations using computational fluid dynamics (CFD) for patient-specific geometry to investigate the influence of severe angular (≥ 60°) neck on AAA's hemodynamic and wall shear stress. The patient's geometrical characteristics were obtained from a computed tomography images database of AAA patients. The AAA geometry was reconstructed using Mimics software. In computational method, blood was assumed Newtonian fluid and an inlet varying velocity waveform in a cardiac cycle was assigned. The CFD study was performed with ANSYS software. The results of flow behaviours indicated that the blood flow through severe bending of angular neck leads to high turbulence and asymmetry of flows within the aneurysm sac resulting in blood recirculation. The high wall shear stress (WSS) occurred near the AAA neck and on surface of aneurysm sac. This study explained and showed flow behaviours and WSS progression within high angular neck AAA and risk prediction of abdominal aorta rupture. We expect that the visualization of blood flow and hemodynamic changes resulted from CFD simulation could be as an extra tool to assist clinicians during a decision making when estimation the risks of interventional procedures.

Non-invasive Evaluation of Fluid Dynamic of Aortoiliac Atherosclerotic Disease: Impact of Bifurcation Angle and Different Stent Configurations

Objectives

To non-invasively evaluate by computational fluid dynamic (CFD) analysis the physiology and rheology of aortoiliac bifurcation disease at different angles and different stent configurations.

Material and methods

For the analysis, we considered a physiologic model of abdominal aorta with an iliac bifurcation set at 30°, 45° and 70° without stenosis. Subsequently, a bilateral ostial common iliac stenosis of 80% was considered for each type of bifurcation. For the stent simulation, we reconstructed Zilver vascular self-expanding (Zilver; Cook, Bloomington, MN) and Palmaz Genesis Peripheral (Cordis, Miami, FL) stents.

Results

The physiologic model, across the different angles, static pressure, Reynolds number and stream function, were lower for the 30° bifurcation angle with a gradient from 70° to 30° angles, whereas all the other parameters were inversely higher. After stenting, all the fluid parameters decreased homogenously independent of the stent type, maintaining a gradient in favour of 30° compared to 45° and 70° angles. The absolute greater deviation from physiology was observed for low kissing when self-expandable stents were used across all angles; in particular, the wall shear stress was high at at 45° angle.

Conclusion

Bifurcation angle deeply impacts the physiology of aortoiliac bifurcations, which are used to predict the fluid dynamic profile after stenting. CFD, having the potential to be derived both from computed tomography scan or invasive angiography, appears to be an ideal tool to predict fluid dynamic profile before and after stenting in aortoiliac bifurcation.

Non-invasive Evaluation of Fluid Dynamic of Aortoiliac Atherosclerotic Disease: Impact of Bifurcation Angle and Different Stent Configurations

Objectives

To non-invasively evaluate by computational fluid dynamic (CFD) analysis the physiology and rheology of aortoiliac bifurcation disease at different angles and different stent configurations.

Material and methods

For the analysis, we considered a physiologic model of abdominal aorta with an iliac bifurcation set at 30°, 45° and 70° without stenosis. Subsequently, a bilateral ostial common iliac stenosis of 80% was considered for each type of bifurcation. For the stent simulation, we reconstructed Zilver vascular self-expanding (Zilver; Cook, Bloomington, MN) and Palmaz Genesis Peripheral (Cordis, Miami, FL) stents.

Results

The physiologic model, across the different angles, static pressure, Reynolds number and stream function, were lower for the 30° bifurcation angle with a gradient from 70° to 30° angles, whereas all the other parameters were inversely higher. After stenting, all the fluid parameters decreased homogenously independent of the stent type, maintaining a gradient in favour of 30° compared to 45° and 70° angles. The absolute greater deviation from physiology was observed for low kissing when self-expandable stents were used across all angles; in particular, the wall shear stress was high at at 45° angle.

Conclusion

Bifurcation angle deeply impacts the physiology of aortoiliac bifurcations, which are used to predict the fluid dynamic profile after stenting. CFD, having the potential to be derived both from computed tomography scan or invasive angiography, appears to be an ideal tool to predict fluid dynamic profile before and after stenting in aortoiliac bifurcation.

Effects on Aortoiliac Fluid Dynamics After Endovascular Sealing of Abdominal Aneurysms

OBJECTIVES:

To evaluate the effects on aortoiliac fluid dynamics after the implantation of an endograft based on endovascular aneurysm sealing (EVAS) versus endovascular aneurysm repair (EVAR) strategy.

METHODS:

An adaptive geometrical deformable model was used for aortic lumen segmentation in 8 patients before and after the surgery. Abdominal aneurysms were treated with an endograft based on the EVAS system (Nellix, n = 4) and with a device based on an anatomical fixation technology (n = 4). Pressure, blood velocity, and wall shear stress (WSS) were estimated at different aortic regions using computational fluid dynamics methods. Physiologic inlet/outlet flow values at the abdominal aorta, the celiac trunk, and the mesenteric and the renal arteries were set. Pressure references were set at iliac arteries outlet.

RESULTS:

Maximum aneurysm sizes were similar for both groups in the preoperative scans. The lumen area was lower after EVAR ( P < .05) and EVAS ( P < .01) compared to preoperative aortic lumen sizes. Pressure increase was higher in the proximal abdominal aorta after EVAS compared to EVAR (2.3 ± 0.3 mm Hg vs 0.9 ± 0.3 mm Hg, P < .001). Peak blood velocities inside the endografts were 3-fold higher for EVAS compared to EVAR (54 ± 5 cm/s vs 17 ± 4 cm/s, P < .01). Velocities at the iliac arteries also remained higher for EVAS (38 ± 4 cm/s vs 24 ± 4 cm/s, P < .05). Peak WSS at the iliac arteries remained higher for EVAS compared to EVAR group ( P < .05).

CONCLUSION:

The significant modification of the aortic bifurcation anatomy after EVAS alters aortoiliac fluid dynamics, showing a pressure impact at the renal arteries level and an acceleration of the blood velocity at the iliac region with a concomitant increase in peak WSS.

Deformation and dynamic response of abdominal aortic aneurysm sealing

Endovascular sealing is a new technique for the repair of abdominal aortic aneurysms. Commercially available in Europe since 2013, it takes a revolutionary approach to aneurysm repair through minimally invasive techniques. Although aneurysm sealing may be thought as more stable than conventional endovascular stent graft repairs, post-implantation movement of the endoprosthesis has been described, potentially leading to late complications. The paper presents for the first time a model, which explains the nature of forces, in static and dynamic regimes, acting on sealed abdominal aortic aneurysms, with references to real case studies. It is shown that elastic deformation of the aorta and of the endoprosthesis induced by static forces and vibrations during daily activities can potentially promote undesired movements of the endovascular sealing structure.

Deformation and dynamic response of abdominal aortic aneurysm sealing

Endovascular sealing is a new technique for the repair of abdominal aortic aneurysms. Commercially available in Europe since 2013, it takes a revolutionary approach to aneurysm repair through minimally invasive techniques. Although aneurysm sealing may be thought as more stable than conventional endovascular stent graft repairs, post-implantation movement of the endoprosthesis has been described, potentially leading to late complications. The paper presents for the first time a model, which explains the nature of forces, in static and dynamic regimes, acting on sealed abdominal aortic aneurysms, with references to real case studies. It is shown that elastic deformation of the aorta and of the endoprosthesis induced by static forces and vibrations during daily activities can potentially promote undesired movements of the endovascular sealing structure.

New imaging tools in cardiovascular medicine: computational fluid dynamics and 4D flow MRI

Blood flow imaging is a novel technology in cardiovascular medicine and surgery. Today, two types of blood flow imaging tools are available: measurement-based flow visualization including 4D flow MRI (or 3D cine phase-contrast magnetic resonance imaging), or echocardiography flow visualization software, and computer flow simulation modeling based on computational fluid dynamics (CFD). MRI and echocardiography flow visualization provide measured blood flow but have limitations in temporal and spatial resolution, whereas CFD flow calculates the flow according to assumptions instead of flow measurement, and it has sufficiently fine resolution up to the computer memory limit, and it enables even virtual surgery when combined with computer graphics. Blood flow imaging provides profound insight into the pathophysiology of cardiovascular diseases, because it quantifies and visualizes mechanical stress on the vessel walls or heart ventricle. Wall shear stress (WSS) is a stress on the endothelial wall caused by the near wall blood flow, and it is thought to be a predictor of atherosclerosis progression in coronary or aortic diseases. Flow energy loss (EL) is the loss of blood flow energy caused by viscous friction of turbulent diseased flow, and it is expected to be a predictor of ventricular workload on various heart diseases including heart valve disease, cardiomyopathy, and congenital heart diseases. Blood flow imaging can provide useful information for developing predictive medicine in cardiovascular diseases, and may lead to breakthroughs in cardiovascular surgery, especially in the decision-making process.

Influence of endoleak positions on the pressure shielding ability of stent-graft after endovascular aneurysm repair (EVAR) of abdominal aortic aneurysm (AAA)

Background

Abdominal aortic aneurysm (AAA) is a kind of dangerous aortic vascular disease, which is characterized by abdominal aorta partial enlargement. At present, endovascular aneurysm repair (EVAR) is one of the main treatments of abdominal aortic aneurysm. However for some patients after EVAR the aneurysm re-expanded and even ruptured, leading to poor postoperative effect. The stent-graft endoleak after EVAR was realized to influence the AAA in-sac pressure and contribute to the aneurysm re-enlargement.

Methods

In order to analyze the influence of endoleaks positions on the pressure shielding ability of stent-graft after EVAR, type I and type III endoleak models were reconstructed based on computed tomography (CT) scan images, and the hemodynamic environment in AAA was numerically simulated.

Results

When the endoleak was at the proximal position the pressure shielding ability will be obviously weakened. While, the pressure shielding ability was higher in the systole phase than that in diastole phase when the endoleak located at the middle or distal positions. Unfortunately, when the endoleak located at the proximal position, the pressure shielding ability would be relatively weak in the whole cardiac cycle.

Conclusions

The results revealed that the influence of endoleaks on pressure shielding ability of stent-graft was both location and time specific.

Stent graft surface movement after infrarenal abdominal aortic aneurysm repair: comparison of patients with and without a type 2 endoleak

OBJECTIVES

The aim was to compare multidirectional stent graft movement in patients with and without a type 2 endoleak.

METHODS

This was a retrospective case control study of patients being followed up after elective endovascular aneurysm repair of abdominal aortic aneurysms. The post-procedural and final follow up multislice computed tomography (MSCT) of 69 patients with and 74 without a type 2 endoleak were analyzed. Three dimensional (3D) surface models of the stent graft, delimited by landmarks using custom built software, were derived from these MSCT data. The stent graft was segmented in different zones, and the proportion of the total stent graft surface moving >9 mm between the post-procedural and the final follow up MSCT was calculated, given in percentages, and compared between groups. Changes of infrarenal neck, renal artery to stent graft distance, and freedom from stent graft related endoleaks were evaluated.

RESULTS

Overall surface movement was higher in the no endoleak (18.8%, IQR 0.1-45.1%) than in the type 2 endoleak group (5.3%, IQR 0-29.7%; p = .06). Furthermore, significantly higher surface movement in the no endoleak group was found in the proximal anchoring zone (p = .04) and the distal left limb (p = .01), which was the modular limb in 81.1% (p < .01). Neck diameter increase (1.0 mm, IQR 0-3.0 mm; p < .01) and renal artery to stent graft distance difference (0 mm, IQR 0-3.3 mm; p < .01) were significantly higher in the no endoleak group. Five patients in the no endoleak and one patient in the type 2 endoleak group suffered from a stent graft related endoleak (p = .27).

CONCLUSIONS

The presence of a type 2 endoleak is associated with decreased surface movement of the proximal anchoring zone and the distal modular limb of bifurcated stent grafts.

Patient-specific analysis of displacement forces acting on fenestrated stent grafts for endovascular aneurysm repair

Treatment options for abdominal aortic aneurysm (AAA) include highly invasive open surgical repair or minimally invasive endovascular aneurysm repair (EVAR). Despite being minimally invasive, some patients are not suitable for EVAR due to hostile AAA morphology. Fenestrated-EVAR (F-EVAR) was introduced to address these limitations of standard EVAR, where AAA is treated using a Fenestrated Stent Graft (FSG). In order to assess durability of F-EVAR, displacement forces acting on FSGs were analysed in this study, based on patient-specific geometries reconstructed from computed tomography (CT) scans. The magnitude and direction of the resultant displacement forces acting on the FSG were numerically computed using computational fluid dynamics (CFD) with a rigid wall assumption. Although displacement force arises from blood pressure and friction due to blood flow, numerical simulations elucidated that net blood pressure is the dominant contributor to the overall displacement force; as a result, time dependence of the resultant displacement force followed pressure waveform very closely. The magnitude of peak displacement force varied from 1.9N to 14.3N with a median of 7.0N. A strong positive correlation was found between inlet cross-sectional area (CSA), anterior/posterior (A/P) angle and the peak displacement force i.e. as inlet CSA or A/P angle increases, the magnitude of resultant displacement increases. This study manifests that while loads exerted by the pulsatile flow dictates the cyclic variation of the displacement force, its magnitude depends not only on blood pressure but also the FSG morphology, with the latter determining the direction of the displacement force.

In Vitro Validation of Patient-Specific Hemodynamic Simulations in Coronary Aneurysms Caused by Kawasaki Disease

To perform experimental validation of computational fluid dynamics (CFD) applied to patient specific coronary aneurysm anatomy of Kawasaki disease. We quantified hemodynamics in a patient-specific coronary artery aneurysm physical phantom under physiologic rest and exercise flow conditions. Using phase contrast MRI (PCMRI), we acquired 3-component flow velocity at two slice locations in the aneurysms. We then performed numerical simulations with the same geometry and inflow conditions, and performed qualitative and quantitative comparisons of velocities between experimental measurements and simulation results. We observed excellent qualitative agreement in flow pattern features. The quantitative spatially and temporally varying differences in velocity between PCMRI and CFD were proportional to the flow velocity. As a result, the percent discrepancy between simulation and experiment was relatively constant regardless of flow velocity variations. Through 1D and 2D quantitative comparisons, we found a 5-17% difference between measured and simulated velocities. Additional analysis assessed wall shear stress differences between deformable and rigid wall simulations. This study demonstrated that CFD produced good qualitative and quantitative predictions of velocities in a realistic coronary aneurysm anatomy under physiological flow conditions. The results provide insights on factors that may influence the level of agreement, and a set of in vitro experimental data that can be used by others to compare against CFD simulation results. The findings of this study increase confidence in the use of CFD for investigating hemodynamics in the specialized anatomy of coronary aneurysms. This provides a basis for future hemodynamics studies in patient-specific models of Kawasaki disease.

Three-dimensional hemodynamic design optimization of stents for cerebral aneurysms

Flow-diverting stents occlude aneurysms by diverting the blood flow from entering the aneurysm sac. Their effectiveness is determined by the thrombus formation rate, which depends greatly on stent design. The aim of this study was to provide a general framework for efficient stent design using design optimization methods, with a focus on stent hemodynamics as the starting point. Kriging method was used for completing design optimization. Three different cases of idealized stents were considered, and 40–60 samples from each case were evaluated using computational fluid dynamics. Using maximum velocity and vorticity reduction as objective functions, the optimized designs were identified from the samples. A number of optimized stent designs have been found from optimization, which revealed that a combination of high pore density and thin struts is desired. Additionally, distributing struts near the proximal end of aneurysm neck was found to be effective. The success of the methods and framework devised in this study offers a future possibility of incorporating other disciplines to carry out multidisciplinary design optimization.

Studying the Flow Dynamics in an Aortic Endograft with Crossed-limbs

Purpose

To evaluate the flow phenomena within an aortic endograft with crossed-limbs, comparing to an endograft with the ordinary limb bifurcation.

Methods

An endograft model with crossed-limbs was computationally reconstructed based on Computed Tomography patient-specific data, using commercially available software. Accordingly, its analogue model was reconstructed in the ordinary fashion (ordinary bifurcation). Computational fluid dynamics analysis was performed to determine and compare the flow fields, velocity profiles, pressure and shear stress distribution throughout the different parts of both endograft configurations, in different phases of the cardiac cycle.

Results

The flow patterns between the “Ballerina” and the classic endograft were similar, with flow disturbance near the inlet zone at late diastole and smooth flow patterns during the systolic phase. Both configurations presented similar pressure distribution patterns throughout the cardiac cycle. The highest and lowest pressures were demonstrated in the inlet-main body area and the iliac limbs, respectively. Marked differences were observed in the velocity profiles of the proximal limb segments between the two configurations, mostly in the peak- and end-systolic phase. The regions of lower velocities correlated well to low shear values. Differences in the shear stress distribution were noted between the two configurations in the systolic and, predominantly, in the diastolic phase.

Conclusions

There are differences in the velocity profiles and shear distribution between the limbs of the two endograft configurations. The pathophysiologic implication of our findings and their possible association with clinical events, such as thrombus apposition, deserves further investigation.

Preliminary computational hemodynamics study of double aortic aneurysms under multistage surgical procedures: an idealised model study

Double aortic aneurysm (DAA) falls under the category of multiple aortic aneurysms. Repair is generally done through staged surgery due to low invasiveness. In this approach, one aneurysm is cured per operation. Therefore, two operations are required for DAA. However, post-first-surgery rupture cases have been reported. Although the problems involved with managing staged surgery have been discussed for more than 30 years, investigation from a hemodynamic perspective has not been attempted. Hence, this is the first computational fluid dynamics approach to the DAA problem. Three idealized geometries were prepared: presurgery, thoracic aortic aneurysm (TAA) cured, and abdominal aortic aneurysm (AAA) cured. By applying identical boundary conditions for flow rate and pressure, the Navier-Stokes equation and continuity equations were solved under the Newtonian fluid assumption. Average pressure in TAA was increased by AAA repair. On the other hand, average pressure in AAA was decreased after TAA repair. Average wall shear stress was decreased at the peak in post-first-surgery models. However, the wave profile of TAA average wall shear stress was changed in the late systole phase after AAA repair. Since the average wall shear stress in the post-first-surgery models decreased and pressure at TAA after AAA repair increased, the TAA might be treated first to prevent rupture.

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.

Hybrid Endografts Combinations for the Treatment of Endoleak in Endovascular Abdominal Aortic Aneurysm Repair

Hybrid endografting in endovascular abdominal aortic aneurysm repair (EVAR) is defined as the process of placing a series of two or more different types of covered stents, usually to treat a complex abdominal aortic aneurysm (AAA) or a primary or secondary endoleak. We describe the treatment of a type III, a type Ib, and a type Ia endoleak in three patients respectively, using hybrid solutions, assembling components from different manufacturers. An update of the current clinical and experimental evidence on the application of anatomically compatible, hybrid endograft systems in conventional EVAR is also provided.

Abdominal aortic aneurysm: Treatment options, image visualizations and follow-up procedures

Abdominal aortic aneurysm is a common vascular disease that affects elderly population. Open surgical repair is regarded as the gold standard technique for treatment of abdominal aortic aneurysm, however, endovascular aneurysm repair has rapidly expanded since its first introduction in 1990s. As a less invasive technique, endovascular aneurysm repair has been confirmed to be an effective alternative to open surgical repair, especially in patients with co-morbid conditions. Computed tomography (CT) angiography is currently the preferred imaging modality for both preoperative planning and post-operative follow-up. 2D CT images are complemented by a number of 3D reconstructions which enhance the diagnostic applications of CT angiography in both planning and follow-up of endovascular repair. CT has the disadvantage of high cummulative radiation dose, of particular concern in younger patients, since patients require regular imaging follow-ups after endovascular repair, thus, exposing patients to repeated radiation exposure for life. There is a trend to change from CT to ultrasound surveillance of endovascular aneurysm repair. Medical image visualizations demonstrate excellent morphological assessment of aneurysm and stent-grafts, but fail to provide hemodynamic changes caused by the complex stent-graft device that is implanted into the aorta. This article reviews the treatment options of abdominal aortic aneurysm, various image visualization tools, and follow-up procedures with use of different modalities including both imaging and computational fluid dynamics methods. Future directions to improve treatment outcomes in the follow-up of endovascular aneurysm repair are outlined.

Functionality, growth and accelerated aging of tissue engineered living autologous vascular grafts

Living autologous tissue engineered vascular-grafts (TEVGs) with growth-capacity may overcome the limitations of contemporary artificial-prostheses. However, the multi-step in vitro production of TEVGs requires extensive ex vivo cell-manipulations with unknown effects on functionality and quality of TEVGs due to an accelerated biological age of the cells. Here, the impact of biological cell-age and tissue-remodeling capacity of TEVGs in relation to their clinical long-term functionality are investigated. TEVGs were implanted as pulmonary-artery (PA) replacements in juvenile sheep and followed for up to 240 weeks (∼4.5years). Telomere length and telomerase activity were compared amongst TEVGs and adjacent native tissue. Telomerase-activity of in vitro expanded autologous vascular-cells prior to seeding was <5% as compared to a leukemic cell line, indicating biological-aging associated with decreasing telomere-length with each cellular-doubling. Up to 100 weeks, the cells in the TEVGs had consistently shorter telomeres compared to the native counterpart, whereas no significant differences were detectable at 240 weeks. Computed tomography (CT) analysis demonstrated physiological wall-pressures, shear-stresses, and flow-pattern comparable to the native PA. There were no signs of degeneration detectable and continuous native-analogous growth was confirmed by vessel-volumetry. TEVGs exhibit a higher biological age compared to their native counterparts. However, despite of this tissue engineering technology related accelerated biological-aging, growth-capacity and long-term functionality was not compromised. To the contrary, extensive in-vivo remodeling processes with substantial endogenous cellular turnover appears to result in "TEVG rejuvenation" and excellent clinical performance. As these large-animal results can be extrapolated to approximately 20 human years, this study suggests long-term clinical-safety of cardiovascular in vitro tissue engineering and may contribute to safety-criteria as to first-in-man clinical-trials.

Haemodynamic imaging of thoracic stent-grafts by computational fluid dynamics (CFD): presentation of a patient-specific method combining magnetic resonance imaging and numerical simulations

OBJECTIVES

In the last decade, there was been increasing interest in finding imaging techniques able to provide a functional vascular imaging of the thoracic aorta. The purpose of this paper is to present an imaging method combining magnetic resonance imaging (MRI) and computational fluid dynamics (CFD) to obtain a patient-specific haemodynamic analysis of patients treated by thoracic endovascular aortic repair (TEVAR).

METHODS

MRI was used to obtain boundary conditions. MR angiography (MRA) was followed by cardiac-gated cine sequences which covered the whole thoracic aorta. Phase contrast imaging provided the inlet and outlet profiles. A CFD mesh generator was used to model the arterial morphology, and wall movements were imposed according to the cine imaging. CFD runs were processed using the finite volume (FV) method assuming blood as a homogeneous Newtonian fluid.

RESULTS

Twenty patients (14 men; mean age 62.2 years) with different aortic lesions were evaluated. Four-dimensional mapping of velocity and wall shear stress were obtained, depicting different patterns of flow (laminar, turbulent, stenosis-like) and local alterations of parietal stress in-stent and along the native aorta.

CONCLUSIONS

A computational method using a combined approach with MRI appears feasible and seems promising to provide detailed functional analysis of thoracic aorta after stent-graft implantation.

KEY POINTS

• Functional vascular imaging of the thoracic aorta offers new diagnostic opportunities • CFD can model vascular haemodynamics for clinical aortic problems • Combining CFD with MRI offers patient specific method of aortic analysis • Haemodynamic analysis of stent-grafts could improve clinical management and follow-up.

Impact of an aortic nitinol stent graft on flow measurements by time-resolved three-dimensional velocity-encoded MRI

RATIONALE AND OBJECTIVES

Three-dimensional (3D) velocity-encoded cine (VEC) magnetic resonance imaging (MRI) has the potential to quantify 3D hemodynamic aspects known from computational fluid dynamics and to be used to identify hemodynamic risk factors for complications of endovascular aortic repair. The purpose of this study was to investigate the impact of an aortic nickel-titanium (nitinol) stent graft on the accuracy of flow measurements by 3D VEC MRI.

MATERIALS AND METHODS

A pump generated pulsatile aortic flow in an elastic tube phantom mimicking the aorta. Stacked two-dimensional three-directional VEC MRI (stacked-2D-3dir-MRI), 3D three-directional VEC MRI (3D-3dir-MRI), and gold-standard 2D through-plane VEC MRI were applied before and after the insertion of an aortic nitinol stent graft. Six equidistant levels were analyzed twice by the same reader. The percentage difference of the measured flow rate from the gold standard was defined as the parameter of accuracy.

RESULTS

The overall accuracy of in-stent flow measurements related to the gold standard was -5.4% for stacked-2D-3dir-MRI and -4.1% for 3D-3dir-MRI, demonstrating significant overall underestimation compared to the gold standard (P = .016 and P = .013). However, flow measurements with the stent graft were significantly overestimated by 4.1% using stacked-2D-3dir-MRI (P < .001) and by 5.4% using 3D-3dir-MRI (P = .003) compared to identical measurements without the stent graft. In stacked-2D-3dir-MRI, this positive bias was significantly greater at the proximal and distal ends of the stent graft (P = .025). In 3D-3dir-MRI, measurements along the whole length of the stent graft were affected (P = .006). Intraobserver agreement was excellent, with intraclass correlation coefficients of 0.94 for stacked-2D-3dir-MRI (P < .001) and 0.90 for 3D-3dir-MRI (P < .001).

CONCLUSIONS

Flow measurements within an aortic nitinol stent graft by 3D VEC MRI are feasible, but stent grafts may cause a significant positive bias.

Computational fluid dynamics analysis of the effect of plaques in the left coronary artery

This study was to investigate the hemodynamic effect of simulated plaques in left coronary artery models, which were generated from a sample patient's data. Plaques were simulated and placed at the left main stem and the left anterior descending (LAD) to produce at least 60% coronary stenosis. Computational fluid dynamics analysis was performed to simulate realistic physiological conditions that reflect the in vivo cardiac hemodynamics, and comparison of wall shear stress (WSS) between Newtonian and non-Newtonian fluid models was performed. The pressure gradient (PSG) and flow velocities in the left coronary artery were measured and compared in the left coronary models with and without presence of plaques during cardiac cycle. Our results showed that the highest PSG was observed in stenotic regions caused by the plaques. Low flow velocity areas were found at postplaque locations in the left circumflex, LAD, and bifurcation. WSS at the stenotic locations was similar between the non-Newtonian and Newtonian models although some more details were observed with non-Newtonian model. There is a direct correlation between coronary plaques and subsequent hemodynamic changes, based on the simulation of plaques in the realistic coronary models.

Computational fluid dynamics in cardiovascular disease

Computational fluid dynamics (CFD) is a mechanical engineering field for analyzing fluid flow, heat transfer, and associated phenomena, using computer-based simulation. CFD is a widely adopted methodology for solving complex problems in many modern engineering fields. The merit of CFD is developing new and improved devices and system designs, and optimization is conducted on existing equipment through computational simulations, resulting in enhanced efficiency and lower operating costs. However, in the biomedical field, CFD is still emerging. The main reason why CFD in the biomedical field has lagged behind is the tremendous complexity of human body fluid behavior. Recently, CFD biomedical research is more accessible, because high performance hardware and software are easily available with advances in computer science. All CFD processes contain three main components to provide useful information, such as pre-processing, solving mathematical equations, and post-processing. Initial accurate geometric modeling and boundary conditions are essential to achieve adequate results. Medical imaging, such as ultrasound imaging, computed tomography, and magnetic resonance imaging can be used for modeling, and Doppler ultrasound, pressure wire, and non-invasive pressure measurements are used for flow velocity and pressure as a boundary condition. Many simulations and clinical results have been used to study congenital heart disease, heart failure, ventricle function, aortic disease, and carotid and intra-cranial cerebrovascular diseases. With decreasing hardware costs and rapid computing times, researchers and medical scientists may increasingly use this reliable CFD tool to deliver accurate results. A realistic, multidisciplinary approach is essential to accomplish these tasks. Indefinite collaborations between mechanical engineers and clinical and medical scientists are essential. CFD may be an important methodology to understand the pathophysiology of the development and progression of disease and for establishing and creating treatment modalities in the cardiovascular field.

In vitro validation of finite-element model of AAA hemodynamics incorporating realistic outlet boundary conditions

The purpose of this study is to validate numerical simulations of flow and pressure in an abdominal aortic aneurysm (AAA) using phase-contrast magnetic resonance imaging (PCMRI) and an in vitro phantom under physiological flow and pressure conditions. We constructed a two-outlet physical flow phantom based on patient imaging data of an AAA and developed a physical Windkessel model to use as outlet boundary conditions. We then acquired PCMRI data in the phantom while it operated under conditions mimicking a resting and a light exercise physiological state. Next, we performed in silico numerical simulations and compared experimentally measured velocities, flows, and pressures in the in vitro phantom to those computed in the in silico simulations. There was a high degree of agreement in all of the pressure and flow waveform shapes and magnitudes between the experimental measurements and simulated results. The average pressures and flow split difference between experiment and simulation were all within 2%. Velocity patterns showed good agreement between experimental measurements and simulated results, especially in the case of whole-cycle averaged comparisons. We demonstrated methods to perform in vitro phantom experiments with physiological flows and pressures, showing good agreement between numerically simulated and experimentally measured velocity fields and pressure waveforms in a complex patient-specific AAA geometry.

Computation of hemodynamics in the left coronary artery with variable angulations

The purpose of this study was to investigate the hemodynamic effect of variations in the angulations of the left coronary artery, based on simulated and realistic coronary artery models. Twelve models consisting of four realistic and eight simulated coronary artery geometries were generated with the inclusion of left main stem, left anterior descending and left circumflex branches. The simulated models included various coronary artery angulations, namely, 15°, 30°, 45°, 60°, 75°, 90°, 105° and 120°. The realistic coronary angulations were based on selected patient's data with angles ranging from narrow angles of 58° and 73° to wide angles of 110° and 120°. Computational fluid dynamics analysis was performed to simulate realistic physiological conditions that reflect the in vivo cardiac hemodynamics. The wall shear stress, wall shear stress gradient, velocity flow patterns and wall pressure were measured in simulated and realistic models during the cardiac cycle. Our results showed that a disturbed flow pattern was observed in models with wider angulations, and wall pressure was found to reduce when the flow changed from the left main stem to the bifurcated regions, based on simulated and realistic models. A low wall shear stress gradient was demonstrated at left bifurcations with wide angles. There is a direct correlation between coronary angulations and subsequent hemodynamic changes, based on realistic and simulated models. Further studies based on patients with different severities of coronary artery disease are required to verify our results.

A relation between near-wall particle-hemodynamics and onset of thrombus formation in abdominal aortic aneurysms

A novel computational particle-hemodynamics analysis of key criteria for the onset of an intraluminal thrombus (ILT) in a patient-specific abdominal aortic aneurysm (AAA) is presented. The focus is on enhanced platelet and white blood cell residence times as well as their elevated surface-shear loads in near-wall regions of the AAA sac. The generalized results support the hypothesis that a patient's AAA geometry and associated particle-hemodynamics have the potential to entrap activated blood particles, which will play a role in the onset of ILT. Although the ILT history of only a single patient was considered, the modeling and simulation methodology provided allow for the development of an efficient computational tool to predict the onset of ILT formation in complex patient-specific cases.

In vitro validation of flow measurements in an aortic nitinol stent graft by velocity-encoded MRI

PURPOSE

To validate flow measurements within an aortic nickel-titanium (nitinol) stent graft using velocity-encoded cine magnetic resonance imaging (VEC MRI) and to assess intraobserver agreement of repeated flow measurements.

MATERIALS AND METHODS

An elastic tube phantom mimicking the descending aorta was developed with the possibility to insert an aortic nitinol stent graft. Different flow patterns (constant, sinusoidal and pulsatile aortic flow) were applied by a gear pump. A two-dimensional phase-contrast sequence was used to acquire VEC perpendicular cross-sections at six equidistant levels along the phantom. Each acquisition was performed twice with and without stent graft, and each dataset was analysed twice by the same reader. The percental difference of the measured flow volume to the gold standard (pump setting) was defined as the parameter for accuracy. Furthermore, the intraobserver agreement was assessed.

RESULTS

Mean accuracy of flow volume measurements was -0.45±1.63% without stent graft and -0.18±1.45% with stent graft. Slightly lower accuracy was obtained for aortic flow both without (-2.31%) and with (-1.29%) stent graft. Accuracy was neither influenced by the measurement position nor by repeated acquisitions. There was significant intraobserver agreement with an intraclass correlation coefficient of 0.87 (without stent graft, p<0.001) and 0.80 (with stent graft, p<0.001). The coefficient of variance was 0.25% without stent graft and 0.28% with stent graft.

CONCLUSION

This study demonstrated high accuracy and excellent intraobserver agreement of flow measurements within an aortic nitinol stent graft using VEC MRI. VEC MRI may give new insights into the haemodynamic consequences of endovascular aortic repair.