A novel mathematical model of activation and sensitization of platelets subjected to dynamic stress histories

Fluid-structure interaction simulation of mechanical aortic valves: a narrative review exploring its role in total product life cycle

Over the last years computer modelling and simulation has emerged as an effective tool to support the total product life cycle of cardiovascular devices, particularly in the device preclinical evaluation and post-market assessment. Computational modelling is particularly relevant for heart valve prostheses, which require an extensive assessment of their hydrodynamic performance and of risks of hemolysis and thromboembolic complications associated with mechanically-induced blood damage. These biomechanical aspects are typically evaluated through a fluid-structure interaction (FSI) approach, which enables valve fluid dynamics evaluation accounting for leaflets movement. In this context, the present narrative review focuses on the computational modelling of bileaflet mechanical aortic valves through FSI approach, aiming to foster and guide the use of simulations in device total product life cycle. The state of the art of FSI simulation of heart valve prostheses is reviewed to highlight the variety of modelling strategies adopted in the literature. Furthermore, the integration of FSI simulations in the total product life cycle of bileaflet aortic valves is discussed, with particular emphasis on the role of simulations in complementing and potentially replacing the experimental tests suggested by international standards. Simulations credibility assessment is also discussed in the light of recently published guidelines, thus paving the way for a broader inclusion of in silico evidence in regulatory submissions. The present narrative review highlights that FSI simulations can be successfully framed within the total product life cycle of bileaflet mechanical aortic valves, emphasizing that credible in silico models evaluating the performance of implantable devices can (at least) partially replace preclinical in vitro experimentation and support post-market biomechanical evaluation, leading to a reduction in both time and cost required for device development.

Optimization design of semi-open impeller based on thrombogenicity in a blood pump

BACKGROUND

Blood clots are composed of aggregated fibrin and platelets, and thrombosis is the body's natural response to repairing injured blood vessels or stopping bleeding. However, when this process is activated abnormally, such as in a mechanical blood pump, it can lead to excessive thrombus formation. Therefore, how to avoid or reduce the probability of thrombus formation is an important indicator of the stable operation of a blood pump.

METHODS

In this paper, Lagrangian particle tracking trajectories are simulated to study platelet transport in a blood pump. The design of the thrombus blood pump was optimized using an orthogonal design method based on three factors: inlet angle, outlet angle, and blade number. The effect of blood pump pressure, rotational speed, impeller outlet angle, inlet angle, and number of blades on thrombus formation was analysed using Fluent software. The thrombogenic potential was derived by analyzing the trajectory and flow parameters of platelet particles in the blood pump, as well as the statistical parameters of residence time and stress accumulation thrombus in the platelet pump.

RESULTS

When the impeller inlet angle is 30°, the outlet angle is 20°, and the number of blades is 6, the probability of thrombus formation is minimized in the orthogonal design method, aligning with the requirements for blood pump performance.

CONCLUSIONS

These design parameters serve as a numerical guideline for optimizing the geometry of the semi-open impeller in blood pumps and provide a theoretical foundation for subsequent in vitro experiments.

Numerical methods for hemolysis and thrombus evaluation in the percutaneous ventricular assist device

BACKGROUND

A percutaneous ventricular assist device (pVAD) is an effective method to treat heart failure, but its complications, mainly hemolysis and thrombus formation, cannot be ignored. Accurate evaluation of hemolysis and thrombus formation in pVAD is essential to guide the development of pVAD and reduce the incidence of complications.

METHODS

This study optimized the numerical model to predict hemolysis and thrombus formation in pVAD. The hemolysis model is based on the power law function, and the multi-component thrombus prediction model is improved by introducing the von Willebrand factor.

RESULTS

The error between the numerical simulation and the hydraulic performance experiment is within 5%. The numerical results of hemolysis are in good agreement with those of in vitro experiments. Meanwhile, the thrombus location predicted by the numerical model is the same as that found in the in vivo experiment.

CONCLUSION

The numerical model suggested in this study may therefore accurately assess the possible hemolytic and thrombotic dangers in pVAD, making it an effective tool to support the development of pVAD.

Shear Histories Alter Local Shear Effects on Thrombus Nucleation and Growth

Our goal was to determine the impact of physiological and pathological shear histories on platelet nucleation and thrombus growth at various local shear rates. We designed and characterized a microfluidic device capable of subjecting platelets to shear histories reaching as high as 6700 s- 1 in a single passage. Time-lapse videos of platelets and thrombi are captured using fluorescence microscopy. Thrombi are tracked, and the degree of thrombosis is evaluated through surface coverage, platelet nucleation maps, and ensemble-averaged aggregate areas and intensities. Surface coverage rates were the lowest when platelets deposited at high shear rates following a pathological shear history and were highest at low shear rates following a pathological shear history. Early aggregate area growth rates were significantly larger for thrombi developing at high shear following physiological shear history than at high shear following a pathological shear history. Aggregate vertical growth was restricted when depositing at low shear following a pathological shear history. In contrast, thrombi grew faster vertically following physiological shear histories. These results show that physiological shear histories pose thrombotic risks via volumetric growth, and pathological shear histories drastically promote nucleation. These findings may inform region-based geometries for biomedical devices and refine thrombosis simulations.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSION

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

A multi-constituent model for assessing the effect of impeller shroud on the thrombosis potential of a centrifugal blood pump

Centrifugal blood pumps can be used for treating heart failure patients. However, pump thrombosis has remained one of the complications that trouble clinical treatment. This study analyzed the effect of impeller shroud on the thrombosis risk of the blood pump, and predicted areas prone to thrombosis. Multi-constituent transport equations were presented, considering mechanical activation and biochemical activation. It was found that activated platelets concentration can increase with shear stress and adenosine diphosphate(ADP) concentration increasing, and the highest risk of thrombosis inside the blood pump was under extracorporeal membrane oxygenation (ECMO) mode. Under the same condition, ADP concentration and thrombosis index of semi-shroud impeller can increase by 7.3% and 7.2% compared to the closed-shroud impeller. The main reason for the increase in thrombosis risk was owing to elevated scalar shear stress and more coagulation promoting factor-ADP released. The regions with higher thrombosis potential were in the center hole, top and bottom clearance. As a novelty, the findings revealed that impeller shroud can influence mechanical and biochemical activation factors. It is useful for identifying potential risk regions of thrombus formation based on relative comparisons.

Vibro-Acoustic Platelet Activation: An Additive Mechanism of Prothrombosis with Applicability to Snoring and Obstructive Sleep Apnea

Introduction: Obstructive sleep apnea (OSA) and loud snoring are conditions with increased cardiovascular risk and notably an association with stroke. Central in stroke are thrombosis and thromboembolism, all related to and initiaing with platelet activation. Platelet activation in OSA has been felt to be driven by biochemical and inflammatory means, including intermittent catecholamine exposure and transient hypoxia. We hypothesized that snore-associated acoustic vibration (SAAV) is an activator of platelets that synergizes with catecholamines and hypoxia to further amplify platelet activation. Methods: Gel-filtered human platelets were exposed to snoring utilizing a designed vibro-acoustic exposure device, varying the time and intensity of exposure and frequency content. Platelet activation was assessed via thrombin generation using the Platelet Activity State assay and scanning electron microscopy. Comparative activation induced by epinephrine and hypoxia were assessed individually as well as additively with SAAV, as well as the inhibitory effect of aspirin. Results: We demonstrate that snore-associated acoustic vibration is an independent activator of platelets, which is dependent upon the dose of exposure, i.e., intensity x time. In snoring, acoustic vibrations associated with low-frequency sound content (200 Hz) are more activating than those associated with high frequencies (900 Hz) (53.05% vs. 22.08%, p = 0.001). Furthermore, SAAV is additive to both catecholamines and hypoxia-mediated activation, inducing synergistic activation. Finally, aspirin, a known inhibitor of platelet activation, has no significant effect in limiting SAAV platelet activation. Conclusion: Snore-associated acoustic vibration is a mechanical means of platelet activation, which may drive prothrombosis and thrombotic risk clinically observed in loud snoring and OSA.

Toward modeling thrombosis and thromboembolism in laminar and turbulent flow regimes

Thrombosis and thromboembolism are deadly risk factors in blood-contacting biomedical devices, and in-silico models of thrombosis are attractive tools to understand the mechanics of these processes, though the simulation of thromboembolism remains underdeveloped. The purpose of this study is to modify an existing computational thrombosis model to allow for thromboembolism and to investigate the behavior of the modified model at a range of flow rates. The new and existing models are observed to lead to similar predictions of thrombosis in a canonical backward-facing step geometry across flow rates, and neither model predicts thrombosis in a turbulent flow. Simulations are performed by increasing flow rates in the case of a clot formed at lower flow to induce embolization. While embolization is observed, most of the clot breakdown is by shear rather than by breakup and subsequent transport of clotted material, and further work is required in the formulation and validation of embolization. This model provides a framework to further investigate thromboembolization.

The Effects of Implantation Orientation of a Bileaflet Mechanical Heart Valve in an Anatomic Left Ventricle-Aorta Configuration

We have performed three-dimensional high-resolution numerical simulations of a bi-leaflet mechanical heart valve (BMHV) implanted at different orientations in an anatomic left ventricle-aorta obtained from magnetic resonance imaging (MRI) of a volunteer. The thoroughly validated overset curvilinear-immersed boundary (overset-CURVIB) fluid-structure interaction (FSI) flow solver is used in which the aorta and LV are discretized with boundary-conforming and non-conforming curvilinear grids, respectively. The motion of the LV wall is prescribed based on a lumped parameter model while the motion of the leaflets are calculated using a strong coupled FSI algorithm enhanced with Aitken convergence technique. We carried out simulations for three valve orientations, which differ from each other by 45 degrees and compared the leaflet motion and flow field for multiple cycles. Our results show reproducible and relatively symmetrical opening for all valve orientations. The presence of small-scale vortical structures after peak systole, cause significant cycle-to-cycle variations in valve kinematics during the closing phase for all valve orientations. Furthermore, our results show that valve orientation does not have a significant effect on the distribution of viscous shear stress in the ascending aorta. Additionally, two different mathematical activation models including linear level of activation and Soares model are used to quantify the platelet activation in the ascending aorta. The results show that the valve orientation does not significantly affect (less than 8%) the total platelet activation in the ascending aorta.

Models of Shear-Induced Platelet Activation and Numerical Implementation With Computational Fluid Dynamics Approaches

Shear-induced platelet activation is one of the critical outcomes when blood is exposed to elevated shear stress. Excessively activated platelets in the circulation can lead to thrombus formation and platelet consumption, resulting in serious adverse events such as thromboembolism and bleeding. While experimental observations reveal that it is related to the shear stress level and exposure time, the underlying mechanism of shear-induced platelet activation is not fully understood. Various models have been proposed to relate shear stress levels to platelet activation, yet most are modified from the empirically calibrated power-law model. Newly developed multiscale platelet models are tested as a promising approach to capture a single platelet's dynamic shape during activation, but it would be computationally expensive to employ it for a large-scale analysis. This paper summarizes the current numerical models used to study the shear-induced platelet activation and their computational applications in the risk assessment of a particular flow pattern and clot formation prediction.

Ramifications of Vorticity on Aggregation and Activation of Platelets in Bi-Leaflet Mechanical Heart Valve: Fluid-Structure-Interaction Study

Bileaflet Mechanical Heart Valves (BMHV) are widely implanted to replace diseased heart valves. Despite many improvements in design, these valves still suffer from various complications, such as valve dysfunction, tissue overgrowth, hemolysis, and thromboembolism. Thrombosis and thromboembolism are believed to be initiated by platelet activation due to contact with foreign surfaces and non-physiological flow patterns. The implantation of the valve causes non-physiological patterns of vortex shedding behind the leaflets. The present study signifies the importance of vorticity in platelet activation and aggregation in BMHV implants. A two-phase model with the first Eulerian phase for blood and the second Discrete phase for platelets are used here. The generalized cross model of viscosity has been used to simulate the non-Newtonian viscosity of blood. A Fluid-Structure-Interaction model has been used to simulate the motion of leaflets. The present study has also estimated Platelet Activation State (PAS), which is the mathematical estimation of the degree of activation of platelets due to flow-induced shear stresses that cause thrombus formation. The regions in the fluid domain with a higher vorticity field have been found to contain platelets with relatively higher PAS than regions with relatively lower vorticity fields. Also, this study has quantitatively reported the effect of vorticity on platelet aggregation. The densities of platelets in the fluid areas with higher vorticity fields are higher than densities in the fluid regions with relatively lower vorticity fields, which indicate aggregation of highly activated platelets in areas with somewhat higher vorticity.

An Accelerated Thrombosis Model for Computational Fluid Dynamics Simulations in Rotary Blood Pumps

Purpose

Thrombosis ranks among the major complications in blood-carrying medical devices and a better understanding to influence the design related contribution to thrombosis is desirable. Over the past years many computational models of thrombosis have been developed. However, numerically cheap models able to predict localized thrombus risk in complex geometries are still lacking. The aim of the study was to develop and test a computationally efficient model for thrombus risk prediction in rotary blood pumps.

Methods

We used a two-stage approach to calculate thrombus risk. The first stage involves the computation of velocity and pressure fields by computational fluid dynamic simulations. At the second stage, platelet activation by mechanical and chemical stimuli was determined through species transport with an Eulerian approach. The model was compared with existing clinical data on thrombus deposition within the HeartMate II. Furthermore, an operating point and model parameter sensitivity analysis was performed.

Results

Our model shows good correlation (R² > 0.93) with clinical data and identifies the bearing and outlet stator region of the HeartMate II as the location most prone to thrombus formation. The calculation of thrombus risk requires an additional 10–20 core hours of computation time.

Conclusion

The concentration of activated platelets can be used as a surrogate and computationally low-cost marker to determine potential risk regions of thrombus deposition in a blood pump. Relative comparisons of thrombus risk are possible even considering the intrinsic uncertainty in model parameters and operating conditions.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13239-021-00606-y.

Blood Pumps for Extracorporeal Membrane Oxygenation: Platelet Activation During Different Operating Conditions

Extracorporeal membrane oxygenation (ECMO) is a therapy used in severe cardiopulmonary failure. Blood is pumped through an artificial circuit exposing it to nonphysiologic conditions, which promote platelet activation and coagulation. Centrifugal pumps used at lower flow rates than their design point may lose pump efficiency and increase the risk of hemolysis. In this study, thrombogenic properties of two ECMO pumps designed for adult and neonatal use were evaluated using simulations in different flow scenarios. Three scenarios, adult pump in adult mode (4 L/min), adult pump in baby mode (300 ml/min), and neonatal pump used in its design point (300 ml/min), were simulated using computational fluid dynamics. The flow was numerically seeded with platelets, whose activation state was computed considering the stress history that acted along their respective path lines. Statistical distributions of activation state and residence time were drawn. The results showed that using the adult pump in baby mode increased the fraction of platelets with higher activation state confirming that low-pump flow rate impacts thrombogenicity. The neonatal pump showed a backflow at the inlet, which carried platelets in a retrograde motion contributing to an increased thrombogenic potential compared with the adult mode scenario.

Discrete responses of erythrocytes, platelets, and von Willebrand factor to shear

Despite decades of technological advancements in blood-contacting medical devices, complications related to shear flow-induced blood trauma are still frequently observed in clinic. Blood trauma includes haemolysis, platelet activation, and degradation of High Molecular Weight von Willebrand Factor (HMW vWF) multimers, all of which are dependent on the exposure time and magnitude of shear stress. Specifically, accumulating evidence supports that when blood is exposed to shear stresses above a certain threshold, blood trauma ensues; however, it remains unclear how various constituents of blood are affected by discrete shears experimentally. The aim of this study was to expose blood to discrete shear stresses and evaluate blood trauma indices that reflect red cell, platelet, and vWF structure. Citrated human whole blood (n = 6) was collected and its haematocrit was adjusted to 30 ± 2% by adding either phosphate buffered saline (PBS) or polyvinylpyrrolidone (PVP). Viscosity of whole blood was adjusted to 3.0, 12.5, 22.5 and 37.5 mPa·s to yield stresses of 3, 6, 9, 12, 50, 90 and 150 Pa in a custom-developed shearing system. Blood samples were exposed to shear for 0, 300, 600 and 900 s. Haemolysis was measured using spectrophotometry, platelet activation using flow cytometry, and HMW vWF multimer degradation was quantified with gel electrophoresis and immunoblotting. For tolerance to 300, 600 and 900 s of exposure time, the critical threshold of haemolysis was reached after blood was exposed to 90 Pa for 600 s (P < 0.05), platelet activation and HMW vWF multimer degradation were 50 Pa for 600 s and 12 Pa for 300 s respectively (P < 0.05). Our experimental results provide simultaneous comparison of blood trauma indices and thus also the relation between shear duration and magnitude required to induce damage to red cells, platelets, and vWF. Our results also demonstrate that near-physiological shear stress (<12 Pa) is needed in order to completely avoid any form of blood trauma. Therefore, there is an urgent need to design low shear-flow medical devices in order to avoid blood trauma in this blood-contacting medical device field.

The effect of turbulence modelling on the assessment of platelet activation

Pathological platelet activation by abnormal shear stresses is regarded as a main clinical complication in recipients of cardiovascular mechanical devices. In order to improve their performance computational fluid dynamics (CFD) are used to evaluate flow fields and related shear stresses. CFD models are coupled with mathematical models that describe the relation between fluid dynamics variables, and in particular shear stresses, and the platelet activation state (PAS). These models typically use a Lagrangian approach to compute the shear stresses along possible platelet trajectories. However, in the case of turbulent flow, the choice of the proper turbulence closure is still debated for both concerning its effect on shear stress calculation and Lagrangian statistics. In this study different numerical simulations of the flow through a mechanical heart valve were performed and then compared in terms of Eulerian and Lagrangian quantities: a direct numerical simulation (DNS), a large eddy simulation (LES), two Reynolds-averaged Navier-Stokes (RANS) simulations (SST k-ω and RSM) and a "laminar" (no turbulence modelling) simulation. Results exhibit a large variability in the PAS assessment depending on the turbulence model adopted. "Laminar" and RSM estimates of platelet activation are about 60% below DNS, while LES is 16% less. Surprisingly, PAS estimated from the SST k- ω velocity field is only 8% less than from DNS data. This appears more artificial than physical as can be inferred after comparing frequency distributions of PAS and of the different Lagrangian variables of the mechano-biological model of platelet activation. Our study indicates how much turbulence closures may affect platelet activation estimates, in comparison to an accurate DNS, when assessing blood damage in blood contacting devices.

Influence of shear rate and surface chemistry on thrombus formation in micro-crevice

Thromboembolic complications remain a central issue in management of patients on mechanical circulatory support. Despite the best practices employed in design and manufacturing of modern ventricular assist devices, complexity and modular nature of these systems often introduces internal steps and crevices in the flow path which can serve as nidus for thrombus formation. Thrombotic potential is influenced by multiple factors including the characteristics of the flow and surface chemistry of the biomaterial. This study explored these elements in the setting of blood flow over a micro-crevice using a multi-constituent numerical model of thrombosis. The simulations reproduced the platelet deposition patterns observed experimentally and elucidated the role of flow, shear rate, and surface chemistry in shaping the deposition. The results offer insights for design and operation of blood-contacting devices.

Advances in extracorporeal membrane oxygenator design for artificial placenta technology

Extreme prematurity, defined as a gestational age of fewer than 28 weeks, is a significant health problem worldwide. It carries a high burden of mortality and morbidity, in large part due to the immaturity of the lungs at this stage of development. The standard of care for these patients includes support with mechanical ventilation, which exacerbates lung pathology. Extracorporeal life support (ECLS), also called artificial placenta technology when applied to extremely preterm (EPT) infants, offers an intriguing solution. ECLS involves providing gas exchange via an extracorporeal device, thereby doing the work of the lungs and allowing them to develop without being subjected to injurious mechanical ventilation. While ECLS has been successfully used in respiratory failure in full-term neonates, children, and adults, it has not been applied effectively to the EPT patient population. In this review, we discuss the unique aspects of EPT infants and the challenges of applying ECLS to these patients. In addition, we review recent progress in artificial placenta technology development. We then offer analysis on design considerations for successful engineering of a membrane oxygenator for an artificial placenta circuit. Finally, we examine next-generation oxygenators that might advance the development of artificial placenta devices.

Fluid dynamics characterization and thrombogenicity assessment of a levitating centrifugal pump with different impeller designs

Technical guidelines nowadays recommend and regulate the use Computational Fluid Dynamics (CFD) to assess the performance of medical devices. CFD coupled to blood damage models has emerged as a powerful tool to evaluate the hemocompatibility of blood recirculating devices. The present study is aimed at evaluating the hydrodynamic performance and the thrombogenic potential of two prototypes of magnetically levitating centrifugal pumps. The two devices differ in the impeller configuration - 6-blades vs. 12-blades - and have been designed to be used in Cardiopulmonary Bypass (CPB) circuits during open heart surgery and in Extracorporeal Membrane Oxygenation (ECMO) to support patients with severe cardiac or respiratory failure. The pumps have been modelled using Direct Numerical Simulation coupled to Lagrangian analysis to predict platelet activation due to abnormal shear stress histories. Numerical results have been compared with experimental data in terms of head generation for different working points. Results show that the 6-blades pump has i) smaller stagnation areas, ii) lower stress levels and iii) higher strain rate, resulting in a lower thrombogenic potential, whereas the 12-blade impeller guarantees a more stable performance at high flow rates, suggesting its preferential use for more demanding applications, such as CPB.

Evaluating medical device and material thrombosis under flow: current and emerging technologies

This review highlights the importance of flow in medical device thrombosis and explores current and emerging technologies to evaluate dynamic biomaterial Thrombosis in vitro.

Numerical investigation on the effect of bileaflet mechanical heart valve's implantation tilting angle and aortic root geometry on intermittent regurgitation and platelet activation

Platelet activation induced by shear stresses and non-physiological flow field generated by bileaflet mechanical heart valves (BMHVs) leads to thromboembolism, which can cause fatal consequences. One of the causes of platelet activation could be intermittent regurgitation, which arises due to asynchronous movement and rebound of BMHV leaflets during the valve closing phase. In this numerical study, the effect of intermittent regurgitation on the platelet activation potential of BMHVs was quantified by modeling a BMHV in the straight and anatomic aorta at implantation tilt angles 0°, 5°, 10°, and 20°. A fully implicit Arbitrary Lagrangian-Eulerian-based Fluid-Structure Interaction formulation was adopted with blood modeled as a multiphase, non-Newtonian fluid. Results showed that the intermittent regurgitation and consequently the platelet activation level increases with the increasing implantation tilt of BMHV. For the straight aorta, the leaflet of the 20° tilted BMHV underwent a rebound of approximately 20° after initially closing, whereas the leaflet of the 10°, 5°, and 0° tilted BMHVs underwent a rebound of 8.5°, 3°, and 0°, respectively. For the anatomic aorta, the leaflet of the 20° tilted BMHV underwent a rebound of approximately 24° after initially closing, whereas the leaflet of the 10°, 5°, and 0° tilted BMHVs underwent a rebound of 14°, 10°, and 7°, respectively. For all the implantation orientations of BMHVs, intermittent regurgitation and platelet activation were always higher in the anatomic aorta than in the straight aorta. The study concludes that the pivot axis of BMHV must be implanted parallel to the aortic root's curvature to minimize intermittent regurgitation and platelet activation.

Modeling sensitivity and uncertainties in platelet activation models applied on centrifugal pumps for extracorporeal life support

Two platelet activation models were studied with respect to uncertainties of model parameters and variables. The sensitivity was assessed using two direct/deterministic approaches as well as the statistical Monte Carlo method. The first two, are linear in character whereas the latter is non-linear. The platelet activation models were applied on platelets moving within an extracorporeal centrifugal blood pump. The phenomenological, Lagrangian stress- and time-based power law-based models under consideration, have experimentally calibrated parameters and the stress expressed in a scalar form. The sensitivity of the model with respect to model parameters and the expression of the scalar stress was examined focusing on a smaller group of platelets associated with an elevated risk of activation. The results showed a high disparity between the models in terms of platelet activation state, found to depend on the platelets’ trajectory in the pump and the expression used for the scalar stress. Monte Carlo statistics was applied to the platelets at risk for activation and not to the entire platelet population. The method reveals the non-linear sensitivity of the activation models. The results imply that power-law based models have a restricted range of validity. The conclusions of this study apply to both platelet activation and hemolysis models.

Flow-induced platelet activation in components of the extracorporeal membrane oxygenation circuit

Extracorporeal membrane oxygenation (ECMO) is used for rescue in severe respiratory and/or circulatory failure. The patient’s blood is pumped over artificial surfaces in the ECMO circuit. A platelet activation model was applied to study the potential thrombogenicity of ECMO circuit components: the centrifugal blood pump, cannulae, and tubing connectors. Based on the accumulated effect of the scalar form of the stress acting on the platelet over time, the activation model enables assessment of platelet activation and pinpoints regions of elevated activation risk in a component. Numerical simulations of the flow in different components of the ECMO circuit was carried out where the activation level is a function of the impact of local stress and its history along the path that the platelets follow. The results showed that the pump carried the largest risk for platelet activation followed by the reinfusion cannula and lastly the connectors. Pump thrombogenicity was mainly due to long residence time and high shear-rate while the connector showed a high level of non-stationary shear-rate that in turn may contribute to the formation of aggregates through direct platelet activation or through high shear-rate modulation of the vWF multimers.

Methods for determination of stagnation in pneumatic ventricular assist devices

Background:

A pneumatic paediatric ventricular assist device developed at the Foundation of Cardiac Surgery Development, Zabrze, equipped with valves based on J. Moll’s design, with later modifications introduced at the Institute of Turbomachinery, Lodz University of Technology, was tested numerically and experimentally. The main aim of those investigations was to detect stagnation zones within the ventricular assist device and indicate advantages and limitations of both approaches.

Methods:

In the numerical transient test, a motion of the diaphragm and discs was simulated. Two different methods were used to illustrate stagnation zones in the ventricular assist device. The flow pattern inside the chamber was represented by velocity contours and vectors to validate the results using images obtained in the laser particle image velocimetry experiment.

Results:

The experimental light-based method implied problems with proper illumination of regions in the wall vicinity. High-resolution flow data and other important parameters as stagnation regions or flow patterns in regions not accessible for light in the particle image velocimetry method are derived in the numerical solution. However, computations of a single case are much more time-consuming if compared to a single experiment conducted on a well-calibrated stand.

Conclusion:

The resulting main vortexes in the central part of the pump chamber and the velocity magnitudes are correlated in both methods, which are complementary and when used together offer better insight into the flow structure inside the ventricular assist device and enable a deeper analysis of the results.

Microfluidic flow-based platforms for induction and analysis of dynamic shear-mediated platelet activation-Initial validation versus the standardized hemodynamic shearing device

A microfluidic flow-based platform (μFP), able to stimulate platelets via exposure of shear stress patterns pertinent to cardiovascular devices and prostheses, was compared to the Hemodynamic Shearing Device (HSD)-a state-of-the-art bench-top system for exposure of platelets to defined levels and patterns of shear. Platelets were exposed to time-varying shear stress patterns in the two systems; in detail, platelets were recirculated in the μFP or stimulated in the HSD to replicate comparable exposure time. Shear-mediated platelet activation was evaluated via (i) the platelet activity state assay, allowing the measurement of platelet-mediated thrombin generation and associated prothrombotic tendencies, (ii) scanning electron microscopy to evaluate morphological changes of sheared platelets, and (iii) flow cytometry for the determination of platelet phosphatidylserine exposure as a marker of shear activation. The results revealed good matching and comparability between the two systems, with similar trends of platelet activation, formation of microaggregates, and analogous trends of activation marker exposure for both the HSD and microfluidic-stimulated samples. These findings support future translation of the microfluidic platform as a Point-of-Care facsimile system for the diagnosis of thrombotic risk in patients implanted with cardiovascular devices.

High Frequency Components of Hemodynamic Shear Stress Profiles are a Major Determinant of Shear-Mediated Platelet Activation in Therapeutic Blood Recirculating Devices

We systematically analyzed the relative contributions of frequency component elements of hemodynamic shear stress waveforms encountered in cardiovascular blood recirculating devices as to overall platelet activation over time. We demonstrated that high frequency oscillations are the major determinants for priming, triggering and yielding activated “prothrombotic behavior” for stimulated platelets, even if the imparted shear stress has low magnitude and brief exposure time. Conversely, the low frequency components of the stress signal, with limited oscillations over time, did not induce significant activation, despite being of high magnitude and/or exposure time. In vitro data were compared with numerical predictions computed according to a recently proposed numerical model of shear-mediated platelet activation. The numerical model effectively resolved the correlation between platelet activation and the various frequency components examined. However, numerical predictions exhibited a different activation trend compared to experimental results for different time points of a stress activation sequence. With this study we provide a more fundamental understanding for the mechanobiological responsiveness of circulating platelets to the hemodynamic environment of cardiovascular devices, and the importance of these environments in mediating life-threatening thromboembolic complications associated with shear-mediated platelet activation. Experimental data will guide further optimization of the thromboresistance of cardiovascular implantable therapeutic devices.

Prediction of Thrombus Growth: Effect of Stenosis and Reynolds Number

Shear stresses play a major role in platelet-substrate interactions and thrombus formation and growth in blood flow, where under both pathological and physiological conditions platelet adhesion and accumulation occur. In this study, a shear-dependent continuum model for platelet activation, adhesion and aggregation is presented. The model was first verified under three different shear conditions and at two heparin levels. Three-dimensional simulations were then carried out to evaluate the performance of the model for severely damaged (stripped) aortas with mild and severe stenosis degrees in laminar flow regime. For these cases, linear shear-dependent functions were developed for platelet-surface and platelet-platelet adhesion rates. It was confirmed that the platelet adhesion rate is not only a function of Reynolds number (or wall shear rate) but also the stenosis severity of the vessel. General correlations for adhesion rates of platelets as functions of stenosis and Reynolds number were obtained based on these cases. Finally using the new platelet adhesion rates, the model was applied to different experimental systems and shown to agree well with measured platelet deposition.

Development of a platelet adhesion transport equation for a computational thrombosis model

Thrombosis is a significant issue for cardiovascular device development and use. While thrombosis models are available, very few are device-related and none have been thoroughly validated experimentally. Here, we introduce a surface adherent platelet transport equation into a continuum model to account for the biomaterial interface/blood interaction. Using a rotating disc system and polyurethane-urea material, we characterize steady and pulsatile flow fields using laser Doppler velocimetry. In vitro measurements of platelet adhesion are used in combination with the LDV data to provide further experimental validation. The rotating disc system is computationally studied using the device-induced thrombosis model with the surface platelet adherent transport equation. The results indicate that the flow field is in excellent agreement to the experimental LDV data and that the platelet adhesion simulations are in good agreement with the in vitro platelet data. These results provide good evidence that this transport equation can be used to express the relationship between blood and a biomaterial if the correct platelet adhesion characteristics are known for the biomaterial. Further validation is necessary with other materials.

Extending the Power-Law Hemolysis Model to Complex Flows

Hemolysis (damage to red blood cells) is a long-standing problem in blood contacting devices, and its prediction has been the goal of considerable research. The most popular model relating hemolysis to fluid stresses is the power-law model, which was developed from experiments in pure shear only. In the absence of better data, this model has been extended to more complex flows by replacing the shear stress in the power-law equation with a von Mises-like scalar stress. While the validity of the scalar stress also remains to be confirmed, inconsistencies exist in its application, in particular, two forms that vary by a factor of 2 have been used. This article will clarify the proper extension of the power law to complex flows in a way that maintains correct results in the limit of pure shear.

Micro and nanotechnologies in heart valve tissue engineering

Due to the increased morbidity and mortality resulting from heart valve diseases, there is a growing demand for off-the-shelf implantable tissue engineered heart valves (TEHVs). Despite the significant progress in recent years in improving the design and performance of TEHV constructs, viable and functional human implantable TEHV constructs have remained elusive. The recent advances in micro and nanoscale technologies including the microfabrication, nano-microfiber based scaffolds preparation, 3D cell encapsulated hydrogels preparation, microfluidic, micro-bioreactors, nano-microscale biosensors as well as the computational methods and models for simulation of biological tissues have increased the potential for realizing viable, functional and implantable TEHV constructs. In this review, we aim to present an overview of the importance and recent advances in micro and nano-scale technologies for the development of TEHV constructs.

Development of a computational model for macroscopic predictions of device-induced thrombosis

While cardiovascular device-induced thrombosis is associated with negative patient outcomes, the convoluted nature of the processes resulting in a thrombus makes the full thrombotic network too computationally expensive to simulate in the complex geometries and flow fields associated with devices. A macroscopic, continuum computational model is developed based on a simplified network, which includes terms for platelet activation (chemical and mechanical) and thrombus deposition and growth in regions of low wall shear stress (WSS). Laminar simulations are performed in a two-dimensional asymmetric sudden expansion geometry and compared with in vitro thrombus size data collected using whole bovine blood. Additionally, the predictive power of the model is tested in a flow cell containing a series of symmetric sudden expansions and contractions. Thrombi form in the low WSS area downstream of the asymmetric expansion and grow into the nearby recirculation region, and thrombus height and length largely remain within 95 % confidence intervals calculated from the in vitro data for 30 min of blood flow. After 30 min, predicted thrombus height and length are 0.94 and 4.32 (normalized by the 2.5 mm step height). Importantly, the model also correctly predicts locations of thrombus deposition observed in the in vitro flow cell of expansions and contractions. As the simulation results, which rely on a greatly reduced model of the thrombotic network, are still able to capture the macroscopic behavior of the full network, the model shows promise for timely predictions of device-induced thrombosis toward optimizing and expediting the device development process.

Computational simulation of platelet interactions in the initiation of stent thrombosis due to stent malapposition

Coronary stenting is one of the most commonly used approaches to open coronary arteries blocked due to atherosclerosis. Stent malapposition can induce thrombosis but the microscopic process is poorly understood. The objective of this study was to determine the platelet-level process by which different extents of stent malapposition affect the initiation of stent thrombosis. We utilized a discrete element model to computationally simulate the transport, adhesion, and activation of thousands of individual platelets and red blood cells during thrombus initiation in stented coronary arteries. Simulated arteries contained a malapposed stent with a specified gap distance (0, 10, 25, 50, or 200 μm) between the struts and endothelium. Platelet-level details of thrombus formation near the proximal-most strut were measured during the simulations. The relationship between gap distance and amount of thrombus in the artery varied depending on different conditions (e.g., amount of dysfunctional endothelium, shear-induced activation of platelets, and thrombogenicity of the strut). Without considering shear-induced platelet activation, the largest gap distance (200 μm) produced no recirculation and less thrombus than the smallest two gap distances (0 and 10 μm) that created recirculation downstream of the strut. However, with the occurrence of shear-induced platelet activation, the largest gap distance produced more thrombus than the two smallest gap distances, but less thrombus than an intermediate gap distance (25 μm). A large gap distance was not necessarily the most thrombogenic, in contrast to implications of some computational fluid dynamics studies. The severity of stent malapposition affected initial stent thrombosis differently depending on various factors related to fluid recirculation, platelet trajectories, shear stress, and endothelial condition.

Repetitive Hypershear Activates and Sensitizes Platelets in a Dose-Dependent Manner

Implantation of mechanical circulatory support (MCS) devices-ventricular assist devices and the total artificial heart-has emerged as a vital therapy for advanced and end-stage heart failure. Unfortunately, MCS patients face the requirement of life-long antiplatelet and anticoagulant therapy to combat thrombotic complications resulting from the dynamic and supraphysiologic shear stress conditions associated with such devices, whose effect on platelet activation is poorly understood. We developed a syringe-capillary viscometer-the "platelet hammer"-that repeatedly exposed platelets to average shear stresses up to 1000 dyne/cm(2) for as short as 25 ms. Platelet activation state was measured using a modified prothrombinase assay, with morphological changes analyzed using scanning electron microscopy. We observed an increase in platelet activation state and post-high shear platelet activation rate, or sensitization, with an increase in stress accumulation (SA), the product of shear stress and exposure time. A significant increase in platelet activation state was observed beyond an SA of 1500 dyne-s/cm(2) , with a marked increase in pseudopod length visible beyond an SA of 1000 dyne-s/cm(2) . Utility of the platelet hammer extends to studies of other shear-dependent pathologies, and may assist development of approaches to enhance the safety and effectiveness of MCS devices and objective antithrombotic pharmacotherapy management.

Physical Characterization and Platelet Interactions under Shear Flows of a Novel Thermoset Polyisobutylene-based Co-polymer

Over the years, several polymers have been developed for use in prosthetic heart valves as alternatives to xenografts. However, most of these materials are beset with a variety of issues, including low material strength, biodegradation, high dynamic creep, calcification, and poor hemocompatibility. We studied the mechanical, surface, and flow-mediated thrombogenic response of poly(styrene-coblock-4-vinylbenzocyclobutene)-polyisobutylene-poly(styrene-coblock-4-vinylbenzocylcobutene) (xSIBS), a thermoset version of the thermoplastic elastomeric polyolefin poly(styrene-block-isobutylene-block-styrene) (SIBS), which has been shown to be resistant to in vivo hydrolysis, oxidation, and enzymolysis. Uniaxial tensile testing yielded an ultimate tensile strength of 35 MPa, 24.5 times greater than that of SIBS. Surface analysis yielded a mean contact angle of 82.05° and surface roughness of 144 nm, which was greater than for poly(ε-caprolactone) (PCL) and poly(methyl methacrylate) (PMMA). However, the change in platelet activation state, a predictor of thrombogenicity, was not significantly different from PCL and PMMA after fluid exposure to 1 dyn/cm(2) and 20 dyn/cm(2). In addition, the number of adherent platelets after 10 dyn/cm(2) flow exposure was on the same order of magnitude as PCL and PMMA. The mechanical strength and low thrombogenicity of xSIBS therefore suggest it as a viable polymeric substrate for fabrication of prosthetic heart valves and other cardiovascular devices.

Hemodynamic and thrombogenic analysis of a trileaflet polymeric valve using a fluid-structure interaction approach

Surgical valve replacement in patients with severe calcific aortic valve disease using either bioprosthetic or mechanical heart valves is still limited by structural valve deterioration for the former and thrombosis risk mandating anticoagulant therapy for the latter. Prosthetic polymeric heart valves have the potential to overcome the inherent material and design limitations of these valves, but their development is still ongoing. The aim of this study was to characterize the hemodynamics and thrombogenic potential of the Polynova polymeric trileaflet valve prototype using a fluid-structure interaction (FSI) approach. The FSI model replicated experimental conditions of the valve as tested in a left heart simulator. Hemodynamic parameters (transvalvular pressure gradient, flow rate, maximum velocity, and effective orifice area) were compared to assess the validity of the FSI model. The thrombogenic footprint of the polymeric valve was evaluated using a Lagrangian approach to calculate the stress accumulation (SA) values along multiple platelet trajectories and their statistical distribution. In the commissural regions, platelets were exposed to the highest SA values because of highest stress levels combined with local reverse flow patterns and vortices. Stress-loading waveforms from representative trajectories in regions of interest were emulated in our hemodynamic shearing device (HSD). Platelet activity was measured using our platelet activation state (PAS) assay and the results confirmed the higher thrombogenic potential of the commissural hotspots. In conclusion, the proposed method provides an in depth analysis of the hemodynamic and thrombogenic performance of the polymer valve prototype and identifies locations for further design optimization.

Mechanical platelet activation potential in abdominal aortic aneurysms

Intraluminal thrombus (ILT) in abdominal aortic aneurysms (AAA) has potential implications to aneurysm growth and rupture risk; yet, the mechanisms underlying its development remain poorly understood. Some researchers have proposed that ILT development may be driven by biomechanical platelet activation within the AAA, followed by adhesion in regions of low wall shear stress. Studies have investigated wall shear stress levels within AAA, but platelet activation potential (AP) has not been quantified. In this study, patient-specific computational fluid dynamic (CFD) models were used to analyze stress-induced AP within AAA under rest and exercise flow conditions. The analysis was conducted using Lagrangian particle-based and Eulerian continuum-based approaches, and the results were compared. Results indicated that biomechanical platelet activation is unlikely to play a significant role for the conditions considered. No consistent trend was observed in comparing rest and exercise conditions, but the functional dependence of AP on stress magnitude and exposure time can have a large impact on absolute levels of anticipated platelet AP. The Lagrangian method obtained higher peak AP values, although this difference was limited to a small percentage of particles that falls below reported levels of physiologic background platelet activation.

Micro-scale blood particulate dynamics using a non-uniform rational B-spline-based isogeometric analysis

The current research presents a novel method in which blood particulates - biconcave red blood cells (RBCs) and spherical cells are modeled using isogeometric analysis, specifically Non-Uniform Rational B-Splines (NURBS) in 3-D. The use of NURBS ensures that even with a coarse representation, the geometry of the blood particulates maintains an accurate description when subjected to large deformations. The fundamental advantage of this method is the coupling of the geometrical description and the stress analysis of the cell membrane into a single, unified framework. Details on the modeling approach, implementation of boundary conditions and the membrane mechanics analysis using isogeometric modeling are presented, along with validation cases for spherical and biconcave cells. Using NURBS - based isogeometric analysis, the behavior of individual cells in fluid flow is presented and analyzed in different flow regimes using as few as 176 elements for a spherical cell and 220 elements for a biconcave RBC. This work provides a framework for modeling a large number of 3-D deformable biological cells, each with its own geometric description and membrane properties. To the best knowledge of the authors, this is the first application of the NURBS - based isogeometric analysis to model and simulate blood particulates in flow in 3D.

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

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