Search

Blood flowing through microvascular bifurcations has been an active research topic for many decades, while the partitioning pattern of nanoscale solutes in the blood remains relatively unexplored. Here we demonstrate a multiscale computational framework for direct numerical simulation of the nanoparticle (NP) partitioning through physiologically relevant vascular bifurcations in the presence of red blood cells (RBCs). The computational framework is established by embedding a particulate suspension inflow-outflow boundary condition into a multiscale blood flow solver. The computational framework is verified by recovering a tubular blood flow without a bifurcation and validated against the experimental measurement of an intravital bifurcation flow. The classic Zweifach-Fung (ZF) effect is shown to be well captured by the method. Moreover, we observe that NPs exhibit a ZF-like heterogeneous partition in response to the heterogeneous partition of the RBC phase. The NP partitioning prioritizes the high-flow-rate daughter branch except for extreme (large or small) suspension flow partition ratios under which the complete phase separation tends to occur. By analyzing the flow field and the particle trajectories, we show that the ZF-like heterogeneity in the NP partition can be explained by the RBC-entrainment effect caused by the deviation of the flow separatrix preceded by the tank treading of RBCs near the bifurcation junction. The recovery of homogeneity in the NP partition under extreme flow partition ratios is due to the plasma skimming of NPs in the cell-free layer. These findings, based on the multiscale computational framework, provide biophysical insights to the heterogeneous distribution of NPs in microvascular beds that are observed pathophysiologically.

A theory is presented for the transition between the coflowing and the drop-generation regimes observed in microfluidic channels with a rectangular cross section. This transition is characterized by a critical ratio of the dispersed- to continuous-phase volume flow rates. At flow-rate ratios higher than this critical value, drop generation is suppressed. The critical ratio corresponds to the fluid cross section where the dispersed-phase fluid is just tangent to the channel walls. The transition criterion is a function of the ratio of the fluid viscosities, the three-phase contact angle formed between the fluid phases and the channel walls, and the aspect ratio of the channel cross section; the transition is independent of interfacial tension. Hysteretic behavior of drop generation with respect to the flow-rate ratio is predicted for partially wetting dispersed-phase fluids. Experimental data are consistent with this theory.

A thin flow-focusing microfluidic channel is evaluated for generating monodisperse liquid droplets. The microfluidic device is used in its native state, which is hydrophilic, or treated with OTS to make it hydrophobic. Having both hydrophilic and hydrophobic surfaces allows for creation of both oil-in-water and water-in-oil emulsions, facilitating a large parameter study of viscosity ratios (droplet fluid/continuous fluid) ranging from 0.05 to 96 and flow rate ratios (droplet fluid/continuous fluid) ranging from 0.01 to 2 in one geometry. The hydrophilic chip provides a partially-wetting surface (contact angle less than 90°) for the inner fluid. This surface, combined with the unusually thin channel height, promotes a flow regime where the inner fluid wets the top and bottom of the channel in the orifice and a stable jet is formed. Through confocal microscopy, this fluid stabilization is shown to be highly influenced by the contact angle of the liquids in the channel. Non-wetting jets undergo breakup and produce drops when the jet is comparable to or smaller than the channel thickness. In contrast, partially-wetting jets undergo breakup only when they are much smaller than the channel thickness. Drop sizes are found to scale with a modified capillary number based on the total flow rate regardless of wetting behavior.