Potential fluid mechanic pathways of platelet activation

Platelet activation is a precursor for blood clotting, which plays leading roles in many vascular complications and causes of death. Platelets can be activated by chemical or mechanical stimuli. Mechanically, platelet activation has been shown to be a function of elevated shear stress and exposure time. These contributions can be combined by considering the cumulative stress or strain on a platelet as it is transported. Here, we develop a framework for computing a hemodynamic-based activation potential that is derived from a Lagrangian integral of strain rate magnitude. We demonstrate that such a measure is generally maximized along, and near to, distinguished material surfaces in the flow. The connections between activation potential and these structures are illustrated through stenotic flow computations. We uncover two distinct structures that may explain observed thrombus formation at the apex and downstream of stenoses. More broadly, these findings suggest fundamental relationships may exist between potential fluid mechanic pathways for mechanical platelet activation and the mechanisms governing their transport.     

Computational modeling of blood component transport related to coronary artery thrombosis in Kawasaki disease

Coronary artery thrombosis is the major risk associated with Kawasaki disease (KD). Long-term management of KD patients with persistent aneurysms requires a thrombotic risk assessment and clinical decisions regarding the administration of anticoagulation therapy. Computational fluid dynamics has demonstrated that abnormal KD coronary artery hemodynamics can be associated with thrombosis. However, the underlying mechanisms of clot formation are not yet fully understood. Here we present a new model incorporating data from patient-specific simulated velocity fields to track platelet activation and accumulation. We use a system of Reaction-Advection-Diffusion equations solved with a stabilized finite element method to describe the evolution of non-activated platelets and activated platelet concentrations [AP], local concentrations of adenosine diphosphate (ADP) and poly-phosphate (PolyP). The activation of platelets is modeled as a function of shear-rate exposure and local concentration of agonists. We compared the distribution of activated platelets in a healthy coronary case and six cases with coronary artery aneurysms caused by KD, including three with confirmed thrombosis. Results show spatial correlation between regions of higher concentration of activated platelets and the reported location of the clot, suggesting predictive capabilities of this model towards identifying regions at high risk for thrombosis. Also, the concentration levels of ADP and PolyP in cases with confirmed thrombosis are higher than the reported critical values associated with platelet aggregation (ADP) and activation of the intrinsic coagulation pathway (PolyP). These findings suggest the potential initiation of a coagulation pathway even in the absence of an extrinsic factor. Finally, computational simulations show that in regions of flow stagnation, biochemical activation, as a result of local agonist concentration, is dominant. Identifying the leading factors to a pro-coagulant environment in each case—mechanical or biochemical—could help define improved strategies for thrombosis prevention tailored for each patient.     

Inhibition of Hb Binding to GP1bα Abrogates Hb-Mediated Thrombus Formation on Immobilized VWF and Collagen under Physiological Shear Stress

BackgroundReports including our own describe that intravascular hemolysis increases the risk of thrombosis in hemolytic disorders. Our recent study shows that plasma Hb concentrations correlate directly with platelet activation in patients with paroxysmal nocturnal hemoglobinuria (PNH). The binding of Hb to glycoprotein1bα (GP1bα) increases platelet activation. A peptide AA1-50, designed from N-terminal amino acid sequence of GP1bα significantly inhibits the Hb binding to GP1bα as well as Hb-induced platelet activation. This study further examined if the Hb-mediated platelet activation plays any significant role in thrombus formation on subendothelium matrix under physiological flow shear stresses and the inhibition of Hb-platelet interaction can abrogate the above effects of Hb.

Methods and Results

Study performed thrombus formation assay in vitro by perfusing whole blood over immobilized VWF or collagen type I in presence of Hb under shear stresses simulating arterial or venous flow. The Hb concentrations ranging from 5 to 10 μM, commonly observed level in plasma of the hemolytic patients including PNH, dose-dependently increased thrombus formation on immobilized VWF under higher shear stress of 25 dyne/cm², but not at 5 dyne/cm². The above Hb concentrations also increased thrombus formation on immobilized collagen under both shear stresses of 5 and 25 dyne/cm². The peptide AA1-50 abrogated invariably the above effects of Hb on thrombus formation.

Conclusions and Significance

This study therefore indicates that the Hb-induced platelet activation plays a crucial role in thrombus formation on immobilized VWF or collagen under physiological flow shear stresses. Thus suggesting a probable role of this mechanism in facilitating thrombosis under hemolytic conditions.     

Simulation of Platelet,Thrombus and Erythrocyte Hydrodynamic Interactions in a 3D Arteriole with In Vivo Comparison

Cylindrical blood vessels, ellipsoid platelets and biconcave-shaped deformable erythrocytes (RBCs) are important participants in hemostasis and thrombosis. However, due to the challenge of combining these components in simulation tools, few simulation studies have included all of them in realistic three-dimensional models. In the present study, we apply a recently developed simulation model to incorporate these components and analyze the flow in a thrombotic tubular arteriole, particularly the detailed hydrodynamic interactions between the thrombus shape, RBCs and platelets. It was found that at certain azimuth positions, the velocity drops in the proximity of both the upstream and downstream edge of the thrombus, which is accompanied by a rapid velocity increase in the narrowed region. The RBCs alter the flow profiles significantly from the typical low Reynolds (Re) number flow, and also enhance the deposition of free flowing platelets onto the thrombus. By evaluating the platelet-thrombus interaction and platelet-RBC interaction together, several mechanisms of platelet deposition augmentation are identified. With in vivo data comparison, our model illustrates the potential of future thrombosis studies that incorporate detailed receptor-ligand adhesion modules.     

Computational haemodynamics in two idealised cerebral wide-necked aneurysms after stent placement

Endovascular stents are being commonly used to treat cerebral wide-necked aneurysms recently. The effect of a stent placed in the parent artery is not only to protect the parent artery from occlusion, due to extension of coils and thrombosis, but also to act as flow diverter to vary the haemodynamics in the aneurysm. In this article, two idealised cerebral wide-necked aneurysms were created, one was sidewall aneurysm with curved parent vessel and the other was terminal aneurysm with the bifurcated parent vessel. The plexiglass models of the two aneurysms were 'treated' with commercial porous intravascular stents. The stented physical models were scanned by Micro-CT and the numerical models of the two idealised cerebral wide-necked aneurysms after stent placement were constructed from the scanned image files. The pulsatile flow of non-Newtonian fluid inside the models was simulated by using computational fluid dynamics package. From the simulated flow dynamics, various haemodynamic characteristics such as velocity contours, wall shear stress and oscillatory shear index (OSI) were computed. The velocity of the jet entering the sacs reduced after stent was deployed across the necks of both sidewall and terminal aneurysms; the wall shear stress on the distal neck of sidewall aneurysm reduced, the wall shear stress on the dome of the terminal aneurysm increased and the OSI on the dome of the terminal aneurysm reduced. Therefore, stent placement not only promotes thrombus formation in both aneurysm models but also reduces the regrowth risk of the sidewall aneurysm and the rupture risk of the terminal aneurysm.     

Computational haemodynamics in two idealised cerebral wide-necked aneurysms after stent placement

Endovascular stents are being commonly used to treat cerebral wide-necked aneurysms recently. The effect of a stent placed in the parent artery is not only to protect the parent artery from occlusion, due to extension of coils and thrombosis, but also to act as flow diverter to vary the haemodynamics in the aneurysm. In this article, two idealised cerebral wide-necked aneurysms were created, one was sidewall aneurysm with curved parent vessel and the other was terminal aneurysm with the bifurcated parent vessel. The plexiglass models of the two aneurysms were ‘treated’ with commercial porous intravascular stents. The stented physical models were scanned by Micro-CT and the numerical models of the two idealised cerebral wide-necked aneurysms after stent placement were constructed from the scanned image files. The pulsatile flow of non-Newtonian fluid inside the models was simulated by using computational fluid dynamics package. From the simulated flow dynamics, various haemodynamic characteristics such as velocity contours, wall shear stress and oscillatory shear index (OSI) were computed. The velocity of the jet entering the sacs reduced after stent was deployed across the necks of both sidewall and terminal aneurysms; the wall shear stress on the distal neck of sidewall aneurysm reduced, the wall shear stress on the dome of the terminal aneurysm increased and the OSI on the dome of the terminal aneurysm reduced. Therefore, stent placement not only promotes thrombus formation in both aneurysm models but also reduces the regrowth risk of the sidewall aneurysm and the rupture risk of the terminal aneurysm.     

A fatigue microcrack alters fluid velocities in a computational model of interstitial fluid flow in cortical bone

Targeted remodeling is activated by fatigue microcracks and plays an important role in maintaining bone integrity. It is widely believed that fluid flow-induced shear stress plays a major role in modulating the mechanotransduction process. Therefore, it is likely that fluid flow-induced shear stress plays a major role in the initiation of the repair of fatigue damage. Since no in vivo measurements of fluid flow within bone exist, computational and mathematical models must be employed to investigate the fluid flow field and the shear stress occurring within cortical bone. We developed a computational fluid dynamic model of cortical bone to examine the effect of a fatigue microcrack on the fluid flow field. Our results indicate that there are alterations in the fluid flow field as far as 150 microm away from the crack, and that at distances farther than this, the fluid flow field is similar to the fluid flow field of intact bone. Through the crack and immediately above and below it, the fluid velocity is higher, while at the lateral edges it is lower than that calculated for the intact model, with a maximum change of 29%. Our results suggest that the presence of a fatigue microcrack can alter the shear stress in regions near the crack. These alterations in shear stress have the potential to significantly alter mechanotransduction and may play a role in the initiation of the repair of fatigue microcracks.     

Developing pulsatile flow in a deployed coronary stent

A major consequence of stent implantation is restenosis that occurs due to neointimal formation. This patho-physiologic process of tissue growth may not be completely eliminated. Recent evidence suggests that there are several factors such as geometry and size of vessel, and stent design that alter hemodynamic parameters, including local wall shear stress distributions, all of which influence the restenosis process. The present three-dimensional analysis of developing pulsatile flow in a deployed coronary stent quantifies hemodynamic parameters and illustrates the changes in local wall shear stress distributions and their impact on restenosis. The present model evaluates the effect of entrance flow, where the stent is placed at the entrance region of a branched coronary artery. Stent geometry showed a complex three-dimensional variation of wall shear stress distributions within the stented region. Higher order of magnitude of wall shear stress of 530 dyn/cm2 is observed on the surface of cross-link intersections at the entrance of the stent. A low positive wall shear stress of 10 dyn/cm2 and a negative wall shear stress of -10 dyn/cm2 are seen at the immediate upstream and downstream regions of strut intersections, respectively. Modified oscillatory shear index is calculated which showed persistent recirculation at the downstream region of each strut intersection. The portions of the vessel where there is low and negative wall shear stress may represent locations of thrombus formation and platelet accumulation. The present results indicate that the immediate downstream regions of strut intersections are areas highly susceptible to restenosis, whereas a high shear stress at the strut intersection may cause platelet activation and free emboli formation.     

Activation and extinction models for platelet adhesion

Adherent platelets are an important part of both thrombus formation and in certain stages of atherogenesis. Platelets can be activated by potent chemicals released from adherent platelets and adhere far more readily than unactivated ones. An analytical and numerical model is presented utilising high Peclet number for the activation and adhesion of platelets in shear flows. The model uses a similarity transformation, which characterises the relationship between convective, diffusive transport and the bulk platelet activating reaction mechanism. A first order surface reaction mechanism is used to model platelet adhesion at the wall (cell) surface. The reduced Damk?hler number, M, characterises the importance of the bulk reaction and includes both convective and diffusive terms. For a high rate of blood flow (M-->0) the activation of platelets can effectively be terminated. In contrast, for (M-->infinity) an inner layer of activated platelets exists with an infinitesimally thin reaction sheet separating activated and non-activated platelets. This characterisation by the Damk?hler number highlights results found clinically, in that thrombus forms in areas of low shear (high M) and in some cases an increased blood flow (low M) can inhibit the activation of platelets completely. The model shows the critical balance that exists between convection, diffusion and reaction.     

Mathematical analysis of mural thrombogenesis. Concentration profiles of platelet-activating agents and effects of viscous shear flow.

The concentration profiles of adenosine diphosphate (ADP), thromboxane A2 (TxA2), thrombin, and von Willebrand factor (vWF) released extracellularly from the platelet granules or produced metabolically on the platelet membrane during thrombus growth, were estimated using finite element simulation of blood flow over model thrombi of various shapes and dimensions. The wall fluxes of these platelet-activating agents were estimated for each model thrombus at three different wall shear rates (100 s-1, 800 s-1, and 1,500 s-1), employing experimental data on thrombus growth rates and sizes. For that purpose, whole human blood was perfused in a parallel-plate flow chamber coated with type l fibrillar human collagen, and the kinetic data collected and analyzed by an EPl-fluorescence video microscopy system and a digital image processor. It was found that thrombin concentrations were large enough to cause irreversible platelet aggregation. Although heparin significantly accelerated thrombin inhibition by antithrombin lll, the remaining thrombin levels were still significantly above the minimum threshold required for irreversible platelet aggregation. While ADP concentrations were large enough to cause irreversible platelet aggregation at low shear rates and for small aggregate sizes, TxA2 concentrations were only sufficient to induce platelet shape change over the entire range of wall shear rates and thrombi dimensions studied. Our results also indicated that the local concentration of vWF multimers released from the platelet alpha-granules could be sufficient to modulate platelet aggregation at low and intermediate wall shear rates (less than 1,000 s-1). The sizes of standing vortices formed adjacent to a growing aggregate and the embolizing stresses and the torque, acting at the aggregate surface, were also estimated in this simulation. It was found that standing vortices developed on both sides of the thrombus even at low wall shear rates. Their sizes increased with thrombus size and wall shear rate, and were largely dependent upon thrombus geometry. The experimental observation that platelet aggregation occurred predominantly in the spaces between adjacent thrombi, confirmed the numerical prediction that those standing vortices are regions of reduced fluid velocities and high concentrations of platelet-activating substances, capable of trapping and stimulating platelets for aggregation. The average shear stress and normal stress, as well as the torque, acting to detach the thrombus, increased with increasing wall shear rate. Both stresses were found to be nearly independent of thrombus size and only weekly dependent upon thrombus geometry. Although both stresses had similar values at low wall shear rates, the average shear stress became the predominant embolizing stress at high wall shear rates.     

Clot Permeability,Agonist Transport,and Platelet Binding Kinetics in Arterial Thrombosis

The formation of wall-adherent platelet aggregates is a critical process in arterial thrombosis. A growing aggregate experiences frictional drag forces exerted on it by fluid moving over or through the aggregate. The magnitude of these forces is strongly influenced by the permeability of the developing aggregate; the permeability depends on the aggregate’s porosity. Aggregation is mediated by formation of ensembles of molecular bonds; each bond involves a plasma protein bridging the gap between specific receptors on the surfaces of two different platelets. The ability of the bonds existing at any time to sustain the drag forces on the aggregate determines whether it remains intact or sheds individual platelets or larger fragments (emboli). We investigate platelet aggregation in coronary-sized arteries using both computational simulations and in vitro experiments. The computational model tracks the formation and breaking of bonds between platelets and treats the thrombus as an evolving porous, viscoelastic material, which moves differently from the background fluid. This relative motion generates drag forces which the fluid and thrombus exert on one another. These forces are computed from a permeability-porosity relation parameterized from experimental measurements. Basing this relation on measurements from occlusive thrombi formed in our flow chamber experiments, along with other physiological parameter values, the model produced stable dense thrombi on a similar timescale to the experiments. When we parameterized the permeability-porosity relation using lower permeabilities reported by others, bond formation was insufficient to balance drag forces on an early thrombus and keep it intact. Under high shear flow, soluble agonist released by platelets was limited to the thrombus and a boundary layer downstream, thus restricting thrombus growth into the vessel lumen. Adding to the model binding and activation of unactivated platelets through von Willebrand-factor-mediated processes allowed greater growth and made agonist-induced activation more effective.     

An integrated fluid–structure interaction and thrombosis model for type B aortic dissection

False lumen thrombosis (FLT) in type B aortic dissection has been associated with the progression of dissection and treatment outcome. Existing computational models mostly assume rigid wall behavior which ignores the effect of flap motion on flow and thrombus formation within the FL. In this study, we have combined a fully coupled fluid–structure interaction (FSI) approach with a shear-driven thrombosis model described by a series of convection–diffusion reaction equations. The integrated FSI-thrombosis model has been applied to an idealized dissection geometry to investigate the interaction between vessel wall motion and growing thrombus. Our simulation results show that wall compliance and flap motion can influence the progression of FLT. The main difference between the rigid and FSI models is the continuous development of vortices near the tears caused by drastic flap motion up to 4.45 mm. Flap-induced high shear stress and shear rates around tears help to transport activated platelets further to the neighboring region, thus speeding up thrombus formation during the accelerated phase in the FSI models. Reducing flap mobility by increasing the Young’s modulus of the flap slows down the thrombus growth. Compared to the rigid model, the predicted thrombus volume is 25% larger using the FSI-thrombosis model with a relatively mobile flap. Furthermore, our FSI-thrombosis model can capture the gradual effect of thrombus growth on the flow field, leading to flow obstruction in the FL, increased blood viscosity and reduced flap motion. This model is a step closer toward simulating realistic thrombus growth in aortic dissection, by taking into account the effect of intimal flap and vessel wall motion.

     

Plant food delphinidin-3-glucoside significantly inhibits platelet activation and thrombosis: novel protective roles against cardiovascular diseases

Delphinidin-3-glucoside (Dp-3-g) is one of the predominant bioactive compounds of anthocyanins in many plant foods. Although several anthocyanin compounds have been reported to be protective against cardiovascular diseases (CVDs), the direct effect of anthocyanins on platelets, the key players in atherothrombosis, has not been studied. The roles of Dp-3-g in platelet function are completely unknown. The present study investigated the effects of Dp-3-g on platelet activation and several thrombosis models in vitro and in vivo. We found that Dp-3-g significantly inhibited human and murine platelet aggregation in both platelet-rich plasma and purified platelets. It also markedly reduced thrombus growth in human and murine blood in perfusion chambers at both low and high shear rates. Using intravital microscopy, we observed that Dp-3-g decreased platelet deposition, destabilized thrombi, and prolonged the time required for vessel occlusion. Dp-3-g also significantly inhibited thrombus growth in a carotid artery thrombosis model. To elucidate the mechanisms, we examined platelet activation markers via flow cytometry and found that Dp-3-g significantly inhibited the expression of P-selectin, CD63, CD40L, which reflect platelet α- and δ-granule release, and cytosol protein secretion, respectively. We further demonstrated that Dp-3-g downregulated the expression of active integrin αIIbβ3 on platelets, and attenuated fibrinogen binding to platelets following agonist treatment, without interfering with the direct interaction between fibrinogen and integrin αIIbβ3. We found that Dp-3-g reduced phosphorylation of adenosine monophosphate-activated protein kinase, which may contribute to the observed inhibitory effects on platelet activation. Thus, Dp-3-g significantly inhibits platelet activation and attenuates thrombus growth at both arterial and venous shear stresses, which likely contributes to its protective roles against thrombosis and CVDs.     

Introducing the pro-coagulant contact system in the numerical assessment of device-related thrombosis

Thrombosis is a major concern in blood-coated medical devices. Contact activation, which is the initial part of the coagulation cascade in device-related thrombosis, is not considered in current thrombus formation models. In the present study, pro-coagulant reactions including the contact activation system are coupled with a fluid solver in order to evaluate the potential of the contact system to initiate thrombin production. The biochemical/fluid model is applied to a backward-facing step configuration, a flow configuration that frequently appears in medical devices. In contrast to the in vivo thrombosis models in which a specific thrombotic zone (injury region) is set a priori by the user to initiate the coagulation reaction, a reactive surface boundary condition is applied to the whole device wall. Simulation results show large thrombin concentration in regions related to recirculation zones without the need of an a priori knowledge of the thrombus location. The numerical results align well with the regions prone to thrombosis observed in experimental results reported in the literature. This approach could complement thrombus formation models that take into account platelet activity and thrombus growth to optimize a wide range of medical devices.     

The hemodynamic forces acting on thrombi, from incipient attachment of single cells to maturity and embolization

We consider the steady fluid forces acting on a thrombus from the time of first contact of a single cell with a natural or artificial surface, through the attachment process and growth to embolization. For a hemi-spherical or cylindrical attached cell of height less than 1/100-1/20th of the channel width, shear and tensile stresses are solely dependent on viscosity and on the ratio of average fluid velocity to channel width vt/Dt (shear rate). Large values of this ratio reduce adhesion and increase embolization. The average shear stress on such cells is approximately 1-10 Pa (10-100 dyn cm2), the average tensile stress about three times higher. For other shapes and larger protrusions, stress varies with protrusion height as well. Maturing thrombi composed of cell aggregates embedded in a fibrin mesh do not appear to allow significant fluid flow through their porous structure. The interior forces are then due solely to hydrostatic pressure and initially vary directly with vt/Dt and inversely with thrombus height Hp, thus favouring embolization at an early stage and in arterial systems. Rough surfaces are identified as causing an increase in dwell-time and possibly immobilizing an unattached cell due to 'negative lift'.     

Pathologic high shear stress induces apoptosis events in human platelets

Recently, it has been discovered that apoptosis of anucleate platelets can be induced by chemical agonists. Other studies demonstrated that mechanical forces (shear stresses) stimulate platelet activation and signaling in the absence of exogenous chemical stimuli. We analyzed whether shear stresses can trigger platelet apoptosis, a question that has not yet been studied. Using a cone-and-plate viscometer, we exposed human platelet-rich plasma to different shear stresses, ranging from physiologic arterial and arteriole levels (10-44 dyn/cm2) to pathologic high levels (117-388 dyn/cm2) occurring in stenotic vessels. We found that pathologic shear stresses induce not only platelet activation (P-selectin upregulation and GPIbalpha downregulation) but also trigger apoptosis events, including mitochondrial transmembrane potential depolarization, caspase 3 activation, phosphatidylserine exposure, and platelet shrinkage and fragmentation, whereas physiological shear stresses are not effective.This novel finding suggests that shear-induced platelet apoptosis can be mediated by mechanoreceptors, does not require nuclear participation, and may affect platelet clearance.     

Continuous Modeling of Arterial Platelet Thrombus Formation Using a Spatial Adsorption Equation

In this study, we considered a continuous model of platelet thrombus growth in an arteriole. A special model describing the adhesion of platelets in terms of their concentration was derived. The applications of the derived model are not restricted to only describing arterial platelet thrombus formation; the model can also be applied to other similar adhesion processes. The model reproduces an auto-wave solution in the one-dimensional case; in the two-dimensional case, in which the surrounding flow is taken into account, the typical torch-like thrombus is reproduced. The thrombus shape and the growth velocity are determined by the model parameters. We demonstrate that the model captures the main properties of the thrombus growth behavior and provides us a better understanding of which mechanisms are important in the mechanical nature of the arterial thrombus growth.     

Cytometric analysis of high shear-induced platelet microparticles and effect of cytokines on microparticle generation

BACKGROUND: Microparticles released from platelets may play a role in the normal hemostatic response to vascular injury, because they exhibit prothrombinase activity. Microparticles are generated by high shear stress and may be formed in diseased small arteries and arterioles in various clinical settings. However, the surface composition of high shear-induced platelet microparticles is unknown. It was recently shown that some cytokines modulate platelet activation. However, no reports are available concerning the effect of cytokines on high shear-induced platelet aggregation (SIPA) microparticle generation. MATERIALS AND METHODS: Measurement of SIPA was performed with a cone-plate viscometer. The conformational characteristics of high shear (108 dynes/cm(2))-induced platelet microparticles were analyzed by flow cytometry and confocal laser scanning microscopy. Effects of cytokines for high SIPA microparticle generation were also analyzed using flow cytometry. RESULTS: The overall pattern of monoclonal antibody binding in high shear-induced microparticles was almost the same as that in activated platelets under high shear stress. Microparticles exhibited markedly increased Annexin V binding. In fluorescent confocal images, small and fine regions of fluorescence (microparticles) were recognized separate from platelet fluorescence. Thrombopoietin not only induced platelet activation, as demonstrated by CD62P expression, but also increased the number of microparticles. Erythropoietin and interleukin-6 enhanced only microparticle generation. CONCLUSIONS: These results suggest that microparticles possessing procoagulant activity are released by platelet activation when levels of certain cytokines increase under high shear stress in various clinical settings.     

Comparison of the hemodynamic and thrombogenic performance of two bileaflet mechanical heart valves using a CFD/FSI model

The hemodynamic and the thrombogenic performance of two commercially available bileaflet mechanical heart valves (MHVs)--the ATS Open Pivot Valve (ATS) and the St. Jude Regent Valve (SJM), was compared using a state of the art computational fluid dynamics-fluid structure interaction (CFD-FSI) methodology. A transient simulation of the ATS and SJM valves was conducted in a three-dimensional model geometry of a straight conduit with sudden expansion distal the valves, including the valve housing and detailed hinge geometry. An aortic flow waveform (60 beats/min, cardiac output 4 l/min) was applied at the inlet. The FSI formulation utilized a fully implicit coupling procedure using a separate solver for the fluid problem (FLUENT) and for the structural problem. Valve leaflet excursion and pressure differences were calculated, as well as shear stress on the leaflets and accumulated shear stress on particles released during both forward and backward flow phases through the open and closed valve, respectively. In contrast to the SJM, the ATS valve opened to less than maximal opening angle. Nevertheless, maximal and mean pressure gradients and velocity patterns through the valve orifices were comparable. Platelet stress accumulation during forward flow indicated that no platelets experienced a stress accumulation higher than 35 dyne x s/cm2, the threshold for platelet activation (Hellums criterion). However, during the regurgitation flow phase, 0.81% of the platelets in the SJM valve experienced a stress accumulation higher than 35 dyne x s/cm2, compared with 0.63% for the ATS valve. The numerical results indicate that the designs of the ATS and SJM valves, which differ mostly in their hinge mechanism, lead to different potential for platelet activation, especially during the regurgitation phase. This numerical methodology can be used to assess the effects of design parameters on the flow induced thrombogenic potential of blood recirculating devices.     

Slowness of blood flow and resultant thrombus formation in cerebral aneurysms

It is not yet fully understood what causes cerebral aneurysms to rupture. Although no definitive conclusion has been reached, it is considered that there are haemodynamic, biochemical and physiological factors contributing to rupture. Numerical techniques seem promising for investigation of this multi-physical phenomenon. In fact, recent intensive numerical studies with computational fluid dynamics have revealed detailed haemodynamic features of the flow in cerebral aneurysms such as velocity, pressure and wall shear stress distributions. It is, therefore, expected that biochemical and physiological aspects of aneurysmal rupture will also be actively investigated using numerical approaches. Considering this background, the authors have been working on modelling of thrombus formation in cerebral aneurysms caused by stagnant blood flow, because many studies have suggested that slow blood flow and resulting low wall shear stress are connected with rupture. Firstly, in this review paper, slowness of the intra-aneurysmal flow is reviewed with an energy balance theory, and secondly, thrombus formation in cerebral bifurcation aneurysms is discussed from the viewpoint of numerical modelling. A computational result obtained by application of the authors’ platelet aggregation–adhesion model is also provided.