Fluid Dynamic Study of the Penn State Pediatric Total Artificial Heart.

Penn State University is developing a pediatric total artificial heart (TAH) as a bridge-to-transplant device that supports infants and small children with single ventricle anomalies or biventricular heart failure to address high waitlist mortality rates for pediatric patients with severe congenital heart disease (CHD). Two issues with mechanical circulatory support devices are thrombus formation and thromboembolic events. This in vitro study characterizes flow within Penn State's pediatric total artificial heart under physiological operating conditions. Particle image velocimetry (PIV) is used to quantify flow within the pump and to calculate wall shear rates (WSRs) along the internal pump surface to identify potential thrombogenic regions. Results show that the diastolic inflow jets produce sufficient wall shear rates to reduce thrombus deposition potential along the inlet side of the left and right pumps. The inlet jet transitions to rotational flow, which promotes wall washing along the apex of the pumps, prevents flow stasis, and aligns flow with the outlet valve prior to systolic ejection. However, inconsistent high wall shear rates near the pump apex cause increased thrombogenic potential. Strong systolic outflow jets produce high wall shear rates near the outlet valve to reduce thrombus deposition risk. The right pump, which has a modified outlet port angle to improve anatomical fit, produces lower wall shear rates and higher thrombus susceptibility potential (TSP) compared to the left pump. In summary, this study provides a fluid dynamic understanding of a new pediatric total artificial heart and indicates thrombus susceptibility is primarily confined to the apex, consistent with similar pulsatile heart pumps.

The Accelerated Transcatheter Heart Valve Testing Environment: Loading, Motion, and Fluid Dynamics.

Transcatheter aortic valve replacements (TAVRs) are an increasingly common treatment for aortic valve disease due to their minimally invasive delivery. As TAVR designs require thinner leaflets to facilitate catheter-based delivery, they experience greater leaflet operational stresses and potentially greater durability issues than conventional surgical valves. Yet, our understanding of TAVR durability remains largely unexplored. Currently, preclinical TAVR durability is evaluated within an ISO:5840 compliant accelerated wear tester (AWT) up to a required 200 × 106 cycles, corresponding to approximately five years in vivo. While AWTs use high cycle frequencies (10-20 Hz) to achieve realistic timeframes, the resulting valve loading behaviors and fluid dynamics are not representative of the in vivo environment and thus may not accurately predict failure mechanisms. Despite the importance of fatigue and failure predictions for replacement heart valves, surprisingly, little quantitative information exists on the dynamic AWT environment. To better understand this environment, we examined frequency and diastolic period effects for the first time using high-speed enface imaging and particle image velocimetry to quantify valve motion and flow, respectively, using a Durapulse™ AWT at frequencies of 10, 15, and 20 Hz. Regardless of operating condition, no waveform achieved a physiologically relevant transvalvular loading pressure, despite having an ISO compliant geometric orifice area (GOA) and waveform. General fluid mechanics were consistent with in vivo but the AWT geometry developed secondary flow structures, which could impact mechanical loading. Therefore, the nonphysiologic loading and variability induced by changes in operating condition must be carefully regulated to ensure physiologically relevant fatigue.

Heparin-Peptide Nanogranules for Thrombosis-Actuated Anticoagulation.

Despite nearly a century of clinical use as a blood thinner, heparin's rapid serum clearance and potential to induce severe bleeding events continue to urge the development of more effective controlled delivery strategies. Subcutaneous depots that steadily release the anticoagulant into circulation represent a promising approach to reducing overdose frequency, sustaining therapeutic concentrations of heparin in plasma, and prolonging anticoagulant activity in a safe and effective manner. Subcutaneously deliverable heparin-peptide nanogranules that allow for long-lasting heparin bioavailability in the circulatory system, while enabling on-demand activation of heparin's anticoagulant effects in the thrombus microenvironment, are reported. Biophysical studies demonstrate this responsive behavior is due to the sequestration of heparin within self-assembling peptide nanofibrils and its mechanically actuated decoupling to elicit antithrombotic effects at the clotting site. In vivo studies show these unique properties converge to allow subcutaneous nanogranule depots to extend heparin serum concentrations for an order of magnitude longer than standard dosing regimens while enabling prolonged and controlled anticoagulant activity. This biohybrid delivery system demonstrates a potentially scalable platform for the development of safer, easier to administer, and more effective antithrombotic nanotechnologies.

Experimental Hemodynamics Within the Penn State Fontan Circulatory Assist Device.

For children born with a single functional ventricle, the Fontan operation bypasses the right ventricle by forming a four-way total cavopulmonary connection and adapts the existing ventricle for the systemic circulation. However, upon reaching adulthood, many Fontan patients exhibit low cardiac output and elevated venous pressure, eventually requiring a heart transplantation. Despite efforts in developing a new device or using an existing device for failing Fontan support, there is still no Food and Drug Administration-approved device for subpulmonary support. Penn State University is developing a hydrodynamically levitated Fontan circulatory assist device (FCAD) for bridge-to-transplant or destination therapy. The hemodynamics within the FCAD, at both steady and patient averaged pulsatile conditions for three physiological pump operating conditions, were quantified using particle image velocimetry (PIV) to determine the velocity magnitudes and Reynolds normal and shear stresses within the device. Data were acquired at three planes (0 mm and ±25% of the radius) for the inferior and superior vena cavae inlets and the pulmonary artery outlet. The inlets had a blunt velocity profile that became skewed toward the collecting volute as fluid approached the rotor. At the outlet, regardless of the flow condition, a high-velocity jet exited the volute and moved downstream in a helical pattern. Turbulent stresses observed at the volute exit were influenced by the rotor's rotation. Regardless of inlet conditions, the pump demonstrated advantageous behavior for clinical use with a predictable flow field and a low risk of platelet adhesion and hemolysis based on calculated wall shear rates and turbulent stresses, respectively.

In vitro real-time magnetic resonance imaging for quantification of thrombosis.

OBJECTIVES: Thrombosis is a leading cause of failure for cardiovascular devices. While computational simulations are a powerful tool to predict thrombosis and evaluate the risk for medical devices, limited experimental data are available to validate the simulations. The aim of the current study is to provide experimental data of a growing thrombus for device-induced thrombosis. MATERIALS AND METHODS: Thrombosis within a backward-facing step (BFS), or sudden expansion was investigated, using bovine and human blood circulated through the BFS model for 30 min, with a constant inflow rate of 0.76 L/min. Real-time three-dimensional flow-compensated magnetic resonance imaging (MRI), supported with Magnevist, a contrast agent improving thrombus delineation, was applied to quantify thrombus deposition and growth within the model. RESULTS: The study showed that the BFS model induced a flow recirculation region, which facilitated thrombosis. By 30 min, in comparison to bovine blood, human blood resulted in smaller thrombus formation, in terms of the length (13.3 ± 0.6 vs. 18.1 ± 1.3 mm), height (2.3 ± 0.1 vs. 2.6 ± 0.04 mm), surface area exposed to blood (0.67 ± 0.03 vs 1.05 ± 0.08 cm2), and volume (0.069 ± 0.004 vs. 0.093 ± 0.007 cm3), with p < 0.01. Normalization of the thrombus measurements, which excluded the flow recirculation effects, suggested that the thrombus sizes increased during the first 15 min and stabilized after 20 min. Blood properties, including viscosity, hematocrit, and platelet count affected thrombosis. CONCLUSION: For the first time, contrast agent-supported real-time MRI was performed to investigate thrombus deposition and growth within a sudden expansion. This study provides experimental data for device-induced thrombosis, which is valuable for validation of computational thrombosis simulations.

Computational modeling of the Food and Drug Administration's benchmark centrifugal blood pump.

In order to simulate hemodynamics within centrifugal blood pumps and to predict pump hemolysis, CFD simulations must be thoroughly validated against experimental data. They must also account for and accurately model the specific working fluid in the pump, whether that is a blood-analog solution to match an experimental PIV study or animal blood in a hemolysis experiment. Therefore, the Food and Drug Administration (FDA) benchmark centrifugal blood pump and its database of experimental PIV and hemolysis data were used to thoroughly validate CFD simulations of the same blood pump. A Newtonian blood model was first used to compare to the PIV data with a blood analog fluid while hemolysis data were compared using a power-law hemolysis model fit to porcine blood data. A viscoelastic blood model was then incorporated into the CFD solver to investigate the importance of modeling blood's viscoelasticity in centrifugal pumps. The established computational framework, including a dynamic rotating mesh, animal blood-specific fluid properties and hemolysis modeling, and a k-ω SST turbulence model, was shown to more accurately predict pump pressure heads, velocity fields, and hemolysis compared to previously published CFD studies of the FDA centrifugal pump. The CFD simulations were able to match the FDA pressure and hemolysis data for multiple pump operating conditions, with the CFD results being within the standard deviations of the experimental results. While CFD radial velocity profiles between the impeller blades also compared well to the PIV velocity results, more work is still needed to address the large variability among both experimental and computational predictions of velocity in the diffuser outlet jet. Small differences were observed between the Newtonian and viscoelastic blood models in pressure head and hemolysis at the higher flow rate cases (FDA Conditions 4 and 5) but were more significant at lower flow rate and pump impeller speeds (FDA Condition 1). These results suggest that the importance of accounting for blood's viscoelasticity may be dependent on the specific blood pump operating conditions. This detailed computational framework with improved modeling techniques and an extensive validation procedure will be used in future CFD studies of centrifugal blood pumps to aid in device design and predictions of their biological responses.

Characterizing the HeartMate II Left Ventricular Assist Device Outflow Using Particle Image Velocimetry.

Ventricular assist devices (VADs) are implanted in patients with a diseased ventricle to maintain peripheral perfusion as a bridge-to-transplant or as destination therapy. However, some patients with continuous flow VADs (e.g., HeartMate II (HMII)) have experienced gastrointestinal (GI) bleeding, in part caused by the proteolytic cleavage or mechanical destruction of von Willebrand factor (vWF), a clotting glycoprotein. in vitro studies were performed to measure the flow located within the HMII outlet cannula under both steady and physiological conditions using particle image velocimetry (PIV). Under steady flow, a mock flow loop was used with the HMII producing a flow rate of 3.2 L/min. The physiological experiment included a pulsatile pump operated at 105 BPM with a ventricle filling volume of 50 mL and in conjunction with the HMII producing a total flow rate of 5.0 L/min. Velocity fields, Reynolds normal stresses (RNSs), and Reynolds shear stresses (RSSs) were analyzed to quantify the outlet flow's potential contribution to vWF degradation. Under both flow conditions, the HMII generated principal Reynolds stresses that are, at times, orders of magnitude higher than those needed to unfurl vWF, potentially impacting its physiological function. Under steady flow, principal RNSs were calculated to be approximately 500 Pa in the outlet cannula. Elevated Reynolds stresses were observed throughout every phase of the cardiac cycle under physiological flow with principal RNSs approaching 1500 Pa during peak systole. Prolonged exposure to these conditions may lead to acquired von Willebrand syndrome (AvWS), which is accompanied by uncontrollable bleeding episodes.

Toward the Virtual Benchmarking of Pneumatic Ventricular Assist Devices: Application of a Novel Fluid-Structure Interaction-Based Strategy to the Penn State 12 cc Device.

The pediatric use of pneumatic ventricular assist devices (VADs) as a bridge to heart transplant still suffers for short-term major complications such as bleeding and thromboembolism. Although numerical techniques are increasingly exploited to support the process of device optimization, an effective virtual benchmark is still lacking. Focusing on the 12 cc Penn State pneumatic VAD, we developed a novel fluid-structure interaction (FSI) model able to capture the device functioning, reproducing the mechanical interplay between the diaphragm, the blood chamber, and the pneumatic actuation. The FSI model included the diaphragm mechanical response from uniaxial tensile tests, realistic VAD pressure operative conditions from a dedicated mock loop system, and the behavior of VAD valves. Our FSI-based benchmark effectively captured the complexity of the diaphragm dynamics. During diastole, the initial slow diaphragm retraction in the air chamber was followed by a more rapid phase; asymmetries were noticed in the diaphragm configuration during its systolic inflation in the blood chamber. The FSI model also captured the major features of the device fluid dynamics. In particular, during diastole, a rotational wall washing pattern is promoted by the penetrating inlet jet with a low-velocity region located in the center of the device. Our numerical analysis of the 12 cc Penn State VAD points out the potential of the proposed FSI approach well resembling previous experimental evidences; if further tested and validated, it could be exploited as a virtual benchmark to deepen VAD-related complications and to support the ongoing optimization of pediatric devices.

Platelet adhesion to polyurethane urea under pulsatile flow conditions.

Platelet adhesion to a polyurethane urea surface is a precursor to thrombus formation within blood-contacting cardiovascular devices, and platelets have been found to adhere strongly to polyurethane surfaces below a shear rate of approximately 500 s(-1). The aim of the current work is to determine the properties of platelet adhesion to the polyurethane urea surface as a function of time-varying shear exposure. A rotating disk system was used to study the influence of steady and pulsatile flow conditions (e.g., cardiac inflow and sawtooth waveforms) for platelet adhesion to the biomaterial surface. All experiments were conducted with the same root mean square angular rotation velocity (29.63 rad/s) and waveform period. The disk was rotated in platelet-rich bovine plasma for 2 h, with adhesion quantified by confocal microscopy measurements of immunofluorescently labeled bovine platelets. Platelet adhesion under pulsating flow was found to decay exponentially with increasing shear rate. Adhesion levels were found to depend upon peak platelet flux and shear rate, regardless of rotational waveform. In combination with flow measurements, these results may be useful for predicting regions susceptible to thrombus formation within ventricular assist devices.

A computational method for predicting inferior vena cava filter performance on a patient-specific basis.

A computational methodology for simulating virtual inferior vena cava (IVC) filter placement and IVC hemodynamics was developed and demonstrated in two patient-specific IVC geometries: a left-sided IVC and an IVC with a retroaortic left renal vein. An inverse analysis was performed to obtain the approximate in vivo stress state for each patient vein using nonlinear finite element analysis (FEA). Contact modeling was then used to simulate IVC filter placement. Contact area, contact normal force, and maximum vein displacements were higher in the retroaortic IVC than in the left-sided IVC (144 mm(2), 0.47 N, and 1.49 mm versus 68 mm(2), 0.22 N, and 1.01 mm, respectively). Hemodynamics were simulated using computational fluid dynamics (CFD), with four cases for each patient-specific vein: (1) IVC only, (2) IVC with a placed filter, (3) IVC with a placed filter and model embolus, all at resting flow conditions, and (4) IVC with a placed filter and model embolus at exercise flow conditions. Significant hemodynamic differences were observed between the two patient IVCs, with the development of a right-sided jet, larger flow recirculation regions, and lower maximum flow velocities in the left-sided IVC. These results support further investigation of IVC filter placement and hemodynamics on a patient-specific basis.

In vitro quantification of time dependent thrombus size using magnetic resonance imaging and computational simulations of thrombus surface shear stresses.

Thrombosis and thromboembolization remain large obstacles in the design of cardiovascular devices. In this study, the temporal behavior of thrombus size within a backward-facing step (BFS) model is investigated, as this geometry can mimic the flow separation which has been found to contribute to thrombosis in cardiac devices. Magnetic resonance imaging (MRI) is used to quantify thrombus size and collect topographic data of thrombi formed by circulating bovine blood through a BFS model for times ranging between 10 and 90 min at a constant upstream Reynolds number of 490. Thrombus height, length, exposed surface area, and volume are measured, and asymptotic behavior is observed for each as the blood circulation time is increased. Velocity patterns near, and wall shear stress (WSS) distributions on, the exposed thrombus surfaces are calculated using computational fluid dynamics (CFD). Both the mean and maximum WSS on the exposed thrombus surfaces are much more dependent on thrombus topography than thrombus size, and the best predictors for asymptotic thrombus length and volume are the reattachment length and volume of reversed flow, respectively, from the region of separated flow downstream of the BFS.

Current Approaches and Methods to Understand Acute Ischemic Stroke Treatment Using Aspiration Thrombectomy.

Acute ischemic stroke occurs when a blood clot occludes a cerebral artery. Mechanical interventions, primarily stent retrievers and aspiration thrombectomy, are used currently for removing the occluding clot and restoring blood flow. Aspiration involves using a long catheter to traverse the cerebral vasculature to reach the blood clot, followed by application of suction through the catheter bore. Aspiration is also used in conjunction with other techniques such as stent retrievers and balloon guide catheters. Despite the wide use of aspiration, our physical understanding of the process and the causes of the failure of aspiration to retrieve cerebral clots in certain scenarios is not well understood. Experimental and computational studies can help develop the capability to provide deeper insights into the procedure and enable development of new devices and more effective treatment methods. We recapitulate the aspiration-based thrombectomy techniques in clinical practice and provide a perspective of existing engineering methods for aspiration. We articulate the current knowledge gap in the understanding of aspiration and highlight possible directions for future engineering studies to bridge this gap, help clinical translation of engineering studies, and develop new patient-specific stroke therapy.

Influence of Hematocrit Level and Integrin αIIbßIII Function on vWF-Mediated Platelet Adhesion and Shear-Induced Platelet Aggregation in a Sudden Expansion.

Purpose: Shear-mediated thrombosis is a clinically relevant phenomenon that underlies excessive arterial thrombosis and device-induced thrombosis. Red blood cells are known to mechanically contribute to physiological hemostasis through margination of platelets and vWF, facilitating the unfurling of vWF multimers, and increasing the fraction of thrombus-contacting platelets. Shear also plays a role in this phenomenon, increasing both the degree of margination and the near-wall forces experienced by vWF and platelets leading to unfurling and activation. Despite this, the contribution of red blood cells in shear-induced platelet aggregation has not been fully investigated-specifically the effect of elevated hematocrit has not yet been demonstrated. Methods: Here, a microfluidic model of a sudden expansion is presented as a platform for investigating platelet adhesion at hematocrits ranging from 0 to 60% and shear rates ranging from 1000 to 10,000 s-1. The sudden expansion geometry models nonphysiological flow separation characteristic to mechanical circulatory support devices, and the validatory framework of the FDA benchmark nozzle. PDMS microchannels were fabricated and coated with human collagen. Platelets were fluorescently tagged, and blood was reconstituted at variable hematocrit prior to perfusion experiments. Integrin function of selected blood samples was inhibited by a blocking antibody, and platelet adhesion and aggregation over the course of perfusion was monitored. Results: Increasing shear rates at physiological and elevated hematocrit levels facilitate robust platelet adhesion and formation of large aggregates. Shear-induced platelet aggregation is demonstrated to be dependent on both αIIbßIII function and the presence of red blood cells. Inhibition of αIIbßIII results in an 86.4% reduction in overall platelet adhesion and an 85.7% reduction in thrombus size at 20-60% hematocrit. Hematocrit levels of 20% are inadequate for effective platelet margination and subsequent vWF tethering, resulting in notable decreases in platelet adhesion at 5000 and 10,000 s-1 compared to 40% and 60%. Inhibition of αIIbßIII triggered dramatic reductions in overall thrombus coverage and large aggregate formation. Stability of platelets tethered by vWF are demonstrated to be αIIbßIII-dependent, as adhesion of single platelets treated with A2A9, an anti-αIIbßIII blocking antibody, is transient and did not lead to sustained thrombus formation. Conclusions: This study highlights driving factors in vWF-mediated platelet adhesion that are relevant to clinical suppression of shear-induced thrombosis and in vitro assays of platelet adhesion. Primarily, increasing hematocrit promotes platelet margination, permitting shear-induced platelet aggregation through αIIbßIII-mediated adhesion at supraphysiological shear rates. Supplementary Information: The online version contains supplementary material available at 10.1007/s12195-024-00796-0.

Venous Thromboembolism: Review of Clinical Challenges, Biology, Assessment, Treatment, and Modeling.

Venous thromboembolism (VTE) is a massive clinical challenge, annually affecting millions of patients globally. VTE is a particularly consequential pathology, as incidence is correlated with extremely common risk factors, and a large cohort of patients experience recurrent VTE after initial intervention. Altered hemodynamics, hypercoagulability, and damaged vascular tissue cause deep-vein thrombosis and pulmonary embolism, the two permutations of VTE. Venous valves have been identified as likely locations for initial blood clot formation, but the exact pathway by which thrombosis occurs in this environment is not entirely clear. Several risk factors are known to increase the likelihood of VTE, particularly those that increase inflammation and coagulability, increase venous resistance, and damage the endothelial lining. While these risk factors are useful as predictive tools, VTE diagnosis prior to presentation of outward symptoms is difficult, chiefly due to challenges in successfully imaging deep-vein thrombi. Clinically, VTE can be managed by anticoagulants or mechanical intervention. Recently, direct oral anticoagulants and catheter-directed thrombolysis have emerged as leading tools in resolution of venous thrombosis. While a satisfactory VTE model has yet to be developed, recent strides have been made in advancing in silico models of venous hemodynamics, hemorheology, fluid-structure interaction, and clot growth. These models are often guided by imaging-informed boundary conditions or inspired by benchtop animal models. These gaps in knowledge are critical targets to address necessary improvements in prediction and diagnosis, clinical management, and VTE experimental and computational models.

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 analysis of effects of clot length on Acute ischemic stroke recanalization under different cyclic aspiration loading conditions.

Acute ischemic stroke, the second leading cause of death worldwide, results from occlusion of a cerebral artery by a blood clot. Application of cyclic aspiration using an aspiration catheter is a current therapy for the removal of lodged clots. In this study, we perform finite element simulations to analyze deformation of long clots, having length to radius ratio of 2-10, which corresponds to clot-length of 2.85-14.25 mm, under peak-to-peak cyclic aspiration pressures of 10-50 mmHg, and frequencies of 0.5, 1, and 2 Hz. Our computational system comprises of a nonlinear viscoelastic solid clot, a hyperelastic artery, and a nonlinear viscoelastic cohesive zone, the latter modeling the clot-artery interface. We observe that clots having length-to-radius ratio approximately greater than two separate from the inner arterial surface somewhere between the axial and distal ends, irrespective of the cyclic aspiration loading conditions. The stress distribution within the clot shows large tensile stresses in the clot interior, indicating the possibility of simultaneous fragmentation of the clot. Thus, this study shows us the various failure mechanisms simultaneously present in the clot during cyclic aspiration. Similarly, the stress distribution within the artery implies a possibility of endothelial damage to the arterial wall near the end where the aspiration pressure is applied. This framework provides a foundation for further investigation to clot fracture and adhesion characterization.

Examining the universality of the hemolysis power law model from simulations of the FDA nozzle using calibrated model coefficients.

Computational fluid dynamics (CFD) is widely used to predict mechanical hemolysis in medical devices. The most popular hemolysis model is the stress-based power law model that is based on an empirical correlation between hemoglobin release from red blood cells (RBCs) and the magnitude of flow-induced stress and exposure time. Empirical coefficients are traditionally calibrated using data from experiments in simplified Couette-type blood-shearing devices with uniform-shear laminar flow and well-defined exposure times. Use of such idealized coefficients in simulations of real medical devices with complex hemodynamics is thought to be a primary reason for the historical inaccuracy of absolute hemolysis predictions using the power law model. Craven et al. (Biomech Model Mechanobiol 18:1005-1030, 2019) recently developed a CFD-based Kriging surrogate modeling approach for calibrating empirical coefficients in real devices that could potentially be used to more accurately predict absolute hemolysis. In this study, we use the FDA benchmark nozzle to investigate whether utilizing such calibrated coefficients improves the predictive accuracy of the standard Eulerian power law model. We first demonstrate the credibility of our CFD flow simulations by comparing with particle image velocimetry measurements. We then perform hemolysis simulations and compare the results with in vitro experiments. Importantly, the simulations use coefficients calibrated for the flow of a suspension of bovine RBCs through a small capillary tube, which is relatively comparable to the flow of bovine blood through the FDA nozzle. The results show that the CFD predictions of relative hemolysis in the FDA nozzle are reasonably accurate. The absolute predictions are, however, highly inaccurate with modified index of hemolysis values from CFD in error by roughly three orders of magnitude compared with the experiments, despite using calibrated model coefficients from a relatively similar geometry. We rigorously examine the reasons for the inaccuracy that include differences in the flow conditions in the hemolytic regions of each device and the lack of universality of the hemolysis power law model that is entirely empirical. Thus, while the capability to predict relative hemolysis is valuable for product development, further improvements are needed before the power law model can be relied upon to accurately predict the absolute hemolytic potential of a medical device.

Clot embolization studies and computational framework for embolization in a canonical tube model.

Despite recent advances in the development of computational methods of modeling thrombosis, relatively little effort has been made in developing methods of modeling blood clot embolization. Such a model would provide substantially greater understanding of the mechanics of embolization, as in-vitro and in-vivo characterization of embolization is difficult. Here, a method of computationally simulating embolization is developed. Experiments are performed of blood clots formed in a polycarbonate tube, where phosphate-buffered saline is run through the tube at increasing flow rates until the clot embolizes. The experiments revealed embolization can be initiated by leading edge and trailing edge detachment or by non-uniform detachment. Stress-relaxation experiments are also performed to establish values of constitutive parameters for subsequent simulations. The embolization in the tube is reproduced in silico using a multiphase volume-of-fluid approach, where the clot is modeled as viscoelastic. By varying the constitutive parameters at the wall, embolization can be reproduced in-silico at varying flow rates, and a range of constitutive parameters fitting the experiments is reported. Here, the leading edge embolization is simulated at flow rates consistent with the experiments demonstrating excellent agreement in this specific behavior.

Modeling Flow in an In Vitro Anatomical Cerebrovascular Model with Experimental Validation.

Acute ischemic stroke (AIS) is a leading cause of mortality that occurs when an embolus becomes lodged in the cerebral vasculature and obstructs blood flow in the brain. The severity of AIS is determined by the location and how extensively emboli become lodged, which are dictated in large part by the cerebral flow and the dynamics of embolus migration which are difficult to measure in vivo in AIS patients. Computational fluid dynamics (CFD) can be used to predict the patient-specific hemodynamics and embolus migration and lodging in the cerebral vasculature to better understand the underlying mechanics of AIS. To be relied upon, however, the computational simulations must be verified and validated. In this study, a realistic in vitro experimental model and a corresponding computational model of the cerebral vasculature are established that can be used to investigate flow and embolus migration and lodging in the brain. First, the in vitro anatomical model is described, including how the flow distribution in the model is tuned to match physiological measurements from the literature. Measurements of pressure and flow rate for both normal and stroke conditions were acquired and corresponding CFD simulations were performed and compared with the experiments to validate the flow predictions. Overall, the CFD simulations were in relatively close agreement with the experiments, to within ±7% of the mean experimental data with many of the CFD predictions within the uncertainty of the experimental measurement. This work provides an in vitro benchmark data set for flow in a realistic cerebrovascular model and is a first step towards validating a computational model of AIS.

Performance of a Hydrogel Coated Nitinol with Oligonucleotide-Modified Nanoparticles Within Turbulent Conditions of Blood-Contacting Devices.

INTRODUCTION: Hydrogels offer a wide range of applications in the antithrombotic modification of biomedical devices. The functionalization of these hydrogels with potentially drug-laden nanoparticles in the context of deviceassociated turbulence is critically under-studied. Thus, the purpose of this study was to use a hydrogel-coating nitinol surface as a model to understand the functions of hydrogels and the capture of nanoparticles under clinically relevant flow conditions. METHODS: Nitinol was coated by an oligonucleotide (ON) functionalized hydrogel. Nanoparticles were functionalized with complementary oligonucleotides (CONs). The capture of CONfunctionalized nanoparticles by the ON-functionalized hydrogel surfaces was studied under both static and dynamic attachment conditions. Fluorescent-labelling of nanoparticles was utilized to assess capture efficacy and resistance to removal by device-relevant flow conditions. RESULTS: The specificity of the ON-CON bond was verified, exhibiting a dose-dependent attachment response. The hydrogel coating was resistant to stripping by flow, retaining >95% after exposure to one hour of turbulent flow. Attachment of nanoparticles to the hydrogel was higher in the static condition than under laminar flow (p < 0.01), but comparable to that of attachment under turbulent flow. Modified nitinol samples underwent one hour of flow treatment under both laminar and turbulent regimes and demonstrated decreased nanoparticle loss following static conjugation rather than turbulent conjugation (36.1% vs 53.8%, p < 0.05). There was no significant difference in nanoparticle functionalization by upstream injection between laminar and turbulent flow. CONCLUSION: The results demonstrate promising potential of hydrogelfunctionalized nitinol for capturing nanoparticles using nucleic acid hybridization. The hydrogel structure and ONCON bond integrity both demonstrated a resistance to mechanical damage and loss of biomolecular functionalization by exposure to turbulence. Further investigation is warranted to highlight drug delivery and antithrombogenic modification applications of nanoparticle-functionalized hydrogels.