A Computational Investigation of the Effects of Temporal Synchronization of Left Ventricular Assist Device Speed Modulation with the Cardiac Cycle on Intraventricular Hemodynamics

Patients with advanced heart failure are implanted with a left ventricular assist device (LVAD) as a bridge-to-transplantation or destination therapy. Despite advances in pump design, the risk of stroke remains high. LVAD implantation significantly alters intraventricular hemodynamics, where regions of stagnation or elevated shear stresses promote thrombus formation. Third generation pumps incorporate a pulsatility mode that modulates rotational speed of the pump to enhance in-pump washout. We investigated how the timing of the pulsatility mode with the cardiac cycle affects intraventricular hemodynamic factors linked to thrombus formation. Computational fluid dynamics simulations with Lagrangian particle tracking to model platelet behavior in a patient-specific left ventricle captured altered intraventricular hemodynamics due to LVAD implantation. HeartMate 3 incorporates a pulsatility mode that modulates the speed of the pump every two seconds. Four different timings of this pulsatility mode with respect to the cardiac cycle were investigated. A strong jet formed between the mitral valve and inflow cannula. Blood stagnated in the left ventricular outflow tract beneath a closed aortic valve, in the near-wall regions off-axis of the jet, and in a large counterrotating vortex near the anterior wall. Computational results showed good agreement with particle image velocimetry results. Synchronization of the pulsatility mode with peak systole decreased stasis, reflected in the intraventricular washout of virtual contrast and Lagrangian particles over time. Temporal synchronization of HeartMate 3 pulsatility with the cardiac cycle reduces intraventricular stasis and could be beneficial for decreasing thrombogenicity.

Numerical Modeling of Flow in the Cerebral Vasculature: Understanding Changes in Collateral Flow Directions in the Circle of Willis for a Cohort of Vasospasm Patients Through Image-Based Computational Fluid Dynamics

The Circle of Willis (CoW) is a ring-like network of blood vessels that perfuses the brain. Flow in the collateral pathways that connect major arterial inputs in the CoW change dynamically in response to vessel narrowing or occlusion. Vasospasm is an involuntary constriction of blood vessels following subarachnoid hemorrhage (SAH), which can lead to stroke. This study investigated interactions between localization of vasospasm in the CoW, vasospasm severity, anatomical variations, and changes in collateral flow directions. Patient-specific computational fluid dynamics (CFD) simulations were created for 25 vasospasm patients. Computed tomographic angiography scans were segmented capturing the anatomical variation and stenosis due to vasospasm. Transcranial Doppler ultrasound measurements of velocity were used to define boundary conditions. Digital subtraction angiography was analyzed to determine the directions and magnitudes of collateral flows as well as vasospasm severity in each vessel. Percent changes in resistance and viscous dissipation were analyzed to quantify vasospasm severity and localization of vasospasm in a specific region of the CoW. Angiographic severity correlated well with percent changes in resistance and viscous dissipation across all cerebral vessels. Changes in flow direction were observed in collateral pathways of some patients with localized vasospasm, while no significant changes in flow direction were observed in others. CFD simulations can be leveraged to quantify the localization and severity of vasospasm in SAH patients. These factors as well as anatomical variation may lead to changes in collateral flow directions. Future work could relate localization and vasospasm severity to clinical outcomes like the development of infarct.

Atrial fibrillation increases thrombogenicity of LVAD therapy

BACKGROUND

This study investigates the hypothesis that presence of atrial fibrillation (AF) in LVAD patients increases thrombogenicity in the left ventricle (LV) and exacerbates stroke risk.

METHODS

Using an anatomical LV model implanted with an LVAD inflow cannula, we analyze thrombogenic risk and blood flow patterns in either AF or sinus rhythm (SR) using unsteady computational fluid dynamics (CFD). To analyze platelet activation and thrombogenesis in the LV, hundreds of thousands of platelets are individually tracked to quantify platelet residence time (RT) and shear stress accumulation history (SH).

RESULTS

The irregular and chaotic mitral inflow associated with AF results in markedly different intraventricular flow patterns, with profoundly negative impact on blood flow-induced stimuli experienced by platelets as they traverse the LV. Twice as many platelets accumulated very high SH in the LVAD + AF case, resulting in a 36% increase in thrombogenic potential score, relative to the LVAD + SR case.

CONCLUSIONS

This supports the hypothesis that AF results in unfavorable blood flow patterns in the LV adding to an increased stroke risk for LVAD + AF patients. Quantification of thrombogenic risk associated with AF for LVAD patients may help guide clinical decision-making on interventions to mitigate the increased risk of thromboembolic events.

Surface Reconstruction of the Pediatric Larynx via Structure from Motion Photogrammetry: A Pilot Study

Endoscopy is the gold standard for characterizing pediatric airway disorders, however, it is limited for quantitative analysis due to lack of three-dimensional (3D) vision and poor stereotactic depth perception. We utilize structure from motion (SfM) photogrammetry, to reconstruct 3D surfaces of pathologic and healthy pediatric larynges from monocular two-dimensional (2D) endoscopy. Models of pediatric subglottic stenosis were 3D printed and airway endoscopies were simulated. 3D surfaces were successfully reconstructed from endoscopic videos of all models using an SfM analysis toolkit. Average subglottic surface error between SfM reconstructed surfaces and 3D printed models was 0.65 mm as measured by Modified Hausdorff Distance. Average volumetric similarity between SfM surfaces and printed models was 0.82 as measured by Jaccard Index. SfM can be used to accurately reconstruct 3D surface renderings of the larynx from 2D endoscopy video. This technique has immense potential for use in quantitative analysis of airway geometry and virtual surgical planning.

Hemodynamics of Saline Flushing in Endoscopic Imaging of Partially Occluded Coronary Arteries

PURPOSE

Intravascular endoscopy can aid in the diagnosis of coronary atherosclerosis by providing direct color images of coronary plaques. The procedure requires a blood-free optical path between the catheter and plaque, and achieving clearance safely remains an engineering challenge. In this study, we investigate the hemodynamics of saline flushing in partially occluded coronary arteries to advance the development of intravascular forward-imaging catheters that do not require balloon occlusion.

METHODS

In-vitro experiments and CFD simulations are used to quantify the influence of plaque size, catheter stand-off distance, saline injection flowrate, and injection orientation on the time required to achieve blood clearance.

RESULTS

Experiments and simulation of saline injection from a dual-lumen catheter demonstrated that flushing times increase both as injection flow rate (Reynolds number) decreases and as the catheter moves distally away from the plaque. CFD simulations demonstrated that successful flushing was achieved regardless of lumen axial orientation in a 95% occluded artery. Flushing time was also found to increase as plaque size decreases for a set injection flowrate, and a lower limit for injection flowrate was found to exist for each plaques size, below which clearance was not achieved. For the three occlusion sizes investigated (90, 95, 97% by area), successful occlusion was achieved in less than 1.2 s. Investigation of the pressure fields developed during injection, highlight that rapid clearance can be achieved while keeping the arterial overpressure to < 1 mmHg.

CONCLUSIONS

A dual lumen saline injection catheter was shown to produce clearance safely and effectively in models of partially occluded coronary arteries. Clearance was achieved across a range of engineering and clinical parameters without the use of a balloon occlusion, providing development guideposts for a fluid injection system in forward-imaging coronary endoscopes.

Standardized viscosity as a source of error in computational fluid dynamic simulations of cerebral aneurysms

BACKGROUND

Computational fluid dynamics (CFD) simulations are a powerful tool for studying cerebral aneurysms, capable of evaluating hemodynamics in a way that is infeasible with imaging alone. However, the difficulty of incorporating patient-specific information and inherent obstacles of in vivo validation have limited the clinical usefulness of CFD of cerebral aneurysms. In this work we investigate the effect of using standardized blood viscosity values in CFD simulations of cerebral aneurysms when compared to simulations of the same aneurysms using patient-specific viscosity values derived from hematocrit measurements.

PURPOSE

The objective of this work is to determine the level of error, on average, that is, caused by using standardized values of viscosity in CFD simulations of cerebral aneurysms. By quantifying this error, we demonstrate the need for incorporating patient-specific viscosity in future CFD investigations of cerebral aneurysms.

METHODS

CFD simulations of forty-one cerebral aneurysms were conducted using patient-specific boundary conditions. For each aneurysm two simulations were conducted, one utilizing patient-specific blood viscosity derived from hematocrit measurements and another using a standardized value for blood viscosity. Hemodynamic parameters such as wall shear stress (WSS), wall shear stress gradient (WSSG), and the oscillatory shear index (OSI) were calculated for each of the simulations for each aneurysm. Paired t-tests for differences in the time-averaged maps of these hemodynamic parameters between standardized and patient-specific viscosity simulations were conducted for each aneurysm. Bland-Altman analysis was used to examine the cohort-wide changes in the hemodynamic parameters. Subjects were broken into two groups, those with higher than standard viscosity and those with lower than standard viscosity. An unpaired t-test was used to compare the percent change in WSS, WSSG, and OSI between patient-specific and standardized viscosity simulations for the two cohorts. The percent changes in hemodynamic parameters were correlated against the direction and magnitude of percent change in viscosity, aneurysm size, and aneurysm location. For all t-tests, a Bonferroni-corrected significance level of 0.0167 was used.

RESULTS

63.2%, 41.5%, and 48.7% of aneurysms showed statistically significant differences between patient-specific and standardized viscosity simulations for WSS, WSSG, and OSI respectively. No statistically significant difference was found in the percent changes in WSS, WSSG, and OSI between the group with higher than standard viscosity and those with lower than standard viscosity, indicating an increase in viscosity can cause either an increase or decrease in each of the hemodynamic parameters. On a study-wide level no significant bias was found in either direction for WSS, WSSG, or OSI between the simulation groups due to the bidirectional effect of changing viscosity. No correlation was found between percent change of viscosity and percent change of WSS, WSSG, or OSI, meaning an after-the-fact correction for patient-specific viscosity is not feasible.

CONCLUSION

Standardizing viscosity values in CFD of cerebral aneurysms has a large and unpredictable impact on the calculated WSS, WSSG, and OSI when compared to CFD simulations of the same aneurysms using a patient-specific viscosity. We recommend implementing hematocrit-based patient-specific blood viscosity values for all CFD simulations of cerebral aneurysms.

Fusiform versus Saccular Intracranial Aneurysms-Hemodynamic Evaluation of the Pre-Aneurysmal, Pathological, and Post-Interventional State

Minimally-invasive therapies are well-established treatment methods for saccular intracranial aneurysms (SIAs). Knowledge concerning fusiform IAs (FIAs) is low, due to their wide and alternating lumen and their infrequent occurrence. However, FIAs carry risks like ischemia and thus require further in-depth investigation. Six patient-specific IAs, comprising three position-identical FIAs and SIAs, with the FIAs showing a non-typical FIA shape, were compared, respectively. For each model, a healthy counterpart and a treated version with a flow diverting stent were created. Eighteen time-dependent simulations were performed to analyze morphological and hemodynamic parameters focusing on the treatment effect (TE). The stent expansion is higher for FIAs than SIAs. For FIAs, the reduction in vorticity is higher (Δ35-75% case 2/3) and the reduction in the oscillatory velocity index is lower (Δ15-68% case 2/3). Velocity is reduced equally for FIAs and SIAs with a TE of 37-60% in FIAs and of 41-72% in SIAs. Time-averaged wall shear stress (TAWSS) is less reduced within FIAs than SIAs (Δ30-105%). Within this study, the positive TE of FDS deployed in FIAs is shown and a similarity in parameters found due to the non-typical FIA shape. Despite the higher stent expansion, velocity and vorticity are equally reduced compared to identically located SIAs.

Towards Prediction of Blood Flow in Coiled Aneurysms Before Treatment: A Porous Media Approach

Modeling blood flow in aneurysms treated with coils could be used to understand the complete embolization of the aneurysm, through thrombus formation that fills the entire sac. Modeling of the endovascular coil mass as a porous medium is a technique that allows for study of aneurysm hemodynamics, efficiently for patient-specific treatment outcome predictions. Models in the literature use mean porosity of coils in the aneurysmal volume, proving inadequate for outcome prediction. However, models that consider heterogeneous porosity distribution have shown more accurate hemodynamics. We recently published the porous crown model, considering the heterogeneous coil mass distribution, validated on two patients. This study aims (i) to validate the porous crown model for a larger cohort (eight patients), and (ii) to propose a porous medium model translatable to clinical practice in treatment planning. We analyzed the porosity distribution of the endovascular coils deployed inside the cerebral aneurysm phantoms of eight patients using 3D x-ray synchrotron images. The permeability and inertial factor of the porous crown model are calculated using previously published methodology. We propose a new "bilinear" porous model, that uses the same hypothesis, but the permeability and inertial factor can be defined from just basic information available in the neuro-suite, i.e., the aneurysmal sac volume and the coil volume fraction targeted by the neurosurgeon. These two models are compared to the coil-resolved simulations, considered as the gold standard. The results show that both the porous crown model and the bilinear model produce similarly accurate hemodynamics in the aneurysm. The error in the standard (mean porosity) porous model is 66%, whereas the error of the bilinear model is 26%, compared to the coil-resolved. The bilinear model is promising as a means of treatment outcome prediction at time of intervention.

Solar-powered shape-changing origami microfliers

Using wind to disperse microfliers that fall like seeds and leaves can help automate large-scale sensor deployments. Here, we present battery-free microfliers that can change shape in mid-air to vary their dispersal distance. We designed origami microfliers using bistable leaf-out structures and uncovered an important property: A simple change in the shape of these origami structures causes two dramatically different falling behaviors. When unfolded and flat, the microfliers exhibit a tumbling behavior that increases lateral displacement in the wind. When folded inward, their orientation is stabilized, resulting in a downward descent that is less influenced by wind. To electronically transition between these two shapes, we designed a low-power electromagnetic actuator that produces peak forces of up to 200 millinewtons within 25 milliseconds while powered by solar cells. We fabricated a circuit directly on the folded origami structure that includes a programmable microcontroller, a Bluetooth radio, a solar power-harvesting circuit, a pressure sensor to estimate altitude, and a temperature sensor. Outdoor evaluations show that our 414-milligram origami microfliers were able to electronically change their shape mid-air, travel up to 98 meters in a light breeze, and wirelessly transmit data via Bluetooth up to 60 meters away, using only power collected from the sun.

Improving the accuracy of computational fluid dynamics simulations of coiled cerebral aneurysms using finite element modeling

Cerebral aneurysms are a serious clinical challenge, with ∼half resulting in death or disability. Treatment via endovascular coiling significantly reduces the chances of rupture, but the techniquehas failure rates of ∼20 %. This presents a pressing need to develop a method fordetermining optimal coildeploymentstrategies. Quantification of the hemodynamics of coiled aneurysms using computational fluid dynamics (CFD) has the potential to predict post-treatment outcomes, but representing the coil mass in CFD simulations remains a challenge. We use the Finite Element Method (FEM) for simulating patient-specific coil deployment for n = 4 ICA aneurysms for which 3D printed in vitro models were also generated, coiled, and scanned using ultra-high resolution synchrotron micro-CT. The physical and virtual coil geometries were voxelized onto a binary structured grid and porosity maps were generated for geometric comparison. The average binary accuracy score is 0.8623 and the average error in porosity map is 4.94 %. We then conduct patient-specific CFD simulations of the aneurysm hemodynamics using virtual coils geometries, micro-CT generated oil geometries, and using the porous medium method to represent the coil mass. Hemodynamic parameters including Neck Inflow Rate (Qneck) and Wall Shear Stress (WSS) were calculated for each of the CFD simulations. The average relative error in Qneck and WSS from CFD using FEM geometry were 6.6 % and 21.8 % respectively, while the error from CFD using a porous media approximation resulted in errors of 55.1 % and 36.3 % respectively; demonstrating a marked improvement in the accuracy of CFD simulations using FEM generated coil geometries.

In Vitro Investigation of the Effect of the Timing of Left Ventricular Assist Device Speed Modulation on Intraventricular Flow Patterns

Thromboembolic events remain a common complication for left ventricular assist device (LVAD) patients. To prevent in-pump thrombosis, third-generation LVADs use speed modulation, which is not synchronized with the native left ventricle (LV) contractility. This study aims to investigate the effect of speed modulation on intraventricular flow patterns, and specifically, the impact of timing relative to pressure variations in the LV. Stereo-particle image velocimetry measurements were performed in a patient-derived LV implanted with an LVAD, for different timings of the speed modulation and speed. Speed modulation has a strong effect on instantaneous afterload and flowrate (-16% and +20%). The different timings of the speed modulation resulted in different flowrate waveforms, exhibiting different maxima (5.3-5.9 L/min, at constant average flowrate). Moreover, the timing of the speed modulation was found to strongly influence intraventricular flow patterns, specifically, stagnation areas within the LV. These experiments highlight, once more, the complex relationship between LVAD speed, hemodynamic resistance, and intraventricular pressure. Overall, this study demonstrates the importance of considering native LV contractility in future LVAD controls, to improve hemocompatibility and reduce the risk of thromboembolic complications.

A Novel Patient-Specific Computational Fluid Dynamics Study of the Activation of Primary Collateral Pathways in the Circle of Willis During Vasospasm

The Circle of Willis (CoW) is a redundant network of blood vessels that perfuses the brain. The ringlike anatomy mitigates the negative effects of stroke by activating collateral pathways that help maintain physiological perfusion. Previous studies have investigated the activation of these pathways during embolic stroke and internal carotid artery occlusion. However, the role of collateral pathways during cerebral vasospasm-an involuntary constriction of blood vessels after subarachnoid hemorrhage-is not well-documented. This study presents a novel technique to create patient-specific computational fluid dynamics (CFD) simulations of the Circle of Willis before and during vasospasm. Computed tomographic angiography (CTA) scans are segmented to model the vasculature, and transcranial Doppler ultrasound (TCD) measurements of blood flow velocity are applied as boundary conditions. Bayesian analysis leverages information about the uncertainty in the measurements of vessel diameters and velocities to find an optimized parameter set that satisfies mass conservation and that is applied in the final simulation. With this optimized parameter set, the diameters, velocities, and flow rates fall within typical literature values. Virtual angiograms modeled using passive scalar transport agree closely with clinical angiography. A sensitivity analysis quantifies the changes in collateral flow rates with respect to changes in the inlet and outlet flow rates. This analysis can be applied in the future to a cohort of patients to investigate the relationship between the locations and severities of vasospasm, the patient-to-patient anatomical variability in the Circle of Willis, and the activation of collateral pathways.

Modeling the Mechanical Microenvironment of Coiled Cerebral Aneurysms

Successful occlusion of cerebral aneurysms using coil embolization is contingent upon stable thrombus formation, and the quality of the thrombus depends upon the biomechanical environment. The goal of this study was to investigate how coil embolization alters the mechanical micro-environment within the aneurysm dome. Inertialess particles were injected in three-dimensional, computational simulations of flow inside patient aneurysms using patient-specific boundary conditions. Coil embolization was simulated as a homogenous porous medium of known permeability and inertial constant. Lagrangian particle tracking was used to calculate the residence time and shear stress history for particles in the flow before and after treatment. The percentage of particles entering the aneurysm dome correlated with the neck surface area before and after treatment (pretreatment: R2 = 0.831, P < 0.001; post-treatment: R2 = 0.638, P < 0.001). There was an inverse relationship between the change in particles entering the dome and coil packing density (R2 = 0.600, P < 0.001). Following treatment, the particles with the longest residence times tended to remain within the dome even longer while accumulating lower shear stress. A significant correlation was observed between the treatment effect on residence time and the ratio of the neck surface area to porosity (R2 = 0.390, P = 0.007). The results of this study suggest that coil embolization triggers clot formation within the aneurysm dome via a low shear stress-mediated pathway. This hypothesis links independently observed findings from several benchtop and clinical studies, furthering our understanding of this treatment strategy.

Dialysate Regeneration via Urea Photodecomposition with TiO2 Nanowires at Therapeutic Rates

BACKGROUND

The standard weekly treatment for end-stage renal disease patients is three 4 h long hemodialysis sessions with each session consuming over 120 L of clean dialysate, which prevents the development of portable or continuous ambulatory dialysis treatments. The regeneration of a small (~1 L) amount of dialysate would enable treatments that give conditions close to continuous hemostasis and improve patient quality of life through mobility.

METHODS

Small scale studies have shown that nanowires of TiO₂ are highly efficient at photodecomposing urea into CO₂ and N₂ when using an applied bias and an air permeable cathode. To enable the demonstration of a dialysate regeneration system at therapeutically useful rates, a scalable microwave hydrothermal synthesis of single crystal TiO₂ nanowires grown directly from conductive substrates was developed. These were incorporated into an 1,810 cm² flow channel arrays. Regenerated dialysate samples were treated with activated carbon (2 min at 0.2 g/mL).

RESULTS

The photodecomposition system achieved the therapeutic target of 14.2 g urea removal in 24 h. TiO₂ electrode had a high urea removal photocurrent efficiency of 91%, with less than 1% of the decomposed urea generating NH₄ ⁺ (1.04 μg/h/cm² ), 3% generating NO₃ ^- and 0.5% generating chlorine species. Activated carbon treatment could reduce total chlorine concentration from 0.15 mg/L to <0.02 mg/L. The regenerated dialysate showed significant cytotoxicity which could be removed by treatment with activated carbon. Additionally, a forward osmosis membrane with sufficient urea flux can cut off the mass transfer of the by-products back into the dialysate.^1,2 CONCLUSION: Urea could be removed from spent dialysate at a therapeutic rate using a TiO₂ based photooxidation unit, which can enable portable dialysis systems.

Modeling Flow in Cerebral Aneurysm After Coils Embolization Treatment: A Realistic Patient-Specific Porous Model Approach

PURPOSE

Computational fluid dynamics (CFD) has been used to evaluate the efficiency of endovascular treatment in coiled cerebral aneurysms. The explicit geometry of the coil mass cannot typically be incorporated into CFD simulations since the coil mass cannot be reconstructed from clinical images due to its small size and beam hardening artifacts. The existing methods use imprecise porous medium representations. We propose a new porous model taking into account the porosity heterogeneity of the coils deployed in the aneurysm.

METHODS

The porosity heterogeneity of the coil mass deployed inside two patients' cerebral aneurysm phantoms is first quantified based on 3D X-ray synchrotron images. These images are also used to compute the permeability and the inertial factor arising in porous models. A new homogeneous porous model (porous crowns model), considering the coil's heterogeneity, is proposed to recreate the flow within the coiled aneurysm. Finally, the validity of the model is assessed through comparisons with coil-resolved simulations.

RESULTS

The strong porosity gradient of the coil measured close to the aneurysmal wall is well captured by the porous crowns model. The permeability and the inertial factor values involved in this model are closed to the ideal homogeneous porous model leading to a mean velocity in the aneurysmal sac similar as in the coil-resolved model.

CONCLUSION

The porous crowns model allows for an accurate description of the mean flow within the coiled cerebral aneurysm.

Endothelial Responses to Curvature-Induced Flow Patterns in Engineered Cerebral Aneurysms

Hemodynamic factors have long been associated with clinical outcomes in the treatment of cerebral aneurysms. Computational studies of cerebral aneurysm hemodynamics have provided valuable estimates of the mechanical environment experienced by the endothelium in both the parent vessel and aneurysmal dome walls and have correlated them with disease state. These computational-clinical studies have recently been correlated with the response of endothelial cells (EC) using either idealized or patient-specific models. Here, we present a robust workflow for generating anatomic-scale aneurysm models, establishing luminal cultures of ECs at physiological relevant flow profiles, and comparing EC responses to curvature mediated flow. We show that flow patterns induced by parent vessel curvature produce changes in wall shear stress (WSS) and wall shear stress gradients (WSSG) that are correlated with differences in cell morphology and cellular protein localization. Cells in higher WSS regions align better with the flow and display strong Notch1-extracellular domain (ECD) polarization, while, under low WSS, differences in WSSG due to curvature change were associated with less alignment and attenuation of Notch1-ECD polarization in ECs of the corresponding regions. These proof-of-concept results highlight the use of engineered cellularized aneurysm models for connecting computational fluid dynamics to the underlying endothelial biology that mediates disease.

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.

The effect of Dean, Reynolds, and Womersley number on the flow in a spherical cavity on a curved round pipe. Part 1. Fluid mechanics in the cavity as a canonical flow representing intracranial aneurysms

Flow in side-wall cerebral aneurysms can be ideally modelled as the combination of flow over a spherical cavity and flow in a curved circular pipe, two canonical flows. Flow in a curved pipe is known to depend on the Dean number De, combining the effects of Reynolds number, Re, and of the curvature along the pipe centreline, κ. Pulsatility in the flow introduces a dependency on the Womersley number Wo. Using stereo PIV measurements, this study investigated the effect of these three key non-dimensional parameters, by modifying pipe curvature (De), flow-rate (Re), and pulsatility frequency (Wo), on the flow patterns in a spherical cavity. A single counter-rotating vortex was observed in the cavity for all values of pipe curvature κ and Re, for both steady and pulsatile inflow conditions. Increasing the pipe curvature impacted both the flow patterns in the pipe and the cavity, by shifting the velocity profile towards the cavity opening and increasing the flow rate into the cavity. The circulation in the cavity was found to collapse well with only the Dean number, for both steady and pulsatile inflows. For pulsatile inflow, the counter-rotating vortex was unstable and the location of its centre over time was impacted by the curvature of the pipe, as well as the Re and the Wo in the freestream. The circulation in the cavity was higher for steady inflow than for the equivalent average Reynolds and Dean number pulsatile inflow, with very limited impact of the Womersley in the range studied.

The effect of Dean, Reynolds and Womersley numbers on the flow in a spherical cavity on a curved round pipe. Part 2. The haemodynamics of intracranial aneurysms treated with flow-diverting stents

The flow in a spherical cavity on a curved round pipe is a canonical flow that describes well the flow inside a sidewall aneurysm on an intracranial artery. Intracranial aneurysms are often treated with a flow-diverting stent (FDS), a low-porosity metal mesh that covers the entrance to the cavity, to reduce blood flow into the aneurysm sac and exclude it from mechanical stresses imposed by the blood flow. Successful treatment is highly dependent on the degree of reduction of flow inside the cavity, and the resulting altered fluid mechanics inside the aneurysm following treatment. Using stereoscopic particle image velocimetry, we characterize the fluid mechanics in a canonical configuration representative of an intracranial aneurysm treated with a FDS: a spherical cavity on the side of a curved round pipe covered with a metal mesh formed by an actual medical FDS. This porous mesh coverage is the focus of Part 2 of the paper, characterizing the effects of parent vessel Re, De and pulsatility, Wo, on the fluid dynamics, compared with the canonical configuration with no impediments to flow into the cavity that is described in Part 1 (Chassagne et al., J. Fluid Mech., vol. 915, 2021, A123). Coverage with a FDS markedly reduces the flow Re in the aneurysmal cavity, creating a viscous-dominated flow environment despite the parent vessel Re > 100. Under steady flow conditions, the topology that forms inside the cavity is shown to be a function of the parent vessel De. At low values of De, flow enters the cavity at the leading edge and remains attached to the wall before exiting at the trailing edge, a novel behaviour that was not found under any conditions of the high-Re, unimpeded cavity flow described in Part 1. Under these conditions, flow in the cavity co-rotates with the direction of the free-stream flow, similar to Stokes flow in a cavity. As De increases, the flow along the leading edge begins to separate, and the recirculation zone grows with increasing De, until, above De ≈ 180, the flow inside the cavity is fully recirculating, counter-rotating with respect to the free-stream flow. Under pulsatile flow conditions, the vortex inside the cavity progresses through the same cycle - switching from attached and co-rotating with the free-stream flow at the beginning of the cycle (low velocity and positive acceleration) to separated and counter-rotating as De reaches a critical value. The location of separation within the harmonic cycle is shown to be a function of both De and Wo. The values of aneurysmal cavity Re based on both the average velocity and the circulation inside the cavity are shown to increase with increasing values of De, while Wo is shown to have little influence on the time-averaged metrics. As De increases, the strength of the secondary flow in the parent vessel grows, due to the inertial instability in the curved pipe, and the flow rate entering the cavity increases. Thus, the effectiveness of FDS treatment to exclude the aneurysmal cavity from the haemodynamic stresses is compromised for aneurysms located on high-curvature arteries, i.e. vessels with high De, and this can be a fluid mechanics criterion to guide treatment selection.

Lagrangian Trajectory Simulation of Platelets and Synchrotron Microtomography Augment Hemodynamic Analysis of Intracranial Aneurysms Treated With Embolic Coils

As frequency of endovascular treatments for intracranial aneurysms increases, there is a growing need to understand the mechanisms for coil embolization failure. Computational fluid dynamics (CFD) modeling often simplifies modeling the endovascular coils as a homogeneous porous medium (PM), and focuses on the vascular wall endothelium, not considering the biomechanical environment of platelets. These assumptions limit the accuracy of computations for treatment predictions. We present a rigorous analysis using X-ray microtomographic imaging of the coils and a combination of Lagrangian (platelet) and Eulerian (endothelium) metrics. Four patient-specific, anatomically accurate in vitro flow phantoms of aneurysms are treated with the same patient-specific endovascular coils. Synchrotron tomography scans of the coil mass morphology are obtained. Aneurysmal hemodynamics are computationally simulated before and after coiling, using patient-specific velocity/pressure measurements. For each patient, we analyze the trajectories of thousands of platelets during several cardiac cycles, and calculate residence times (RTs) and shear exposure, relevant to thrombus formation. We quantify the inconsistencies of the PM approach, comparing them with coil-resolved (CR) simulations, showing the under- or overestimation of key hemodynamic metrics used to predict treatment outcomes. We fully characterize aneurysmal hemodynamics with converged statistics of platelet RT and shear stress history (SH), to augment the traditional wall shear stress (WSS) on the vascular endothelium. Incorporating microtomographic scans of coil morphology into hemodynamic analysis of coiled intracranial aneurysms, and augmenting traditional analysis with Lagrangian platelet metrics improves CFD predictions, and raises the potential for understanding and clinical translation of computational hemodynamics for intracranial aneurysm treatment outcomes.

Shear stress associated with cardiopulmonary bypass induces expression of inflammatory cytokines and necroptosis in monocytes

Cardiopulmonary bypass (CPB) is required during most cardiac surgeries. CBP drives systemic inflammation and multiorgan dysfunction that is especially severe in neonatal patients. Limited understanding of molecular mechanisms underlying CPB-associated inflammation presents a significant barrier to improve clinical outcomes. To better understand these clinical issues, we performed mRNA sequencing on total circulating leukocytes from neonatal patients undergoing CPB. Our data identify myeloid cells, particularly monocytes, as the major cell type driving transcriptional responses to CPB. Furthermore, IL-8 and TNF-α were inflammatory cytokines robustly upregulated in leukocytes from both patients and piglets exposed to CPB. To delineate the molecular mechanism, we exposed THP-1 human monocytic cells to CPB-like conditions, including artificial surfaces, high shear stress, and cooling/rewarming. Shear stress was found to drive cytokine upregulation via calcium-dependent signaling pathways. We also observed that a subpopulation of THP-1 cells died via TNF-α-mediated necroptosis, which we hypothesize contributes to post-CPB inflammation. Our study identifies a shear stress-modulated molecular mechanism that drives systemic inflammation in pediatric CPB patients. These are also the first data to our knowledge to demonstrate that shear stress causes necroptosis. Finally, we observe that calcium and TNF-α signaling are potentially novel targets to ameliorate post-CPB inflammation.

Computational fluid dynamics modeling aiding surgical planning in a toddler with Parkes Weber syndrome

Parkes Weber syndrome is a fast-flow and slow-flow vascular anomaly with limb overgrowth that can lead to congestive heart failure and limb ischemia. Current management strategies have focused on symptom management with focal embolization. A pediatric case with early onset heart failure is reported. We discuss the use of computational fluid dynamics (CFD) modeling to guide a surgical management strategy in a toddler with an MAP2K1 mutation.

Accuracy of Doppler blood pressure measurement in HeartMate 3 ventricular assist device patients

Aims

Optimal blood pressure (BP) control is imperative to reduce complications, especially strokes, in continuous flow ventricular assist device (VAD) patients. Doppler BP has been shown to be an accurate and reliable non‐invasive BP measurement method in HeartMate II and HVAD patients. We examined whether Doppler BP is also accurate in patients with the HeartMate 3 VAD.

Methods and results

In a prospective, longitudinal cohort of HeartMate 3 patients, arterial line BP and simultaneously measured Doppler opening pressure were obtained. Correlation and agreement between Doppler opening pressure and arterial line mean arterial pressure (MAP) versus systolic blood pressure (SBP) were analysed, as well as the effect of pulse pressure on the accuracy of Doppler opening pressure.

A total of 589 pairs of simultaneous Doppler opening pressure and arterial line pressure readings were obtained in 43 patients. Doppler opening pressure had good correlation with intra‐arterial MAP (r = 0.754) and more closely approximated MAP than SBP (mean error 2.0 vs. −8.6 mmHg). Pulse pressure did not have a clinically significant impact on the accuracy of the Doppler BP method. These results in HeartMate 3 patients are very similar to previous results in HeartMate II and HVAD patients.

Conclusions

Doppler BP method should be the default non‐invasive BP measurement method in continuous flow VAD patients including patients implanted with the HeartMate 3.

Small Left Ventricular Size Is an Independent Risk Factor for Ventricular Assist Device Thrombosis

The prevalence of ventricular assist device (VAD) therapy has continued to increase due to a stagnant donor supply and growing advanced heart failure (HF) population. We hypothesize that left ventricular (LV) size strongly influences biocompatibility and risk of thrombosis. Unsteady computational fluid dynamics (CFD) was used in conjunction with patient-derived computational modeling and virtual surgery with a standard, apically implanted inflow cannula. A dual-focus approach of evaluating thrombogenicity was employed: platelet-based metrics to characterize the platelet environment and flow-based metrics to investigate hemodynamics. Left ventricular end-diastolic dimensions (LVEDds) ranging from 4.5 to 6.5 cm were studied and ranked according to relative thrombogenic potential. Over 150,000 platelets were individually tracked in each LV model over 15 cardiac cycles. As LV size decreased, platelets experienced markedly increased shear stress histories (SHs), whereas platelet residence time (RT) in the LV increased with size. The complex interplay between increased SH and longer RT has profound implications on thrombogenicity, with a significantly higher proportion of platelets in small LVs having long RT times and being subjected to high SH, contributing to thrombus formation. Our data suggest that small LV size, rather than decreased VAD speed, is the primary pathologic mechanism responsible for the increased incidence of thrombosis observed in VAD patients with small LVs.

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.

Left Ventricular Assist Device Inflow Cannula Angle and Thrombosis Risk

BACKGROUND

As heart failure prevalence continues to increase in the setting of a static donor supply, left ventricular assist device (LVAD) therapy for end-stage heart failure continues to grow. Anecdotal evidence suggests that malalignment of the LVAD inflow cannula may increase thrombosis risk, but this effect has not been explored mechanistically or quantified statistically. Our objective is to elucidate the impact of surgical angulation of the inflow cannula on thrombogenicity.

METHODS AND RESULTS

Unsteady computational fluid dynamics is used in conjunction with computational modeling and virtual surgery to model flow through the left ventricle for 5 different inflow cannula angulations. We use a holistic approach to evaluate thrombogenicity: platelet-based (Lagrangian) metrics to evaluate the platelet mechanical environment, combined with flow-based (Eulerian) metrics to investigate intraventricular hemodynamics. The thrombogenic potential of each LVAD inflow cannula angulation is quantitatively evaluated based on platelet shear stress history and residence time. Intraventricular hemodynamics are strongly influenced by LVAD inflow cannula angulation. Platelet behavior indicates elevated thrombogenic potential for certain inflow cannula angles, potentially leading to platelet activation. Our analysis demonstrates that the optimal range of inflow angulation is within 0±7° of the left ventricular apical axis.

CONCLUSIONS

Angulation of the inflow cannula >7° from the apical axis (axis connecting mitral valve and ventricular apex) leads to markedly unfavorable hemodynamics as determined by computational fluid dynamics. Computational hemodynamic simulations incorporating Lagrangian and Eulerian metrics are a powerful tool for studying optimization of LVAD implantation strategies, with the long-term potential of improving outcomes.

Analysis of Factors Influencing Accuracy of Volume Flow Measurement in Dialysis Access Fistulas Based on Duplex Ultrasound Simulation

OBJECTIVE

We developed a duplex ultrasound simulator and used it to assess accuracy of volume flow measurements in dialysis access fistula (DAF) models.

METHODS

The simulator consists of a mannequin, computer, and mock transducer. Each case is built from a patient's B-mode images that are used to create a 3-dimensional surface model of the DAF. Computational fluid dynamics is used to determine blood flow velocities based on model vessel geometry. The simulator displays real-time B-mode and color-flow images, and Doppler spectral waveforms are generated according to user-defined settings. Accuracy was assessed by scanning each case and measuring volume flow in the inflow artery and outflow vein for comparison with true volume flow values.

RESULTS

Four examiners made 96 volume flow measurements on four DAF models. Measured volume flow deviated from the true value by 35 ± 36%. Mean absolute deviation from true volume flow was lower for arteries than veins (22 ± 19%, N = 48 vs. 58 ± 33%, N = 48, p < 0.0001). This finding is attributed to eccentricity of outflow veins which resulted in underestimating true cross-sectional area. Regression analysis indicated that error in measuring cross-sectional area was a predictor of error in volume flow measurement (β = 0.948, p < 0.001). Volume flow error was reduced from 35 ± 36% to 9 ± 8% (p < 0.000001) by calculating vessel area as an ellipse.

CONCLUSIONS

Duplex volume flow measurements are based on a circular vessel shape. DAF inflow arteries are circular, but outflow veins can be elliptical. Simulation-based analysis showed that error in measuring volume flow is mainly due to assumption of a circular vessel.

Accuracy of Doppler blood pressure measurement in continuous-flow left ventricular assist device patients

Aims

Accurate blood pressure (BP) measurement in continuous‐flow ventricular assist device (CF‐VAD) patients is imperative to reduce stroke risk. This study assesses the accuracy of the Doppler opening pressure method compared with the gold standard arterial line method in CF‐VAD patients.

Methods and results

In a longitudinal cohort of HeartMate II and HVAD patients, arterial line BP and simultaneously measured Doppler opening pressure were obtained. Overall correlation, agreement between Doppler opening pressure and arterial line mean vs. systolic pressure, and the effect of arterial pulsatility on the accuracy of Doppler opening pressure were analysed. A total of 1933 pairs of Doppler opening pressure and arterial line pressure readings within 1 min of each other were identified in 154 patients (20% women, mean age 55 ± 15, 50% HeartMate II and 50% HVAD). Doppler opening pressure had good correlation with invasive mean arterial pressure (r = 0.742, P < 0.0001) and more closely approximated mean than systolic BP (mean error 2.4 vs. −8.4 mmHg). Arterial pulsatility did not have a clinically significant effect on the accuracy of the Doppler opening pressure method.

Conclusions

Doppler opening pressure should be the standard non‐invasive method of BP measurement in CF‐VAD patients.

Genetic correlates of wall shear stress in a patient-specific 3D-printed cerebral aneurysm model

OBJECTIVES

To study the correlation between wall shear stress and endothelial cell expression in a patient-specific, three-dimensional (3D)-printed model of a cerebral aneurysm.

MATERIALS AND METHODS

A 3D-printed model of a cerebral aneurysm was created from a patient's angiogram. After populating the model with human endothelial cells, it was exposed to media under flow for 24 hours. Endothelial cell morphology was characterized in five regions of the 3D-printed model using confocal microscopy. Endothelial cells were then harvested from distinct regions of the 3D-printed model for mRNA collection and gene analysis via quantitative polymerase chain reaction (qPCR.) Cell morphology and mRNA measurement were correlated with computational fluid dynamics simulations.

RESULTS

The model was successfully populated with endothelial cells, which survived under flow for 24 hours. Endothelial morphology showed alignment with flow in the proximal and distal parent vessel and aneurysm neck, but disorganization in the aneurysm dome. Genetic analysis of endothelial mRNA expression in the aneurysm dome and distal parent vessel was compared with the proximal parent vessels. ADAMTS-1 and NOS3 were downregulated in the aneurysm dome, while GJA4 was upregulated in the distal parent vessel. Disorganized morphology and decreased ADAMTS-1 and NOS3 expression correlated with areas of substantially lower wall shear stress and wall shear stress gradient in computational fluid dynamics simulations.

CONCLUSIONS

Creating 3D-printed models of patient-specific cerebral aneurysms populated with human endothelial cells is feasible. Analysis of these cells after exposure to flow demonstrates differences in both cell morphology and genetic expression, which correlate with areas of differential hemodynamic stress.

A simplified and versatile calibration method for multi-camera optical systems in 3D particle imaging

This article describes a stereoscopic multi-camera calibration method that does not require any optical model. It is based on a measure of the light propagation within the measurement volume only instead of modeling its entire path up to the sensors. The calibration uses simple plane by plane transformations which allow us to directly link pixel coordinates to light rays. The appeal of the proposed method relies on the combination of its simplicity of implementation (it is particularly easy to apply in any sophisticated optical imaging setup), its versatility (it can easily handle index-of-refraction gradients, as well as complex optical arrangements), and its accuracy {we show that the proposed method gives better accuracy than commonly used techniques, based on Tsai's simple pinhole camera model [R. Tsai, J. Rob. Autom. 3, 323 (1987)], while its numerical implementation remains extremely simple}. Based on ideas that have been available in the fluid mechanics community, this method is a compact turn-key algorithm that can be implemented with open-source routines.

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.

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.

Analytical model of electro-hydrodynamic flow in corona discharge

We present an analytical model for electro-hydrodynamic flow that describes the relationship between the corona voltage, electric field, and ion charge density. The interaction between the accelerated ions and the neutral gas molecules is modeled as an external body force in the Navier-Stokes equation. The gas flow characteristics are solved from conservation principles with spectral methods. This multiphysics model is shown to match experimental data for a point-to-ring corona configuration, shedding new insights into mass, charge, and momentum transport phenomena, and can be readily implemented in any numerical simulation.

Evaluation of Examiner Performance Using a Duplex Ultrasound Simulator. Flow Velocity Measurements in Dialysis Access Fistula Models

We developed a duplex ultrasound simulator for training and assessment of scanning skills. We used the simulator to test examiner performance in the measurement of flow velocities in dialysis access fistulas. Test cases were created from 3-D ultrasound scans of two dialysis access fistulas by reconstructing 3-D blood vessel models and simulating blood flow velocity fields within the lumens. The simulator displays a 2-D B-mode or color Doppler image corresponding to transducer position on a mannequin; a spectral waveform is generated according to Doppler sample volume location and system settings. Examiner performance was assessed by comparing the measured peak systolic velocity (PSV) with the true PSV provided by the computational flow model. The PSV measured by four expert examiners deviated from the true value by 7.8 ± 6.1%. The results indicate the ability of the simulator to objectively assess an examiner's measurement accuracy in complex vascular targets.

Intermittent Aortic Valve Opening and Risk of Thrombosis in Ventricular Assist Device Patients

The current study evaluates quantitatively the impact that intermittent aortic valve (AV) opening has on the thrombogenicity in the aortic arch region for patients under left ventricular assist device (LVAD) therapy. The influence of flow through the AV, opening once every five cardiac cycles, on the flow patterns in the ascending aortic is measured in a patient-derived computed tomography image-based model, after LVAD implantation. The mechanical environment of flowing platelets is investigated, by statistical treatment of outliers in Lagrangian particle tracking, and thrombogenesis metrics (platelet residence times and activation state characterized by shear stress accumulation) are compared for the cases of closed AV versus intermittent AV opening. All hemodynamics metrics are improved by AV opening, even at a reduced frequency and flow rate. Residence times of platelets or microthrombi are reduced significantly by transvalvular flow, as are the shear stress history experienced and the shear stress magnitude and gradients on the aortic root endothelium. The findings of this device-neutral study support the multiple advantages of management that enables AV opening, providing a rationale for establishing this as a standard in long-term treatment and care for advanced heart failure patients.

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.

Computational fluid dynamics of cerebral aneurysm coiling using high-resolution and high-energy synchrotron X-ray microtomography: comparison with the homogeneous porous medium approach

BACKGROUND

Computational modeling of intracranial aneurysms provides insights into the influence of hemodynamics on aneurysm growth, rupture, and treatment outcome. Standard modeling of coiled aneurysms simplifies the complex geometry of the coil mass into a homogeneous porous medium that fills the aneurysmal sac. We compare hemodynamics of coiled aneurysms modeled from high-resolution imaging with those from the same aneurysms modeled following the standard technique, in an effort to characterize sources of error from the simplified model.

MATERIALS

Physical models of two unruptured aneurysms were created using three-dimensional printing. The models were treated with coil embolization using the same coils as those used in actual patient treatment and then scanned by synchrotron X-ray microtomography to obtain high-resolution imaging of the coil mass. Computational modeling of each aneurysm was performed using patient-specific boundary conditions. The coils were modeled using the simplified porous medium or by incorporating the X-ray imaged coil surface, and the differences in hemodynamic variables were assessed.

RESULTS

X-ray microtomographic imaging of coils and incorporation into computational models were successful for both aneurysms. Porous medium calculations of coiled aneurysm hemodynamics overestimated intra-aneurysmal flow, underestimated oscillatory shear index and viscous dissipation, and over- or underpredicted wall shear stress (WSS) and WSS gradient compared with X-ray-based coiled computational fluid dynamics models.

CONCLUSIONS

Computational modeling of coiled intracranial aneurysms using the porous medium approach may inaccurately estimate key hemodynamic variables compared with models incorporating high-resolution synchrotron X-ray microtomographic imaging of complex aneurysm coil geometry.

Surface acoustic wave nebulization device with dual interdigitated transducers improves SAWN-MS performance

We compared mass spectrometric (MS) performance of surface acoustic wave nebulization (SAWN) generated by a single interdigitated transducer (IDT) designed to produce a progressive wave (PW) to one with a dual IDT that can in theory generate standing waves (SW). Given that devices using dual IDTs had been shown to produce fewer large size droplets on average, we hypothesized they would improve MS performance by improving the efficiency of desolvation. Indeed, the SW-SAWN chip provided an improved limit of detection of 1 femtomole of peptide placed on chip making it 100× more sensitive than the PW design. However, as measured by high-speed image recording and phase Doppler particle analyzer measurements, there was only a 26% increase in the small diameter (1-10 µm) droplets produced from the new device, precluding a conclusion that the decrease in droplet size was solely responsible for the improvement in MS signal/noise. Given that the dual IDT design produced a more instantaneous plume than the PW design, the more likely contributor to improved MS signal/noise was concluded to be a higher ion flux entering the mass spectrometer for the dual IDT designs. Notably, the dual IDT device allowed production of much higher quality protein mass spectra up to about 20 kDa, compared with the single IDT device. Copyright © 2016 John Wiley & Sons, Ltd.

Development of a Duplex Ultrasound Simulator and Preliminary Validation of Velocity Measurements in Carotid Artery Models

OBJECTIVE

Duplex ultrasound scanning with B-mode imaging and both color Doppler and Doppler spectral waveforms is relied upon for diagnosis of vascular pathology and selection of patients for further evaluation and treatment. In most duplex ultrasound applications, classification of disease severity is based primarily on alterations in blood flow velocities, particularly the peak systolic velocity (PSV) obtained from Doppler spectral waveforms. We developed a duplex ultrasound simulator for training and assessment of scanning skills.

METHODS

Duplex ultrasound cases were prepared from 2-dimensional (2D) images of normal and stenotic carotid arteries by reconstructing the common carotid, internal carotid, and external carotid arteries in 3 dimensions and computationally simulating blood flow velocity fields within the lumen. The simulator displays a 2D B-mode image corresponding to transducer position on a mannequin, overlaid by color coding of velocity data. A spectral waveform is generated according to examiner-defined settings (depth and size of the Doppler sample volume, beam steering, Doppler beam angle, and pulse repetition frequency or scale). The accuracy of the simulator was assessed by comparing the PSV measured from the spectral waveforms with the true PSV which was derived from the computational flow model based on the size and location of the sample volume within the artery.

RESULTS

Three expert examiners made a total of 36 carotid artery PSV measurements based on the simulated cases. The PSV measured by the examiners deviated from true PSV by 8% ± 5% (N = 36). The deviation in PSV did not differ significantly between artery segments, normal and stenotic arteries, or examiners.

CONCLUSION

To our knowledge, this is the first simulation of duplex ultrasound that can create and display real-time color Doppler images and Doppler spectral waveforms. The results demonstrate that an examiner can measure PSV from the spectral waveforms using the settings on the simulator with a mean absolute error in the velocity measurement of less than 10%. With the addition of cases with a range of pathologies, this duplex ultrasound simulator will be a useful tool for training health-care providers in vascular ultrasound applications and for assessing their skills in an objective and quantitative manner.

Computational Modeling of Venous Sinus Stenosis in Idiopathic Intracranial Hypertension

BACKGROUND AND PURPOSE

Idiopathic intracranial hypertension has been associated with dural venous sinus stenosis in some patients, but the hemodynamic environment of the dural venous sinuses has not been quantitatively described. Here, we present the first such computational fluid dynamics model by using patient-specific blood pressure measurements.

MATERIALS AND METHODS

Six patients with idiopathic intracranial hypertension and at least 1 stenosis or atresia at the transverse/sigmoid sinus junction underwent MR venography followed by cerebral venography and manometry throughout the dural venous sinuses. Patient-specific computational fluid dynamics models were created by using MR venography anatomy, with venous pressure measurements as boundary conditions. Blood flow and wall shear stress were calculated for each patient.

RESULTS

Computational models of the dural venous sinuses were successfully reconstructed in all 6 patients with patient-specific boundary conditions. Three patients demonstrated a pathologic pressure gradient (≥8 mm Hg) across 4 dural venous sinus stenoses. Small sample size precludes statistical comparisons, but average overall flow throughout the dural venous sinuses of patients with pathologic pressure gradients was higher than in those without them (1041.00 ± 506.52 mL/min versus 358.00 ± 190.95 mL/min). Wall shear stress was also higher across stenoses in patients with pathologic pressure gradients (37.66 ± 48.39 Pa versus 7.02 ± 13.60 Pa).

CONCLUSIONS

The hemodynamic environment of the dural venous sinuses can be computationally modeled by using patient-specific anatomy and physiologic measurements in patients with idiopathic intracranial hypertension. There was substantially higher blood flow and wall shear stress in patients with pathologic pressure gradients.

Convective Leakage Makes Heparin Locking of Central Venous Catheters Ineffective Within Seconds: Experimental Measurements in a Model Superior Vena Cava

Central venous catheters (CVCs), placed in the superior vena cava (SVC) for hemodialysis or chemotherapy, are routinely filled while not in use with heparin, an anticoagulant, to maintain patency and prevent thrombus formation at the catheter tip. The heparin-locking procedure, however, places the patient at risk for systemic bleeding, as heparin is known to leak from the catheter into the blood stream. We provide evidence from detailed in vitro experiments that shows the driving mechanism behind heparin leakage to be convective-diffusive transport due to the pulsatile flow surrounding the catheter. This novel mechanism is supported by experimental planar laser-induced fluorescence (PLIF) and particle image velocimetry (PIV) measurements of flow velocity and heparin transport from a CVC placed inside a model SVC inside a pulsatile flow loop. The results predict an initial, fast (<10 s), convection-dominated phase that rapidly depletes the concentration of heparin in the near-tip region, the region of the catheter with side holes. This is followed by a slow, diffusion-limited phase inside the catheter lumen, where the concentration is still high, that is insufficient at replenishing the lost heparin concentration in the near-tip region. The results presented here, which are consistent with previous in vivo estimates of 24 hour leakage rates, predict that the concentration of heparin in the near-tip region is essentially zero for the majority of the interdialytic phase, rendering the heparin locking procedure ineffective.

In vitro validation of endovascular Doppler-derived flow rates in models of the cerebral circulation

This study presents validation of endovascular Doppler velocimetry-based volumetric flow rate measurements conducted in a pulsatile flow loop simulating conditions in both the internal carotid and basilar artery. In vitro models of cerebral vessels, each containing an aneurysm, were fabricated from patient anatomies extracted from 3D rotational angiography. Flow velocity measurements were collected with three different experimental techniques: an endovascular Doppler wire, Particle Image Velocimetry, and a time-resolved ultrasonic flow meter. Womersley's theory of pulsatile flow in a cylindrical vessel was used to compute time-resolved volumetric flow rates from the endovascular Doppler velocity. The volumetric flow rates computed from the Doppler measurements were compared to those from the Particle Image Velocimetry profile measurements, and the direct measurements from the ultrasonic flow meter. The study establishes confidence intervals for any systematic or random errors associated with the wire-derived flow rates as benchmarked to the other two modalities. There is an approximately 10% random error in the Doppler-derived peak and time-averaged flow rates. There is a measurable uniform bias, about 15% too low, in the time-averaged Doppler-derived flow rates. There is also a small proportional bias in the peak systolic Doppler-derived flow rates. Potential sources of error are also discussed.

Incomplete restoration of homeostatic shear stress within arteriovenous fistulae

Arteriovenous fistulae are surgically created to provide adequate access for dialysis patients suffering from end-stage renal disease. It has long been hypothesized that the rapid blood vessel remodeling occurring after fistula creation is, in part, a process to restore the mechanical stresses to some preferred level, i.e., mechanical homeostasis. We present computational hemodynamic simulations in four patient-specific models of mature arteriovenous fistulae reconstructed from 3D ultrasound scans. Our results suggest that these mature fistulae have remodeled to return to ''normal'' shear stresses away from the anastomoses: about 1.0 Pa in the outflow veins and about 2.5 Pa in the inflow arteries. Large parts of the anastomoses were found to be under very high shear stresses >15 Pa, over most of the cardiac cycle. These results suggest that the remodeling process works toward restoring mechanical homeostasis in the fistulae, but that the process is limited or incomplete, even in mature fistulae, as evidenced by the elevated shear at or near the anastomoses. Based on the long term clinical viability of these dialysis accesses, we hypothesize that the elevated nonhomeostatic shear stresses in some portions of the vessels were not detrimental to fistula patency.

Accuracy of computational cerebral aneurysm hemodynamics using patient-specific endovascular measurements

Computational hemodynamic simulations of cerebral aneurysms have traditionally relied on stereotypical boundary conditions (such as blood flow velocity and blood pressure) derived from published values as patient-specific measurements are unavailable or difficult to collect. However, controversy persists over the necessity of incorporating such patient-specific conditions into computational analyses. We perform simulations using both endovascularly-derived patient-specific and typical literature-derived inflow and outflow boundary conditions. Detailed three-dimensional anatomical models of the cerebral vasculature are developed from rotational angiography data, and blood flow velocity and pressure are measured in situ by a dual-sensor pressure and velocity endovascular guidewire at multiple peri-aneurysmal locations in 10 unruptured cerebral aneurysms. These measurements are used to define inflow and outflow boundary conditions for computational hemodynamic models of the aneurysms. The additional in situ measurements which are not prescribed in the simulation are then used to assess the accuracy of the simulated flow velocity and pressure drop. Simulated velocities using patient-specific boundary conditions show good agreement with the guidewire measurements at measurement locations inside the domain, with no bias in the agreement and a random scatter of ≈25%. Simulated velocities using the simplified, literature-derived values show a systematic bias and over-predicted velocity by ≈30% with a random scatter of ≈40%. Computational hemodynamics using endovascularly measured patient-specific boundary conditions have the potential to improve treatment predictions as they provide more accurate and precise results of the aneurysmal hemodynamics than those based on commonly accepted reference values for boundary conditions.

Effects of wall distensibility in hemodynamic simulations of an arteriovenous fistula

Arteriovenous fistulae are created surgically to provide adequate access for dialysis patients suffering from end-stage renal disease. It has long been hypothesized that the rapid blood vessel remodeling occurring after fistula creation is in part a process to restore the mechanical stresses to some preferred level, i.e., mechanical homeostasis. The current study presents fluid-structure interaction (FSI) simulations of a patient-specific model of a mature arteriovenous fistula reconstructed from 3D ultrasound scans. The FSI results are compared with previously published data of the same model but with rigid walls. Ultrasound-derived wall motion measurements are also used to validate the FSI simulations of the wall motion. Very large time-averaged shear stresses, 10-15 Pa, are calculated at the fistula anastomosis in the FSI simulations, values which are much larger than what is typically thought to be the normal homeostatic shear stress in the peripheral vasculature. Although this result is systematically lower by as much as 50% compared to the analogous rigid-walled simulations, the inclusion of distensible vessel walls in hemodynamic simulations does not reduce the high anastomotic shear stresses to "normal" values. Therefore, rigid-walled analyses may be acceptable for identifying high shear regions of arteriovenous fistulae.