Computational fluid dynamics analysis of surgical adjustment of left ventricular assist device implantation to minimise stroke risk

Influence of the Outflow Graft Angular Position on the Outcomes in Patients With a Left Ventricular Assist Device

This study aimed to explore the potential impact of the angular position of the outflow graft on thromboembolic events and aortic valve regurgitation in people with a left ventricular assist device (LVAD). We analyzed contrast computed tomography (CT) data of patients with LVAD implantation between 2016 and 2021. Three-dimensional reconstructions of the outflow graft and aortic arch were performed to calculate the horizontal (azimuth) angle and vertical (polar) angle, as well as the relative distance between the outflow graft, aortic valve, and brachiocephalic artery. Among 59 patients (median age 57, 68% male), a vertical angle ≥107° correlated significantly with increased cerebrovascular accidents (hazard ratio [HR]: 5.8, 95% confidence interval [CI]: 1.3-26.3, p = 0.022) and gastrointestinal bleeding (HR: 3.4, 95% CI: 1.0-11.2, p = 0.049) during a median 25 month follow-up. No significant differences were found between the vertical angle and aortic valve regurgitation or survival. The horizontal angle and relative distance did not show differences regarding clinical adverse events. This study emphasizes the importance of the LVAD outflow graft angular position to prevent life-threatening thromboembolic events. This study suggests the need for prospective research to further validate these findings.

Development of idealized human aortic models for in vitro and in silico hemodynamic studies

Background

The aorta, a central component of the cardiovascular system, plays a pivotal role in ensuring blood circulation. Despite its importance, there is a notable lack of idealized models for experimental and computational studies.

Objective

This study aims to develop computer-aided design (CAD) models for the idealized human aorta, intended for studying hemodynamics or solid mechanics in both in vitro and in silico settings.

Methods

Various parameters were extracted from comprehensive literature sources to evaluate major anatomical characteristics of the aorta in healthy adults, including variations in aortic arch branches and corresponding dimensions. The idealized models were generated based on averages weighted by the cohort size of each study for several morphological parameters collected and compiled from image-based or cadaveric studies, as well as data from four recruited subjects. The models were used for hemodynamics assessment using particle image velocimetry (PIV) measurements and computational fluid dynamics (CFD) simulations.

Results

Two CAD models for the idealized human aorta were developed, focusing on the healthy population. The CFD simulations, which align closely with the PIV measurements, capture the main global flow features and wall shear stress patterns observed in patient-specific cases, demonstrating the capabilities of the designed models.

Conclusions

The collected statistical data on the aorta and the two idealized aorta models, covering prevalent arch variants known as Normal and Bovine types, are shown to be useful for examining the hemodynamics of the aorta. They also hold promise for applications in designing medical devices where anatomical statistics are needed.

Computational investigation of outflow graft variation impact on hemocompatibility profile in LVADs

BACKGROUND

Hemocompatibility-related adverse events (HRAE) occur commonly in patients with left ventricular assist devices (LVADs) and add to morbidity and mortality. It is unclear whether the outflow graft orientation can impact flow conditions leading to HRAE. This study presents a simulation-based approach using exact patient anatomy from medical images to investigate the influence of outflow cannula orientation in modulating flow conditions leading to HRAEs.

METHODS

A 3D model of a proximal aorta and outflow graft was reconstructed from a computed tomography (CT) scan of an LVAD patient and virtually modified to model multiple cannula orientations (n = 10) by varying polar (cranio-caudal) (n = 5) and off-set (anterior-posterior) (n = 2) angles. Time-dependent computational flow simulations were then performed for each anatomical orientation. Qualitative and quantitative hemodynamics metrics of thrombogenicity including time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), endothelial cell platelet activation potential (ECAP), particle residence time (PRT), and platelet activation potential (PLAP) were analyzed.

RESULTS

Within the simulations performed, endothelial cell activation potential (ECAP) and particle residence time (PRT) were found to be lowest with a polar angle of 85°, regardless of offset angle. However, polar angles that produced parameters at levels least associated with thrombosis varied when the offset angle was changed from 0° to 12°. For offset angles of 0° and 12° respectively, flow shear was lowest at 65° and 75°, time averaged wall shear stress (TAWSS) was highest at 85° and 35°, and platelet activation potential (PLAP) was lowest at 65° and 45°.

CONCLUSION

This study suggests that computational fluid dynamic modeling based on patient-specific anatomy can be a powerful analytical tool when identifying optimal positioning of an LVAD. Contrary to previous work, our findings suggest that there may be an "ideal" outflow cannula for each individual patient based on a CFD-based hemocompatibility profile.

Tailoring left ventricular assist device cannula implantation using coupled multi-scale multi-objective optimization

BACKGROUND

The frequent occurrence of thromboembolic cerebral events continues to limit the widespread implementation of Ventricular Assist Devices (VAD) despite continued advancements in VAD design and anti-coagulation treatments. Recent studies point to the optimal positioning of the outflow graft (OG) as a potential mitigator of post implantation thromboembolism.

OBJECTIVE

This study aims to examine the tailoring of the OG implantation orientation with the goal of minimizing the number of thrombi reaching the cerebral vessels by means of a formal shape optimization scheme incorporated into a multi-scale hemodynamics analysis.

METHODS

A 3-D patient-specific computational fluid dynamics model is loosely coupled in a two-way manner to a 0-D lumped parameter model of the peripheral circulation. A Lagrangian particle-tracking scheme models and tracks thrombi as non-interacting solid spheres. The loose coupling between CFD and LPM is integrated into a geometric shape optimization scheme which aims to optimize an objective function that targets a drop in cerebral embolization, and an overall reduction in particle residence times.

RESULTS

The results elucidate the importance of OG anastomosis orientation and placement particularly in the case that studied particle release from the OG, as a fivefold decrease in cerebral embolization was observed between the optimal and non-optimal implantations. Another case considered particle release from the ventricle and aortic root walls, in which optimal implantation was achieved with a shallow insertion angle. Particle release from all three origins was investigated in the third case, demonstrating that the optimal configurations were generally characterized by VAD flow directed along the central lumen of the aortic arch. Because optimal configurations depended on the anatomic origin of the thrombus, it is important to determine, in clinical studies, the most likely sites of thrombus formation in VAD patients.

Quantitative Assessment of Aortic Hemodynamics for Varying Left Ventricular Assist Device Outflow Graft Angles and Flow Pulsation

Left ventricular assist devices (LVADs) comprise a primary treatment choice for advanced heart failure patients. Treatment with LVAD is commonly associated with complications like stroke and gastro-intestinal (GI) bleeding, which adversely impacts treatment outcomes, and causes fatalities. The etiology and mechanisms of these complications can be linked to the fact that LVAD outflow jet leads to an altered state of hemodynamics in the aorta as compared to baseline flow driven by aortic jet during ventricular systole. Here, we present a framework for quantitative assessment of aortic hemodynamics in LVAD flows realistic human vasculature, with a focus on quantifying the differences between flow driven by LVAD jet and the physiological aortic jet when no LVAD is present. We model hemodynamics in the aortic arch proximal to the LVAD outflow graft, as well as in the abdominal aorta away from the LVAD region. We characterize hemodynamics using quantitative descriptors of flow velocity, stasis, helicity, vorticity and mixing, and wall shear stress. These are used on a set of 27 LVAD scenarios obtained by parametrically varying LVAD outflow graft anastomosis angles, and LVAD flow pulse modulation. Computed descriptors for each of these scenarios are compared against the baseline flow, and a detailed quantitative characterization of the altered state of hemodynamics due to LVAD operation (when compared to baseline aortic flow) is compiled. These are interpreted using a conceptual model for LVAD flow that distinguishes between flow originating from the LVAD outflow jet (and its impingement on the aorta wall), and flow originating from aortic jet during aortic valve opening in normal physiological state.

A Computational Hemodynamics Approach to Left Ventricular Assist Device (LVAD) Optimization Validated in a Large Patient Cohort

With increasing use of left ventricular assist devices (LVAD) it is critical to devise strategies to optimize LVAD speed while controlling mean arterial pressure (MAP) and flow according to patient physiology. The complex interdependency between LVAD speed, MAP, and flow frequently makes optimization difficult under clinical conditions. We propose a method to guide this procedure in silico, narrowing the conditions to test clinically. A computational model of the circulatory network that simulates HF and LVAD support, incorporating LVAD pressure-flow curves was applied retrospectively to anonymized patient hemodynamics data from the University of Washington Medical Center. MAP management on 61 patient-specific computational models with a target of 70 mm Hg, resulting flow for a given LVAD speed was analyzed, and compared to a target output of 5 L/min. Before performing virtual MAP management, 51% had a MAP>70 mm Hg and CO>5 L/min, and 33% had a MAP>70 mm Hg and CO<5 L/min. After changing systemic resistance to meet the MAP target (without adjusting LVAD speed), 84% of cases resulted in CO higher than 5 L/min, with a median CO of 6.79 L/min, using the computational predictive model. Blood pressure management alone is insufficient in meeting both MAP and CO targets, due to the risk of hypervolemia, and requires appropriate LVAD speed optimization to achieve both targets, while preserving right heart health. Such computational tools can narrow down conditions to be tested for each patient, providing significant insight into the pump-patient interplay. LVAD hemodynamic optimization has the potential to reduce complications and improve outcomes.

In-silico analysis of outflow graft implantation orientation and cerebral thromboembolism incidence for full LVAD support

We investigate tailoring cannula implantation angles of left ventricle assist devices (LVAD) to reduce cerebral embolism risk for full LVAD support. We resolve pulsatile hemodynamics with a multi-scale computational fluid dynamics model coupled to a Lagrangian scheme tracking 2-5 mm particles for three cannula implantations. Blood is modeled as non-Newtonian. Cerebral flow distribution is altered depending on anastomosis angle and comparison of means embolization rates between steady and unsteady flow models show that unsteady modeling is more accurate even in the full LVAD support case. Intermediate angle implantation yields lowest cerebral embolization incidence of 11% vs 29% for normal and 36% for shallow implantation.

Computational Fluid Dynamics Study of Cerebral Thromboembolism Risk in Ventricular Assist Device Patients: Effects of Pulsatility and Thrombus Origin

This study investigates the hypothesis that by surgically manipulating the outflow graft (OG) implantation during ventricle assist device placement, it may be possible to reduce the risk of cerebral embolism. We investigate this hypothesis using a computational approach on a patient-specific basis under fully pulsatile hemodynamics with a multiscale computational fluid dynamics model incorporating a coupled Eulerian-Lagrangian scheme that effectively tracks emboli in the fluid domain. Blood is modeled as a non-Newtonian fluid based on the hematocrit level. Preliminary flow analysis shows that depending on the anastomosis angle the left ventricular assist device (LVAD) can enhance the flow to the cerebral circulation by nearly 31%. Z-test results suggest that unsteady-flow modeling ought to be an integral part of any cardiovascular simulation with residual ventricular function. Assuming unsteady-flow conditions, a shallow LVAD outflow graft anastomosis angle is the most optimal if thrombi are released from the aortic-root reducing cerebral embolization incidence to 15.5% and from the ventricle to 17%, while a more pronounced anastomosis angle becomes advantageous when particles originate from the LVAD with an embolization rate of 16.9%. Overall, computations suggest that a pronounced LVAD anastomosis angle is the better implementation. Unsteady modeling is shown to be necessary for the presence of significant antegrade aortic-root flow which induces cyclical flow patterns due to residual pulsatility. On the other hand, depending on thrombus origin and ventricular assist devices (VAD) anastomosis angle there is a strong tradeoff in embolization rates.

A Comprehensive Review of Risk Factor, Mechanism, and Management of Left Ventricular Assist Device-Associated Stroke

The use of left ventricular assist devices (LVADs) has been increasing in the last decade, along with the number of patients with advanced heart failure refractory to medical therapy. Ischemic stroke and intracranial hemorrhage remain the leading causes of morbidity and mortality in LVAD patients. Despite the common occurrence and the significant outcome impact, underlying mechanisms and management strategies of stroke in LVAD patients are controversial. In this article, we review our current knowledge on pathophysiology and risk factors of LVAD-associated stroke, outline the diagnostic approach, and discuss treatment strategies.

Identification of the haemodynamic environment permissive for plaque erosion

Endothelial erosion of atherosclerotic plaques is the underlying cause of approximately 30% of acute coronary syndromes (ACS). As the vascular endothelium is profoundly affected by the haemodynamic environment to which it is exposed, we employed computational fluid dynamic (CFD) analysis of the luminal geometry from 17 patients with optical coherence tomography (OCT)-defined plaque erosion, to determine the flow environment permissive for plaque erosion. Our results demonstrate that 15 of the 17 cases analysed occurred on stenotic plaques with median 31% diameter stenosis (interquartile range 28–52%), where all but one of the adherent thrombi located proximal to, or within the region of maximum stenosis. Consequently, all flow metrics related to elevated flow were significantly increased (time averaged wall shear stress, maximum wall shear stress, time averaged wall shear stress gradient) with a reduction in relative residence time, compared to a non-diseased reference segment. We also identified two cases that did not exhibit an elevation of flow, but occurred in a region exposed to elevated oscillatory flow. Our study demonstrates that the majority of OCT-defined erosions occur where the endothelium is exposed to elevated flow, a haemodynamic environment known to evoke a distinctive phenotypic response in endothelial cells.

Magnetic Resonance Imaging Based Flow Field and Lagrangian Particle Tracking From a Left Ventricular Assist Device

This study explores the optimal left ventricular assist device (LVAD) cannula outflow configuration in a patient-specific replica of the aorta. The volumetric velocity field is measured using phase-contrast magnetic resonance imaging (PC-MRI) under a physiologically relevant steady flow. The effect of the LVAD outflow graft insertion site and anastomosis angle on the transport of embolic particles to cranial vessels is studied by solving the particle equation of motion for spheres in the range of 0.1-1.0 mm using the measured three-dimensional (3D) velocity field. Results show that for a given aorta anatomy, it is possible to design the cannula graft location and terminal curvature so that the probability of embolic transport to the cranial vessels is significantly minimized. This is particularly important since the complex flow pattern in each cannula case affects the embolic trajectories differently, and hence the common assumption that particles distribute by the volumetric flow division does not hold.

Platelet Dynamics and Hemodynamics of Cerebral Aneurysms Treated with Flow-Diverting Stents

Flow-diverting stents (FDS) are used to treat cerebral aneurysms. They promote the formation of a stable thrombus within the aneurysmal sac and, if successful, isolate the aneurysmal dome from mechanical stresses to prevent rupture. Platelet activation, a mechanism necessary for thrombus formation, is known to respond to biomechanical stimuli, particularly to the platelets' residence time and shear stress exposure. Currently, there is no reliable method for predicting FDS treatment outcomes, either a priori or after the procedure. Eulerian computational fluid dynamic (CFD) studies of aneurysmal flow have searched for predictors of endovascular treatment outcome; however, the hemodynamics of thrombus formation cannot be fully understood without considering the platelets' trajectories and their mechanics-triggered activation. Lagrangian analysis of the fluid mechanics in the aneurysmal vasculature provides novel metrics by tracking the platelets' residence time (RT) and shear history (SH). Eulerian and Lagrangian parameters are compared for 19 patient-specific cases, both pre- and post-treatment, to assess the degree of change caused by the FDS and subsequent treatment efficacy.

Embolus Analog Trajectory Paths Under Physiological Flowrates Through Patient-Specific Aortic Arch Models

Atrial fibrillation (AF) is the most common irregular heartbeat among the world's population and is a major contributor to cardiogenic embolisms and acute ischemic stroke (AIS). However, the role AF flow plays in the trajectory paths of cardiogenic emboli has not been experimentally investigated. A physiological simulation system was designed to analyze the trajectory patterns of bovine embolus analogs (EAs) (n = 720) through four patient-specific models, under three flow conditions: steady flow, normal pulsatile flow, and AF pulsatile flow. It was seen that EA trajectory paths were proportional to the percentage flowrate split of 25–31% along the branching vessels. Overall, AF flow conditions increased trajectories through the left- (LCCA) and right (RCCA)-common carotid artery by 25% with respect to normal pulsatile flow. There was no statistical difference in the distribution of clot trajectories when the clot was released from the right, left, or anterior positions. Significantly, more EAs traveled through the brachiocephalic trunk (BCT) than through the LCCA or the left subclavian. Yet of the EAs that traveled through the common carotid arteries, there was a greater affiliation toward the LCCA compared to the RCCA (p < 0.05).

Innovative Modeling Techniques and 3D Printing in Patients with Left Ventricular Assist Devices: A Bridge from Bench to Clinical Practice

Left ventricular assist devices (LVAD) cause altered flow dynamics that may result in complications such as stroke, pump thrombosis, bleeding, or aortic regurgitation. Understanding altered flow dynamics is important in order to develop more efficient and durable pump configurations. In patients with LVAD, hemodynamic assessment is limited to imaging techniques such as echocardiography which precludes detailed assessment of fluid dynamics. In this review article, we present some innovative modeling techniques that are often used in device development or for research purposes, but have not been utilized clinically. Computational fluid dynamic (CFD) modeling is based on computer simulations and particle image velocimetry (PIV) employs ex vivo models that helps study fluid characteristics such as pressure, shear stress, and velocity. Both techniques may help elaborate our understanding of complications that occur with LVAD and could be potentially used in the future to troubleshoot LVAD-related alarms. These techniques coupled with 3D printing may also allow for patient-specific device implants, lowering the risk of complications increasing device durability.

Blood damage in Left Ventricular Assist Devices: Pump thrombosis or system thrombosis?

Introduction:

Despite significant technical advancements in the design and manufacture of Left Ventricular Assist Devices, post-implant thrombotic and thromboembolic complications continue to affect long-term outcomes. Previous efforts, aimed at optimizing pump design as a means of reducing supraphysiologic shear stresses generated within the pump and associated prothrombotic shear-mediated platelet injury, have only partially altered the device hemocompatibility.

Methods:

We examined hemodynamic mechanisms that synergize with hypershear within the pump to contribute to the thrombogenic potential of the overall Left Ventricular Assist Device system.

Results:

Numerical simulations of blood flow in differing regions of the Left Ventricular Assist Device system, that is the diseased native left ventricle, the pump inflow cannula, the impeller, the outflow graft and the anastomosed downstream aorta, reveal that prothrombotic hemodynamic conditions might occur at these specific sites. Furthermore, we show that beyond hypershear, additional hemodynamic abnormalities exist within the pump, which may elicit platelet activation, such as recirculation zones and stagnant platelet trajectories. We also provide evidences that particular Left Ventricular Assist Device implantation configurations and specific post-implant patient management strategies, such as those allowing aortic valve opening, are more hemodynamically favorable and reduce the thrombotic risk.

Conclusion:

We extend the perspective of pump thrombosis secondary to the supraphysiologic shear stress environment of the pump to one of Left Ventricular Assist Device system thrombosis, raising the importance of comprehensive characterization of the different prothrombotic risk factors of the total system as the target to achieve enhanced hemocompatibility and improved clinical outcomes.

Impact of LVAD Implantation Site on Ventricular Blood Stagnation

Treatment of end-stage heart failure includes cardiac transplantation or ventricular assist device (VAD) therapy. Although increasingly prevalent, current VAD therapy has inherent complications, including thrombosis. Studies have demonstrated that VAD implantation alters intracardiac blood flow, creating areas of stagnation that predispose to thrombus formation. Two potential surgical configurations exist for VAD implantation: through the apical or diaphragmatic surfaces of the heart. We hypothesized that diaphragmatic implantation causes more stagnation than apical implantation. We also hypothesized that intermittent aortic valve (AV) opening reduces stagnation of blood inside the left ventricle (LV) when compared with a closed AV. To test these hypotheses, a human LV geometry was recreated in silico and a VAD inflow cannula was virtually implanted in each configuration. A computational indicator-dilution study was conducted where "virtually dyed blood" was washed out of the LV by injecting blood with no dye. Simulations demonstrated a substantial reduction in stagnation with intermittent AV opening. In addition, virtual dye was cleared slightly faster in the apical configuration. Simulations from our study demonstrate the clinical importance of VAD management to allow intermittent opening of the AV to prevent subvalvular stagnation, and also suggests that apical configuration might be more hemodynamically favorable.

Prevention and Treatment of Thrombotic and Hemorrhagic Complications in Patients Supported by Continuous-Flow Left Ventricular Assist Devices

PURPOSE OF REVIEW

The purpose of this review is to describe the current knowledge in prevention and treatment of thrombotic (pump thrombosis and ischemic stroke) and bleeding (gastrointestinal and hemorrhagic stroke) complications in patients supported by continuous-flow left ventricular assist devices (CF-LVAD).

RECENT FINDINGS

Left ventricular assist devices (LVADs) are now widely used for the management of end-stage heart failure. Unfortunately, in spite of the indisputable positive impact LVADs have on patients, the frequency and severity of complications are limitations of this therapy. Stroke, pump thrombosis, and gastrointestinal bleeding are among the most serious and frequent complications in these patients. The balance between hemorrhagic and thrombotic complications in patients supported with CF-LVAD is difficult as most patients do not necessarily fit a "bleeder" or a "clotter" profile but rather move from one side to the other of the thrombotic/bleeding spectrum. Further research is necessary to better understand the risk factors and mechanisms involved in the development of these complications.

LVAD Outflow Graft Angle and Thrombosis Risk

This study quantifies thrombogenic potential (TP) of a wide range of left ventricular assist device (LVAD) outflow graft anastomosis angles through state-of-the-art techniques: 3D imaged-based patient-specific models created via virtual surgery and unsteady computational fluid dynamics with Lagrangian particle tracking. This study aims at clarifying the influence of a single parameter (outflow graft angle) on the thrombogenesis associated with flow patterns in the aortic root after LVAD implantation. This is an important and poorly-understood aspect of LVAD therapy, because several studies have shown strong inter and intrapatient thrombogenic variability and current LVAD implantation strategies do not incorporate outflow graft angle optimization. Accurate platelet-level investigation, enabled by statistical treatment of outliers in Lagrangian particle tracking, demonstrates a strong influence of outflow graft anastomoses angle on thrombogenicity (platelet residence times and activation state characterized by shear stress accumulation) with significantly reduced TP for acutely-angled anastomosed outflow grafts. The methodology presented in this study provides a device-neutral platform for conducting comprehensive thrombogenicity evaluation of LVAD surgical configurations, empowering optimal patient-focused surgical strategies for long-term treatment and care for advanced heart failure patients.

Role of Aortic Geometry on Stroke Propensity based on Simulations of Patient-Specific Models

Stroke is a life threatening event that is expected to more than double over the next 40 years. Atrial fibrillation (AF) has been reported as a strong independent risk factor for stroke. We have previously shown that a hemodynamic perturbation by AF or reduced cardiac output and cycle length may have a significant impact on clot trajectory and thus embolic stroke propensity through the left common carotid artery using an idealized aortic arch model. Here, we show the dependence of flow patterns and hence stroke propensity on geometry of patient-specific aortas. We performed computational fluid dynamics (CFD) simulations to determine the variations of AF-induced stroke propensity over various image-based patient-dependent aorta models. The results demonstrated that curvature pattern of aorta can play a determinant role in AF-induced stroke propensity alteration. Specifically, it was shown that the hemodynamic perturbation by AF considered led to substantial increase in stroke propensity (i.e., 2.5~3.8 fold elevation) for lower curvature angle <90° while the changes in stroke propensity by AF are negligible for higher curvature angle >90°. The present simulations suggest that aortic arch curvature is an important risk factor for embolic stroke which should be tested in future clinical trials.

Effect of Outflow Graft Size on Flow in the Aortic Arch and Cerebral Blood Flow in Continuous Flow Pumps: Possible Relevance to Strokes

One of the most devastating complications of continuous flow left ventricular devices (CFLVADS) is stroke, with a higher incidence in HeartWare Ventricular Assist Device (HVAD) as compared with HEARTMATE II. The reason for the observed difference in stroke rates is unclear. Because outflow graft diameters are different, we hypothesized that this could contribute to the difference in stroke rates. A computational fluid-structure interaction model was created from the computed tomography (CT) scan of a patient. Pressures were used as the boundary condition and the flow through the cerebral vessels was derived as outputs. Flow into the innominate artery was very sensitive to the anastomosis angle for a 10 mm as compared with a 14 mm graft, with the net innominate flow severely compromised with a 10 mm graft at 45° angle. Aortic insufficiency seems to affect cerebral blood flow nonlinearly with an 80% decrease at certain angles of outflow graft anastomosis. Arterial return in to the arch through a narrow graft has important jet effects and results in significant flow perturbations in the aortic arch and cerebral vessels and stasis. A 10 mm graft is more sensitive to angle of insertion than a 14 mm graft. Under some conditions, serious hypoperfusion of the innominate artery is possible. Aortic incompetence results in significant decrease of cerebral blood flow. No stasis was found in the pulsatile flow compared with LVAD flow.

The hemodynamic effects of the LVAD outflow cannula location on the thrombi distribution in the aorta: A primary numerical study

Although a growing number of patients undergo LVAD implantation for heart failure treatment, thrombi are still the devastating complication for patients who used LVAD. LVAD outflow cannula location and thrombi generation sources were hypothesized to affect the thrombi distribution in the aorta. To test this hypothesis, numerical studies were conducted by using computational fluid dynamic (CFD) theory. Two anastomotic configurations, in which the LVAD outflow cannula is anastomosed to the anterior and lateral ascending aortic wall (named as anterior configurations and lateral configurations, respectively), are designed. The particles, whose sized are same as those of thrombi, are released at the LVAD output cannula and the aortic valve (named as thrombiP and thrombiL, respectively) to calculate the distribution of thrombi. The simulation results demonstrate that the thrombi distribution in the aorta is significantly affected by the LVAD outflow cannula location. In anterior configuration, the thrombi probability of entering into the three branches is 23.60%, while that in lateral configuration is 36.68%. Similarly, in anterior configuration, the thrombi probabilities of entering into brachiocephalic artery, left common carotid artery and left subclavian artery, is 8.51%, 9.64%, 5.45%, respectively, while that in lateral configuration it is 11.39%, 3.09%, 22.20% respectively. Moreover, the origins of thrombi could affect their distributions in the aorta. In anterior configuration, the thrombiP has a lower probability to enter into the three branches than thrombiL (12% vs. 25%). In contrast, in lateral configuration, the thrombiP has a higher probability to enter into the three branches than thrombiL (47% vs. 35%). In brief, the LVAD outflow cannula location significantly affects the distribution of thrombi in the aorta. Thus, in the clinical practice, the selection of outflow location of LVAD and the risk of thrombi formed in the left ventricle should be paid more attention than before.

Effect of Different Rotational Directions of BJUT-II VAD on Aortic Swirling Flow Characteristics: A Primary Computational Fluid Dynamics Study

Background

The BJUT-II VAD is a novel left ventricular assist device (LVAD), which is thought to have significant effects on the characteristics of aortic swirling flow. However, the precise mechanism of the rotational direction of BJTU-II VAD in the aortic swirling flow is unclear.

Material/Methods

A patient-specific aortic geometric model was reconstructed based on the CT data. Three pump’s output flow profiles with varied rotational direction, termed “counterclockwise”, “flat profile”, and “clockwise”, were used as the boundary conditions. The helicity density, area-weighted average of helicity density (Ha), localized normalized helicity (LNH), wall shear stress (WSS), and WSS spatial gradient (WSSG) were calculated to evaluate the swirling flow characteristics in the aorta.

Results

The results demonstrated that the swirling flow characteristics in the aorta and 3 branches are directly affected by the output blood flow of BJUT-II VAD. In the aortic arch, the helicity density, supported by the clockwise case, achieved the highest value. In the 3 branches, the flat profile case achieved the highest helicity density, whereas the maximum WSS and WSSG generated by clockwise case were lower than in other cases.

Conclusions

The outflow of the BJUT-II VAD has significant effects on the aortic hemodynamics and swirling flow characteristics. The helical blood profiles can enhance the strength of aortic swirling flow, and reduce the areas of low WSS and WSSG regions. The clockwise case may have a benefit for preventing development of atherosclerosis in the aorta.

Physiological impact of continuous flow on end-organ function: clinical implications in the current era of left ventricular assist devices

The clinical era of continuous-flow left ventricular assist devices has debunked many myths about the dire need of a pulse for human existence. While this therapy has been documented to provide a clear survival benefit in end-stage heart failure patients, we are now faced with certain morbidity challenges that as of yet have no easy mechanistic physiological explanation. The effect of physiological changes on end-organ function in patients supported by continuous-flow ventricular assist devices may offer insight into some of these morbidities. We therefore present a review of current evidence documenting the impact of continuous flow on end-organ function.

Review of recent results using computational fluid dynamics simulations in patients receiving mechanical assist devices for end-stage heart failure

Many end-stage heart failure patients are not eligible to undergo heart transplantation due to organ shortage, and even those under consideration for transplantation might suffer long waiting periods. A better understanding of the hemodynamic impact of left ventricular assist devices (LVAD) on the cardiovascular system is therefore of great interest. Computational fluid dynamics (CFD) simulations give the opportunity to study the hemodynamics in this patient population using clinical imaging data such as computed tomographic angiography. This article reviews a recent study series involving patients with pulsatile and constant-flow LVAD devices in which CFD simulations were used to qualitatively and quantitatively assess blood flow dynamics in the thoracic aorta, demonstrating its potential to enhance the information available from medical imaging.

Computational analysis of pediatric ventricular assist device implantation to decrease cerebral particulate embolization

Stroke is the most devastating complication after ventricular assist device (VAD) implantation with a 19% incidence and 65% mortality in the pediatric population. Current pediatric VAD technology and anticoagulation strategies alone are suboptimal. VAD implantation assisted by computational methods (CFD) may contribute reducing the risk of cerebral embolization. Representative three-dimensional aortic arch models of an infant and a child were generated. An 8 mm VAD outflow-graft (VAD-OG) anastomosed to the aorta was rendered and CFD was applied to study blood flow patterns. Particle tracks, originating in the VAD, were computed with a Lagrangian phase model and the percentage of particles entering the cerebral vessels was calculated. Eight implantation configurations (infant = 5 and child = 3) and 5 particle sizes (0.5, 1, 2, 3, and 4 mm) were considered. For the infant model, percentage of particles entering the cerebral vessels ranged from 15% for a VAD-OG anastomosed at 90° to the aorta, to 31% for 30° VAD-OG anastomosis (overall percentages: X(2) = 10,852, p < 0.0001). For the child model, cerebral embolization ranged from 9% for the 30° VAD-OG anastomosis to 15% for the 60° anastomosis (overall percentages: χ(2) = 10,323, p < 0.0001). Using detailed CFD calculations, we demonstrate that the risk of stroke depends significantly on the VAD implantation geometry. In turn, the risk probably depends on patient-specific anatomy. CFD can be used to optimize VAD implantation geometry to minimize stroke risk.

Thrombus deflector stent for stroke prevention: A simulation study

Atrial fibrillation (AF) is a dysfunction of heart rhythm and represents an increased predisposition to ischemic stroke in AF patients. It has been shown that the AF-induced hemodynamic conditions may contribute to the increased embolic propensity through the carotid arteries. We simulated a stroke-prevention device with a unique strut structure to deflect the trajectory of a blood clot to the carotid artery. We identified the important determinants of functionality in a device design using computational fluid dynamics simulations. Quantitative assessment of deflection efficacy over various clot dimensions was carried out for the device with different strut configurations under AF flow conditions. The simulations demonstrate that the trajectory of a clot destined to the left common carotid artery (LCCA) can be deflected by a strut-structured device at the LCCA inlet with virtually no change in flow resistance. The deflection efficacy of the device is dependent on the clot properties and strut designs of the device. A configuration of 0.75 mm thick and 0.75 mm distant struts with 50% of surface convexity were found to provide maximum deflection efficacy (e.g., 36% greater deflection efficacy than a flat filter) among the strut structures considered. The results suggest that a deflector stent implanted in the aortic branch may be an effective stroke-prevention device. The present simulations motivate pre-clinical animal studies as well as further studies on patient-specific design of the device that maximize the deflection efficacy while minimizing device safety issues.

The study on hemodynamic effect of varied support models of BJUT-II VAD on coronary artery: a primary CFD study

BJUT-II VAD (Beijing University of Technology ventricular assist device II) is a novel left ventricular assist device. Because of the special connection between the pump and native heart, the hemodynamic effects of BJUT-II VAD on coronary artery are still unclear. Hence, numerical simulations have been conducted to clarify changes in hemodynamic effects of different support modes. A patient-specific left coronary arterial geometric model is reconstructed based on the computed tomography (CT) data. Three support modes, "constant speed mode," "co-pulse mode," and "counter pulse mode," are used in this study. The wall shear stress (WSS), wall shear stress gradient (WSSG), cycle-averaged wall shear stress (avWSS), oscillatory shear index (OSI), and the flow pattern are calculated to evaluate the hemodynamic states of coronary artery. The computational results demonstrate that the hemodynamic states of coronary artery are directly affected by the support modes. The co-pulse modes could achieve the highest blood perfusion (constant speed: 153 ml/min vs. co-pulse: 775 ml/min vs. counter pulse: 140 ml/min) and the highest avWSS (constant speed: 18.1 Pa vs. co-pulse: 42.6 Pa vs. counter pulse: 22.6 Pa). In addition, both the WSS and WSSG at the time of peak blood velocity under the constant speed mode are lower than those under other two support modes. In contrast, the counter pulse mode generates the highest OSI value (constant speed: 0.365 vs. co-pulse: 0.379 vs. counter pulse: 0.426). BJUT-II VAD under co-pulse mode may have benefits for improving coronary perfusion and preventing the development of atherosclerosis; however, the constant speed mode may have benefit for preventing the development of plaque vulnerability.

Lagrangian postprocessing of computational hemodynamics

Recent advances in imaging, modeling, and computing have rapidly expanded our capabilities to model hemodynamics in the large vessels (heart, arteries, and veins). This data encodes a wealth of information that is often under-utilized. Modeling (and measuring) blood flow in the large vessels typically amounts to solving for the time-varying velocity field in a region of interest. Flow in the heart and larger arteries is often complex, and velocity field data provides a starting point for investigating the hemodynamics. This data can be used to perform Lagrangian particle tracking, and other Lagrangian-based postprocessing. As described herein, Lagrangian methods are necessary to understand inherently transient hemodynamic conditions from the fluid mechanics perspective, and to properly understand the biomechanical factors that lead to acute and gradual changes of vascular function and health. The goal of the present paper is to review Lagrangian methods that have been used in post-processing velocity data of cardiovascular flows.

Recent advances in computational methodology for simulation of mechanical circulatory assist devices

Ventricular assist devices (VADs) provide mechanical circulatory support to offload the work of one or both ventricles during heart failure. They are used in the clinical setting as destination therapy, as bridge to transplant, or more recently as bridge to recovery to allow for myocardial remodeling. Recent developments in computational simulation allow for detailed assessment of VAD hemodynamics for device design and optimization for both children and adults. Here, we provide a focused review of the recent literature on finite element methods and optimization for VAD simulations. As VAD designs typically fall into two categories, pulsatile and continuous flow devices, we separately address computational challenges of both types of designs, and the interaction with the circulatory system with three representative case studies. In particular, we focus on recent advancements in finite element methodology that have increased the fidelity of VAD simulations. We outline key challenges, which extend to the incorporation of biological response such as thrombosis and hemolysis, as well as shape optimization methods and challenges in computational methodology.

Stroke propensity is increased under atrial fibrillation hemodynamics: a simulation study

Atrial fibrillation (AF) is the most common sustained dysfunction in heart rhythm clinically and has been identified as an independent risk factor for stroke through formation and embolization of thrombi. AF is associated with reduced cardiac output and short and irregular cardiac cycle length. Although the effect of AF on cardiac hemodynamic parameters has been reported, it remains unclear how the hemodynamic perturbations affect the potential embolization of blood clots to the brain that can cause stroke. To understand stroke propensity in AF, we performed computer simulations to describe trajectories of blood clots subject to the aortic flow conditions that represent normal heart rhythm and AF. Quantitative assessment of stroke propensity by blood clot embolism was carried out for a range of clot properties (e.g., 2–6 mm in diameter and 0–0.8 m/s ejection speed) under normal and AF flow conditions. The simulations demonstrate that the trajectory of clot is significantly affected by clot properties as well as hemodynamic waveforms which lead to significant variations in stroke propensity. The predicted maximum difference in stroke propensity in the left common carotid artery was shown to be about 60% between the normal and AF flow conditions examined. The results suggest that the reduced cardiac output and cycle length induced by AF can significantly increase the incidence of carotid embolism. The present simulations motivate further studies on patient-specific risk assessment of stroke in AF.

Aortic outflow cannula tip design and orientation impacts cerebral perfusion during pediatric cardiopulmonary bypass procedures

Poor perfusion of the aortic arch is a suspected cause for peri- and post-operative neurological complications associated with cardiopulmonary bypass (CPB). High-speed jets from 8 to 10FR pediatric/neonatal cannulae delivering ~1 L/min of blood can accrue sub-lethal hemolytic damage while also subjecting the aorta to non-physiologic flow conditions that compromise cerebral perfusion. Therefore, we emphasize the importance of cannulation strategy and hypothesize engineering better CPB perfusion through a redesigned aortic cannula tip. This study employs computational fluid dynamics to investigate novel diffuser-tipped aortic cannulae for shape sensitivity to cerebral perfusion, in an in silico cross-clamped aortic arch model modeled with fixed outflow resistances. 17 parametrically altered configurations of an 8FR end-hole and several diffuser cone angled tips in combination with jet incidence angles toward or away from the head-neck vessels were studied. Experimental pressure-flow characterizations were also conducted on these cannula tip designs. An 8FR end-hole aortic cannula delivering 1 L/min along the transverse aortic arch was found to give rise to backflow from the brachicephalic artery (BCA), irrespective of angular orientation, for the chosen ascending aortic insertion location. Parametric alteration of the cannula tip to include a diffuser cone angle (tested up to 7°) eliminated BCA backflow for any tested angle of jet incidence. Experiments revealed that a 1 cm long 10° diffuser cone tip demonstrated the best pressure-flow performance improvement in contrast with either an end-hole tip or diffuser cone angles greater than 10°. Performance further improved when the diffuser was preceded by an expanded four-lobe swirl inducer attachment-a novel component. In conclusion, aortic cannula orientation is crucial in determining net head-neck perfusion but precise angulations and insertion-depths are difficult to achieve practically. Altering the cannula tip to include a diffuser cone angle has been shown for the first time to have potential in ensuring a net positive outflow at the BCA. Cannula insertion distanced from the BCA inlet may also avoid backflow owing to the Venturi effect, but the diffuser tipped cannula design presents a promising solution to mitigate this issue irrespective of in vivo cannula tip orientation.