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.

Comparison Between Bench-Top and Computational Modelling of Cerebral Thromboembolism in Ventricular Assist Device Circulation.

Despite improvements in ventricular assist devices (VAD) design, VAD-induced stroke rates remain remarkably high at 14-47%. We previously employed computational fluid dynamics (CFD) to propose adjustment of VAD outflow graft (VAD-OG) implantation to reduce stoke. Herein, we present an in-vitro model of cerebral vessel embolization in VAD-assisted circulation, and compare benchtop results to CFD predictions. The benchtop flow-loop consists of a 3D printed aortic bed using Accura 60 polymer driven by a continuous-flow pump. Three hundred spherical particles simulating thrombi of 2, 3.5, and 5 mm diameters were injected at the mock VAD-OG inlet. A water and glycerin mixture (3.8 cP viscosity) synthetically mimicked blood. The flowrate was adjusted to match the CFD Reynolds number. Catch cans were used to capture and count particles reaching cerebral vessels. VAD-OG geometries were evaluated using comparison of means Z-score range of -1.96 ≤ Z ≤ 1.96 to demonstrate overall agreement between computational and in-vitro techniques. Z-scores were: (i) Z = -1.05 for perpendicular (0°), (ii) Z = 0.32 for intermediate (30°), and (iii) Z = -0.52 for shallow (60°) anastomosis and confirmed agreement for all geometries. This study confirmed added benefits of using a left carotid artery bypass-graft with percent embolization reduction: 22.6% for perpendicular, 21.2% for intermediate, and 11.9% for shallow anastomoses. The shallow anastomosis demonstrated lower degrees of aortic arch flow recirculation, consistent with steady-flow computations. Quantitatively and qualitatively, contemporary steady-flow computational models for predicting VAD-induced cerebral embolization can be achieved in-vitro to validate the CFD equivalent.

Mathematical modeling of patient-specific ventricular assist device implantation to reduce particulate embolization rate to cerebral vessels.

Stroke is the most devastating complication after ventricular assist device (VAD) implantation, with an incidence of 14%-47% despite improvements in device design and anticoagulation. This complication continues to limit the widespread implementation of VAD therapy. Patient-specific computational fluid dynamics (CFD) analysis may elucidate ways to reduce this risk. A patient-specific three-dimensional model of the aortic arch was generated from computed tomography. A 12 mm VAD outflow-graft (VAD-OG) "anastomosed" to the aorta was rendered. CFD was applied to study blood flow patterns. Particle tracks, originating from the VAD, were computed with a Lagrangian phase model and percentage of particles entering the cerebral vessels was calculated. Twelve implantation configurations of the VAD-OG and three particle sizes (2, 4, and 5 mm) were considered. Percentage of particles entering the cerebral vessels ranged from 6% for the descending aorta VAD-OG anastomosis, to 14% for the ascending aorta at 90 deg VAD-OG anastomosis. Values were significantly different among all configurations (X(2) = 3925, p < 0.0001). Shallower and more cephalad anastomoses prevented formation of zones of recirculation in the ascending aorta. In this computational model and within the range of anatomic parameters considered, the percentage of particles entering the cerebral vessels from a VAD-OG is reduced by nearly 60% by optimizing outflow-graft configuration. Ascending aorta recirculation zones, which may be thrombogenic, can also be eliminated. CFD methods coupled with patient-specific anatomy may aid in identifying the optimal location and angle for VAD-OG anastomosis to minimize stroke risk.

Mechanics of hip dysplasia reductions in infants using the Pavlik harness: a physics-based computational model.

Biomechanical factors influencing the reduction of dislocated hips with the Pavlik harness in patients of Developmental Dysplasia of the Hip (DDH) were studied using a three-dimensional computer model simulating hip reduction dynamics in (1) subluxated and (2) fully dislocated hip joints. Five hip adductor muscles were identified as key mediators of DDH prognosis, and the non-dimensional force contribution of each in the direction necessary to achieve concentric hip reductions was determined. Results point to the adductor muscles as mediators of subluxated hip reductions, as their mechanical action is a function of the degree of hip dislocation. For subluxated hips in abduction and flexion, the Pectineus, Adductor Brevis, Adductor Longus, and proximal Adductor Magnus contribute positively to reduction, while the rest of the Adductor Magnus contributes negatively. In full dislocations all muscles contribute detrimentally to reduction, elucidating the need for traction to reduce Graf IV type dislocations. Reduction of dysplastic hips was found to occur in two distinct phases: (a) release phase and (b) reduction phase.

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

BACKGROUND: Currently, mechanical support is the most promising alternative to cardiac transplantation. Ventricular assist devices (VADs) were originally used to provide mechanical circulatory support in patients awaiting planned heart transplantation ('bridge-to-transplantation' therapy). The success of short-term bridge devices led to clinical trials evaluating the clinical suitability of long-term support ('destination' therapy) with left ventricular assist devices (LVADs). The first larger scale, randomised trial that tested long-term support with an LVAD reported a 44% reduction in the risk of stroke or death in patients with an LVAD. In spite of the success of LVADs as bridge-to-transplantation and long-term support, patients managed by these devices are still at risk of several adverse events. The most devastating complication is caused by embolisation of thrombi formed within the LVAD or inside the heart into the brain. Prevention of thrombi formation is attempted through anticoagulation management and by improving LVADs design; however, there is still significant occurrence of thromboembolic events in patients. Investigators have reported that the incidence of thromboembolic cerebral events ranges from 14% to 47% over a period of 6-12 months. METHODS AND APPROACH: An alternative method to reduce the incidence of cerebral embolisation is proposed by the co-authors, and the hypothesis is that it is possible to minimise the number of thrombi flowing into the carotid and vertebral arteries by an optimal placement of the LVAD outflow conduit, with or without the addition of aortic bypass connecting the ascending aorta and the innominate artery (IA), or left carotid artery. This paper presents the computational fluid dynamics (CFD) analysis of the aortic arch haemodynamics using a representative geometry of the human aortic arch with or without an alternative aortic bypass. In order to study the trajectory of the thrombi within the aortic arch bed, the CFD code, Fluent 6.3, is utilised to resolve the flow field and to solve the Lagrangian particle tracking of thrombi released randomly at the inlet of the LVAD cannula. RESULTS: Results are presented for simulations of thrombi in the range of 2-5 mm. The percentage of individual diameter as well as aggregate diameter thrombi flowing to the carotid and vertebral arteries as a function of LVAD conduit placement and aortic bypass implantation is reported. The influence of the LVAD conduit implantation and bypass reveals a nearly 50% variation in predicted cerebral embolism rates. CONCLUSIONS: The adjustment of the location of the anastomosis of the LVAD outflow cannula as well as its angle of incidence plays a significant role in the level of thromboembolisms. By proper adjustment in this CFD study of a synthetic model of an aortic arch bed, we found that nearly a 50% reduction in cerebral embolism could be achieved for a configuration consisting of a shallow angle of implantation over a baseline normal incidence of the LVAD cannula. Within the limitations of our model, we have established that the LVAD implantation geometry is an important factor and should be taken into consideration when implanting an LVAD. It is possible that other parameters such as distance of the LVAD outflow cannula to the root of the IA could affect the thrombi embolisation probabilities. However, the results of this study suggest that the risk of stroke may be significantly reduced by as much as 50% by tailoring the VAD implantation by a simple surgical manoeuvre. The results of this line of research may ultimately lead to techniques that can be used to estimate the optimal LVAD configuration in a patient-specific manner by pre-operative imaging.