Assessment of Rheological Models Applied to Blood Flow in Human Thoracic Aorta.

PURPOSE: The purpose of this study is to assess the importance of non-Newtonian rheological models on blood flow in the human thoracic aorta. METHODS: The pulsatile flow in the aorta is simulated using the models of Casson, Quemada and Walburn-Schneck in addition to a case of fixed (Newtonian) viscosity. The impact of the four rheological models (using constant hematocrit) was assessed with respect to (i) magnitude and deviation of the viscosity relative to a reference value (the Newtonian case); (ii) wall shear stress (WSS) and its time derivative; (iii) common WSS-related indicators, OSI, TAWSS and RRT; (iv) relative volume and surface-based retrograde flow; and (v) the impact of rheological models on the transport of small particles in the thoracic aorta. RESULTS: The time-dependent flow in the thoracic aorta implies relatively large variations in the instantaneous WSS, due to variations in the instantaneous viscosity by as much as an order of magnitude. The largest effect was observed for low shear rates (tens s-1). The different viscosity models had a small impact in terms of time- and spaced-averaged quantities. The significance of the rheological models was clearly demonstrated in the instantaneous WSS, for the space-averaged WSS (about 10%) and the corresponding temporal derivative of WSS (up to 20%). The longer-term accumulated effect of the rheological model was observed for the transport of spherical particles of 2 mm and 2 mm in diameter (density of 1200 kg/m3). Large particles' total residence time in the brachiocephalic artery was 60% longer compared to the smaller particles. For the left common carotid artery, the opposite was observed: the smaller particles resided considerably longer than their larger counterparts. CONCLUSIONS: The dependence on the non-Newtonian properties of blood is mostly important at low shear regions (near walls, stagnation regions). Time- and space-averaging parameters of interest reduce the impact of the rheological model and may thereby lead to under-estimation of viscous effects. The rheological model affects the local WSS and its temporal derivative. In addition, the transport of small particles includes the accumulated effect of the blood rheological model as the several forces (e.g., drag, added mass and lift) acting on the particles are viscosity dependent. Mass transport is an essential factor for the development of pathologies in the arterial wall, implying that rheological models are important for assessing such risks.

Red blood cell distribution width is associated with increased interactions of blood cells with vascular wall.

The mechanism underlying the association between elevated red cell distribution width (RDW) and poor prognosis in variety of diseases is unknown although many researchers consider RDW a marker of inflammation. We hypothesized that RDW directly affects intravascular hemodynamics, interactions between circulating cells and vessel wall, inducing local changes predisposing to atherothrombosis. We applied different human and animal models to verify our hypothesis. Carotid plaques harvested from patients with high RDW had increased expression of genes and proteins associated with accelerated atherosclerosis as compared to subjects with low RDW. In microfluidic channels samples of blood from high RDW subjects showed flow pattern facilitating direct interaction with vessel wall. Flow pattern was also dependent on RDW value in mouse carotid arteries analyzed with Magnetic Resonance Imaging. In different mouse models of elevated RDW accelerated development of atherosclerotic lesions in aortas was observed. Therefore, comprehensive biological, fluid physics and optics studies showed that variation of red blood cells size measured by RDW results in increased interactions between vascular wall and circulating morphotic elements which contribute to vascular pathology.

Flows and function of the Infant Flow neonatal continuous positive airway pressure device investigated with computational fluid dynamics.

AIM: The first dedicated neonatal continuous positive airway pressure (CPAP) device using variable flow was the Infant Flow. The system was pressure stable with a low resistance to breathing. The aim of this study was to describe the flow and function of the Infant Flow geometry using simulated breathing and computational fluid dynamics. METHOD: The original Infant Flow geometry was used with simulated term infant breathing at three levels of CPAP. The large eddy simulation methodology was applied in combination with the WALE sub-grid scale model. RESULTS: The simulation fully resolved the flow phenomena in the Infant Flow geometry. The main flow feature during inspiration was support by gas entrainment and mixing. During expiration, the jet deflected towards the outlet with unstable impingement of the jet at the opposing edge. CONCLUSION: The proposed mechanism has previously been based on theoretical reasoning, and our results present the first detailed description of the Infant Flow. The pressure stability was based on a jet supporting inspiration by gas entrainment and then being deflected during exhalation. This confirmed previously assumed principles of function and flows within the geometry and provided a base for further developments.

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.

Detection of Thrombosis in the Extracorporeal Membrane Oxygenation Circuit by Infrasound: Proof of Concept.

As of today, there exist no reliable, objective methods for early detection of thrombi in the extracorporeal membrane oxygenators (ECMO) system. Within the ECMO system, thrombi are not always fixed to a certain component or location in the circuit. Thus, clot fragments of different shapes and consistencies may circulate and give rise to vibrations and sound generation. By bedside sound measurements and additional laboratory experiments (although not detailed herein), we found that the presence of particles (clots or aggregates and fragments of clots) can be detected by analyzing the strength of infra-sound (< 20 Hz) modes of the spectrum near the inlet and outlet of the centrifugal pump in the ECMO circuit. For the few patients that were considered in this study, no clear false positive or negative examples were found when comparing the spectral approach with clinical observations. A laboratory setup provided insight to the flow in and out of the pump, confirming that in the presence of particles a low-amplitude low-frequency signal is strongly amplified, enabling the identification of a clot.