Local Hemodynamics Research Articles

Semi-analytical solution for pressure-driven non-Newtonian flows in concentric annuli

Hemodynamics during cardiac catheterization can be modeled as non-Newtonian fluid flow within concentric annuli. Rapid and precise evaluation of flow rate variations is crucial for understanding the catheter's impact on local hemodynamics. We propose a general semi-analytical solution for flow rates in pressure-driven non-Newtonian flows within concentric annuli. This solution is suitable for commonly used non-Newtonian models. Validation using the Casson and Carreau–Yasuda models across various diameter ratios confirms the accuracy of the proposed solution, with a maximum relative error of 3.55%. It also reliably predicts time-averaged flow rates in pulsatile flows, with a maximum relative error of only 0.53%.

In vivo volumetric analysis of retinal vascular hemodynamics in mice with spatio-temporal optical coherence tomography.

Microcirculation and neurovascular coupling are important parameters to study in neurological and neuro-ophthalmic conditions. As the retina shares many similarities with the cerebral cortex and is optically accessible, a special focus is directed to assessing the chorioretinal structure, microvasculature, and hemodynamics of mice, a vital animal model for vision and neuroscience research. We aim to introduce an optical imaging tool enabling in vivo volumetric mouse retinal monitoring of vascular hemodynamics with high temporal resolution. We translated the spatio-temporal optical coherence tomography (STOC-T) technique into the field of small animal imaging by designing a new optical system that could compensate for the mouse eye refractive error. We also developed post-processing algorithms, notably for the assessment of (i) localized hemodynamics from the analysis of pulse wave-induced Doppler artifact modulation and (ii) retinal tissue displacement from phase-sensitive measurements. We acquired high-quality, in vivo volumetric mouse retina images at a rate of 113Hz over a lateral field of view of . We presented high-resolution en face images of the retinal and choroidal structure and microvasculature from various layers, after digital aberration correction. We were able to measure the pulse wave velocity in capillaries of the outer plexiform layer with a mean speed of 0.35mm/s and identified venous and arterial pulsation frequency and phase delay. We quantified the modulation amplitudes of tissue displacement near major vessels (with peaks of 150nm), potentially carrying information about the biomechanical properties of the retinal layers involved. Last, we identified the delays between retinal displacements due to the passing of venous and arterial pulse waves. The developed STOC-T system provides insights into the hemodynamics of the mouse retina and choroid that could be beneficial in the study of neurovascular coupling and vasculature and flow speed anomalies in neurological and neuro-ophthalmic conditions.

Radiomic profiling of high-risk aneurysms with blebs: an exploratory study

BackgroundBlebs significantly increase rupture risk of intracranial aneurysms. Radiomic analysis offers a robust characterization of the aneurysm wall. However, the unique radiomic profile of various compartments, including blebs, remains unexplored....

A generalized reduced-order model for trans-stenotic pressure drop with and without a guidewire

Guidewire-based pressure measurement is essential for diagnosing coronary artery disease. However, the impact of the guidewire on local hemodynamics and diagnostic outcomes is not fully understood. In this study, we propose a generalized reduced-order model (ROM) to accurately predict the trans-stenotic pressure drop in arteries. A key advantage of this model is that the viscous term does not rely on empirical parameters, making it applicable to both scenarios with and without guidewire insertion, and across varying stenosis severities. The proposed model demonstrates good accuracy compared to 3D idealized numerical models, achieving an average prediction error of 3.61% for cases without a guidewire and 4.53% for cases with a guidewire. Furthermore, when applied to a patient-specific model, it achieves comparable or better results than previously published ROMs. Finally, this ROM is employed to investigate the shifting relative importance of different components of the trans-stenotic pressure drop at various stenosis severities, and to provide further insights into the guidewire’s influence on FFR measurements.

Vascular endothelial cell morphology and alignment regulate VEGF-induced endothelial nitric oxide synthase activation.

Nitric oxide (NO) production by endothelial nitric oxide synthase (eNOS) inhibits platelet and leukocyte adhesion while promoting vasorelaxation in smooth muscle cells. Dysfunctional regulation of eNOS is a hallmark of various vascular pathologies, notably atherosclerosis, often associated with areas of low shear stress on endothelial cells (ECs). While the link between EC morphology and local hemodynamics is acknowledged, the specific impact of EC morphology on eNOS regulation remains unclear. Morphological differences between elongated, aligned ECs and polygonal, randomly oriented ECs correspond to variations in focal adhesion and cytoskeletal organization, suggesting differing levels of cytoskeletal prestress. However, the functional outcomes of cytoskeletal prestress, particularly in the absence of shear stress, are not extensively studied in ECs. Some evidence suggests that elongated ECs exhibit decreased immunogenicity and enhanced NO production. This study aims to elucidate the signaling pathways governing VEGF-stimulated eNOS regulation in the aligned EC phenotype characterized by elongated and aligned cells within a monolayer. Using anisotropic topographic cues, bovine aortic endothelial cells (BAECs) were elongated and aligned, followed by VEGF treatment in the presence or absence of cytoskeletal tension inhibitors. Phosphorylation of eNOS ser1179, AKT ser437 and FAK Tyr397 in response to VEGF challenge were significantly heightened in aligned ECs compared to unaligned ECs. Moreover this response proved to be robustly tied to cytoskeletal tension as evinced by the abrogation of responses in the presence of the myosin II ATPase inhibitor, blebbistatin. Notably, this work demonstrates for the first time the reliance on FAK phosphorylation in VEGF-mediated eNOS activation and the comparatively greater contribution of the cytoskeletal machinery in propagating VEGF-eNOS signaling in aligned and elongated ECs. This research underscores the importance of utilizing appropriate vascular models in drug development and sheds light on potential mechanisms underlying vascular function and pathology that can help inform vascular graft design.

The hemodynamic and geometric characteristics of carotid artery atherosclerotic plaque formation.

Ischemic stroke, which has a high incidence, disability, and mortality rate, is mainly caused by carotid atherosclerotic plaque. The difference in the geometric structures of the carotid arteries inevitably leads to the variability in the local hemodynamics, which plays a key role in the formation of carotid atherosclerosis. At present, the combined mechanisms of hemodynamic and geometric in the formation of carotid atherosclerotic plaque are not clear. Thus, this study characterized the geometric and hemodynamic characteristics of carotid atherosclerotic plaque formation using four-dimensional (4D) flow magnetic resonance imaging (MRI). Ultimately, 122 carotid arteries from 61 patients were examined in this study. According to the presence of plaques at the bifurcation of the carotid artery on cervical vascular ultrasound (US), carotid arteries were placed into a plaque group (N=69) and nonplaque group (N=53). The ratio of the maximum internal carotid artery (ICA) inner diameter to the maximum common carotid artery (CCA) inner diameter (ICA-CCA diameter ratio), bifurcation angle, and tortuosity were measured using neck three-dimensional time-of-flight magnetic resonance angiography (3D TOF-MRA). Meanwhile, 4D flow MRI was used to obtain the following hemodynamic parameters of the carotid arteries: volume flow rate, velocity, wall shear stress (WSS), and pressure gradient (PG). Independent sample t-tests were used to compare carotid artery geometry and hemodynamic changes between the plaque group and nonplaque group. The ICA-CCA diameter ratio between the plaque group and the nonplaque group was not significantly different (P=0.124), while there were significant differences in the bifurcation angle (P=0.005) and tortuosity (P=0.032). The bifurcation angle of the plaque group was greater than that of the nonplaque group (60.70°±20.75° vs. 49.32°±22.90°), and the tortuosity was smaller than that of the nonplaque group (1.07±0.04 vs. 1.09±0.05). There were no significant differences between the two groups in terms of volume flow rate (P=0.351) and the maximum value of velocity (velocitymax) (P=0.388), but the axial, circumferential, and 3D WSS values were all significantly different, including their mean values (all P values <0.001) and the maximum value of 3D WSS (P<0.001), with the mean axial, circumferential, 3D WSS values, along with the maximum 3D WSS value, being lower in the plaque group. The two groups also differed significantly in terms of maximum PG value (P=0.030) and mean PG value (P=0.026), with these values being greater in the nonplaque group than in the plaque group. A large bifurcation angle and a low tortuosity of the carotid artery are geometric risk factors for plaque formation in this area. Low WSS and low PG values are associated with carotid atherosclerotic plaque formation.

Evolution of Natural Myocardial Shear Wave Behavior in Young Hearts: Determinant Factors and Reproducibility Analysis

Myocardial diastolic function assessment in children by conventional echocardiography is challenging. High-frame rate echocardiography facilitates the assessment of myocardial stiffness, a key factor in diastolic function, by measuring the propagation velocities of myocardial shear waves (SWs). However, normal values of natural SWs in children are currently lacking. The aim of this study was to explore the behavior of natural SWs among children and adolescents, their reproducibility, and the factors affecting SW velocities from childhood into adulthood. One hundred six healthy children (2-18years of age) and 62 adults (19-80years of age) were recruited. High-frame rate images were acquired using a modified commercial scanner. An anatomic M-mode line was drawn along the ventricular septum, and propagation velocities of natural SWs after mitral valve closure were measured in the tissue acceleration-coded M-mode display. Throughout life, SW velocities after mitral valve closure exhibited pronounced age dependency (r=0.73; P<.001). Among the pediatric population, SW velocities correlated significantly with measures of cardiac geometry (septal thickness and left ventricular end-diastolic dimension), local hemodynamics (systolic blood pressure), and echocardiographic parameters of systolic and diastolic function (global longitudinal strain, mitral E/e' ratio, isovolumic relaxation time, and mitral deceleration time) (P<.001). In a multivariate analysis including all these factors, the predictors of SW velocities were age, mitral E/e', and global longitudinal strain (r=0.81). Natural myocardial SW velocities in children can be detected and measured. SW velocities showed significant dependence on age and diastolic function. Natural SWs could be a promising additive tool for the assessment of diastolic function among children.

Clinical Impact Of Assessing Plaque Quantification And Local Hemodynamics Of Intermediate Lesions In The Prediction Of Acute Coronary Syndrome

Clinical Impact Of Assessing Plaque Quantification And Local Hemodynamics Of Intermediate Lesions In The Prediction Of Acute Coronary Syndrome

Clinical Impact Of Assessing Plaque Quantification And Local Hemodynamics Of Intermediate Lesions In The Prediction Of Acute Coronary Syndrome

Clinical Impact Of Assessing Plaque Quantification And Local Hemodynamics Of Intermediate Lesions In The Prediction Of Acute Coronary Syndrome

Role of the vessel morphology on the lenticulostriate arteries hemodynamics during atrial fibrillation: A CFD-based multivariate regression analysis

Background and objectiveAtrial fibrillation (AF) is the most common cardiac arrhythmia, inducing accelerated and irregular beating. Beside well-known disabling symptoms - such as palpitations, reduced exercise tolerance, and chest discomfort - there is growing evidence that an alteration of deep cerebral hemodynamics due to AF increases the risk of vascular dementia and cognitive impairment, even in the absence of clinical strokes. The alteration of deep cerebral circulation in AF represents one of the least investigated among the possible mechanisms. Lenticulostriate arteries (LSAs) are small perforating arteries mainly departing from the middle cerebral artery (MCA) and susceptible to small vessel disease, which is one of the mechanisms of subcortical vascular dementia development. The purpose of this study is to investigate the impact of different LSAs morphologies on the cerebral hemodynamics during AF. MethodsBy combining a computational fluid dynamics (CFD) analysis of LSAs with 7T high-resolution magnetic resonance imaging (MRI), we performed different CFD-based multivariate regression analyses to detect which geometrical and morphological vessel features mostly affect AF hemodynamics in terms of wall shear stress. We exploited 17 cerebral 7T-MRI derived LSA vascular geometries extracted from 10 subjects and internal carotid artery data from validated 0D cardiovascular-cerebral modeling as inflow conditions. ResultsOur results revealed that few geometrical variables - namely the size of the MCA and the bifurcation angles between MCA and LSA - are able to satisfactorily predict the AF impact. In particular, the present study indicates that LSA morphologies exhibiting markedly obtuse LSA-MCA inlet angles and small MCA size downstream of the LSA-MCA bifurcation may be more prone to vascular damage induced by AF. ConclusionsThe present MRI-based computational study has been able for the first time to: (i) investigate the net impact of LSAs vascular morphologies on cerebral hemodynamics during AF events; (ii) detect which combination of morphological features worsens the hemodynamic response in the presence of AF. Awaiting necessary clinical confirmation, our analysis suggests that the local hemodynamics of LSAs is affected by their geometrical features and some LSA morphologies undergo greater hemodynamic alterations in the presence of AF.

Distribution of rupture sites and blebs on intracranial aneurysm walls suggests distinct rupture patterns in ACom and MCA aneurysms.

The mechanisms behind intracranial aneurysm formation and rupture are not fully understood, with factors such as location, patient demographics, and hemodynamics playing a role. Additionally, the significance of anatomical features like blebs in ruptures is debated. This highlights the necessity for comprehensive research that combines patient-specific risk factors with a detailed analysis of local hemodynamic characteristics at bleb and rupture sites. Our study analyzed 359 intracranial aneurysms from 268 patients, reconstructing patient-specific models for hemodynamic simulations based on 3D rotational angiographic images and intraoperative videos. We identified aneurysm subregions and delineated rupture sites, characterizing blebs and their regional overlap, employing statistical comparisons across demographics, and other risk factors. This work identifies patterns in aneurysm rupture sites, predominantly at the dome, with variations across patient demographics. Hypertensive and anterior communicating artery (ACom) aneurysms showed specific rupture patterns and bleb associations, indicating two pathways: high-flow in ACom with thin blebs at impingement sites and low-flow, oscillatory conditions in middle cerebral artery (MCA) aneurysms fostering thick blebs. Bleb characteristics varied with gender, age, and smoking, linking rupture risks to hemodynamic factors and patient profiles. These insights enhance understanding of the hemodynamic mechanisms leading to rupture events. This analysis elucidates the role of localized hemodynamics in intracranial aneurysm rupture, challenging the emphasis on location by revealing how flow variations influence stability and risk. We identify two pathways to wall failure-high-flow and low-flow conditions-highlighting the complexity of aneurysm behavior. Additionally, this research advances our knowledge of how inherent patient-specific characteristics impact these processes, which need further investigation.

Linking wall shear stress and vorticity topologies: Toward a unified theory of cardiovascular flow disturbances

Broadening current knowledge about the complex relationship at the blood-vessel wall interface is a main challenge in hemodynamics research. Moving from the consideration that wall shear stress (WSS) provides a signature for the near-wall velocity dynamics and vorticity is considered the skeleton of fluid motion, here we present a unified theory demonstrating the existing link between surface vorticity (SV) and WSS topological skeletons, the latter recently emerged as a predictor of vascular disease. The analysis focused on WSS and SV fixed points, i.e., points where the fields vanish, as they play a major role in shaping the main vector field features. The theoretical analysis proves that: (i) all SV fixed points on the surface must necessarily be WSS fixed points, although with differences in nature and stability and (ii) a WSS fixed point is not necessarily a SV fixed point. In the former case, WSS fixed points are the consequence of flow patterns where only shear contributes to vorticity; in the latter case, WSS fixed points are the consequence of flow impingement to/emanation from the vessel wall. Moreover, fluid structures interacting with the wall characterized by zero or non-zero rotational momentum generate WSS fixed points of different nature/stability. High-fidelity computational fluid dynamics simulations in intracranial aneurysm models confirmed the applicability of the theoretical considerations. The presented unified theory unambiguously explains the mechanistic link between near-wall flow disturbances and the underlying intravascular flow features expressed in terms of vorticity, ultimately facilitating a clearer interpretation of the role of local hemodynamics in vascular pathophysiology.

Local hemodynamics effects of ischemic preconditioning on non-fatiguing voluntary contractions

Local hemodynamics effects of ischemic preconditioning on non-fatiguing voluntary contractions

Validation of ultrasound velocimetry and computational fluid dynamics for flow assessment in femoral artery stenotic disease.

To investigate the accuracy of high-framerate echo particle image velocimetry (ePIV) and computational fluid dynamics (CFD) for determining velocity vectors in femoral bifurcation models through comparison with optical particle image velocimetry (oPIV). Separate femoral bifurcation models were built for oPIV and ePIV measurements of a non-stenosed (control) and a 75%-area stenosed common femoral artery. A flow loop was used to create triphasic pulsatile flow. In-plane velocity vectors were measured with oPIV and ePIV. Flow was simulated with CFD using boundary conditions from ePIV and additional duplex-ultrasound (DUS) measurements. Mean differences and 95%-limits of agreement (1.96*SD) of the velocity magnitudes in space and time were compared, and the similarity of vector complexity (VC) and time-averaged wall shear stress (TAWSS) was assessed. Similar flow features were observed between modalities with velocities up to 110 and in the control and the stenosed model, respectively. Relative to oPIV, ePIV and CFD-ePIV showed negligible mean differences in velocity (), with limits of agreement of (control) and (stenosed). CFD-DUS overestimated velocities with limits of agreements of and for the control and stenosed model, respectively. VC showed good agreement, whereas TAWSS showed similar trends but with higher values for ePIV, CFD-DUS, and CFD-ePIV compared to oPIV. EPIV and CFD-ePIV can accurately measure complex flow features in the femoral bifurcation and around a stenosis. CFD-DUS showed larger deviations in velocities making it a less robust technique for hemodynamical assessment. The applied ePIV and CFD techniques enable two- and three-dimensional assessment of local hemodynamics with high spatiotemporal resolution and thereby overcome key limitations of current clinical modalities making them an attractive and cost-effective alternative for hemodynamical assessment in clinical practice.

Subjects with carotid webs demonstrate pro-thrombotic hemodynamics compared to subjects with carotid atherosclerosis

Carotid artery webs (CaW) are non-atherosclerotic projections into the vascular lumen and have been linked to up to one-third of cryptogenic strokes in younger patients. Determining how CaW affects local hemodynamics is essential for understanding clot formation and stroke risk. Computational fluid dynamics simulations were used to investigate patient-specific hemodynamics in carotid artery bifurcations with CaW, bifurcations with atherosclerotic lesions having a similar degree of lumen narrowing, and with healthy carotid bifurcations. Simulations were conducted using segmented computed tomography angiography geometries with inlet boundary conditions extracted from 2D phase contrast MRI scans. The study included carotid bifurcations with CaW (n = 13), mild atherosclerosis (n = 7), and healthy bifurcation geometries (n = 6). Hemodynamic parameters associated with vascular dysfunction and clot formation, including shear rate, oscillatory shear index (OSI), low velocity, and flow stasis were calculated and compared between the subject groups. Patients with CaW had significantly larger regions containing low shear rate, high OSI, low velocity, and flow stasis in comparison to subjects with mild atherosclerosis or normal bifurcations. These abnormal hemodynamic metrics in patients with CaW are associated with clot formation and vascular dysfunction and suggest that hemodynamic assessment may be a tool to assess stroke risk in these patients.

Impact of Age-Related Change in Caval Flow Ratio on Hepatic Flow Distribution in the Fontan Circulation.

The Fontan operation is a palliative technique for patients born with single ventricle heart disease. The superior vena cava (SVC), inferior vena cava (IVC), and hepatic veins are connected to the pulmonary arteries in a total cavopulmonary connection by an extracardiac conduit or a lateral tunnel connection. A balanced hepatic flow distribution (HFD) to both lungs is essential to prevent pulmonary arteriovenous malformations and cyanosis. HFD is highly dependent on the local hemodynamics. The effect of age-related changes in caval inflows on HFD was evaluated using cardiac magnetic resonance data and patient-specific computational fluid dynamics modeling. SVC and IVC flow from 414 patients with Fontan were collected to establish a relationship between SVC:IVC flow ratio and age. Computational fluid dynamics modeling was performed in 60 (30 extracardiac and 30 lateral tunnel) patient models to quantify the HFD that corresponded to patient ages of 3, 8, and 15 years, respectively. SVC:IVC flow ratio inverted at ≈8 years of age, indicating a clear shift to lower body flow predominance. Our data showed that variation of HFD in response to age-related changes in caval inflows (SVC:IVC, 2, 1, and 0.5 corresponded to ages, 3, 8, and 15+, respectively) was not significant for extracardiac but statistically significant for lateral tunnel cohorts. For all 3 caval inflow ratios, a positive correlation existed between the IVC flow distribution to both the lungs and the HFD. However, as the SVC:IVC ratio changed from 2 to 0.5 (age, 3-15+) years, the correlation's strength decreased from 0.87 to 0.64, due to potential flow perturbation as IVC flow momentum increased. Our analysis provided quantitative insights into the impact of the changing caval inflows on Fontan's long-term HFD, highlighting the importance of SVC:IVC variations over time on Fontan's long-term hemodynamics. These findings broaden our understanding of Fontan hemodynamics and patient outcomes.

A parametric study of the effect of 3D plaque shape on local hemodynamics and implications for plaque instability

The vast majority of heart attacks occur when vulnerable plaques rupture, releasing their lipid content into the blood stream leading to thrombus formation and blockage of a coronary artery. Detection of these unstable plaques before they rupture remains a challenge. Hemodynamic features including wall shear stress (WSS) and wall shear stress gradient (WSSG) near the vulnerable plaque and local inflammation are known to affect plaque instability. In this work, a computational workflow has been developed to enable a comprehensive parametric study detailing the effects of 3D plaque shape on local hemodynamics and their implications for plaque instability. Parameterized geometric 3D plaque models are created within a patient-specific coronary artery tree using a NURBS (non-uniform rational B-splines)-based vascular modeling pipeline. Realistic blood flow features are simulated by using a Navier–Stokes solver within an isogeometric finite-element analysis framework. Near wall hemodynamic quantities such as WSS and WSSG are quantified, and vascular distribution of an inflammatory marker (VCAM-1) is estimated. Results show that proximally skewed eccentric plaques have the most vulnerable combination of high WSS and high positive spatial WSSG, and the presence of multiple lesions increases risk of rupture. The computational tool developed in this work, in conjunction with clinical data, -could help identify surrogate markers of plaque instability, potentially leading to a noninvasive clinical procedure for the detection of vulnerable plaques before rupture.

Design of experiment approach to identify the dominant geometrical feature of left coronary artery influencing atherosclerosis

Background and Objective. Coronary artery geometry heavily influences local hemodynamics, potentially leading to atherosclerosis. Consequently, the unique geometrical configuration of an individual by birth can be associated with future risk of atherosclerosis. Although current researches focus on exploring the relationship between local hemodynamics and coronary artery geometry, this study aims to identify the order of influence of the geometrical features through systematic experiments, which can reveal the dominant geometrical feature for future risk assessment. Methods. According to Taguchi’s method of design of experiment (DoE), the left main stem (LMS) length (l LMS), curvature (k LMS), diameter (d LMS) and the bifurcation angle between left anterior descending (LAD) and left circumflex (LCx) artery (α LAD-LCx) of two reconstructed patient-specific left coronary arteries (LCA) were varied in three levels to create L9 orthogonal array. Computational fluid dynamic (CFD) simulations with physiological boundary conditions were performed on the resulting eighteen LCA models. Average helicity intensity (h 2) and relative atheroprone area (RAA) of near-wall hemodynamic descriptors were analyzed. Results. The proximal LAD (LAD proximal) was identified to be the most atheroprone region of the left coronary artery due to higher h 2, large RAA of time averaged wall shear stress (TAWSS < 0.4 Pa), oscillatory shear index (OSI ∼ 0.5) and relative residence time (RRT > 4.17 Pa−1). In both patient-specific cases, based on h 2 and TAWSS, d lms is the dominant geometric parameter while based on OSI and RRT, α LAD-LCx is the dominant one influencing hemodynamic condition in proximal LAD (p < 0.05). Based on RRT, the rank of the geometrical factors is: α LAD-LCx > d LMS > l LMS > k LMS, indicating that α LAD-LCx is the most dominant geometrical factor affecting hemodynamics at proximal LAD which may influence atherosclerosis. Conclusion. The proposed identification of the rank of geometrical features of LCA and the dominant feature may assist clinicians in predicting the possibility of atherosclerosis, of an individual, long before it will occur. This study can further be translated to be used to rank the influence of several arterial geometrical features at different arterial locations to explore detailed relationships between the arterial geometrical features and local hemodynamics.

Modeling Anisotropic Electrical Conductivity of Blood: Translating Microscale Effects of Red Blood Cell Motion into a Macroscale Property of Blood.

Cardiovascular diseases are a leading global cause of mortality. The current standard diagnostic methods, such as imaging and invasive procedures, are relatively expensive and partly connected with risks to the patient. Bioimpedance measurements hold the promise to offer rapid, safe, and low-cost alternative diagnostic methods. In the realm of cardiovascular diseases, bioimpedance methods rely on the changing electrical conductivity of blood, which depends on the local hemodynamics. However, the exact dependence of blood conductivity on the hemodynamic parameters is not yet fully understood, and the existing models for this dependence are limited to rather academic flow fields in straight pipes or channels. In this work, we suggest two closely connected anisotropic electrical conductivity models for blood in general three-dimensional flows, which consider the orientation and alignment of red blood cells (RBCs) in shear flows. In shear flows, RBCs adopt preferred orientations through a rotation of their membrane known as tank-treading motion. The two models are built on two different assumptions as to which hemodynamic characteristic determines the preferred orientation. The models are evaluated in two example simulations of blood flow. In a straight rigid vessel, the models coincide and are in accordance with experimental observations. In a simplified aorta geometry, the models yield different results. These differences are analyzed quantitatively, but a validation of the models with experiments is yet outstanding.

Partitioning of dense RBC suspensions in single microfluidic bifurcations: role of cell deformability and bifurcation angle

Red blood cells (RBCs) are a key determinant of human physiology and their behaviour becomes extremely heterogeneous as they navigate in narrow, bifurcating vessels in the microvasculature, affecting local haemodynamics. This is due to partitioning in bifurcations which is dependent on the biomechanical properties of RBCs, especially deformability. We examine the effect of deformability on the haematocrit distributions of dense RBC suspensions flowing in a single, asymmetric Y-shaped bifurcation, experimentally. Human RBC suspensions (healthy and artificially hardened) at 20% haematocrit (Ht) were perfused through the microchannels at different flow ratios between the outlet branches, and negligible inertia, and imaged to infer cell distributions. Notable differences in the shape of the haematocrit distributions were observed between healthy and hardened RBCs near the bifurcation apex. These lead to more asymmetric distributions for healthy RBCs in the daughter and outlet branches with cells accumulating near the inner channel walls, exhibiting distinct hematocrit peaks which are sharper for healthy RBCs. Although the hematocrit distributions differed locally, similar partitioning characteristics were observed for both suspensions. Comparisons with RBC distributions measured in a T-shaped bifurcation showed that the bifurcation angle affects the haematocrit characteristics of the healthy RBCs and not the hardened ones. The extent of RBC partitioning was found similar in both geometries and suspensions. The study highlights the differences between local and global characteristics which impact RBC distribution in more complex, multi-bifurcation networks.