Application of multiphase computational fluid dynamics to analyze monocyte adhesion.

Study of the mechanisms of monocyte adhesion initiating atheroslerotic lesions has engaged investigators for decades. Single-phase computational fluid dynamics (CFD) analyses fail to account for particulate migration. Consequently, inconsistencies arise when correlating adhesion with wall shear stress (WSS). The purpose of this paper is to present, to our knowledge, the first computational analysis of in vitro U937 monocyte-like human cell adhesion data using a coupled multiphase CFD-population balance adhesion model. The CFD model incorporates multiphase non-Newtonian hemodynamic models to compute the spatial distributions of freely flowing monocytes and WSSs in control volumes adjacent to the wall. Measurements of monocyte adhesion onto an E-selectin-coated flow model that included an idealized stenosis and an abrupt expansion were available from the literature. In this study, we develop a new monolayer population balance adhesion model, based on the widely accepted mechanism of ligand-receptor binding, coupled to the CFD results. The monolayer population balance model accounts for the interactions of freely flowing, rolling, and adhering monocytes with surfaces via first-order reactions, transport of rolling cells in the monolayer, and the concept of a WSS detachment threshold, clearly evident in the adhesion experiments. The new paradigm of coupling the multiphase hemodynamic CFD model with the proposed adhesion model is illustrated by determining and interpreting the model parameters for experimental datasets having Reynolds numbers of 100 and 140. The coupled multiphase CFD adhesion model is able to simultaneously predict the spatial variations in freely flowing monocytes, their adherent number density, and carrier fluid WSSs adjacent to ligand-coated flow cell surfaces.

Pulsatile flow in a coronary artery using multiphase kinetic theory.

Pulsatile flow in a model of a right coronary artery (RCA) was previously modeled as a single-phase fluid and as a two-phase fluid using experimental rheological data for blood as a function of hematocrit and shear rate. Here we present a multiphase kinetic theory model which has been shown to compute correctly the viscosity of red blood cells (RBCs) and their migration away from vessel walls: the Fahraeus-Lindqvist effect. The computed RBC viscosity decreases with shear rate and vessel size, consistent with measurements. The pulsatile computations were performed using a typical cardiac waveform until a limit cycle was well established. The RBC volume fractions, shear stresses, shear stress gradients, granular temperatures, viscosities, and phase velocities varied with time and position during each cardiac cycle. Steady-state computations were also performed and were found to compare well with time-averaged transient results. The wall shear stress and wall shear stress gradients (both spatial and temporal) were found to be highest on the inside area of maximum curvature. Potential atherosclerosis sites are identified using these computational results.

Hemodynamic computation using multiphase flow dynamics in a right coronary artery.

Hemodynamic data on the roles of physiologically critical blood particulates are needed to better understand cardiovascular diseases. The blood flow patterns and particulate buildup were numerically simulated using the multiphase non-Newtonian theory of dense suspension hemodynamics in a realistic right coronary artery (RCA) having various cross sections. The local hemodynamic factors, such as wall shear stress (WSS), red blood cell (RBC) buildup, viscosity, and velocity, varied with the spatially nonuniform vessel structures and temporal cardiac cycles. The model generally predicted higher RBC buildup on the inside radius of curvature. A low WSS region was found in the high RBC buildup region, in particular, on the area of maximum curvature of a realistic human RCA. The complex recirculation patterns, the oscillatory flow with flow reversal, and vessel geometry resulted in RBC buildup due to the prolonged particulate residence time, specifically, at the end of the diastole cycle. The increase of the initial plasma viscosity caused the lower WSS. These predictions have significant implications for understanding the local hemodynamic phenomena that may contribute to the earliest stage of atherosclerosis, as clinically observed on the inside curvatures and torsion of coronary arteries.

Multiphase hemodynamic simulation of pulsatile flow in a coronary artery.

A multiphase transient non-Newtonian three-dimensional (3-D) computational fluid dynamics (CFD) simulation has been performed for pulsatile hemodynamics in an idealized curved section of a human coronary artery. We present the first prediction, to the authors' knowledge, of particulate buildup on the inside curvature using the multiphase theory of dense suspension hemodynamics. In this study, the particulates are red blood cells (RBCs). The location of RBC buildup on the inside curvature correlates with lower wall shear stress (WSS) relative to the outside curvature. These predictions provide insight into how blood-borne particulates interact with artery walls and hence, have relevance for understanding atherogenesis since clinical observations show that atherosclerotic plaques generally form on the inside curvatures of arteries. The buildup of RBCs on the inside curvature is driven by the secondary flow and higher residence times. The higher viscosity in the central portion of the curved vessel tends to block their flow, causing them to migrate preferentially through the boundary layer. The reason for this is the nearly neutrally buoyant nature of the dense two-phase hemodynamic flow. The two-phase non-Newtonian viscosity model predicts greater shear thinning than the single-phase non-Newtonian model. Consequently, the secondary flow induced in the curvature is weaker. The waveforms for computed hemodynamic parameters, such as hematocrit, WSS, and viscosity, follow the prescribed inlet velocity waveforms. The lower oscillatory WSS produced on the inside curvature has implications for understanding thickening of the intimal layer.