Development of a computational fluid dynamic model to investigate the hemodynamic impact of REBOA

INVESTIGATING THE RELATIONSHIP BETWEEN BLEEDING, CLOTTING, AND COAGULOPATHY DURING AUTOMATED PARTIAL REBOA STRATEGIES IN A HIGHLY LETHAL PORCINE HEMORRHAGE MODEL

ABSTRACT

Background: Death due to hemorrhagic shock, particularly, noncompressible truncal hemorrhage, remains one of the leading causes of potentially preventable deaths. Automated partial and intermittent resuscitative endovascular balloon occlusion of the aorta (i.e., pREBOA and iREBOA, respectively) are lifesaving endovascular strategies aimed to achieve quick hemostatic control while mitigating distal ischemia. In iREBOA, the balloon is titrated from full occlusion to no occlusion intermittently, whereas in pREBOA, a partial occlusion is maintained. Therefore, these two interventions impose different hemodynamic conditions, which may impact coagulation and the endothelial glycocalyx layer. In this study, we aimed to characterize the clotting kinetics and coagulopathy associated with iREBOA and pREBOA, using thromboelastography (TEG). We hypothesized that iREBOA would be associated with a more hypercoagulopathic response compared with pREBOA due to more oscillatory flow. Methods: Yorkshire swine (n = 8/group) were subjected to an uncontrolled hemorrhage by liver transection, followed by 90 min of automated pREBOA, iREBOA, or no balloon support (control). Hemodynamic parameters were continuously recorded, and blood samples were serially collected during the experiment (i.e., eight key time points: baseline (BL), T0, T10, T30, T60, T90, T120, T210 min). Citrated kaolin heparinase assays were run on a TEG 5000 (Haemonetics, Niles, IL). General linear mixed models were employed to compare differences in TEG parameters between groups and over time using STATA (v17; College Station, TX), while adjusting for sex and weight. Results: As expected, iREBOA was associated with more oscillations in proximal pressure (and greater magnitudes of peak pressure) because of the intermittent periods of full aortic occlusion and complete balloon deflation, compared to pREBOA. Despite these differences in acute hemodynamics, there were no significant differences in any of the TEG parameters between the iREBOA and pREBOA groups. However, animals in both groups experienced a significant reduction in clotting times (R time: P < 0.001; K time: P < 0.001) and clot strength (MA: P = 0.01; G: P = 0.02) over the duration of the experiment. Conclusions: Despite observing acute differences in peak proximal pressures between the iREBOA and pREBOA groups, we did not observe any significant differences in TEG parameters between iREBOA and pREBOA. The changes in TEG profiles were significant over time, indicating that a severe hemorrhage followed by both pREBOA and iREBOA can result in faster clotting reaction times (i.e., R times). Nevertheless, when considering the significant reduction in transfusion requirements and more stable hemodynamic response in the pREBOA group, there may be some evidence favoring pREBOA usage over iREBOA.

Numerical Study on the Impact of Central Venous Catheter Placement on Blood Flow in the Cavo-Atrial Junction

An in silico study is performed to investigate fluid dynamic effects of central venous catheter (CVC) placement within patient-specific cavo-atrial junctions. Prior studies show the CVC infusing a liquid, but this study focuses on the placement without any liquid emerging from the CVC. A 7 or 15-French double-lumen CVC is placed virtually in two patient-specific models; the CVC tip location is altered to understand its effect on the venous flow field. Results show that the CVC impact is trivial on flow in the superior vena cava when the catheter-to-vein ratio ranges from 0.15 to 0.33. Results further demonstrate that when the CVC tip is directly in the right atrium, flow vortices in the right atrium result in elevated wall shear stress near the tip hole. A recirculation region characterizes a spatially variable flow field inside the CVC side hole. Furthermore, flow stagnation is present near the internal side hole corners but an elevated wall shear stress near the curvature of the side hole's exit. These results suggest that optimal CVC tip location is within the superior vena cava, so as to lower the potential for platelet activation due to elevated shear stresses and that CVC geometry and location depth in the central vein significantly influences the local CVC fluid dynamics. A thrombosis model also shows thrombus formation at the side hole and tip hole. After modifying the catheter design, the hemodynamics change, which alter thrombus formation. Future studies are warranted to study CVC design and placement location in an effort to minimize CVC-induced thrombosis incidence.

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.