A numerical study of blood flow using mixture theory

A homogenized two-phase computational framework for meso- and macroscale blood flow simulations

BACKGROUND AND OBJECTIVE

Owing to the complexity of physics linked with blood flow and its associated phenomena, appropriate modeling of the multi-constituent rheology of blood is of primary importance. To this effect, various kinds of computational fluid dynamic models have been developed, each with merits and limitations. However, when additional physics like thrombosis and embolization is included within the framework of these models, computationally efficient scalable translation becomes very difficult. Therefore, this paper presents a homogenized two-phase blood flow framework with similar characteristics to a single fluid model but retains the flow resolution of a classical two-fluid model. The presented framework is validated against four different sets of experiments.

METHODS

The two-phase model of blood presented here is based on the classical diffusion-flux framework. Diffusion flux models are known to be less computationally expensive than two-fluid multiphase models since the numerical implementation resembles single-phase flow models. Diffusion flux models typically use empirical slip velocity correlations to resolve the motion between phases. However, such correlations do not exist for blood. Therefore, a modified slip velocity equation is proposed, derived rigorously from the two-fluid governing equations. An additional drag law for red blood cells (RBCs) as a function of volume fraction is evaluated using a previously published cell-resolved solver. A new hematocrit-dependent expression for lift force on RBCs is proposed. The final governing equations are discretized and solved using the open-source software OpenFOAM.

RESULTS

The framework is validated against four sets of experiments: (i) flow through a rectangular microchannel to validate RBC velocity profiles against experimental measurements and compare computed hematocrit distributions against previously reported simulation results (ii) flow through a sudden expansion microchannel for comparing experimentally obtained contours of hematocrit distributions and normalized cell-free region length obtained at different flowrates and inlet hematocrits, (iii) flow through two hyperbolic channels to evaluate model predictions of cell-free layer thickness, and (iv) flow through a microchannel that mimics crevices of a left ventricular assist device to predict hematocrit distributions observed experimentally. The simulation results exhibit good agreement with the results of all four experiments.

CONCLUSION

The computational framework presented in this paper has the advantage of resolving the multiscale physics of blood flow while still leveraging numerical techniques used for solving single-phase flows. Therefore, it becomes an excellent candidate for addressing more complicated problems related to blood flow, such as modeling mechanical entrapment of RBCs within blood clots, predicting thrombus composition, and visualizing clot embolization.

Effects of cardiac function alterations on the risk of postoperative thrombotic complications in patients receiving endovascular aortic repair

Introduction: Chronic heart disease (CHD) is a common comorbidity of patients receiving endovascular aneurysm repair (EVAR) for abdominal aortic aneurysms (AAA). The explicit relationship between ventricular systolic function and EVAR complication of thrombotic events is unknown. Methods: In this study, we proposed a three-dimensional numerical model coupled with the lumped-elements heart model, which is capable of simulating thrombus formation in diverse systolic functions. The relation of cardiac functions and the predicted risk of thrombus formation in the aorta and/or endograft of 4 patients who underwent EVAR was investigated. Relative risks for thrombus formation were identified using machine-learning algorithms. Results: The computational results demonstrate that thrombus tended to form on the interior side of the aorta arch and iliac branches, and cardiac function can affect blood flow field and affect thrombus formation, which is consistent with the four patients' post-operative imaging follow-up. We also found that RRT, OSI, TAWSS in thrombosis area are lower than whole average. In addition, we found that the thrombus formation has negative correlations with the maximum ventricular contractile force (r = -.281 ± .101) and positive correlations with the minimum ventricular contractile force (r = .238 ± .074), whereas the effect of heart rate (r = -.015 ± .121) on thrombus formation is not significant. Conclusion: In conclusion, changes in ventricular systolic function may alter the risk of thrombotic events after EVAR repair, which could provide insight into the selection of adjuvant therapy strategies for AAA patients with CHD.

von Willebrand factor unfolding mediates platelet deposition in a model of high-shear thrombosis

Thrombosis under high-shear conditions is mediated by the mechanosensitive blood glycoprotein von Willebrand factor (vWF). vWF unfolds in response to strong flow gradients and facilitates rapid recruitment of platelets in flowing blood. While the thrombogenic effect of vWF is well recognized, its conformational response in complex flows has largely been omitted from numerical models of thrombosis. We recently presented a continuum model for the unfolding of vWF, where we represented vWF transport and its flow-induced conformational change using convection-diffusion-reaction equations. Here, we incorporate the vWF component into our multi-constituent model of thrombosis, where the local concentration of stretched vWF amplifies the deposition rate of free-flowing platelets and reduces the shear cleaning of deposited platelets. We validate the model using three benchmarks: in vitro model of atherothrombosis, a stagnation point flow, and the PFA-100, a clinical blood test commonly used for screening for von Willebrand disease (vWD). The simulations reproduced the key aspects of vWF-mediated thrombosis observed in these experiments, such as the thrombus location, thrombus growth dynamics, and the effect of blocking platelet-vWF interactions. The PFA-100 simulations closely matched the reported occlusion times for normal blood and several hemostatic deficiencies, namely, thrombocytopenia, vWD type 1, and vWD type 3. Overall, this multi-constituent model of thrombosis enables macro-scale 3D simulations of thrombus formation in complex geometries over a wide range of shear rates and accounts for qualitative and quantitative hemostatic deficiencies in patient blood.

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.

Suitable CO2 Solubility Models for Determination of the CO2 Removal Performance of Oxygenators

CO₂ removal via membrane oxygenators during lung protective ventilation has become a reliable clinical technique. For further optimization of oxygenators, accurate prediction of the CO₂ removal rate is necessary. It can either be determined by measuring the CO₂ content in the exhaust gas of the oxygenator (sweep flow-based) or using blood gas analyzer data and a CO₂ solubility model (blood-based). In this study, we determined the CO₂ removal rate of a prototype oxygenator utilizing both methods in in vitro trials with bovine and in vivo trials with porcine blood. While the sweep flow-based method is reliably accurate, the blood-based method depends on the accuracy of the solubility model. In this work, we quantified performances of four different solubility models by calculating the deviation of the CO₂ removal rates determined by both methods. Obtained data suggest that the simplest model (Loeppky) performs better than the more complex ones (May, Siggaard-Anderson, and Zierenberg). The models of May, Siggaard-Anderson, and Zierenberg show a significantly better performance for in vitro bovine blood data than for in vivo porcine blood data. Furthermore, the suitability of the Loeppky model parameters for bovine blood (in vitro) and porcine blood (in vivo) is evaluated.

Hemodynamic study in 3D printed stenotic coronary artery models: experimental validation and transient simulation

Atherosclerosis is a progressive disease that can significantly reduce blood supply to vital organs, being one of the main causes of death worldwide. In this work, a numerical and experimental study in 3D printed stenotic coronary arteries, considering both steady and pulsatile blood flow conditions, is presented. The results revealed that a degree of stenosis superior to 50% creates disturbed flows downstream of the contraction, with an accented increase in the wall shear stress measurements at the stenosis throat. Finally, the multiphase mixture was investigated and compared with a single-phase modelling, and only slight differences were observed right after the stenosis throat.

Simulation of blood flow in a sudden expansion channel and a coronary artery

In this paper, we numerically simulate the flow of blood in two benchmark problems: the flow in a sudden expansion channel and the flow through an idealized curved coronary artery with pulsatile inlet velocity. Blood is modeled as a suspension (a non-linear complex fluid) and the movement of the red blood cell (RBCs) is modeled by using a concentration flux equation. The viscosity of blood is obtained from experimental data. In the sudden expansion flow, the predicted velocity profiles for two different Reynolds numbers (based on the inlet velocity) agree well with the available experiments; furthermore, the numerical results also show that after the sudden expansion there exists a RBCs depletion region. For the second problem, the idealized curved coronary artery, it is found that the RBCs move towards and concentrate near the inner surface where the viscosity is higher and the shear stress lower; this phenomenon may be related to the atherosclerotic plaque formation which usually occurs on the inside surface of the arteries.

Simulation of thrombosis in a stenotic microchannel: The effects of vWF-enhanced shear activation of platelets

This study was undertaken to develop a numerical/computational simulation of von Willebrand Factor (vWF) - mediated platelet shear activation and deposition in an idealized stenosis. Blood is treated as a multi-constituent mixture comprised of a linear fluid component and a porous solid component (thrombus). Chemical and biological species involved in coagulation are modeled using a system of coupled convection-reaction-diffusion (CRD) equations. This study considers the cumulative effect of shear stress (history) on platelet activation. The vWF activity is modeled as an enhancement function for the shear stress accumulation and is related to the experimentally-observed unfolding rate of vWF. A series of simulations were performed in an idealized stenosis in which the predicted platelets deposition agreed well with previous experimental observations spatially and temporally, including the reduction of platelet deposition with decreasing expansion angle. Further simulation indicated a direct relationship between vWF-mediated platelet deposition and degree of stenosis. Based on the success with these benchmark simulations, it is hoped that the model presented here may provide additional insight into vWF-mediated thrombosis and prove useful for the development of more hemo-compatible blood-wetted devices in the future.

A Review of Rheological Modeling of Cement Slurry in Oil Well Applications

The rheological behavior of cement slurries is important in trying to prevent and eliminate gas-migration related problems in oil well applications. In this paper, we review the constitutive modeling of cement slurries/pastes. Cement slurries, in general, behave as complex non-linear fluids with the possibility of exhibiting viscoelasticity, thixotropy, yield stress, shear-thinning effects, etc. The shear viscosity and the yield stress are two of the most important rheological characteristics of cement; these have been studied extensively and a review of these studies is provided in this paper. We discuss the importance of changing the concentration of cement particles, water-to-cement ratio, additives/admixtures, shear rate, temperature and pressure, mixing methods, and the thixotropic behavior of cement on the stress tensor. In the concluding remarks, we propose a new constitutive model for cement slurry, considering the basic non-Newtonian nature of the different models.

A non-linear fluid suspension model for blood flow

Motivated by the complex rheological behaviors observed in small/micro scale blood vessels, such as the Fahraeus effect, plasma-skimming, shear-thinning, etc., we develop a non-linear suspension model for blood. The viscosity is assumed to depend on the volume fraction (hematocrit) and the shear rate. The migration of the red blood cells (RBCs) is studied using a concentration flux equation. A parametric study with two representative problems, namely simple shear flow and a pressure driven flow demonstrate the ability of this reduced-order model to reproduce several key features of the two-fluid model (mixture theory approach), with much lower computational cost.

Multi-Constituent Simulation of Thrombus Deposition

In this paper, we present a spatio-temporal mathematical model for simulating the formation and growth of a thrombus. Blood is treated as a multi-constituent mixture comprised of a linear fluid phase and a thrombus (solid) phase. The transport and reactions of 10 chemical and biological species are incorporated using a system of coupled convection-reaction-diffusion (CRD) equations to represent three processes in thrombus formation: initiation, propagation and stabilization. Computational fluid dynamic (CFD) simulations using the libraries of OpenFOAM were performed for two illustrative benchmark problems: in vivo thrombus growth in an injured blood vessel and in vitro thrombus deposition in micro-channels (1.5 mm × 1.6 mm × 0.1 mm) with small crevices (125 μm × 75 μm and 125 μm × 137 μm). For both problems, the simulated thrombus deposition agreed very well with experimental observations, both spatially and temporally. Based on the success with these two benchmark problems, which have very different flow conditions and biological environments, we believe that the current model will provide useful insight into the genesis of thrombosis in blood-wetted devices, and provide a tool for the design of less thrombogenic devices.

Hemodynamic modeling of leukocyte and erythrocyte transport and interactions in intracranial aneurysms by a multiphase approach

Hemodynamics has been recognized as an important factor in the development, growth, and rupture of cerebral aneurysms, and investigated by computational fluid dynamics techniques using a single phase approach. However, flow-dependent cell transport and interactions are usually ignored in single phase models, in which blood is usually treated as a single phase Newtonian fluid. For getting better insight into the underlying pathology of intracranial aneurysm, cell transport and interactions should be covered in hemodynamic studies. In the present study, a multiphase hemodynamic model incorporating cell transport and interactions was developed, in which blood was modeled as multiphase fluid having a continuous phase (plasma) and two particulate phases (erythrocytes and leukocytes). The model showed good agreement with experimental data and observations in the literature, and was applied to four patient-specific aneurysms in a pulsatile manner. Leukocyte accumulations were predicted at locations with flow disturbance and low wall shear stress. The concentrations of leukocyte at accumulation sites were found to exceed 200 to 500% of normal physiological level on three unstable aneurysms, including two ruptured aneurysms and a growing aneurysm where accumulation was observed near a daughter sac and a secondary aneurysm. This suggested that aneurysms with complex secondary flow patterns could be prone to leukocyte accumulation on the wall. As this is the first study to characterize cell transport and interactions in aneurysm hemodynamics, our model can serve as a foundation for future intracranial aneurysm models.

Design of microfluidic channels for magnetic separation of malaria-infected red blood cells

This study is motivated by the development of a blood cell filtration device for removal of malaria-infected, parasitized red blood cells (pRBCs). The blood was modeled as a multi-component fluid using the computational fluid dynamics discrete element method (CFD-DEM), wherein plasma was treated as a Newtonian fluid and the red blood cells (RBCs) were modeled as soft-sphere solid particles which move under the influence of drag, collisions with other RBCs, and a magnetic force. The CFD-DEM model was first validated by a comparison with experimental data from Han et al. 2006 (Han and Frazier 2006) involving a microfluidic magnetophoretic separator for paramagnetic deoxygenated blood cells. The computational model was then applied to a parametric study of a parallel-plate separator having hematocrit of 40% with a 10% of the RBCs as pRBCs. Specifically, we investigated the hypothesis of introducing an upstream constriction to the channel to divert the magnetic cells within the near-wall layer where the magnetic force is greatest. Simulations compared the efficacy of various geometries upon the stratification efficiency of the pRBCs. For a channel with nominal height of 100 µm, the addition of an upstream constriction of 80% improved the proportion of pRBCs retained adjacent to the magnetic wall (separation efficiency) by almost 2 fold, from 26% to 49%. Further addition of a downstream diffuser reduced remixing, hence improved separation efficiency to 72%. The constriction introduced a greater pressure drop (from 17 to 495 Pa), which should be considered when scaling-up this design for a clinical-sized system. Overall, the advantages of this design include its ability to accommodate physiological hematocrit and high throughput - which is critical for clinical implementation as a blood-filtration system.

Study of blood flow in several benchmark micro-channels using a two-fluid approach

It is known that in a vessel whose characteristic dimension (e.g., its diameter) is in the range of 20 to 500 microns, blood behaves as a non-Newtonian fluid, exhibiting complex phenomena, such as shear-thinning, stress relaxation, and also multi-component behaviors, such as the Fahraeus effect, plasma-skimming, etc. For describing these non-Newtonian and multi-component characteristics of blood, using the framework of mixture theory, a two-fluid model is applied, where the plasma is treated as a Newtonian fluid and the red blood cells (RBCs) are treated as shear-thinning fluid. A computational fluid dynamic (CFD) simulation incorporating the constitutive model was implemented using OpenFOAM® in which benchmark problems including a sudden expansion and various driven slots and crevices were studied numerically. The numerical results exhibited good agreement with the experimental observations with respect to both the velocity field and the volume fraction distribution of RBCs.